Mg2+-dependent compaction and folding of yeast tRNAPhe and the catalytic domain of the B. subtilis RNase P RNA determined by small-angle X-ray scattering

Pressure pushes tRNALys3 into excited conformational states

Conformational dynamics play essential roles in RNA function. However, detailed structural characterization of excited states of RNA remains challenging. Here, we apply high hydrostatic pressure (HP) to populate excited conformational states of tRNALys3, and structurally characterize them using a combination of HP 2D-NMR, HP-SAXS (HP-small-angle X-ray scattering), and computational modeling. HP-NMR revealed that pressure disrupts the interactions of the imino protons of the uridine and guanosine U-A and G-C base pairs of tRNALys3. HP-SAXS profiles showed a change in shape, but no change in overall extension of the transfer RNA (tRNA) at HP. Configurations extracted from computational ensemble modeling of HP-SAXS profiles were consistent with the NMR results, exhibiting significant disruptions to the acceptor stem, the anticodon stem, and the D-stem regions at HP. We propose that initiation of reverse transcription of HIV RNA could make use of one or more of these excited states.

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.

GapR binds DNA through dynamic opening of its tetrameric interface

GapR is a nucleoid-associated protein that is an essential regulator of chromosome replication in the cell cycle model Caulobacter crescentus. Here, we demonstrate that free GapR is a homotetramer, but not a dimer as previously reported (Guo et al., Cell 175: 583–597, 2018). We have determined the crystal structure of GapR in complex with a 10-bp A-tract DNA, which has an open tetrameric conformation, different from the closed clamp conformation in the previously reported crystal structure of GapR/DNA complex. The free GapR adopts multiple conformations in dynamic exchange equilibrium, with the major conformation resembling the closed tetrameric conformation, while the open tetrameric conformation is a representative of minor conformers. As it is impossible for the circular genomic DNA to get into the central DNA binding tunnel of the major conformation, we propose that GapR initially binds DNA through the open conformation, and then undergoes structural rearrangement to form the closed conformation which fully encircles the DNA. GapR prefers to bind DNA with 10-bp consecutive A/T base pairs nonselectively (K_d ∼12 nM), while it can also bind GC-rich DNA sequence with a reasonable affinity of about 120 nM. Besides, our results suggest that GapR binding results in widening the minor groove of AT-rich DNA, instead of overtwisting DNA.

Mg2+ Sensing by an RNA Fragment: Role of Mg2+-Coordinated Water Molecules

RNA molecules selectively bind to specific metal ions to populate their functional active states, making it important to understand their source of ion selectivity. In large RNA systems, metal ions interact with the RNA at multiple locations, making it difficult to decipher the precise role of ions in folding. To overcome this complexity, we studied the role of different metal ions (Mg²⁺, Ca²⁺, and K⁺) in the folding of a small RNA hairpin motif (5'-ucCAAAga-3') using unbiased all-atom molecular dynamics simulations. The advantage of studying this system is that it requires specific binding of a single metal ion to fold to its native state. We find that even for this small RNA, the folding free energy surface (FES) is multidimensional as different metal ions present in the solution can simultaneously facilitate folding. The FES shows that specific binding of a metal ion is indispensable for its folding. We further show that in addition to the negatively charged phosphate groups, the spatial organization of electronegative nucleobase atoms drives the site-specific binding of the metal ions. Even though the binding site cannot discriminate between different metal ions, RNA folds efficiently only in a Mg²⁺ solution. We show that the rigid network of Mg²⁺-coordinated water molecules facilitates the formation of important interactions in the transition state. The other metal ions such as K⁺ and Ca²⁺ cannot facilitate the formation of such interactions. These results allow us to hypothesize possible metal-sensing mechanisms in large metalloriboswitches and also provide useful insights into the design of appropriate collective variables for studying large RNA molecules using enhanced sampling methods.

Small-Angle Scattering as a Structural Probe for Nucleic Acid Nanoparticles (NANPs) in a Dynamic Solution Environment

Nucleic acid-based technologies are an emerging research focus area for pharmacological and biological studies because they are biocompatible and can be designed to produce a variety of scaffolds at the nanometer scale. The use of nucleic acids (ribonucleic acid (RNA) and/or deoxyribonucleic acid (DNA)) as building materials in programming the assemblies and their further functionalization has recently established a new exciting field of RNA and DNA nanotechnology, which have both already produced a variety of different functional nanostructures and nanodevices. It is evident that the resultant architectures require detailed structural and functional characterization and that a variety of technical approaches must be employed to promote the development of the emerging fields. Small-angle X-ray and neutron scattering (SAS) are structural characterization techniques that are well placed to determine the conformation of nucleic acid nanoparticles (NANPs) under varying solution conditions, thus allowing for the optimization of their design. SAS experiments provide information on the overall shapes and particle dimensions of macromolecules and are ideal for following conformational changes of the molecular ensemble as it behaves in solution. In addition, the inherent differences in the neutron scattering of nucleic acids, lipids, and proteins, as well as the different neutron scattering properties of the isotopes of hydrogen, combined with the ability to uniformly label biological macromolecules with deuterium, allow one to characterize the conformations and relative dispositions of the individual components within an assembly of biomolecules. This article will review the application of SAS methods and provide a summary of their successful utilization in the emerging field of NANP technology to date, as well as share our vision on its use in complementing a broad suite of structural characterization tools with some simulated results that have never been shared before.

Conformational Changes and Flexibility of DNA Devices Observed by Small-Angle X-ray Scattering

