Affinities and selectivities of divalent cation binding sites within an RNA tertiary structure

TPP Riboswitch Populates Holo-Form-like Structure Even in the Absence of Cognate Ligand at High Mg2+ Concentration

Riboswitches are noncoding RNA that regulate gene expression by folding into specific three-dimensional structures (holo-form) upon binding by their cognate ligand in the presence of Mg²⁺. Riboswitch functioning is also hypothesized to be under kinetic control requiring large cognate ligand concentrations. We ask the question under thermodynamic conditions, can the riboswitches populate structures similar to the holo-form only in the presence of Mg²⁺ and absence of cognate ligand binding. We addressed this question using thiamine pyrophosphate (TPP) riboswitch as a model system and computer simulations using a coarse-grained model for RNA. The folding free energy surface (FES) shows that with the initial increase in Mg²⁺ concentration ([Mg²⁺]), the aptamer domain (AD) of TPP riboswitch undergoes a barrierless collapse in its dimensions. On further increase in [Mg²⁺], intermediates separated by barriers appear on the FES, and one of the intermediates has a TPP ligand-binding competent structure. We show that site-specific binding of the Mg²⁺ aids in the formation of tertiary contacts. For [Mg²⁺] greater than physiological concentration, AD folds into a structure similar to the crystal structure of the TPP holo-form even in the absence of the TPP ligand. The folding kinetics shows that TPP AD populates an intermediate due to the misalignment of two arms present in the structure, which acts as a kinetic trap, leading to larger folding timescales. The predictions of the intermediate structures from the simulations are amenable for experimental verification.

The structural plasticity of nucleic acid duplexes revealed by WAXS and MD

Double-stranded DNA (dsDNA) and RNA (dsRNA) helices display an unusual structural diversity. Some structural variations are linked to sequence and may serve as signaling units for protein-binding partners. Therefore, elucidating the mechanisms and factors that modulate these variations is of fundamental importance. While the structural diversity of dsDNA has been extensively studied, similar studies have not been performed for dsRNA. Because of the increasing awareness of RNA's diverse biological roles, such studies are timely and increasingly important. We integrate solution x-ray scattering at wide angles (WAXS) with all-atom molecular dynamics simulations to explore the conformational ensemble of duplex topologies for different sequences and salt conditions. These tightly coordinated studies identify robust correlations between features in the WAXS profiles and duplex geometry and enable atomic-level insights into the structural diversity of DNA and RNA duplexes. Notably, dsRNA displays a marked sensitivity to the valence and identity of its associated cations.

Salt dependent mesoscopic model for RNA at multiple strand concentrations

Mesoscopic models can be used for the description of the thermodynamic properties of RNA duplexes. With the use of experimental melting temperatures, its parametrization can provide important insights into its hydrogen bonds and stacking interactions as has been done for high sodium concentrations. However, the RNA parametrization for lower salt concentrations is still missing due to the limited amount of published melting temperature data. While the Peyrard-Bishop (PB) parametrization was found to be largely independent of strand concentrations, it requires that all temperatures are provided at the same strand concentrations. Here we adapted the PB model to handle multiple strand concentrations and in this way we were able to make use of an experimental set of temperatures to model the hydrogen bond and stacking interactions at low and intermediate sodium concentrations. For the parametrizations we make a distinction between terminal and internal base pairs, and the resulting potentials were qualitatively similar as we obtained previously for DNA. The main difference from DNA parameters, was the Morse potentials at low sodium concentrations for terminal r(AU) which is stronger than d(AT), suggesting higher hydrogen bond strength.

Ribosomal Protein L11 Selectively Stabilizes a Tertiary Structure of the GTPase Center rRNA Domain

The GTPase Center (GAC) RNA domain in bacterial 23S rRNA is directly bound by ribosomal protein L11, and this complex is essential to ribosome function. Previous cocrystal structures of the 58-nucleotide GAC RNA bound to L11 revealed the intricate tertiary fold of the RNA domain, with one monovalent and several divalent ions located in specific sites within the structure. Here, we report a new crystal structure of the free GAC that is essentially identical to the L11-bound structure, which retains many common sites of divalent ion occupation. This new structure demonstrates that RNA alone folds into its tertiary structure with bound divalent ions. In solution, we find that this tertiary structure is not static, but rather is best described as an ensemble of states. While L11 protein cannot bind to the GAC until the RNA has adopted its tertiary structure, new experimental data show that L11 binds to Mg²⁺-dependent folded states, which we suggest lie along the folding pathway of the RNA. We propose that L11 stabilizes a specific GAC RNA tertiary state, corresponding to the crystal structure, and that this structure reflects the functionally critical conformation of the rRNA domain in the fully assembled ribosome.

Competitive Binding of Mg2+ and Na+ Ions to Nucleic Acids: From Helices to Tertiary Structures

Nucleic acids generally reside in cellular aqueous solutions with mixed divalent/monovalent ions, and the competitive binding of divalent and monovalent ions is critical to the structures of nucleic acids because of their polyanionic nature. In this work, we first proposed a general and effective method for simulating a nucleic acid in mixed divalent/monovalent ion solutions with desired bulk ion concentrations via molecular dynamics (MD) simulations and investigated the competitive binding of Mg²⁺/Na⁺ ions to various nucleic acids by all-atom MD simulations. The extensive MD-based examinations show that single MD simulations conducted using the proposed method can yield desired bulk divalent/monovalent ion concentrations for various nucleic acids, including RNA tertiary structures. Our comprehensive analyses show that the global binding of Mg²⁺/Na⁺ to a nucleic acid is mainly dependent on its structure compactness, as well as Mg²⁺/Na⁺ concentrations, rather than the specific structure of the nucleic acid. Specifically, the relative global binding of Mg²⁺ over Na⁺ is stronger for a nucleic acid with higher effective surface charge density and higher relative Mg²⁺/Na⁺ concentrations. Furthermore, the local binding of Mg²⁺/Na⁺ to a phosphate of a nucleic acid mainly depends on the local phosphate density in addition to Mg²⁺/Na⁺ concentrations.

Divalent ions tune the kinetics of a bacterial GTPase center rRNA folding transition from secondary to tertiary structure

Folding of an RNA from secondary to tertiary structure often depends on divalent ions for efficient electrostatic charge screening (nonspecific association) or binding (specific association). To measure how different divalent cations modify folding kinetics of the 60 nucleotide Ecoli rRNA GTPase center, we combined stopped-flow fluorescence in the presence of Mg²⁺, Ca²⁺, or Sr²⁺ together with time-resolved small angle X-ray scattering (SAXS) in the presence of Mg²⁺ to observe the folding process. Immediately upon addition of each divalent ion, the RNA undergoes a transition from an extended state with secondary structure to a more compact structure. Subsequently, specific divalent ions modulate populations of intermediates in conformational ensembles along the folding pathway with transition times longer than 10 msec. Rate constants for the five folding transitions act on timescales from submillisecond to tens of seconds. The sensitivity of RNA tertiary structure to divalent cation identity affects all but the fastest events in RNA folding, and allowed us to identify those states that prefer Mg²⁺ The GTPase center RNA appears to have optimized its folding trajectory to specifically utilize this most abundant intracellular divalent ion.

