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

We apply synchrotron-based small-angle X-ray scattering to investigate the relationship between compaction, metal binding, and structure formation of two RNAs at 37 degrees C: the 76 nucleotide yeast tRNA(Phe) and the 255 nucleotide catalytic domain of the Bacillus subtilis RNase P RNA. For both RNAs, this method provides direct evidence for the population of a distinct folding intermediate. The relative compaction between the intermediate and the native state does not correlate with the size of the RNA but does correlate well with the amount of surface burial as quantified previously by the urea-dependent m-value. The total compaction process can be described in two major stages. Starting from a completely unfolded state (4-8 M urea, no Mg(2+)), the major amount of compaction occurs upon the dilution of the denaturant and the addition of micromolar amounts of Mg(2+) to form the intermediate. The native state forms in a single transition from the intermediate state upon cooperative binding of three to four Mg(2+) ions. The characterization of this intermediate by small-angle X-ray scattering lends strong support for the cooperative Mg(2+)-binding model to describe the stability of a tertiary RNA.     

The effect of long-range loop-loop interactions on folding of the Tetrahymena self-splicing RNA

Bas?e-pairing between the terminal loops of helices P2.1 and P9.1a (P13) and P2 and P5c (P14) stabilize the folded structure of the Tetrahymena group I intron. Using native gel electrophoresis to analyze the folding kinetics of a natural pre-RNA containing the Tetrahymena intron, we show that P13 and P14 are the only native loop-loop interactions among six possible combinations. Other base-pairing interactions of the loop sequences stabilize misfolded and inactive pre-RNAs. Mismatches in P13 or P14 raised the midpoints and decreased the cooperativity of the Mg(2+)-dependent eqXuilibrium folding transitions. Although some mutations in P13 resulted in slightly higher folding rates, others led to slower folding compared to the wild-type, suggesting that P13 promotes formation of P3 and P7. In contrast, mismatches in P14 increased the rate of folding, suggesting that base-pairing between P5c and P2 stabilizes intermediates in which the catalytic core is misfolded. Although the peripheral helices stabilize the native structure of the catalytic core, our results show that formation of long-range interactions, and competition between correct and incorrect loop-loop base-pairs, decrease the rate at which the active pre-RNA structure is assembled.     

RNA tertiary interactions mediate native collapse of a bacterial group I ribozyme

Large RNAs collapse into compact intermediates in the presence of counterions before folding to the native state. We previously found that collapse of a bacterial group I ribozyme correlates with the formation of helices within the ribozyme core, but occurs at Mg2+ concentrations too low to support stable tertiary structure and catalytic activity. Here, using small-angle X-ray scattering, we show that Mg2+-induced collapse is a cooperative folding transition that can be fit by a two-state model. The Mg2+ dependence of collapse is similar to the Mg2+ dependence of helix assembly measured by partial ribonuclease T1 digestion and of an unfolding transition measured by UV hypochromicity. The correspondence between multiple probes of RNA structure further supports a two-state model. A mutation that disrupts tertiary contacts between the L9 tetraloop and its helical receptor destabilized the compact state by 0.8 kcal/mol, while mutations in the central triplex were less destabilizing. These results show that native tertiary interactions stabilize the compact folding intermediates under conditions in which the RNA backbone remains accessible to solvent.     

Multiple conformations of SAM-II riboswitch detected with SAXS and NMR spectroscopy

Riboswitches are a newly discovered large family of structured functional RNA elements that specifically bind small molecule targets out of a myriad of cellular metabolites to modulate gene expression. Structural studies of ligand-bound riboswitches by X-ray crystallography and nuclear magnetic resonance (NMR) spectroscopy have provided insights into detailed RNA-ligand recognition and interactions. However, the structures of ligand-free riboswitches remain poorly characterized. In this study, we have used a variety of biochemical, biophysical and computational techniques including small-angle X-ray scattering and NMR spectroscopy to characterize the ligand-free and ligand-bound forms of SAM-II riboswitch. Our data demonstrate that the RNA adopts multiple conformations along its folding pathway and suggest that the RNA undergoes marked conformational changes upon Mg(2+) compaction and S-adenosylmethionine (SAM) metabolite binding. Further studies indicated that Mg(2+) ion is not essential for the ligand binding but can stabilize the complex by facilitating loop/stem interactions. In the presence of millimolar concentration of Mg(2+) ion, the RNA samples a more compact conformation. This conformation is near to, but distinct from, the native fold and competent to bind the metabolite. We conclude that the formation of various secondary and tertiary structural elements, including a pseudoknot, occur to sequester the putative Shine-Dalgarno sequence of the RNA only after metabolite binding.     