Self-assembled DNA origami nanostructures enable the creation of precisely defined shapes at the molecular scale. Dynamic DNA devices that are capable of switching between defined conformations could afford completely novel functionalities for diagnostic, therapeutic, or engineering applications. Developing such objects benefits strongly from experimental feedback about conformational changes and 3D structures, ideally in solution, free of potential biases from surface attachment or labeling. Here, we demonstrate that small-angle X-ray scattering (SAXS) can quantitatively resolve the conformational changes of a DNA origami two-state switch device as a function of the ionic strength of the solution. In addition, we show how SAXS data allow for refinement of the predicted idealized three-dimensional structure of the DNA object using a normal mode approach based on an elastic network model. The results reveal deviations from the idealized design geometries that are otherwise difficult to resolve. Our results establish SAXS as a powerful tool to investigate conformational changes and solution structures of DNA origami and we anticipate our methodology to be broadly applicable to increasingly complex DNA and RNA devices.

Footprinting analysis of interactions between the largest eukaryotic RNase P/MRP protein Pop1 and RNase P/MRP RNA components

Ribonuclease (RNase) P and RNase MRP are closely related catalytic ribonucleoproteins involved in the metabolism of a wide range of RNA molecules, including tRNA, rRNA, and some mRNAs. The catalytic RNA component of eukaryotic RNase P retains the core elements of the bacterial RNase P ribozyme; however, the peripheral RNA elements responsible for the stabilization of the global architecture are largely absent in the eukaryotic enzyme. At the same time, the protein makeup of eukaryotic RNase P is considerably more complex than that of the bacterial RNase P. RNase MRP, an essential and ubiquitous eukaryotic enzyme, has a structural organization resembling that of eukaryotic RNase P, and the two enzymes share most of their protein components. Here, we present the results of the analysis of interactions between the largest protein component of yeast RNases P/MRP, Pop1, and the RNA moieties of the enzymes, discuss structural implications of the results, and suggest that Pop1 plays the role of a scaffold for the stabilization of the global architecture of eukaryotic RNase P RNA, substituting for the network of RNA-RNA tertiary interactions that maintain the global RNA structure in bacterial RNase P.

Molecular crowding overcomes the destabilizing effects of mutations in a bacterial ribozyme

The native structure of the Azoarcus group I ribozyme is stabilized by the cooperative formation of tertiary interactions between double helical domains. Thus, even single mutations that break this network of tertiary interactions reduce ribozyme activity in physiological Mg(2+) concentrations. Here, we report that molecular crowding comparable to that in the cell compensates for destabilizing mutations in the Azoarcus ribozyme. Small angle X-ray scattering, native polyacrylamide gel electrophoresis and activity assays were used to compare folding free energies in dilute and crowded solutions containing 18% PEG1000. Crowder molecules allowed the wild-type and mutant ribozymes to fold at similarly low Mg(2+) concentrations and stabilized the active structure of the mutant ribozymes under physiological conditions. This compensation helps explains why ribozyme mutations are often less deleterious in the cell than in the test tube. Nevertheless, crowding did not rescue the high fraction of folded but less active structures formed by double and triple mutants. We conclude that crowding broadens the fitness landscape by stabilizing compact RNA structures without improving the specificity of self-assembly.

Applications of synchrotron-based spectroscopic techniques in studying nucleic acids and nucleic acid-functionalized nanomaterials

In this review, we summarize recent progress in the application of synchrotron-based spectroscopic techniques for nucleic acid research that takes advantage of high-flux and high-brilliance electromagnetic radiation from synchrotron sources. The first section of the review focuses on the characterization of the structure and folding processes of nucleic acids using different types of synchrotron-based spectroscopies, such as X-ray absorption spectroscopy, X-ray emission spectroscopy, X-ray photoelectron spectroscopy, synchrotron radiation circular dichroism, X-ray footprinting and small-angle X-ray scattering. In the second section, the characterization of nucleic acid-based nanostructures, nucleic acid-functionalized nanomaterials and nucleic acid-lipid interactions using these spectroscopic techniques is summarized. Insights gained from these studies are described and future directions of this field are also discussed.

Structural studies of a double-stranded RNA from trypanosome RNA editing by small-angle X-ray scattering

We used small-angle X-ray scattering (SAXS) to evaluate the solution structure of a double-stranded RNA with 32 base pairs. We wanted to compare the solution structure to the crystal structure to assess the impact of the crystal lattice on the overall conformation of the RNA. The RNA was designed to self-anneal and form a head-to-head fusion of two identical mRNA/oligo(U) tail domains (the U-helix) from a trypanosome RNA editing substrate formed by the annealing of a guide RNA to a pre-edited mRNA. This substrate is from the U insertion/deletion RNA editing system of trypanosomes. Each strand in the fusion RNA had 16 purines from the pre-mRNA followed by 16 uracils (Us) from the U-tail at the 3' end of the guide RNA. The strands were designed to form a double helix with blunt ends, but each strand had the potential to form hairpins and single-stranded RNA helices. Hairpins could form by the 3' oligouridylate tract folding back to hybridize with the 5' oligopurine tract and forming an intervening loop. Single-stranded helices could form by the stacking of bases in the polypurine tract. Some of the 16 Us 3' to the polypurine tract may have been unstacked and in random coils. Our SAXS studies showed that the RNA formed a mix of single-stranded structures in the absence of MgCl2. In the presence of MgCl2 at concentrations similar to those in the crystal, the solution structure was consistent with the double-stranded, blunt-ended structure, in agreement with the crystal structure. Here we describe the preparation of RNA samples, data collection with an in-house SAXS instrument designed for biological samples, and the processing and modeling of the scattering data.

Global analysis of riboswitches by small-angle X-ray scattering and calorimetry