RNA as a Complex Polymer with Coupled Dynamics of Ions and Water in the Outer Solvation Sphere

We unravel the internal and collective modes of a widely studied 58-nucleotide rRNA fragment in solvent using atomically detailed molecular dynamics simulations. The variation of lifetimes for water hydrogen bonds with nucleotide groups indicates heterogeneity of water dynamics on the RNA surface. The time scales of interactions of the discrete water molecules with RNA nucleotides extend from several hundred picoseconds to a few nanoseconds. We determine all of the association sites and the spatial distribution of residence times for Mg²⁺, K⁺, and water molecules in those sites. We provide insights into the population of Mg²⁺ and K⁺ ions and water molecules in the outer sphere and how their fluctuations are intricately linked with the kinetics of the 58-mer. We find that many of the long-lived Mg²⁺ sites identified from the simulations agree with the locations of ions in the X-ray structure. We determine the excess ion atmosphere around the rRNA and compare it with experimental data. We investigate the collective behavior of RNA, ions, and water, by performing a joint principle component analysis for the Cartesian coordinates of the RNA phosphorus atoms and for the occupation counts of the association sites. Our results indicate that the 58-mer system is a complex polymer, composed of RNA that is encased by a fluctuating network of associated counterions, co-ions, and water.

Divalent ion competition reveals reorganization of an RNA ion atmosphere upon folding

Although RNA interactions with K+ and Mg2+ have been studied extensively, much less is known about the third most abundant cation in bacterial cells, putrescine2+, and how RNA folding might be influenced by the three ions in combination. In a new approach, we have observed the competition between Mg2+ and putrescine2+ (in a background of K+) with native, partially unfolded and highly extended conformations of an adenine riboswitch aptamer. With the native state, putrescine2+ is a weak competitor when the ratio of the excess Mg2+ (which neutralizes phosphate charge) to RNA is very low, but becomes much more effective at replacing Mg2+ as the excess Mg2+ in the RNA ion atmosphere increases. Putrescine2+ is even more effective in competing Mg2+ from the extended conformation, independent of the Mg2+ excess. To account for these and other results, we propose that both ions closely approach the surface of RNA secondary structure, but the completely folded RNA tertiary structure develops small pockets of very negative electrostatic potential that are more accessible to the compact charge of Mg2+. The sensitivity of RNA folding to the combination of Mg2+ and putrescine2+ found in vivo depends on the architectures of both the unfolded and native conformations.

All-Atom MD Predicts Magnesium-Induced Hairpin in Chemically Perturbed RNA Analog of F10 Therapeutic

Given their increasingly frequent usage, understanding the chemical and structural properties which allow therapeutic nucleic acids to promote the death of cancer cells is critical for medical advancement. One molecule of interest is a 10-mer of FdUMP (5-fluoro-2'-deoxyuridine-5'-O-monophosphate) also called F10. To investigate causes of structural stability, we have computationally restored the 2' oxygen on each ribose sugar of the phosphodiester backbone, creating FUMP[10]. Microsecond time-scale, all-atom, simulations of FUMP[10] in the presence of 150 mM MgCl₂ predict that the strand has a 45% probability of folding into a stable hairpin-like secondary structure. Analysis of 16 μs of data reveals phosphate interactions as likely contributors to the stability of this folded state. Comparison with polydT and polyU simulations predicts that FUMP[10]'s lowest order structures last for one to 2 orders of magnitude longer than similar nucleic acid strands. Here we provide a brief structural and conformational analysis of the predicted structures of FUMP[10], and suggest insights into its stability via comparison to F10, polydT, and polyU.

Computational Assessment of Potassium and Magnesium Ion Binding to a Buried Pocket in GTPase-Associating Center RNA

An experimentally well-studied model of RNA tertiary structures is a 58mer rRNA fragment, known as GTPase-associating center (GAC) RNA, in which a highly negative pocket walled by phosphate oxygen atoms is stabilized by a chelated cation. Although such deep pockets with more than one direct phosphate to ion chelation site normally include magnesium, as shown in one GAC crystal structure, another GAC crystal structure and solution experiments suggest potassium at this site. Both crystal structures also depict two magnesium ions directly bound to the phosphate groups comprising this controversial pocket. Here, we used classical molecular dynamics simulations as well as umbrella sampling to investigate the possibility of binding of potassium versus magnesium inside the pocket and to better characterize the chelation of one of the binding magnesium ions outside the pocket. The results support the preference of the pocket to accommodate potassium rather than magnesium and suggest that one of the closely binding magnesium ions can only bind at high magnesium concentrations, such as might be present during crystallization. This work illustrates the complementary utility of molecular modeling approaches with atomic-level detail in resolving discrepancies between conflicting experimental results.

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.

Measuring the differential stoichiometry and energetics of ligand binding to macromolecules by single-molecule force spectroscopy: an extended theory

Many chemical techniques exist for measuring the stoichiometry of ligand binding to a macromolecule; however, these techniques are often specific to certain ligands or require the presumption of specific binding models. Here, we further develop a previously reported, general, thermodynamic method for extracting the change in number of ligands bound to a macromolecule as that macromolecule undergoes a conformational transition driven by mechanical stretching, for example, by magnetic tweezers or optical trapping. We extend the theory of this method to consider systems with many ligands, experiments conducted in different thermodynamic ensembles (e.g., constant-force, constant-extension), and experiments in which the system is not at equilibrium. Further, we show that analysis of the same single-molecule mechanical manipulation data yields a measure of the differential free energy of stabilization due to ligand binding, that is, the free energy contribution by which ligand binding favors one conformation of the macromolecule over another. We interpret an existing data set measuring ion binding to RNA and DNA in terms of this free energy.

Comparison of interactions of diamine and Mg²⁺ with RNA tertiary structures: similar versus differential effects on the stabilities of diverse RNA folds

Cations play a large role in stabilizing the native state of RNA in vivo. In addition to Mg²⁺, putrescine²⁺ is an abundant divalent cation in bacterial cells, but its effect on the folding of RNA tertiary structure has not been widely explored. In this study, we look at how the stabilities of four structured RNAs, each with a different degree of dependence on K⁺ and Mg²⁺, are affected by putrescine²⁺ relative to Mg²⁺. Through the use of thermal melts, we observe that (i) at a given concentration, putrescine²⁺ is less effective than Mg²⁺ at stabilizing RNA, (ii) the stability imparted to RNA by various diamines is a function of charge density (average separation distance between charges) as well as the flexibility of the counterion, and (iii) when Mg²⁺ is already present in a buffer, further addition of putrescine²⁺ may either destabilize or stabilize RNA structure, depending on whether the native RNA does or does not chelate Mg²⁺ ion, respectively. At ion concentrations likely to be found in vivo, the effect of putrescine²⁺ on the free energy of folding of an RNA tertiary structure is probably quite small compared to that of Mg²⁺, but the ability of mixed Mg²⁺/putrescine²⁺ environments to (in effect) discriminate between different RNA architectures suggests that, in some cells, the evolution of functional RNA structures may have been influenced by the presence of putrescine²⁺.