Cleavage of tRNA by Fe(II)-bleomycin.

We have investigated the action of the chemotherapeutic agent Fe(II)-bleomycin on yeast tRNA(Phe), an RNA of known three-dimensional structure. In the absence of Mg2+ ions, the RNA is cleaved preferentially at two major positions, A31 and G53, both of which are located at the terminal base pairs of hairpin loops, and coincide with the location of tight Mg2+ binding sites. A fragment of the tRNA (residues 47-76) containing the T stem-loop is also cleaved specifically at G53. Cleavage of both the intact tRNA and the tRNA fragment is abolished in the presence of physiological concentrations of Mg2+ (> 0.5 mM). Since Fe(II) is not displaced from bleomycin under these conditions, we infer that tight binding of Mg2+ to tRNA excludes productive interactions between Fe(II)-bleomycin and the RNA. These results also show that loss of cleavage is not due to Mg(2+)-dependent formation of tertiary interactions between the D and T loops. In contrast, cleavage of synthetic DNA analogs of the anticodon and T stem-loops is not detectably inhibited by Mg2+, even at concentrations as high as 50 mM. In addition, the site specificities observed in cleavage of RNA and DNA differ significantly. From these results, and from similar findings with other representative RNA molecules, we suggest that the cleavage of RNA by Fe(II)-bleomycin is unlikely to be important for its therapeutic action.     

Principles of RNA compaction: insights from the equilibrium folding pathway of the P4-P6 RNA domain in monovalent cations

Counterions are required for RNA folding, and divalent metal ions such as Mg(2+) are often critical. To dissect the role of counterions, we have compared global and local folding of wild-type and mutant variants of P4-P6 RNA derived from the Tetrahymena group I ribozyme in monovalent and in divalent metal ions. A remarkably simple picture of the folding thermodynamics emerges. The equilibrium folding pathway in monovalent ions displays two phases. In the first phase, RNA molecules that are initially in an extended conformation enforced by charge-charge repulsion are relaxed by electrostatic screening to a state with increased flexibility but without formation of long-range tertiary contacts. At higher concentrations of monovalent ions, a state that is nearly identical to the native folded state in the presence of Mg(2+) is formed, with tertiary contacts that involve base and backbone interactions but without the subset of interactions that involve specific divalent metal ion-binding sites. The folding model derived from these and previous results provides a robust framework for understanding the equilibrium and kinetic folding of RNA.     

Influence of ground-state structure and Mg2+ binding on folding kinetics of the guanine-sensing riboswitch aptamer domain

Riboswitch RNAs fold into complex tertiary structures upon binding to their cognate ligand. Ligand recognition is accomplished by key residues in the binding pocket. In addition, it often crucially depends on the stability of peripheral structural elements. The ligand-bound complex of the guanine-sensing riboswitch from Bacillus subtilis, for example, is stabilized by extensive interactions between apical loop regions of the aptamer domain. Previously, we have shown that destabilization of this tertiary loop-loop interaction abrogates ligand binding of the G37A/C61U-mutant aptamer domain (Gsw(loop)) in the absence of Mg(2+). However, if Mg(2+) is available, ligand-binding capability is restored by a population shift of the ground-state RNA ensemble toward RNA conformations with pre-formed loop-loop interactions. Here, we characterize the striking influence of long-range tertiary structure on RNA folding kinetics and on ligand-bound complex structure, both by X-ray crystallography and time-resolved NMR. The X-ray structure of the ligand-bound complex reveals that the global architecture is almost identical to the wild-type aptamer domain. The population of ligand-binding competent conformations in the ground-state ensemble of Gsw(loop) is tunable through variation of the Mg(2+) concentration. We quantitatively describe the influence of distinct Mg(2+) concentrations on ligand-induced folding trajectories both by equilibrium and time-resolved NMR spectroscopy at single-residue resolution.     

A thermodynamic framework and cooperativity in the tertiary folding of a Mg2+-dependent ribozyme