Riboswitches are phylogenetically widespread non-coding mRNA domains that directly bind cellular metabolites and regulate transcription, translation, RNA stability or splicing via alternative RNA structures modulated by ligand binding. The details of ligand recognition by many riboswitches have been elucidated using X-ray crystallography and NMR. However, the global dynamics of riboswitch-ligand interactions and their thermodynamic driving forces are less understood. By compiling the work of many laboratories investigating riboswitches using small-angle X-ray scattering (SAXS) and isothermal titration calorimetry (ITC), we uncover general trends and common themes. There is a pressing need for community-wide consensus experimental conditions to allow results of riboswitch studies to be compared rigorously. Nonetheless, our meta-analysis reveals considerable diversity in the extent to which ligand binding reorganizes global riboswitch structures. It also demonstrates a wide spectrum of enthalpy-entropy compensation regimes across riboswitches that bind a diverse set of ligands, giving rise to a relatively narrow range of physiologically relevant free energies and ligand affinities. From the strongly entropy-driven binding of glycine to the predominantly enthalpy-driven binding of c-di-GMP to their respective riboswitches, these distinct thermodynamic signatures reflect the versatile strategies employed by RNA to adapt to the chemical natures of diverse ligands. Riboswitches have evolved to use a combination of long-range tertiary interactions, conformational selection, and induced fit to work with distinct ligand structure, charge, and solvation properties. This article is part of a Special Issue entitled: Riboswitches.

Understanding nucleic acid-ion interactions

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

Molecular crowders and cosolutes promote folding cooperativity of RNA under physiological ionic conditions

Folding mechanisms of functional RNAs under idealized in vitro conditions of dilute solution and high ionic strength have been well studied. Comparatively little is known, however, about mechanisms for folding of RNA in vivo where Mg(2+) ion concentrations are low, K(+) concentrations are modest, and concentrations of macromolecular crowders and low-molecular-weight cosolutes are high. Herein, we apply a combination of biophysical and structure mapping techniques to tRNA to elucidate thermodynamic and functional principles that govern RNA folding under in vivo-like conditions. We show by thermal denaturation and SHAPE studies that tRNA folding cooperativity increases in physiologically low concentrations of Mg(2+) (0.5-2 mM) and K(+) (140 mM) if the solution is supplemented with physiological amounts (∼ 20%) of a water-soluble neutral macromolecular crowding agent such as PEG or dextran. Low-molecular-weight cosolutes show varying effects on tRNA folding cooperativity, increasing or decreasing it based on the identity of the cosolute. For those additives that increase folding cooperativity, the gain is manifested in sharpened two-state-like folding transitions for full-length tRNA over its secondary structural elements. Temperature-dependent SHAPE experiments in the absence and presence of crowders and cosolutes reveal extent of cooperative folding of tRNA on a nucleotide basis and are consistent with the melting studies. Mechanistically, crowding agents appear to promote cooperativity by stabilizing tertiary structure, while those low molecular cosolutes that promote cooperativity stabilize tertiary structure and/or destabilize secondary structure. Cooperative folding of functional RNA under physiological-like conditions parallels the behavior of many proteins and has implications for cellular RNA folding kinetics and evolution.

Advances in the determination of nucleic acid conformational ensembles

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

Molecular crowding favors reactivity of a human ribozyme under physiological ionic conditions

In an effort to relate RNA folding to function under cellular-like conditions, we monitored the self-cleavage reaction of the human hepatitis delta virus-like CPEB3 ribozyme in the background of physiological ionic concentrations and various crowding and cosolute agents. We found that at physiological free Mg(2+) concentrations (∼0.1-0.5 mM), both crowders and cosolutes stimulate the rate of self-cleavage, up to ∼6-fold, but that in 10 mM Mg(2+) (conditions widely used for in vitro ribozyme studies) these same additives have virtually no effect on the self-cleavage rate. We further observe a dependence of the self-cleavage rate on crowder size, wherein the level of rate stimulation is diminished for crowders larger than the size of the unfolded RNA. Monitoring effects of crowding and cosolute agents on rates in biological amounts of urea revealed additive-promoted increases at both low and high Mg(2+) concentrations, with a maximal stimulation of more than 10-fold and a rescue of the rate to its urea-free values. Small-angle X-ray scattering experiments reveal a structural basis for this stimulation in that higher-molecular weight crowding agents favor a more compact form of the ribozyme in 0.5 mM Mg(2+) that is essentially equivalent to the form under standard ribozyme conditions of 10 mM Mg(2+) without a crowder. This finding suggests that at least a portion of the rate enhancement arises from favoring the native RNA tertiary structure. We conclude that cellular-like crowding supports ribozyme reactivity by favoring a compact form of the ribozyme, but only under physiological ionic and cosolute conditions.

Tunable loading of single-stranded DNA on gold nanorods through the displacement of polyvinylpyrrolidone

A quantitative and tunable loading of single-stranded (ss-DNA) molecules onto gold nanorods was achieved through a new method of surfactant exchange. This new method involves the exchange of cetyltrimethylammonium bromide surfactants for an intermediate stabilizing layer of polyvinylpyrrolidone and sodium dodecylsulfate. The intermediate layer of surfactants on the anisotropic gold particles was easily displaced by thiolated ss-DNA, forming a tunable density of single-stranded DNA molecules on the surfaces of the gold nanorods. The success of this ligand exchange process was monitored in part through the combination of extinction, X-ray photoelectron, and infrared absorption spectroscopies. The number of ss-DNA molecules per nanorod for nanorods with a high density of ss-DNA molecules was quantified through a combination of fluorescence measurements and elemental analysis, and the functionality of the nanorods capped with dense monolayers of DNA was assessed using a hybridization assay. Core-satellite assemblies were successfully prepared from spherical particles containing a probe DNA molecule and a nanorod core capped with complementary ss-DNA molecules. The methods demonstrated herein for quantitatively fine tuning and maximizing, or otherwise optimizing, the loading of ss-DNA in monolayers on gold nanorods could be a useful methodology for decorating gold nanoparticles with multiple types of biofunctional molecules.

