Structures of riboswitch RNA reaction states by mix-and-inject XFEL serial crystallography.

Riboswitches are structural RNA elements that are generally located in the 5' untranslated region of messenger RNA. During regulation of gene expression, ligand binding to the aptamer domain of a riboswitch triggers a signal to the downstream expression platform. A complete understanding of the structural basis of this mechanism requires the ability to study structural changes over time. Here we use femtosecond X-ray free electron laser (XFEL) pulses to obtain structural measurements from crystals so small that diffusion of a ligand can be timed to initiate a reaction before diffraction. We demonstrate this approach by determining four structures of the adenine riboswitch aptamer domain during the course of a reaction, involving two unbound apo structures, one ligand-bound intermediate, and the final ligand-bound conformation. These structures support a reaction mechanism model with at least four states and illustrate the structural basis of signal transmission. The three-way junction and the P1 switch helix of the two apo conformers are notably different from those in the ligand-bound conformation. Our time-resolved crystallographic measurements with a 10-second delay captured the structure of an intermediate with changes in the binding pocket that accommodate the ligand. With at least a 10-minute delay, the RNA molecules were fully converted to the ligand-bound state, in which the substantial conformational changes resulted in conversion of the space group. Such notable changes in crystallo highlight the important opportunities that micro- and nanocrystals may offer in these and similar time-resolved diffraction studies. Together, these results demonstrate the potential of 'mix-and-inject' time-resolved serial crystallography to study biochemically important interactions between biomacromolecules and ligands, including those that involve large conformational changes.

Dynamics of biological macromolecules: not a simple slaving by hydration water.

We studied the dynamics of hydrated tRNA using neutron and dielectric spectroscopy techniques. A comparison of our results with earlier data reveals that the dynamics of hydrated tRNA is slower and varies more strongly with temperature than the dynamics of hydrated proteins. At the same time, tRNA appears to have faster dynamics than DNA. We demonstrate that a similar difference appears in the dynamics of hydration water for these biomolecules. The results and analysis contradict the traditional view of slaved dynamics, which assumes that the dynamics of biological macromolecules just follows the dynamics of hydration water. Our results demonstrate that the dynamics of biological macromolecules and their hydration water depends strongly on the chemical and three-dimensional structures of the biomolecules. We conclude that the whole concept of slaving dynamics should be reconsidered, and that the mutual influence of biomolecules and their hydration water must be taken into account.

Dynamics of tRNA at different levels of hydration.

The influence of hydration on the nanosecond timescale dynamics of tRNA is investigated using neutron scattering spectroscopy. Unlike protein dynamics, the dynamics of tRNA is not affected by methyl group rotation. This allows for a simpler analysis of the influence of hydration on the conformational motions in RNA. We find that hydration affects the dynamics of tRNA significantly more than that of lysozyme. Both the characteristic length scale and the timescale of the conformational motions in tRNA depend strongly on hydration. Even the characteristic temperature of the so-called "dynamical transition" appears to be hydration-dependent in tRNA. The amplitude of the conformational motions in fully hydrated tRNA is almost twice as large as in hydrated lysozyme. We ascribe these differences to a more open and flexible structure of hydrated RNA, and to a larger fraction and different nature of hydrophilic sites. The latter leads to a higher density of water that makes the biomolecule more flexible. All-atom molecular-dynamics simulations are used to show that the extent of hydration is greater in tRNA than in lysozyme. We propose that water acts as a "lubricant" in facilitating enhanced motion in solvated RNA molecules.

Persistence length changes dramatically as RNA folds.

We determine the persistence length l(p) for a bacterial group I ribozyme as a function of concentration of monovalent and divalent cations by fitting the distance distribution functions P(r) obtained from small angle x-ray scattering intensity data to the asymptotic form of the calculated P(WLC)(r) for a wormlike chain. The l(p) values change dramatically over a narrow range of Mg(2+) concentration from approximately 21 Angstroms in the unfolded state (U) to approximately 10 Angstroms in the compact (I(C)) and native states. Variations in l(p) with increasing Na(+) concentration are more gradual. In accord with the predictions of polyelectrolyte theory we find l(p) alpha 1/kappa(2) where kappa is the inverse Debye-screening length.