Reverse engineering of monolayers and nanopatterns

In a molecularly decorated surface, the molecular tiles are “glued” to the surface by binding constants and possibly further “glued” to each other by cooperativity factors. At odds with mosaics, these “glues” come with the tiles and cannot be removed or supplemented. Binding polynomials quantify glue amounts from experimental data and may predict molecular self-organization on surfaces that can be exploited in organic (opto-)electronics.

Metal-ion rescue revisited: biochemical detection of site-bound metal ions important for RNA folding

Within the three-dimensional architectures of RNA molecules, divalent metal ions populate specific locations, shedding their water molecules to form chelates. These interactions help the RNA adopt and maintain specific conformations and frequently make essential contributions to function. Defining the locations of these site-bound metal ions remains challenging despite the growing database of RNA structures. Metal-ion rescue experiments have provided a powerful approach to identify and distinguish catalytic metal ions within RNA active sites, but the ability of such experiments to identify metal ions that contribute to tertiary structure acquisition and structural stability is less developed and has been challenged. Herein, we use the well-defined P4-P6 RNA domain of the Tetrahymena group I intron to reevaluate prior evidence against the discriminatory power of metal-ion rescue experiments and to advance thermodynamic descriptions necessary for interpreting these experiments. The approach successfully identifies ligands within the RNA that occupy the inner coordination sphere of divalent metal ions and distinguishes them from ligands that occupy the outer coordination sphere. Our results underscore the importance of obtaining complete folding isotherms and establishing and evaluating thermodynamic models in order to draw conclusions from metal-ion rescue experiments. These results establish metal-ion rescue as a rigorous tool for identifying and dissecting energetically important metal-ion interactions in RNAs that are noncatalytic but critical for RNA tertiary structure.

Characterization of metal ion-nucleic acid interactions in solution

Metal ions are inextricably involved with nucleic acids due to their polyanionic nature. In order to understand the structure and function of RNAs and DNAs, one needs to have detailed pictures on the structural, thermodynamic, and kinetic properties of metal ion interactions with these biomacromolecules. In this review we first compile the physicochemical properties of metal ions found and used in combination with nucleic acids in solution. The main part then describes the various methods developed over the past decades to investigate metal ion binding by nucleic acids in solution. This includes for example hydrolytic and radical cleavage experiments, mutational approaches, as well as kinetic isotope effects. In addition, spectroscopic techniques like EPR, lanthanide(III) luminescence, IR and Raman as well as various NMR methods are summarized. Aside from gaining knowledge about the thermodynamic properties on the metal ion-nucleic acid interactions, especially NMR can be used to extract information on the kinetics of ligand exchange rates of the metal ions applied. The final section deals with the influence of anions, buffers, and the solvent permittivity on the binding equilibria between metal ions and nucleic acids. Little is known on some of these aspects, but it is clear that these three factors have a large influence on the interaction between metal ions and nucleic acids.

Salt contribution to RNA tertiary structure folding stability

Accurate quantification of the ionic contribution to RNA folding stability could greatly enhance our ability to understand and predict RNA functions. Recently, motivated by the potential importance of ion correlation and fluctuation in RNA folding, we developed the tightly bound ion (TBI) model. Extensive experimental tests showed that the TBI model can lead to better treatment of multivalent ions than the Poisson-Boltzmann equation. In this study, we use the model to quantify the contribution of salt (Na(+) and Mg(2+)) to the RNA tertiary structure folding free energy. Folding of the RNA tertiary structure often involves intermediates. We focus on the folding transition from an intermediate state to the native state, and compute the electrostatic folding free energy of the RNA. Based on systematic calculations for a variety of RNA molecules, we derive a set of formulas for the electrostatic free energy for tertiary structural folding as a function of the sequence length and compactness of the RNA and the Na(+) and Mg(2+) concentrations. Extensive comparisons with experimental data suggest that our model and the extracted empirical formulas are quite reliable.

Evidence for a thermodynamically distinct Mg2+ ion associated with formation of an RNA tertiary structure

A folding strategy adopted by some RNAs is to chelate cations in pockets or cavities, where the ions neutralize charge from solvent-inaccessible phosphate. Although such buried Mg(2+)-RNA chelates could be responsible for a significant fraction of the Mg(2+)-dependent stabilization free energy of some RNA tertiary structures, direct measurements have not been feasible because of the difficulty of finding conditions under which the free energy of Mg(2+) chelation is uncoupled from RNA folding and from unfavorable interactions with Mg(2+) ions in other environments. In a 58mer rRNA fragment, we have used a high-affinity thermophilic ribosomal protein to trap the RNA in a structure nearly identical to native; Mg(2+)- and protein-stabilized structures differ in the solvent exposure of a single nucleotide located at the chelation site. Under these conditions, titration of a high affinity chelation site takes place in a micromolar range of Mg(2+) concentration, and is partially resolved from the accumulation of Mg(2+) in the ion atmosphere. From these experiments, we estimate the total and site-specific Mg(2+)-RNA interaction free energies over the range of accessed Mg(2+) concentrations. At 0.1 mM Mg(2+) and 60 mM K(+), specific site binding contributes ∼-3 kcal/mol of the total Mg(2+) interaction free energy of ∼-13 kcal/mol from all sources; at higher Mg(2+) concentrations the site-binding contribution becomes a smaller proportion of the total (-4.5 vs -33 kcal/mol). Under approximately physiological ionic conditions, the specific binding site will be saturated but will provide only a fraction of the total free energy of Mg(2+)-RNA interactions.

The osmolyte TMAO stabilizes native RNA tertiary structures in the absence of Mg2+: evidence for a large barrier to folding from phosphate dehydration

The stabilization of RNA tertiary structures by ions is well known, but the neutral osmolyte trimethylamine oxide (TMAO) can also effectively stabilize RNA tertiary structure. To begin to understand the physical basis for the effects of TMAO on RNA, we have quantitated the TMAO-induced stabilization of five RNAs with known structures. So-called m values, the increment in unfolding free energy per molal of osmolyte at constant KCl activity, are ∼0 for a hairpin secondary structure and between 0.70 and 1.85 kcal mol(-1)m(-1) for four RNA tertiary structures (30-86 nt). Further analysis of two RNAs by small-angle X-ray scattering and hydroxyl radical probing shows that TMAO reduces the radius of gyration of the unfolded ensemble to the same endpoint as seen in titration with Mg(2+) and that the structures stabilized by TMAO and Mg(2+) are indistinguishable. Remarkably, TMAO induces the native conformation of a Mg(2+) ion chelation site formed in part by a buried phosphate, even though Mg(2+) is absent. TMAO interacts weakly, if at all, with KCl, ruling out the possibility that TMAO stabilizes RNA indirectly by increasing salt activity. TMAO is, however, strongly excluded from the vicinity of dimethylphosphate (unfavorable interaction free energy, +211 cal mol(-1)m(-1) for the potassium salt), an ion that mimics the RNA backbone phosphate. We suggest that formation of RNA tertiary structure is accompanied by substantial phosphate dehydration (loss of 66-173 water molecules in the RNA structures studied) and that TMAO works principally by reducing the energetic penalty associated with this dehydration. The strong parallels we find between the effects of TMAO and Mg(2+) suggest that RNA sequence is more important than specific ion interactions in specifying the native structure.