Nucleic acid structure characterization by small angle X-ray scattering (SAXS)

Small angle X-ray scattering (SAXS) is a powerful method for investigating macromolecular structure in solution. SAXS data provide information about the size and shape of a molecule with a resolution of ∼2 to 3 nm. SAXS is particularly useful for the investigation of nucleic acids, which scatter X-rays strongly due to the electron-rich phosphate backbone. Therefore, SAXS has become an increasingly popular method for modeling nucleic acid structures, an endeavor made tractable by the highly regular helical nature of nucleic acid secondary structures. Recently, SAXS was used in combination with NMR to filter and refine all-atom models of a U2/U6 small nuclear RNA complex. In this unit, general protocols for sample preparation, data acquisition, and data analysis and processing are given. Additionally, examples of correctly and incorrectly processed SAXS data and expected results are provided.

A two-stage mechanism of viral RNA compaction revealed by single molecule fluorescence

Long RNAs often exist as multiple conformers in equilibrium. For the genomes of single-stranded RNA viruses, one of these conformers must include a compacted state allowing the RNA to be confined within the virion. We have used single molecule fluorescence correlation spectroscopy to monitor the conformations of viral genomes and sub-fragments in the absence and presence of coat proteins. Cognate RNA-coat protein interactions in two model viruses cause a rapid collapse in the hydrodynamic radii of their respective RNAs. This is caused by protein binding at multiple sites on the RNA that facilitate additional protein-protein contacts. The collapsed species recruit further coat proteins to complete capsid assembly with great efficiency and fidelity. The specificity in RNA-coat protein interactions seen at single-molecule concentrations reflects the packaging selectivity seen for such viruses in vivo. This contrasts with many in vitro reassembly measurements performed at much higher concentrations. RNA compaction by coat protein or polycation binding are distinct processes, implying that defined RNA-coat protein contacts are required for assembly.

Allosteric tertiary interactions preorganize the c-di-GMP riboswitch and accelerate ligand binding

Cyclic diguanylate (c-di-GMP) is a bacterial second messenger important for physiologic adaptation and virulence. Class-I c-di-GMP riboswitches are phylogenetically widespread and thought to mediate pleiotropic genetic responses to the second messenger. Previous studies suggest that the RNA aptamer domain switches from an extended free state to a compact, c-di-GMP-bound conformation in which two helical stacks dock side-by-side. Single molecule fluorescence resonance energy transfer (smFRET) experiments now reveal that the free RNA exists in four distinct populations that differ in dynamics in the extended and docked conformations. In the presence of c-di-GMP and Mg(2+), a stably docked population (>30 min) becomes predominant. smFRET mutant analysis demonstrates that tertiary interactions distal to the c-di-GMP binding site strongly modulate the RNA population structure, even in the absence of c-di-GMP. These allosteric interactions accelerate ligand recognition by preorganizing the RNA, favoring rapid c-di-GMP binding.

Determining RNA three-dimensional structures using low-resolution data

Knowing the 3-D structure of an RNA is fundamental to understand its biological function. Nowadays X-ray crystallography and NMR spectroscopy are systematically applied to newly discovered RNAs. However, the application of these high-resolution techniques is not always possible, and thus scientists must turn to lower resolution alternatives. Here, we introduce a pipeline to systematically generate atomic resolution 3-D structures that are consistent with low-resolution data sets. We compare and evaluate the discriminative power of a number of low-resolution experimental techniques to reproduce the structure of the Escherichia coli tRNA(VAL) and P4-P6 domain of the Tetrahymena thermophila group I intron. We test single and combinations of the most accessible low-resolution techniques, i.e. hydroxyl radical footprinting (OH), methidiumpropyl-EDTA (MPE), multiplexed hydroxyl radical cleavage (MOHCA), and small-angle X-ray scattering (SAXS). We show that OH-derived constraints are accurate to discriminate structures at the atomic level, whereas EDTA-based constraints apply to global shape determination. We provide a guide for choosing which experimental techniques or combination of thereof is best in which context. The pipeline represents an important step towards high-throughput low-resolution RNA structure determination.

SAXS studies of ion-nucleic acid interactions

Positively charged ions, atoms, or molecules compensate the high negative charge of the nucleic acid backbone. Their presence is critical to the biological function of DNA and RNA. This review focuses on experimental studies probing (a) interactions between small ions and nucleic acids and (b) ion-mediated interactions between nucleic acid duplexes. Experimental results on these simple model systems can be compared with specific theoretical models to validate their predictions. Small angle X-ray scattering (SAXS) provides unique insight into these interactions. Anomalous SAXS reports the spatial correlations of condensed (e.g., locally concentrated) counterions to individual DNA or RNA duplexes. SAXS very effectively reports interactions between nucleic acid helices, which range from strongly repulsive to strongly attractive depending on the ionic species present. The sign and strength of interparticle interactions are easily deduced from dramatic changes in the scattering profiles of interacting duplexes.

Thermodynamics and kinetics of the hairpin ribozyme from atomistic folding/unfolding simulations

We report a set of atomistic folding/unfolding simulations for the hairpin ribozyme using a Monte Carlo algorithm. The hairpin ribozyme folds in solution and catalyzes self-cleavage or ligation via a specific two-domain structure. The minimal active ribozyme has been studied extensively, showing stabilization of the active structure by cations and dynamic motion of the active structure. Here, we introduce a simple model of tertiary-structure formation that leads to a phase diagram for the RNA as a function of temperature and tertiary-structure strength. We then employ this model to capture many folding/unfolding events and to examine the transition-state ensemble (TSE) of the RNA during folding to its active "docked" conformation. The TSE is compact but with few tertiary interactions formed, in agreement with single-molecule dynamics experiments. To compare with experimental kinetic parameters, we introduce a novel method to benchmark Monte Carlo kinetic parameters to docking/undocking rates collected over many single molecular trajectories. We find that topology alone, as encoded in a biased potential that discriminates between secondary and tertiary interactions, is sufficient to predict the thermodynamic behavior and kinetic folding pathway of the hairpin ribozyme. This method should be useful in predicting folding transition states for many natural or man-made RNA tertiary structures.