Folding mechanisms of group I ribozymes: role of stability and contact order.

The mechanism by which RNA molecules assemble into unique three-dimensional conformations is important for understanding their function, regulation and interactions with substrates. The Tetrahymena group I ribozyme is an excellent model system for understanding RNA folding mechanisms, because the catalytic activity of the native RNA is easily measured. Folding of the Tetrahymena ribozyme is dominated by intermediates in which the stable P4-P6 domain is correctly formed, but the P3-P9 domain is partially misfolded. The propensity of the RNA to misfold depends on the relative stability of native and non-native interactions. Circular permutation of the Tetrahymena ribozyme shows that the distance in the primary sequence between native interactions also influences the folding pathway.

Role of counterion condensation in folding of the Tetrahymena ribozyme. II. Counterion-dependence of folding kinetics.

Condensed counterions contribute to the stability of compact structures in RNA, largely by reducing electrostatic repulsion among phosphate groups. Varieties of cations induce a collapsed state in the Tetrahymena ribozyme that is readily transformed to the catalytically active structure in the presence of Mg2+. Native gel electrophoresis was used to compare the effects of the valence and size of the counterion on the kinetics of this transition. The rate of folding was found to decrease with the charge of the counterion. Transitions in monovalent ions occur 20- to 40-fold faster than transitions induced by multivalent metal ions. These results suggest that multivalent cations yield stable compact structures, which are slower to reorganize to the native conformation than those induced by monovalent ions. The folding kinetics are 12-fold faster in the presence of spermidine3+ than [Co(NH3)6]3+, consistent with less effective stabilization of long-range RNA interactions by polyamines. Under most conditions, the observed folding rate decreases with increasing counterion concentration. In saturating amounts of counterion, folding is accelerated by addition of urea. These observations indicate that reorganization of compact intermediates involves partial unfolding of the RNA. We find that folding of the ribozyme is most efficient in a mixture of monovalent salt and Mg2+. This is attributed to competition among counterions for binding to the RNA. The counterion dependence of the folding kinetics is discussed in terms of the ability of condensed ions to stabilize compact structures in RNA.

Early events in RNA folding.

We describe a conceptual framework for understanding the way large RNA molecules fold based on the notion that their free-energy landscape is rugged. A key prediction of our theory is that RNA folding can be described by the kinetic partitioning mechanism (KPM). According to KPM a small fraction of molecules folds rapidly to the native state whereas the remaining fraction is kinetically trapped in a low free-energy non-native state. This model provides a unified description of the way RNA and proteins fold. Single-molecule experiments on Tetrahymena ribozyme, which directly validate our theory, are analyzed using KPM. We also describe the earliest events that occur on microsecond time scales in RNA folding. These must involve collapse of RNA molecules that are mediated by counterion-condensation. Estimates of time scales for the initial events in RNA folding are provided for the Tetrahymena ribozyme.

Role of counterion condensation in folding of the Tetrahymena ribozyme. I. Equilibrium stabilization by cations.

Folding of RNA into an ordered, compact structure requires substantial neutralization of the negatively charged backbone by positively charged counterions. Using a native gel electrophoresis assay, we have examined the effects of counterion condensation upon the equilibrium folding of the Tetrahymena ribozyme. Incubation of the ribozyme in the presence of mono-, di- and trivalent ions induces a conformational state that is capable of rapidly forming the native structure upon brief exposure to Mg2+. The cation concentration dependence of this transition is directly correlated with the charge of the counterion used to induce folding. Substrate cleavage assays confirm the rapid onset of catalytic activity under these conditions. These results are discussed in terms of classical counterion condensation theory. A model for folding is proposed which predicts effects of charge, ionic radius and temperature on counterion-induced RNA folding transitions.

Probing RNA folding pathways by RNA fingerprinting.