The folding thermodynamics of the catalytic domain from the Bacillus subtilis RNase P RNA is analyzed using circular dichroism and fluorescence spectroscopies, hydroxyl radical protection, and catalytic activity. Folding of this 255-nucleotide ribozyme can be described with three populated species: unfolded (U), intermediate (I), and native (N) states. The U-to-I transition primarily involves secondary structure formation, whereas the I-to-N transition is dominated by tertiary structure formation. The I-to-N transition is highly cooperative as indicated by the coincidence of the four probes applied here. Two isothermal methods are used to determine the stability of the N state relative to the I state at 10 and 37 degrees C. The first method measures the extent of Mg(2+)-induced folding without urea or at constant urea concentrations. The second method measures the extent of urea-induced unfolding at constant Mg(2+) concentrations. Via application of a cooperative binding analysis, the Mg(2+) transition midpoint (K(Mg)), the Hill constant (n), and the urea-dependent surface burial parameter (m value) determined by both methods are identical, indicating that they report the same, reversible folding event. Three conclusions can be drawn from these results. (i) The folding free energy of a Mg(2+)-dependent tertiary RNA structure can be described by the K(Mg) and n parameters according to a cooperative Mg(2+) binding model. (ii) The Hill constant for this tertiary RNA structure probably represents the differential number of Mg(2+) ions bound in the I-to-N transition. (iii) Under physiological conditions, the stability of this large ribozyme is similar to that of small globular proteins.     

The tertiary structure of the hairpin ribozyme is formed through a slow conformational search

The hairpin ribozyme is a small catalytic RNA comprised of two internal loops carried on two adjacent arms of a four-way helical junction (4WJ). To achieve catalytic activity, the ribozyme folds into a compact conformation that facilitates the formation of tertiary interactions between the two loops. We have investigated the folding kinetics of the natural 4WJ form of the hairpin ribozyme, as well as a minimal construct consisting of just the two loop-containing duplexes, by means of stopped-flow fluorescence resonance energy transfer between donor and acceptor probes attached to the ends of the loop-bearing arms. Folding was initiated by the addition of Mg(2+) ions or a pseudosubstrate strand to the ribozyme, and the ensuing changes in the emission of both donor and acceptor were monitored over time. Both ribozyme constructs exhibited slow, biphasic kinetic behavior, attributed to two parallel folding pathways leading to compact, docked structures. Two distinct folding rates were observed across a range of Mg(2+) concentrations, and increasing amounts of Mg(2+) accelerated both rates. Notably, both rates were essentially independent of temperature, indicating that the corresponding activation enthalpies were negligible, in contrast to the large activation enthalpies generally observed for RNA folding processes. Instead, the slow folding was due to unfavorable entropy changes in reaching the transition state, indicating that the ribozyme tertiary structure forms through a slow conformational search. These features were observed in both forms of the ribozyme, indicating that the conformational search is confined to the two loop regions and is largely independent of the overall ribozyme architecture. Conformational search may be a general mechanism of tertiary structure formation in RNA.     

Applicability of urea in the thermodynamic analysis of secondary and tertiary RNA folding

The equilibrium folding of a series of self-complementary RNA duplexes and the unmodified yeast tRNA(Phe) is studied as a function of urea and Mg(2+) concentration with optical spectroscopies and chemical modification under isothermal conditions. Via application of standard methodologies from protein folding, the folding free energy and its dependence on urea concentration, the m value, are determined. The free energies of the RNA duplexes obtained from the urea titrations are in good agreement with those calculated from thermal melting studies [Freier, S. I., et al. (1986) Proc. Natl. Acad. Sci. U.S.A. 83, 9373]. The m value correlates with the length of the RNA duplex and is not sensitive to ionic conditions and temperature. The folding of the unmodified yeast tRNA(Phe) can be described by two Mg(2+)-dependent transitions, the second of which corresponds to the formation of the native tertiary structure as confirmed by hydroxyl radical protection and partial nuclease digestion. Both transitions are sensitive to urea and have m values of 0.94 and 1.70 kcal mol(-)(1) M(-)(1), respectively. Although the precise chemical basis of urea denaturation of RNA is uncertain, the m values for the duplexes and tRNA(Phe) are proportional to the amount of the surface area buried in the folding transition. This proportionality, 0.099 cal mol(-)(1) M(-)(1) A(-)(2), is very similar to that observed for proteins, 0.11 cal mol(-)(1) M(-)(1) A(-)(2) [Myers, J., Pace, N., and Scholtz, M. (1995) Protein Sci. 4, 2138]. These results indicate that urea titration can be used to measure both the free energy and the magnitude of an RNA folding transition.     

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.     