Compact intermediates in RNA folding

Large noncoding RNAs fold into their biologically functional structures via compact yet disordered intermediates, which couple the stable secondary structure of the RNA with the emerging tertiary fold. The specificity of the collapse transition, which coincides with the assembly of helical domains, depends on RNA sequence and counterions. It determines the specificity of the folding pathways and the magnitude of the free energy barriers to the ensuing search for the native conformation. By coupling helix assembly with nascent tertiary interactions, compact folding intermediates in RNA also play a crucial role in ligand binding and RNA-protein recognition.

Dissecting electrostatic screening, specific ion binding, and ligand binding in an energetic model for glycine riboswitch folding

Riboswitches are gene-regulating RNAs that are usually found in the 5'-untranslated regions of messenger RNA. As the sugar-phosphate backbone of RNA is highly negatively charged, the folding and ligand-binding interactions of riboswitches are strongly dependent on the presence of cations. Using small angle X-ray scattering (SAXS) and hydroxyl radical footprinting, we examined the cation dependence of the different folding stages of the glycine-binding riboswitch from Vibrio cholerae. We found that the partial folding of the tandem aptamer of this riboswitch in the absence of glycine is supported by all tested mono- and divalent ions, suggesting that this transition is mediated by nonspecific electrostatic screening. Poisson-Boltzmann calculations using SAXS-derived low-resolution structural models allowed us to perform an energetic dissection of this process. The results showed that a model with a constant favorable contribution to folding that is opposed by an unfavorable electrostatic term that varies with ion concentration and valency provides a reasonable quantitative description of the observed folding behavior. Glycine binding, on the other hand, requires specific divalent ions binding based on the observation that Mg(2+), Ca(2+), and Mn(2+) facilitated glycine binding, whereas other divalent cations did not. The results provide a case study of how ion-dependent electrostatic relaxation, specific ion binding, and ligand binding can be coupled to shape the energetic landscape of a riboswitch and can begin to be quantitatively dissected.

Dependence of RNA tertiary structural stability on Mg2+ concentration: interpretation of the Hill equation and coefficient

The Mg(2+)-induced folding of RNA tertiary structures is readily observed via titrations of RNA with MgCl(2). Such titrations are commonly analyzed using a site binding formalism that includes a parameter, the Hill coefficient n, which is sometimes deemed the number of Mg(2+) ions bound by the native RNA at specific sites. However, the long-range nature of electrostatic interactions allows ions some distance from the RNA to stabilize an RNA structure. A complete description of all interactions taking place between Mg(2+) and an RNA uses a preferential interaction coefficient, Gamma(2+), which represents the "excess" Mg(2+) neutralizing the RNA charge. The difference between Gamma(2+) for the native and unfolded RNA forms (DeltaGamma(2+)) is the number of Mg(2+) ions "taken up" by an RNA upon folding. Here we determine the conditions under which the Hill coefficient n can be equated to the ion uptake DeltaGamma(2+) and find that two approximations are necessary: (i) the Mg(2+) activity coefficient is independent of concentration during a titration, and (ii) the dependence of DeltaGamma(2+) on Mg(2+) concentration is weak. Titration experiments with a Mg(2+)-binding dye and an adenine-binding riboswitch were designed to test these approximations. Inclusion of a 30-fold excess of KCl over MgCl(2) was sufficient to maintain a constant Mg(2+) activity coefficient. We also observed that Mg(2+) uptake by the RNA varied from near zero to approximately 2.6 as the Mg(2+) concentration increases over an approximately 100-fold range. It is possible to determine DeltaGamma(2+) from Mg(2+)-RNA titrations, but the values are only applicable to a limited range of solution conditions.

Ion-RNA interactions thermodynamic analysis of the effects of mono- and divalent ions on RNA conformational equilibria

RNA secondary and tertiary structures are strongly stabilized by added salts, and a quantitative thermodynamic analysis of the relevant ion-RNA interactions is an important aspect of the RNA folding problem. Because of long-range electrostatic forces, an RNA perturbs the distribution of both cations and anions throughout a large volume. Binding formalisms that require a distinction between "bound" and "free" ions become problematic in such situations. A more fundamental thermodynamic framework is developed here, based on preferential interaction coefficients; linkage equations derived from this framework provide a model-free description of the "uptake" or "release" of cations and anions that accompany an RNA conformational transition. Formulas appropriate for analyzing the dependence of RNA stability on either mono- or divalent salt concentration are presented and their application to experimental data is illustrated. Two example datasets are analyzed with respect to the monovalent salt dependence of tertiary structure formation in different RNAs, and three different experimental methods for quantitating the "uptake" of Mg(2+) ions are applied to the folding of a riboswitch RNA. Advantages and limitations of each method are discussed.

Metal ion dependence of cooperative collapse transitions in RNA

Positively charged counterions drive RNA molecules into compact configurations that lead to their biologically active structures. To understand how the valence and size of the cations influences the collapse transition in RNA, small-angle X-ray scattering was used to follow the decrease in the radius of gyration (R(g)) of the Azoarcus and Tetrahymena ribozymes in different cations. Small, multivalent cations induced the collapse of both ribozymes more efficiently than did monovalent ions. Thus, the cooperativity of the collapse transition depends on the counterion charge density. Singular value decomposition of the scattering curves showed that folding of the smaller and more thermostable Azoarcus ribozyme is well described by two components, whereas collapse of the larger Tetrahymena ribozyme involves at least one intermediate. The ion-dependent persistence length, extracted from the distance distribution of the scattering vectors, shows that the Azoarcus ribozyme is less flexible at the midpoint of transition in low-charge-density ions than in high-charge-density ions. We conclude that the formation of sequence-specific tertiary interactions in the Azoarcus ribozyme overlaps with neutralization of the phosphate charge, while tertiary folding of the Tetrahymena ribozyme requires additional counterions. Thus, the stability of the RNA structure determines its sensitivity to the valence and size of the counterions.

Identification of catalytic metal ion ligands in ribozymes

Site-bound metal ions participate in the catalytic mechanisms of many ribozymes. Understanding these mechanisms therefore requires knowledge of the specific ligands on both substrate and ribozyme that coordinate these catalytic metal ions. A number of different structural and biochemical strategies have been developed and refined for identifying metal ion binding sites within ribozymes, and for assessing the catalytic contributions of the metal ions bound at those sites. We review these approaches and provide examples of their application, focusing in particular on metal ion rescue experiments and their roles in the construction of the transition state models for the Tetrahymena group I and RNase P ribozymes.