This unit provides protocols for using native polyacrylamide gel electrophoresis to distinguish folding and unfolding conformers of RNA. It is useful for studying conformers that can exchange in a period of minutes or seconds, and that are thus difficult to study by solution-based methods. Conformers that have been separated and immobilized in the gel matrix can be used to studying catalytic activity with or without being eluted from the gel. The method can be applied to a wide varied of catalytic RNAs and RNA-protein complexes.

Time-resolved hydroxyl radical footprinting of RNA with X-rays.

RNA tertiary structure and protein interactions can be detected by protection from hydroxyl radical cleavage. Generation of hydroxyl radicals with a synchrotron X-ray beam provides analysis on a short timescale (50 msec to 100 sec), which enables the structures of folding intermediates or other transient conformational states to be determined. This unit provides detailed instructions on the use of the synchrotron beamline for hydroxyl radical footprinting.

Multiple folding pathways for the P4-P6 RNA domain.

We recently described site-specific pyrene labeling of RNA to monitor Mg(2+)-dependent equilibrium formation of tertiary structure. Here we extend these studies to follow the folding kinetics of the 160-nucleotide P4-P6 domain of the Tetrahymena group I intron RNA, using stopped-flow fluorescence with approximately 1 ms time resolution. Pyrene-labeled P4-P6 was prepared using a new phosphoramidite that allows high-yield automated synthesis of oligoribonucleotides with pyrene incorporated at a specific 2'-amino-2'-deoxyuridine residue. P4-P6 forms its higher-order tertiary structure rapidly, with k(obs) = 15-31 s(-1) (t(1/2) approximately 20-50 ms) at 35 degrees C and [Mg(2+)] approximately 10 mM in Tris-borate (TB) buffer. The folding rate increases strongly with temperature from 4 to 45 degrees C, demonstrating a large activation enthalpy DeltaH(double dagger) approximately 26 kcal/mol; the activation entropy DeltaS(double dagger) is large and positive. In low ionic strength 10 mM sodium cacodylate buffer at 35 degrees C, a slow (t(1/2) approximately 1 s) folding component is also observed. The folding kinetics are both ionic strength- and temperature-dependent; the slow phase vanishes upon increasing [Na(+)] in the cacodylate buffer, and the kinetics switch completely from fast at 30 degrees C to slow at 40 degrees C. Using synchrotron hydroxyl radical footprinting, we confirm that fluorescence monitors the same kinetic events as hydroxyl radical cleavage, and we show that the previously reported slow P4-P6 folding kinetics apply only to low ionic strength conditions. One model to explain the fast and slow folding kinetics postulates that some tertiary interactions are present even without Mg(2+) in the initial state. The fast kinetic phase reflects folding that is facilitated by these interactions, whereas the slow kinetics are observed when these interactions are disrupted at lower ionic strength and higher temperature.

Refolding of rRNA exons enhances dissociation of the Tetrahymena intron.

Self-splicing of the Tetrahymena pre-rRNA is inhibited by a conserved rRNA hairpin P(-1) upstream of the 5' splice site. P(-1) inhibits self-splicing by competing with formation of the P1 splice site helix. Here we show that the P(-1) hairpin also enhances dissociation of the spliced products, which was monitored by native gel electrophoresis. Mutations that stabilize the rRNA hairpin increase the rate of dissociation approximately 10-fold, from 0.5 min(-1) for the wild-type RNA to approximately 4 min(-1) at 30 degrees C. Conversely, mutations or oligonucleotides that inhibit refolding of the exons and that stabilize the P1 helix decrease the rate of product release. The results suggest that refolding of products can be used to stimulate the turnover of ribozyme-catalyzed reactions. In the pre-rRNA, this conformational change helps shift the equilibrium of self-splicing toward the mature rRNA.

Folding mechanism of the Tetrahymena ribozyme P4-P6 domain.