An alternative route for the folding of large RNAs: apparent two-state folding by a group II intron ribozyme

Despite a growing literature on the folding of RNA, our understanding of tertiary folding in large RNAs derives from studies on a small set of molecular examples, with primary focus on group I introns and RNase P RNA. To broaden the scope of RNA folding models and to better understand group II intron function, we have examined the tertiary folding of a ribozyme (D135) that is derived from the self-splicing ai5gamma intron from yeast mitochondria. The D135 ribozyme folds homogeneously and cooperatively into a compact, well-defined tertiary structure that includes all regions critical for active-site organization and substrate recognition. When D135 was treated with increasing concentrations of Mg(2+) and then subjected to hydroxyl radical footprinting, similar Mg(2+) dependencies were seen for internalization of all regions of the molecule, suggesting a highly cooperative folding behavior. In this work, we show that global folding and compaction of the molecule have the same magnesium dependence as the local folding previously observed. Furthermore, urea denaturation studies indicate highly cooperative unfolding of the ribozyme that is governed by thermodynamic parameters similar to those for forward folding. In fact, D135 folds homogeneously and cooperatively from the unfolded state to its native, active structure, thereby demonstrating functional reversibility in RNA folding. Taken together, the data are consistent with two-state folding of the D135 ribozyme, which is surprising given the size and multi-domain structure of the RNA. The findings establish that the accumulation of stable intermediates prior to formation of the native state is not a universal feature of RNA folding and that there is an alternative paradigm in which the folding landscape is relatively smooth, lacking rugged features that obstruct folding to the native state.     

Crystal structure of infectious bursal disease virus VP2 subviral particle at 2.6A resolution: implications in virion assembly and immunogenicity

The structural protein VP2 of infectious bursal disease virus (IBDV) spontaneously forms a dodecahedral T=1 subviral particle (SVP), and is a primary immunogen of the virus. In this study, the structure of IBDV SVP was determined in a cubic crystal and refined to 2.6A resolution. It contains 20 independent VP2 subunits in a crystallographic asymmetric unit. Each subunit is folded mainly into a shell domain and a protrusion domain, both with the Swiss-roll topology, plus a small helical base domain. Three VP2 subunits constitute a tight trimer, which is the building block of IBDV (sub)viral particles. The structure revealed a calcium ion bound to three pairs of symmetry-related Asp31 and Asp174 to stabilize the VP2 trimer. Our results of treatment of SVP with EGTA, a Ca(2+)-chelating reagent, indicated that the metal-ion may be important not only in maintaining highly stable quaternary structure but also in regulating the swelling and dissociation of the icosahedral particles. A Ca(2+)-dependent assembly pathway was thus proposed, which involves further interactions between the trimers. The 20 independent subunits showed conformational variations, with the surface loops of the protrusion domain being the most diverse. These loops are targets of the neutralizing antibodies. Several common interactions between the surface loops were clearly observed, suggesting a possible major conformation of the immunogenic epitopes.     

The crystal structure of yeast phenylalanine tRNA at 2.0 A resolution: cleavage by Mg(2+) in 15-year old crystals

We have re-determined the crystal structure of yeast tRNA(Phe) to 2. 0 A resolution using 15 year old crystals. The accuracy of the new structure, due both to higher resolution data and formerly unavailable refinement methods, consolidates the previous structural information, but also reveals novel details. In particular, the water structure around the tightly bound Mg(2+) is now clearly resolved, and hence provides more accurate information on the geometry of the magnesium-binding sites and the role of water molecules in coordinating the metal ions to the tRNA. We have assigned a total of ten magnesium ions and identified a partly conserved geometry for high-affinity Mg(2+ )binding. In the electron density map there is also clear density for a spermine molecule binding in the major groove of the TPsiC arm and also contacting a symmetry-related tRNA molecule. Interestingly, we have also found that two specific regions of the tRNA in the crystals are partially cleaved. The sites of hydrolysis are within the D and anticodon loops in the vicinity of Mg(2+).