The influence of monovalent cation size on the stability of RNA tertiary structures

Many RNA tertiary structures are stable in the presence of monovalent ions alone. To evaluate the degree to which ions at or near the surfaces of such RNAs contribute to stability, the salt-dependent stability of a variety of RNA structures was measured with each of the five group I cations. The stability of hairpin secondary structures and a pseudoknot tertiary structure are insensitive to the ion identity, but the tertiary structures of two other RNAs, an adenine riboswitch and a kissing loop complex, become more stable by 2-3 kcal/mol as ion size decreases. This "default" trend is attributed to the ability of smaller ions to approach the RNA surface more closely. The degree of cation accumulation around the kissing loop complex was also inversely proportional to ion radius, perhaps because of the presence of sterically restricted pockets that can be accessed only by smaller ions. An RNA containing the tetraloop-receptor motif shows a strong (up to approximately 3 kcal/mol) preference for Na(+) or K(+) over other group I ions, consistent with the chelation of K(+) by this motif in some crystal structures. This RNA reverts to the default dependence on ion size when a base forming part of the chelation site is mutated. Lastly, an RNA aptamer for cobinamide, which was originally selected in the presence of high concentrations of LiCl, binds ligand more strongly in the presence of Li(+) than other monovalent ions. On the basis of these trends in RNA stability with group I ion size, it is argued that two features of RNA tertiary structures may promote strong interactions with ions at or near the RNA surface: negative charge densities that are higher than that in secondary structures, and the occasional presence of chelation sites, which are electronegative pockets that selectively bind ions of an optimum size.

Molecular simulation studies of monovalent counterion-mediated interactions in a model RNA kissing loop

A kissing loop is a highly stable complex formed by loop-loop base-pairing between two RNA hairpins. This common structural motif is utilized in a wide variety of RNA-mediated processes, including antisense recognition, substrate recognition in riboswitches, and viral replication. Recent work has shown that the Tar-Tar(*) complex, an archetypal kissing loop, can form without Mg(2+), so long as high concentrations of alkali chloride salts are present. Interestingly, the stability of the complex is found to decrease with increasing cation size. In this work, we used molecular simulations to develop a molecular-level understanding of the origins of the observed counterion specificity. The ionic atmosphere of the Tar-Tar(*) complex was examined in the presence of 800 mm (where m denotes molality) NaCl, KCl, or CsCl. We used spatial free-energy density profiles to analyze differences in counterion accumulation at different spatial extents from the RNA molecule. We found that the lowest free-energy levels, which are situated in the vicinity of the loop-loop interface, can accommodate roughly two counterions, irrespective of counterion type. However, as we moved into higher free-energy levels, which are farther from the loop-loop interface, we observed increased differences in the numbers of accumulated counterions, with this number being largest for Na(+) and smallest for Cs(+). We analyzed the source of these differences and were able to attribute these to two distinct features: The extent of partial dehydration varies based on cation type; the smaller the cation, the greater the degree of dehydration. While smaller ions bind their first-hydration-shell water molecules more tightly than larger ions, they are also able to shed these water molecules for stronger electrostatic interactions with the RNA molecule. Secondly, we observed a distinct asymmetry in the numbers of accumulated cations around each hairpin in the Tar-Tar(*) complex. We were able to ascribe this asymmetry to the presence of a guanine tract in the Tar hairpin, which facilitates partial dehydration of the counterions. However, the smaller ions compensate for this asymmetry by forming a belt around the loop-loop interface in intermediate free-energy levels. As a result, the degree of asymmetry in counterion accumulation around individual hairpins shows an inverse correlation with the experimentally observed cation specificity for the stability of Tar-Tar(*) (i.e., the smaller the asymmetry, the greater the experimentally observed stability). This in turn provides a plausible explanation for why the smaller cations help stabilize the Tar-Tar(*) complex better than the larger cations. These findings suggest that the specific sequence and structural features of the Tar-Tar(*) complex may be the source of the experimentally observed cation specificity in Tar-Tar(*) stability. Our results lead to testable predictions for how changes in sequence might alter the observed counterion specificity in kissing loop stability.

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.

Like polarity ion/ion reactions enable the investigation of specific metal interactions in nucleic acids and their noncovalent assemblies

A rare example of ion/ion reaction between species of like polarity was shown to take place during the transfer of metal cations from nucleic acid substrates to chelating agents in the gas phase. Gaseous anionic reactants were generated from separate solutions of analyte and chelator by using a dual nanospray setup. The respective multiply charged ions shared the same path and were allowed to react for a predetermined interval in an rf-only hexapole before high-resolution analysis by Fourier transform ion cyclotron resonance (FTICR) mass spectrometry. Efficient transfer of sodium and magnesium ions was readily observed with significant reduction of the nonspecific adducts that are typically associated with decreased sensitivity and resolution in the analysis of nucleic acid samples. Metal cations were abstracted from the initial analyte without being replaced by protons, in a process that was clearly dependent on the concentration of chelator in the auxiliary emitter and on the time spent by the reactants in the hexapole element. A survey of the properties of selected anionic chelators showed that their known affinity for a target cation in solution was more critical than their maximum anionic charge in determining the outcome of the transfer process. The analysis of selected assemblies requiring divalent cations to preserve their structural integrity and functional properties demonstrated that ion/ion reactions were clearly capable of discriminating between nonspecific interactions and specific coordination based on transfer susceptibility. These examples demonstrated that the ability to selectively eliminate nonspecific adducts in the gas phase, after the desolvation process is complete, offers a unique opportunity for studying specific metal binding in biological systems without resorting to separation procedures that may adversely affect the position of binding equilibria in solution and disrupt the assemblies under investigation.

Metal ion specificities for folding and cleavage activity in the Schistosoma hammerhead ribozyme

The effects of various metal ions on cleavage activity and global folding have been studied in the extended Schistosoma hammerhead ribozyme. Fluorescence resonance energy transfer was used to probe global folding as a function of various monovalent and divalent metal ions in this ribozyme. The divalent metals ions Ca(2+), Mg(2+), Mn(2+), and Sr(2+) have a relatively small variation (less than sixfold) in their ability to globally fold the hammerhead ribozyme, which contrasts with the very large difference (>10,000-fold) in apparent rate constants for cleavage for these divalent metal ions in single-turnover kinetic experiments. There is still a very large range (>4600-fold) in the apparent rate constants for cleavage for these divalent metal ions measured in high salt (2 M NaCl) conditions where the ribozyme is globally folded. These results demonstrate that the identity of the divalent metal ion has little effect on global folding of the Schistosoma hammerhead ribozyme, whereas it has a very large effect on the cleavage kinetics. Mechanisms by which the identity of the divalent metal ion can have such a large effect on cleavage activity in the Schistosoma hammerhead ribozyme are discussed.

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.