An iron-dependent and transferrin-mediated cellular uptake pathway for plutonium

Plutonium is a toxic synthetic element with no natural biological function, but it is strongly retained by humans when ingested. Using small-angle X-ray scattering, receptor binding assays and synchrotron X-ray fluorescence microscopy, we find that rat adrenal gland (PC12) cells can acquire plutonium in vitro through the major iron acquisition pathway--receptor-mediated endocytosis of the iron transport protein serum transferrin; however, only one form of the plutonium-transferrin complex is active. Low-resolution solution models of plutonium-loaded transferrins derived from small-angle scattering show that only transferrin with plutonium bound in the protein's C-terminal lobe (C-lobe) and iron bound in the N-terminal lobe (N-lobe) (Pu(C)Fe(N)Tf) adopts the proper conformation for recognition by the transferrin receptor protein. Although the metal-binding site in each lobe contains the same donors in the same configuration and both lobes are similar, the differences between transferrin's two lobes act to restrict, but not eliminate, cellular Pu uptake.

Effects of Mg2+ on the free energy landscape for folding a purine riboswitch RNA

There are potentially several ways Mg2+ might promote formation of an RNA tertiary structure: by causing a general "collapse" of the unfolded ensemble to more compact conformations, by favoring a reorganization of structure within a domain to a form with specific tertiary contacts, and by enhancing cooperative linkages between different sets of tertiary contacts. To distinguish these different modes of action, we have studied Mg2+ interactions with the adenine riboswitch, in which a set of tertiary interactions that forms around a purine-binding pocket is thermodynamically linked to the tertiary "docking" of two hairpin loops in another part of the molecule. Each of four RNA forms with different extents of tertiary structure were characterized by small-angle X-ray scattering. The free energy of interconversion between different conformations in the absence of Mg2+ and the free energy of Mg2+ interaction with each form have been estimated, yielding a complete picture of the folding energy landscape as a function of Mg2+ concentration. At 1 mM Mg2+ (50 mM K+), the overall free energy of stabilization by Mg2+ is large, -9.8 kcal/mol, and about equally divided between its effect on RNA collapse to a partially folded structure and on organization of the binding pocket. A strong cooperative linkage between the two sets of tertiary contacts is intrinsic to the RNA. This quantitation of the effects of Mg2+ on an RNA with two distinct sets of tertiary interactions suggests ways that Mg2+ may work to stabilize larger and more complex RNA structures.

Time resolved SAXS and RNA folding

Small angle X-ray scattering provides low resolution structural information about macromolecules in solution. When coupled with rapid mixing methods, SAXS reports time-dependent conformational changes of RNA induced by the addition of Mg(2+) to trigger folding. Thus time-resolved SAXS provides unique information about the global or overall structures of transient intermediates populated during folding. Notably, SAXS provides information about the earliest folding events, which can evade detection by other methods.

Comparisons between chemical mapping and binding to isoenergetic oligonucleotide microarrays reveal unexpected patterns of binding to the Bacillus subtilis RNase P RNA specificity domain

Microarrays with isoenergetic pentamer and hexamer 2'-O-methyl oligonucleotide probes with LNA (locked nucleic acid) and 2,6-diaminopurine substitutions were used to probe the binding sites on the RNase P RNA specificity domain of Bacillus subtilis. Unexpected binding patterns were revealed. Because of their enhanced binding free energies, isoenergetic probes can break short duplexes, merge adjacent loops, and/or induce refolding. This suggests new approaches to the rational design of short oligonucleotide therapeutics but limits the utility of microarrays for providing constraints for RNA structure determination. The microarray results are compared to results from chemical mapping experiments, which do provide constraints. Results from both types of experiments indicate that the RNase P RNA folds similarly in 1 M Na(+) and 10 mM Mg(2+).

Multistage collapse of a bacterial ribozyme observed by time-resolved small-angle X-ray scattering

Ribozymes must fold into compact, native structures to function properly in the cell. The first step in forming the RNA tertiary structure is the neutralization of the phosphate charge by cations, followed by collapse of the unfolded molecules into more compact structures. The specificity of the collapse transition determines the structures of the folding intermediates and the folding time to the native state. However, the forces that enable specific collapse in RNA are not understood. Using time-resolved SAXS, we report that upon addition of 5 mM Mg(2+) to the Azoarcus group I ribozyme up to 80% of chains form compact structures in less than 1 ms. In 1 mM Mg(2+), the collapse transition produces extended structures that slowly approach the folded state, while > or = 1.5 mM Mg(2+) leads to an ensemble of random coils that fold with multistage kinetics. Increased flexibility of molecules in the intermediate ensemble correlates with a Mg(2+)-dependent increase in the fast folding population and a previously unobserved crossover in the collapse kinetics. Partial denaturation of the unfolded RNA with urea also increases the fraction of chains following the fast-folding pathway. These results demonstrate that the preferred collapse mechanism depends on the extent of Mg(2+)-dependent charge neutralization and that non-native interactions within the unfolded ensemble contribute to the heterogeneity of the ribozyme folding pathways at the very earliest stages of tertiary structure formation.

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.

Idiosyncratically tuned switching behavior of riboswitch aptamer domains revealed by comparative small-angle X-ray scattering analysis