Synchrotron X-ray-dependent hydroxyl radical footprinting was used to probe the folding kinetics of the P4-P6 domain of the Tetrahymena group I ribozyme, which forms a stable, closely packed tertiary structure. The 160-nt domain folds independently at a similar rate (approximately 2 s(-1)) as it does in the ribozyme, when folding is measured in 10 mM sodium cacodylate and 10 mM MgCl(2). Surprisingly, tertiary interactions around a three-helix junction (P5abc) within the P4-P6 domain fold at least 25 times more rapidly (k >/= 50 s(-1)) in isolation, than when part of the wild-type P4-P6 RNA. This difference implies that long-range interactions in the P4-P6 domain can interfere with folding of P5abc. P4-P6 was observed to fold much faster at higher ionic strength than in 10 mM sodium cacodylate. Analytical centrifugation was used to measure the sedimentation and diffusion coefficients of the unfolded RNA. The hydrodynamic radius of the RNA decreased from 58 to 46 A over the range of 0-100 mM NaCl. We propose that at low ionic strength, the addition of Mg(2+) causes the domain to collapse to a compact intermediate where P5abc is trapped in a non-native structure. At high ionic strength, the RNA rapidly collapses to the native structure. Faster folding most likely results from a different average initial conformation of the RNA in higher salt conditions.

Recent insights on RNA folding mechanisms from catalytic RNA.

Methods for probing RNA structure in real time have revealed that initial folding steps are complete in less than a second. Refolding of large catalytic RNAs in vitro often results in long-lived intermediates that reach the native structure very slowly. These kinetically trapped intermediates arise from alternative secondary structures that form early in the folding process. In cells, proteins modulate the outcome of RNA folding reactions by stabilizing specific conformations or by accelerating refolding of misfolded intermediates. At the same time, competition between metastable conformations provides a means for regulating the biological activity of transcripts.

Maximizing RNA folding rates: a balancing act.

Large ribozymes typically require very long times to refold into their active conformation in vitro, because the RNA is easily trapped in metastable misfolded structures. Theoretical models show that the probability of misfolding is reduced when local and long-range interactions in the RNA are balanced. Using the folding kinetics of the Tetrahymena ribozyme as an example, we propose that folding rates are maximized when the free energies of forming independent domains are similar to each other. A prediction is that the folding pathway of the ribozyme can be reversed by inverting the relative stability of the tertiary domains. This result suggests strategies for optimizing ribozyme sequences for therapeutics and structural studies.

Fast folding of a ribozyme by stabilizing core interactions: evidence for multiple folding pathways in RNA.

Folding of the Tetrahymena ribozyme under physiological conditions in vitro is limited by slow conversion of long-lived intermediates to the active structure. These intermediates arise because the most stable domain of the ribozyme folds 10-50 times more rapidly than the core region containing helix P3. Native gel electrophoresis and time-resolved X-ray-dependent hydroxyl radical cleavage revealed that mutations that weaken peripheral interactions between domains accelerated folding fivefold, while a point mutation that stabilizes P3 enabled 80 % of the mutant RNA to reach the native conformation within 30 seconds at 22 degrees C. The P3 mutation increased the folding rate of the catalytic core as much as 50-fold, so that both domains of the ribozyme were formed at approximately the same rate. The results show that the ribozyme folds rapidly without significantly populating metastable intermediates when native interactions in the ribozyme core are stabilized relative to peripheral structural elements.

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.

Facilitation of group I splicing in vivo: misfolding of the Tetrahymena IVS and the role of ribosomal RNA exons.

Self-splicing of the group I IVS from Tetrahymena thermophila rDNA is limited by the time required for the RNA to reach its active conformation. In vitro, folding is slow because the pre-rRNA becomes kinetically trapped in inactive structures. In vivo, splicing is 50 times more rapid, implying that misfolding of the pre-rRNA is corrected. Exon mutations that inhibit self-splicing 100-fold in vitro were fully rescued when the pre-rRNA containing the IVS was expressed in E. coli. In contrast, IVS mutations that cause misfolding were only partially suppressed at 42 degrees C, and doubled the activation energy of splicing. These results suggest that intracellular folding of the pre-rRNA involves metastable intermediates similar to those observed in vitro. Precursors with natural rRNA exons were more active and less cold-sensitive than those with non-rRNA exons. This shows that the rRNA reduces misfolding of the IVS, thereby facilitating splicing of the pre-rRNA in vivo.