Low specificity of metal ion binding in the metal ion core of a folded RNA

The structure and activity of nucleic acids depend on their interactions with metal ions. Fundamental to these interactions is the degree of specificity observed between the metal ions and nucleic acids, and a complete description of nucleic acid folding requires that we understand the nature of the interactions with metal ions, including specificity. The prior demonstration that high concentrations of monovalent cations prevent nonspecific association of divalent ions with nucleic acids provides a novel and powerful means to examine site-specific metal ion binding isolated from complicating effects of the ion atmosphere. Using these high monovalent cation solution conditions we have monitored the affinity of a series of divalent metal ions for two site-specific metal ion binding sites in the P4-P6 domain of the Tetrahymena group I intron ribozyme. The metal ion core of this highly structured RNA binds two divalent metal ions under these conditions. Despite multiple metal ion-RNA interactions observed in the X-ray crystallographic structure of P4-P6 RNA at the metal ion binding sites, these sites exhibit low specificity among Mn(2+), Mg(2+), Ca(2+), Ni(2+), and Zn(2+). Nevertheless, the largest divalent metal ions tested, Sr(2+) and Ba(2+), were excluded from binding, exhibiting affinities at least two orders of magnitude weaker than observed for the other metal ions. Thus, a picture emerges of two metal ion binding sites, each with a high tolerance for metal ions with different properties but also with limits to accommodation. These limits presumably arise from steric or electrostatic features of the metal ion binding sites.

Ion atmosphere around nucleic acid

The poyion-ion preferential interaction coefficient Gamma describes the exclusion of coions and accumulations of counterions in the vicinity of a polyion in an aqueous solution. We give tight upper and lower bounds for Gamma when the polyion can be modeled by a cylinder of infinite length but of arbitrary charge density. This case can be used as a model for long strands of DNA or RNA in an aqueous solution containing univalent cations. The salt dependence of Gamma is predicted from low to intermediate and high salt concentrations. We also indicate how the bounds for the infinite polyion can be exploited to place bounds for polyions of length greater than a constant on the order of the inverse Debye screening length.

Gas-Phase Reactions Between Thiourea and Ca2+: New Evidence for the Formation of [Ca(NH3)]2+ and Other Doubly Charged Species

The gas-phase reactions between Ca(2+) and thiourea are investigated by means of electrospray ionization/mass spectrometry experiments. The MS/MS spectra of [Ca(thiourea)](2+) complexes show the appearance of new doubly charged species formed by the loss of NH(3) and HNCS. Other intense peaks at m/z 43, 56, 60, 73, 76 and 98 are also observed, and assigned to monocations produced in different coulomb-explosion processes. The structures and bonding characteristics of the different stationary points of the [Ca(thiourea)](2+) potential energy surface (PES) were theoretically studied by DFT calculations carried out at B3LYP/cc-pWCVTZ level. The analysis of the topology of this PES permits to propose different mechanisms for the loss of ammonia and HNCS, and to identify, the m/z 43, 56, 60, 73, 76 and 98 peaks as H(2)NCNH(+), CaNH(2) (+), H(2)NCS(+), CaSH(+), thiourea(+) and CaNCS(+) ions respectively. There are significant dissimilarities between the reactivity of urea and thiourea, which are related to the lower ionization energy of the latter, and to the fact that thioenols are intrinsically more stable than enols with respect to the corresponding keto forms.

Analysis of biological and synthetic ribonucleic acids by liquid chromatography-mass spectrometry using monolithic capillary columns

Ion-pair reversed-phase high-performance liquid chromatography (IP-RP-HPLC) has been evaluated as a method for the fractionation and desalting of ribonucleic acids prior to their characterization by electrospray ionization mass spectrometry. Monolithic, poly(styrene-divinylbenzene)-based capillary columns allowed the rapid and highly efficient fractionation of both synthetic and biological ribonucleic acids. The common problem of gas-phase cation adduction that is particularly prevalent in the mass spectrometric analysis of ribonucleic acids was tackled through a combination of chromatographic purification and the addition of ethylenediaminetetraacetic acid to the sample at a concentration of 25 mmol/L shortly before on-line analysis. For RNA molecules ranging in size from 10 to 120 nucleotides, the mass accuracies were typically better than 0.02%, which allowed the characterization and identification of failure sequences and byproducts with high confidence. Following injection of a 500 nL sample onto a 60 x 0.2 mm column, the limit of detection for a 120-nucleotide ribosomal RNA transcript from Escherichia coli was in the 50-80 fmol range. The method was applied to the analysis of synthetic oligoribonucleotides, transfer RNAs, and ribosomal RNA. Finally, sequence information was derived for low picomole amounts of a 32-mer RNA upon chromatographic purification and tandem mass spectrometric investigation in an ion trap mass spectrometer. Complete series of fragment ions of the c- and y-types could be assigned in the tandem mass spectrum. In conclusion, IP-RP-HPLC using monolithic capillary columns represents a very useful tool for the structural investigation and quantitative determination of RNAs of synthetic and biological origin.

L11 domain rearrangement upon binding to RNA and thiostrepton studied by NMR spectroscopy

Ribosomal proteins are assumed to stabilize specific RNA structures and promote compact folding of the large rRNA. The conformational dynamics of the protein between the bound and unbound state play an important role in the binding process. We have studied those dynamical changes in detail for the highly conserved complex between the ribosomal protein L11 and the GTPase region of 23S rRNA. The RNA domain is compactly folded into a well defined tertiary structure, which is further stabilized by the association with the C-terminal domain of the L11 protein (L11_ctd). In addition, the N-terminal domain of L11 (L11_ntd) is implicated in the binding of the natural thiazole antibiotic thiostrepton, which disrupts the elongation factor function. We have studied the conformation of the ribosomal protein and its dynamics by NMR in the unbound state, the RNA bound state and in the ternary complex with the RNA and thiostrepton. Our data reveal a rearrangement of the L11_ntd, placing it closer to the RNA after binding of thiostrepton, which may prevent binding of elongation factors. We propose a model for the ternary L11–RNA–thiostrepton complex that is additionally based on interaction data and conformational information of the L11 protein. The model is consistent with earlier findings and provides an explanation for the role of L11_ntd in elongation factor binding.

Mg2+-RNA interaction free energies and their relationship to the folding of RNA tertiary structures