Riboswitches are structured mRNA elements that regulate gene expression upon binding specific cellular metabolites. It is thought that the highly conserved metabolite-binding domains of riboswitches undergo conformational change upon binding their cognate ligands. To investigate the generality of such a mechanism, we employed small-angle X-ray scattering (SAXS). We probed the nature of the global metabolite-induced response of the metabolite-binding domains of four different riboswitches that bind, respectively, thiamine pyrophosphate (TPP), flavin mononucleotide (FMN), lysine, and S-adenosyl methionine (SAM). We find that each RNA is unique in its global structural response to metabolite. Whereas some RNAs exhibit distinct free and bound conformations, others are globally insensitive to the presence of metabolite. Thus, a global conformational change of the metabolite-binding domain is not a requirement for riboswitch function. It is possible that the range of behaviors observed by SAXS, rather than being a biophysical idiosyncrasy, reflects adaptation of riboswitches to the regulatory requirements of their individual genomic context.

Extended structures in RNA folding intermediates are due to nonnative interactions rather than electrostatic repulsion

RNA folding occurs via a series of transitions between metastable intermediate states for Mg(2+) concentrations below those needed to fold the native structure. In general, these folding intermediates are considerably less compact than their respective native states. Our previous work demonstrates that the major equilibrium intermediate of the 154-residue specificity domain (S-domain) of the Bacillus subtilis RNase P RNA is more extended than its native structure. We now investigate two models with falsifiable predictions regarding the origins of the extended intermediate structures in the S-domains of the B. subtilis and the Escherichia coli RNase P RNA that belong to different classes of P RNA and have distinct native structures. The first model explores the contribution of electrostatic repulsion, while the second model probes specific interactions in the core of the folding intermediate. Using small-angle X-ray scattering and Langevin dynamics simulations, we show that electrostatics plays only a minor role, whereas specific interactions largely account for the extended nature of the intermediate. Structural contacts in the core, including a nonnative base pair, help to stabilize the intermediate conformation. We conclude that RNA folding intermediates adopt extended conformations due to short-range, nonnative interactions rather than generic electrostatic repulsion of helical domains. These principles apply to other ribozymes and riboswitches that undergo functionally relevant conformational changes.

The ligand-free state of the TPP riboswitch: a partially folded RNA structure

Riboswitches are elements of mRNA that regulate gene expression by undergoing structural changes upon binding of small ligands. Although the structures of several riboswitches have been solved with their ligands bound, the ligand-free states of only a few riboswitches have been characterized. The ligand-free state is as important for the functionality of the riboswitch as the ligand-bound form, but the ligand-free state is often a partially folded structure of the RNA, with conformational heterogeneity that makes it particularly challenging to study. Here, we present models of the ligand-free state of a thiamine pyrophosphate riboswitch that are derived from a combination of complementary experimental and computational modeling approaches. We obtain a global picture of the molecule using small-angle X-ray scattering data and use an RNA structure modeling software, MC-Sym, to fit local structural details to these data on an atomic scale. We have used two different approaches to obtaining these models. Our first approach develops a model of the RNA from the structures of its constituent junction fragments in isolation. The second approach treats the RNA as a single entity, without bias from the structure of its individual constituents. We find that both approaches give similar models for the ligand-free form, but the ligand-bound models differ for the two approaches, and only the models from the second approach agree with the ligand-bound structure known previously from X-ray crystallography. Our models provide a picture of the conformational changes that may occur in the riboswitch upon binding of its ligand. Our results also demonstrate the power of combining experimental small-angle X-ray scattering data with theoretical structure prediction tools in the determination of RNA structures beyond riboswitches.

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.

Riboswitch conformations revealed by small-angle X-ray scattering

Riboswitches are functional RNA molecules that control gene expression through conformational changes in response to small-molecule ligand binding. In addition, riboswitch 3D structure, like that of other RNA molecules, is dependent on cation-RNA interactions as the RNA backbone is highly negatively charged. Here, we show how small-angle X-ray scattering (SAXS) can be used to probe RNA conformations as a function of ligand and ion concentration. In a recent study of a glycine-binding tandem aptamer from Vibrio cholerae, we have used SAXS data and thermodynamic modeling to investigate how Mg(2+)-dependent folding and glycine binding are energetically coupled. In addition, we have employed ab initio shape reconstruction algorithms to obtain low-resolution models of the riboswitch structure from SAXS data under different solution conditions.

Probing metal binding in the 8-17 DNAzyme by TbIII luminescence spectroscopy

Metal-dependent cleavage activities of the 8-17 DNAzyme were found to be inhibited by Tb(III) ions, and the apparent inhibition constant in the presence of 100 microM of Zn(II) was measured to be 3.3+/-0.3 microM. The apparent inhibition constants increased linearly with increasing Zn(II) concentration, and the inhibition effect could be fully rescued with addition of active metal ions, indicating that Tb(III) is a competitive inhibitor and that the effect is completely reversible. The sensitized Tb(III) luminescence at 543 nm was dramatically enhanced when Tb(III) was added to the DNAzyme-substrate complex. With an inactive DNAzyme in which the GT wobble pair was replaced with a GC Watson-Crick base pair, the luminescence enhancement was slightly decreased. In addition, when the DNAzyme strand was replaced with a complete complementary strand to the substrate, no significant luminescence enhancement was observed. These observations suggest that Tb(III) may bind to an unpaired region of the DNAzyme, with the GT wobble pair playing a role. Luminescence lifetime measurements in D(2)O and H(2)O suggested that Tb(III) bound to DNAzyme is coordinated by 6.7+/-0.2 water molecules and two or three functional groups from the DNAzyme. Divalent metal ions competed for the Tb(III) binding site(s) in the order Co(II)>Zn(II)>Mn(II)>Pb(II)>Ca(II) approximately Mg(II). This order closely follows the order of DNAzyme activity, with the exception of Pb(II). These results indicate that Pb(II), the most active metal ion, competes for Tb(III) binding differently from other metal ions such as Zn(II), suggesting that Pb(II) may bind to a different site from that for the other metal ions including Zn(II) and Tb(III).