Mg2+ ions are very effective at stabilizing tertiary structures in RNAs. In most cases, folding of an RNA is so strongly coupled to its interactions with Mg2+ that it is difficult to separate free energies of Mg2+-RNA interactions from the intrinsic free energy of RNA folding. To devise quantitative models accounting for this phenomenon of Mg2+-induced RNA folding, it is necessary to independently determine Mg2+-RNA interaction free energies for folded and unfolded RNA forms. In this work, the energetics of Mg2+-RNA interactions are derived from an assay that measures the effective concentration of Mg2+ in the presence of RNA. These measurements are used with other measures of RNA stability to develop an overall picture of the energetics of Mg2+-induced RNA folding. Two different RNAs are discussed, a pseudoknot and an rRNA fragment. Both RNAs interact strongly with Mg2+ when partially unfolded, but the two folded RNAs differ dramatically in their inherent stability in the absence of Mg2+ and in the free energy of their interactions with Mg2+. From these results, it appears that any comprehensive framework for understanding Mg2+-induced stabilization of RNA will have to (i) take into account the interactions of ions with the partially unfolded RNAs and (ii) identify factors responsible for the widely different strengths with which folded tertiary structures interact with Mg2+.

Polymorphic Behavior in Protein−Surfactant Mixtures: The Water−Bovine Serum Albumin−Sodium Taurodeoxycholate System

Mixtures containing water, bovine serum albumin (BSA), and sodium taurodeoxycholate (NaTDC), a component of the bile in mammals, have been investigated in a wide range of composition and pH. Depending on the concentration of both solutes and the pH, solutions, precipitates, and gels are formed. Under spontaneous pH conditions, the transport properties in dilute solutions indicate the occurrence of significant interactions between BSA and the surfactant. Conversely, acidic media favor the formation of nonsoluble protein-surfactant complexes, with subsequent precipitation. The nucleation kinetics of the protein-surfactant complexes in solid form and the related precipitation processes can be slow or fast, depending on the overall solute content and the mole ratio. At high concentrations, a gel, extending on both sides of the charge neutralization line, and two-phase regions are observed. Gels shrink in open air and swell in the presence of excess water. Depending on concentration and temperature, the gels transform from an essentially liquidlike behavior to that peculiar to true gels (when G' > or = G''). The thermal gelation threshold, the temperature above which G' > or = G'', depends on BSA and NaTDC content and is concomitant to moderate heat effects, inferred by differential scanning calorimetry (DSC). The above data also indicate that the protein thermal denaturation in the gel is shifted to higher temperatures compared to water. Such a stabilizing effect is presumably related to the occurrence of both electrostatic and hydrophobic interactions with NaTDC. Water self-diffusion in the gels is slightly slower than that in the bulk and poorly sensitive to composition: it is about 65% the value of neat H2O in a wide concentration range, irrespective of the BSA, or NaTDC, concentration. A peculiar behavior is also observed in 23Na longitudinal and transverse relaxation rates. The T1 and T2 values, measured at 105.75 MHz on BSA-NaTDC gels, indicate that the motions determining the NMR relaxation of the sodium ions in the hydration layer of the protein-surfactant aggregates are not slow, having frequencies comparable with the Larmor one. The above properties, especially the rheological and the spectroscopic ones, are important for understanding the behavior of gels based on protein-surfactant mixtures.

Determining the Mg2+ stoichiometry for folding an RNA metal ion core

The folding and catalytic function of RNA molecules depend on their interactions with divalent metal ions, such as magnesium. As with every molecular process, the most basic knowledge required for understanding the close relationship of an RNA with its metal ions is the stoichiometry of the interaction. Unfortunately, inventories of the numbers of divalent ions associated with unfolded and folded RNA states have been unattainable. A common approach has been to interpret Hill coefficients fit to folding equilibria as the number of metal ions bound upon folding. However, this approach is vitiated by the presence of diffusely associated divalent ions in a dynamic ion atmosphere and by the likelihood of multiple transitions along a folding pathway. We demonstrate that the use of molar concentrations of background monovalent salt can alleviate these complications. These simplifying solution conditions allow a precise determination of the stoichiometry of the magnesium ions involved in folding the metal ion core of the P4-P6 domain of the Tetrahymena group I ribozyme. Hill analysis of hydroxyl radical footprinting data suggests that the P4-P6 RNA core folds cooperatively upon the association of two metal ions. This unexpectedly small stoichiometry is strongly supported by counting magnesium ions associated with the P4-P6 RNA via fluorescence titration and atomic emission spectroscopy. By pinpointing the metal ion stoichiometry, these measurements provide a critical but previously missing step in the thermodynamic dissection of the coupling between metal ion binding and RNA folding.

Ions and RNA Folding

The problem of how ions influence the folding of RNA into specific tertiary structures is being addressed from both thermodynamic (by how much do different salts affect the free energy change of folding) and structural (how are ions arranged on or near an RNA and what kinds of environments do they occupy) points of view. The challenge is to link these different approaches in a theoretical framework that relates the energetics of ion-RNA interactions to the spatial distribution of ions. This review distinguishes three different kinds of ion environments that differ in the extent of direct ion-RNA contacts and the degree to which the ion hydration is perturbed, and summarizes the current understanding of the way each environment relates to the overall energetics of RNA folding.

Linkage between proton binding and folding in RNA: a thermodynamic framework and its experimental application for investigating pKa shifting

Perturbation of pKa values can change the favored protonation states of the nucleobases at biological pH and thereby modulate the function of RNA and DNA molecules. In an effort to understand the driving forces for pKa shifting specific to nucleic acids, we developed a thermodynamic framework that relates proton binding to the nucleobases and the helix-coil transition. Key features that emerge from the treatment are a comprehensive description of all the actions of proton binding on RNA folding: acid and alkaline denaturation of the helix and pKa shifting in the folded state. Practical experimental approaches for measuring pKas from thermal denaturation experiments are developed. Microscopic pka values (where ka is the acid dissociation constant) for the unfolded state were determined directly by experiments on unstructured oligonucleotides, which led to a macroscopic pKa for the ensemble of unfolded states shifted toward neutrality. The formalism was then applied to pH-dependent UV melting data for model DNA oligonucleotides. Folded-state pka) values were in good agreement with the outcome of pH titrations, and the acid and alkaline denaturation regions were well described. The formalism developed here is similar to that of Draper and coworkers for Mg2+ binding to RNA, except that the unfolded state is described explicitly owing to the presence of specific proton-binding sites on the bases. A principal conclusion is that it should be possible to attain large pKa shifts by designing RNA molecules that fold cooperatively.

Precipitation of RNA impurities with high salt in a plasmid DNA purification process: use of experimental design to determine reaction conditions

The use of high salt solution to precipitate RNA in a pharmaceutical-grade plasmid DNA purification process was investigated. Five antichaotropic salts were tested for their potential to precipitate RNA. Calcium chloride was by far the best precipitant with high RNA removal in a very short incubation time. Calcium chloride precipitation conditions were investigated at two stages of a plasmid purification process using experimental design techniques. The effect of up to five factors on RNA precipitation and plasmid recovery was assessed by statistical modeling. Optimized conditions for calcium chloride precipitation were then introduced to the plasmid purification process resulting in the efficient removal of most impurities (RNA, chromosomal DNA, proteins, and endotoxins).

Specific binding of ribosome recycling factor (RRF) with the Escherichia coli ribosomes by BIACORE

The direct assays on Biacore with immobilised RRF and purified L11 from E. coli in the flow trough have shown unspecific binding between the both proteins. The interaction of RRF with GTPase domain of E. coli ribosomes, a functionally active complex of L11 with 23S r RNA and L10.(L7/L12)4 was studied by Biacore. In the experiments of binding of RRF with 30S, 50S and 70S ribosomes from E. coli were used the antibiotics thiostrepton, tetracycline and neomycin and factors, influencing the 70S dissociation Mg2+, NH4Cl, EDTA. The binding is strongly dependent from the concentrations of RRF, Mg2+, NH4Cl, EDTA and is inhibited by thiostrepton. The effect is most specific for 50S subunits and indicates that the GTPase centre can be considered as a possible site of interaction of RRF with the ribosome. We can consider an electrostatic character of the interactions with most probable candidate 16S and 23S r RNA at the interface of 30S and 50S ribosomal subunits.

Interactions of the non-coding RNA DsrA and RpoS mRNA with the 30 S ribosomal subunit

Expression of sigma(s), the gene product of rpoS, is controlled translationally in response to many environmental stresses. DsrA, a small 87-nucleotide non-coding RNA molecule, acts to increase translational efficiency of RpoS mRNA under some growth conditions. In this work, we demonstrate that DsrA binds directly to the 30 S ribosomal subunit with an observed equilibrium affinity of 2.8 x 10(7) m(-1). DsrA does not compete with RpoS mRNA or tRNA(f)(Met) for binding to the 30 S subunit. The 5' end of DsrA binds to 30 S subunits with an observed equilibrium association constant of 2.0 x 10(6) m(-1), indicating that the full affinity of the interaction requires the entire DsrA sequence. In order to investigate translational efficiency of RpoS mRNA, we examined both ribosome-binding site accessibility and the binding of RpoS mRNA to 30 S ribosomal subunits. We find that that ribosome-binding site accessibility is modulated as a function of divalent cation concentration during mRNA renaturation and by the presence of an antisense sequence that binds to nucleotides 1-16 of the RpoS mRNA fragment. The ribosome-binding site accessibility correlates with the amount of RpoS mRNA participating in 30 S-mRNA "pre-initiation" translational complex formation and provides evidence that regulation follows a competitive model of regulation.

A thermodynamic framework for the magnesium-dependent folding of RNA

The goal of this review is to present a unified picture of the relationship between ion binding and RNA folding based on recent theoretical and computational advances. In particular, we present a model describing how the association of magnesium ions is coupled to the tertiary structure folding of several well-characterized RNA molecules. This model is developed in terms of the nonlinear Poisson-Boltzmann (NLPB) equation, which provides a rigorous electrostatic description of the interaction between Mg(2+) and specific RNA structures. In our description, most of the ions surrounding an RNA behave as a thermally fluctuating ensemble distributed according to a Boltzmann weighted average of the mean electrostatic potential around the RNA. In some cases, however, individual ions near the RNA may shed some of their surrounding waters to optimize their Coulombic interactions with the negatively charged ligands on the RNA. These chelated ions are energetically distinct from the surrounding ensemble and must be treated explicitly. This model is used to explore several different RNA systems that interact differently with Mg(2+). In each case, the NLPB equation accurately describes the stoichiometric and energetic linkage between Mg(2+) binding and RNA folding without requiring any fitted parameters in the calculation. Based on this model, we present a physical description of how Mg(2+) binds and stabilizes specific RNA structures to promote the folding reaction.

Mechanistic characterization of the HDV genomic ribozyme: classifying the catalytic and structural metal ion sites within a multichannel reaction mechanism

Prior studies of the metal ion dependence of the self-cleavage reaction of the HDV genomic ribozyme led to a mechanistic framework in which the ribozyme can self-cleave by multiple Mg2+ ion-independent and -dependent channels [Nakano et al. (2001) Biochemistry 40, 12022]. In particular, channel 2 involves cleavage in the presence of a structural Mg2+ ion without participation of a catalytic divalent metal ion, while channel 3 involves both structural and catalytic Mg2+ ions. In the present study, experiments were performed to probe the nature of the various divalent ion sites and any specificity for Mg2+. A series of alkaline earth metal ions was tested for the ability to catalyze self-cleavage of the ribozyme under conditions that favor either channel 2 or channel 3. Under conditions that populate primarily channel 3, nearly identical K(d)s were obtained for Mg2+, Ca2+, Ba2+, and Sr2+, with a slight discrimination against Ca2+. In contrast, under conditions that populate primarily channel 2, tighter binding was observed as ion size decreases. Moreover, [Co(NH3)6]3+ was found to be a strong competitive inhibitor of Mg2+ for channel 3 but not for channel 2. The thermal unfolding of the cleaved ribozyme was also examined, and two transitions were found. Urea-dependent studies gave m-values that allowed the lower temperature transition to be assigned to tertiary structure unfolding. The effects of high concentrations of Na+ on the melting temperature for RNA unfolding and the reaction rate revealed ion binding to the folded RNA, with significant competition of Na+ (Hill coefficient of 1.5-1.7) for a structural Mg2+ ion and an unusually high intrinsic affinity of the structural ion for the RNA. Taken together, these data support the existence of two different classes of metal ion sites on the ribozyme: a structural site that is inner sphere with a major electrostatic component and a preference for Mg2+, and a weak catalytic site that is outer sphere with little preference for a particular divalent ion.

Conservation of helical structure contributes to functional metal ion interactions in the catalytic domain of ribonuclease P RNA

Like protein enzymes, catalytic RNAs contain conserved structure motifs important for function. A universal feature of the catalytic domain of ribonuclease P RNA is a bulged-helix motif within the P1-P4 helix junction. Here, we show that changes in bulged nucleotide identity and position within helix P4 affect both catalysis and substrate binding, while a subset of the mutations resulted only in catalytic defects. We find that the proximity of the bulge to sites of metal ion coordination in P4 is important for catalysis; moving the bulge distal to these sites and deleting it had similarly large effects, while moving it proximal to these sites had only a moderate effect on catalysis. To test whether the effects of the mutations are linked to metal ion interactions, we used terbium-dependent cleavage of the phosphate backbone to probe metal ion-binding sites in the wild-type and mutant ribozymes. We detect cleavages at specific sites within the catalytic domain, including helix P4 and J3/4, which have previously been shown to participate directly in metal ion interactions. Mutations introduced into P4 cause local changes in the terbium cleavage pattern due to alternate metal ion-binding configurations with the helix. In addition, a bulge deletion mutation results in a 100-fold decrease in the single turnover cleavage rate constant at saturating magnesium levels, and a reduced affinity for magnesium ions important for catalysis. In light of the alternate terbium cleavage pattern in P4 caused by bulge deletion, this decreased ability to utilize magnesium ions for catalysis appears to be due to localized structural changes in the ribozyme's catalytic core that weaken metal ion interactions in P4 and J3/4. The information reported here, therefore, provides evidence that the universal conservation of the P4 structure is based in part on optimization of metal ion interactions important for catalysis.