Energy barriers, pathways, and dynamics during folding of large, multidomain RNAs

Large, multidomain RNA molecules are generally thought to fold following multiple pathways down rugged landscapes populated with intermediates and traps. A challenge to understanding RNA folding reactions is the complex relationships that exist between the structure of the RNA and its folding landscape. The identification of intermediate species that populate folding landscapes and characterization of elements of their structures are the key components to solving the RNA folding problem. This review explores recent studies that characterize the dominant pathways by which RNA folds, structural and dynamic features of intermediates that populate the folding landscape, and the energy barriers that separate the distinct steps of the folding process.

Reconstructing three-dimensional shape envelopes from time-resolved small-angle X-ray scattering data

Modern computing power has made it possible to reconstruct low-resolution, three-dimensional shapes from solution small-angle X-ray scattering (SAXS) data on biomolecules without a priori knowledge of the structure. In conjunction with rapid mixing techniques, SAXS has been applied to time resolve conformational changes accompanying important biological processes, such as biomolecular folding. In response to the widespread interest in SAXS reconstructions, their value in conjunction with such time-resolved data has been examined. The group I intron from Tetrahymena thermophila and its P4-P6 subdomain are ideal model systems for investigation owing to extensive previous studies, including crystal structures. The goal of this paper is to assay the quality of reconstructions from time-resolved data given the sacrifice in signal-to-noise required to obtain sharp time resolution.

Kinetic barriers and the role of topology in protein and RNA folding

This review compares the folding behavior of proteins and RNAs. Topics covered include the role of topology in the determination of folding rates, major folding events including collapse, properties of denatured states, pathway heterogeneity, and the influence of the mode of initiation on the folding pathway.

X-ray solution scattering (SAXS) combined with crystallography and computation: defining accurate macromolecular structures, conformations and assemblies in solution

Crystallography supplies unparalleled detail on structural information critical for mechanistic analyses; however, it is restricted to describing low energy conformations of macromolecules within crystal lattices. Small angle X-ray scattering (SAXS) offers complementary information about macromolecular folding, unfolding, aggregation, extended conformations, flexibly linked domains, shape, conformation, and assembly state in solution, albeit at the lower resolution range of about 50 A to 10 A resolution, but without the size limitations inherent in NMR and electron microscopy studies. Together these techniques can allow multi-scale modeling to create complete and accurate images of macromolecules for modeling allosteric mechanisms, supramolecular complexes, and dynamic molecular machines acting in diverse processes ranging from eukaryotic DNA replication, recombination and repair to microbial membrane secretion and assembly systems. This review addresses both theoretical and practical concepts, concerns and considerations for using these techniques in conjunction with computational methods to productively combine solution scattering data with high-resolution structures. Detailed aspects of SAXS experimental results are considered with a focus on data interpretation tools suitable to model protein and nucleic acid macromolecular structures, including membrane protein, RNA, DNA, and protein-nucleic acid complexes. The methods discussed provide the basis to examine molecular interactions in solution and to study macromolecular flexibility and conformational changes that have become increasingly relevant for accurate understanding, simulation, and prediction of mechanisms in structural cell biology and nanotechnology.

A large collapsed-state RNA can exhibit simple exponential single-molecule dynamics

The process of large RNA folding is believed to proceed from many collapsed structures to a unique functional structure requiring precise organization of nucleotides. The diversity of possible structures and stabilities of large RNAs could result in non-exponential folding kinetics (e.g. stretched exponential) under conditions where the molecules have not achieved their native state. We describe a single-molecule fluorescence resonance energy transfer (FRET) study of the collapsed-state region of the free energy landscape of the catalytic domain of RNase P RNA from Bacillus stearothermophilus (C(thermo)). Ensemble measurements have shown that this 260 residue RNA folds cooperatively to its native state at >or=1 mM Mg(2+), but little is known about the conformational dynamics at lower ionic strength. Our measurements of equilibrium conformational fluctuations reveal simple exponential kinetics that reflect a small number of discrete states instead of the expected inhomogeneous dynamics. The distribution of discrete dwell times, collected from an "ensemble" of 300 single molecules at each of a series of Mg(2+) concentrations, fit well to a double exponential, which indicates that the RNA conformational changes can be described as a four-state system. This finding is somewhat unexpected under [Mg(2+)] conditions in which this RNA does not achieve its native state. Observation of discrete well-defined conformations in this large RNA that are stable on the seconds timescale at low [Mg(2+)] (<0.1 mM) suggests that even at low ionic strength, with a tremendous number of possible (weak) interactions, a few critical interactions may produce deep energy wells that allow for rapid averaging of motions within each well, and yield kinetics that are relatively simple.

Loop dependence of the stability and dynamics of nucleic acid hairpins

Hairpin loops are critical to the formation of nucleic acid secondary structure, and to their function. Previous studies revealed a steep dependence of single-stranded DNA (ssDNA) hairpin stability with length of the loop (L) as approximately L(8.5 +/- 0.5), in 100 mM NaCl, which was attributed to intraloop stacking interactions. In this article, the loop-size dependence of RNA hairpin stabilities and their folding/unfolding kinetics were monitored with laser temperature-jump spectroscopy. Our results suggest that similar mechanisms stabilize small ssDNA and RNA loops, and show that salt contributes significantly to the dependence of hairpin stability on loop size. In 2.5 mM MgCl2, the stabilities of both ssDNA and RNA hairpins scale as approximately L(4 +/- 0.5), indicating that the intraloop interactions are weaker in the presence of Mg2+. Interestingly, the folding times for ssDNA hairpins (in 100 mM NaCl) and RNA hairpins (in 2.5 mM MgCl2) are similar despite differences in the salt conditions and the stem sequence, and increase similarly with loop size, approximately L(2.2 +/- 0.5) and approximately L(2.6 +/- 0.5), respectively. These results suggest that hairpins with small loops may be specifically stabilized by interactions of the Na+ ions with the loops. The results also reinforce the idea that folding times are dominated by an entropic search for the correct nucleating conformation.

Selective quenching of fluorescence from unbound oligonucleotides by gold nanoparticles as a probe of RNA structure

Binding of small oligonucleotides to the periphery of folded RNA can provide insight into the secondary structure of complex RNA in solution. To discriminate between bound and unbound fluorescein-labeled 2'-O-methyl RNA probes, we use ionically coated gold nanoparticles to selectively adsorb unbound probes and quench their fluorescence. The target is the 3' untranslated region of Bombyx mori R2 RNA. Fluorescence indicates that R2 sequences complementary to some of the probes are accessible for binding in the three-dimensional structure. Hybridization occurs under homogeneous conditions in the absence of the gold nanoparticles so that steric issues associated with chip-based assays are avoided. The assay is compatible with well plate formats, takes less than 5 min, and requires only 2 pmol or less of unlabeled target RNA per probe sequence tested.

Folding of a universal ribozyme: the ribonuclease P RNA

Ribonuclease P is among the first ribozymes discovered, and is the only ubiquitously occurring ribozyme besides the ribosome. The bacterial RNase P RNA is catalytically active without its protein subunit and has been studied for over two decades as a model system for RNA catalysis, structure and folding. This review focuses on the thermodynamic, kinetic and structural frameworks derived from the folding studies of bacterial RNase P RNA.

Evidence that binding of C5 protein to P RNA enhances ribozyme catalysis by influencing active site metal ion affinity

The RNA subunit (P RNA) of the bacterial RNase P ribonucleoprotein is a ribozyme that catalyzes the Mg-dependent hydrolysis of pre-tRNA, but it requires an essential protein cofactor (P protein) in vivo that enhances substrate binding affinities and catalytic rates in a substrate dependent manner. Previous studies of Bacillus subtilis RNase P, containing a Type B RNA subunit, showed that its cognate protein subunit increases the affinity of metal ions important for catalysis, but the functional role of these ions is unknown. Here, we demonstrate that the Mg2+ dependence of the catalytic step for Escherichia coli RNase P, which contains a more common Type A RNA subunit, is also modulated by its cognate protein subunit (C5), indicating that this property is fundamental to P protein. To monitor specifically the binding of active site metal ions, we analyzed quantitatively the rescue by Cd2+ of an inhibitory Rp phosphorothioate modification at the pre-tRNA cleavage site. The results show that binding of C5 protein increases the apparent affinity of the rescuing Cd2+, providing evidence that C5 protein enhances metal ion affinity in the active site, and thus is likely to contribute significantly to rate enhancement at physiological metal ion concentrations.

Effects of osmolytes on RNA secondary and tertiary structure stabilities and RNA-Mg2+ interactions

Osmolytes are small organic molecules accumulated by cells in response to osmotic stress. Although their effects on protein stability have been studied, there has been no systematic documentation of their influence on RNA. Here, the effects of nine osmolytes on the secondary and tertiary structure stabilities of six RNA structures of differing complexity and stability have been surveyed. Using thermal melting analysis, m-values (change in DeltaG degrees of RNA folding per molal concentration of osmolyte) have been measured. All the osmolytes destabilize RNA secondary structure, although to different extents, probably because they favor solubilization of base surfaces. Osmolyte effects on tertiary structure, however, can be either stabilizing or destabilizing. We hypothesize that the stabilizing osmolytes have unfavorable interactions with the RNA backbone, which becomes less accessible to solvent in most tertiary structures. Finally, it was found that as a larger fraction of the negative charge of an RNA tertiary structure is neutralized by hydrated Mg(2+), the RNA becomes less responsive to stabilizing osmolytes and may even be destabilized. The natural selection of osmolytes as protective agents must have been influenced by their effects on the stabilities of functional RNA structures, though in general, the effects of osmolytes on RNA and protein stabilities do not parallel each other. Our results also suggest that some osmolytes can be useful tools for studying intrinsically unstable RNA folds and assessing the mechanisms of Mg(2+)-induced RNA stabilization.

Structural transitions and thermodynamics of a glycine-dependent riboswitch from Vibrio cholerae

Riboswitches are complex folded RNA domains found in noncoding regions of mRNA that regulate gene expression upon small molecule binding. Recently, Breaker and coworkers reported a tandem aptamer riboswitch (VCI-II) that binds glycine cooperatively. Here, we use hydroxyl radical footprinting and small-angle X-ray scattering (SAXS) to study the conformations of this tandem aptamer as a function of Mg(2+) and glycine concentration. We fit a simple three-state thermodynamic model that describes the energetic coupling between magnesium-induced folding and glycine binding. Furthermore, we characterize the structural conformations of each of the three states: In low salt with no magnesium present, the VCI-II construct has an extended overall conformation, presumably representing unfolded structures. Addition of millimolar concentrations of Mg(2+) in the absence of glycine leads to a significant compaction and partial folding as judged by hydroxyl radical protections. In the presence of millimolar Mg(2+) concentrations, the tandem aptamer binds glycine cooperatively. The glycine binding transition involves a further compaction, additional tertiary packing interactions and further uptake of magnesium ions relative to the state in high Mg(2+) but no glycine. Employing density reconstruction algorithms, we obtain low resolution 3-D structures for all three states from the SAXS measurements. These data provide a first glimpse into the structural conformations of the VCI-II aptamer, establish rigorous constraints for further modeling, and provide a framework for future mechanistic studies.