Thermodynamics and folding pathway of tetraloop receptor-mediated RNA helical packing

Lysine-Dependent Entropy Effects in the B. subtilis Lysine Riboswitch: Insights from Single-Molecule Thermodynamic Studies

Riboswitches play an important role in RNA-based sensing/gene regulation control for many bacteria. In particular, the accessibility of multiple conformational states at physiological temperatures allows riboswitches to selectively bind a cognate ligand in the aptamer domain, which triggers secondary structural changes in the expression platform, and thereby "switching" between on or off transcriptional or translational states for the downstream RNA. The present work exploits temperature-controlled, single-molecule total internal reflection fluorescence (TIRF) microscopy to study the thermodynamic landscape of such ligand binding/folding processes, specifically for the Bacillus subtilis lysine riboswitch. The results confirm that the riboswitch folds via an induced-fit (IF) mechanism, in which cognate lysine ligand first binds to the riboswitch before structural rearrangement takes place. The transition state to folding is found to be enthalpically favored (ΔH_fold^‡ < 0), yet with a free-energy barrier that is predominantly entropic (-TΔS_fold^‡ > 0), which results in folding (unfolding) rate constants strongly dependent (independent) of lysine concentration. Analysis of the single-molecule kinetic "trajectories" reveals this rate constant dependence of k_fold on lysine to be predominantly entropic in nature, with the additional lysine conferring preferential advantage to the folding process by the presence of ligands correctly oriented with respect to the riboswitch platform. By way of contrast, van't Hoff analysis reveals enthalpic contributions to the overall folding thermodynamics (ΔH⁰) to be surprisingly constant and robustly independent of lysine concentration. The results demonstrate the crucial role of hydrogen bonding between the ligand and riboswitch platform but with only a relatively modest fraction (45%) of the overall enthalpy change needed to access the transition state and initiate transcriptional switching.

Effects of Molecular Crowders on Single-Molecule Nucleic Acid Folding: Temperature-Dependent Studies Reveal True Crowding vs Enthalpic Interactions

Biomolecular folding in cells can be strongly influenced by spatial overlap/excluded volume interactions (i.e., "crowding") with intracellular solutes. As a result, traditional in vitro experiments with dilute buffers may not accurately recapitulate biomolecule folding behavior in vivo. In order to account for such ubiquitous excluded volume effects, biologically inert polyethylene glycol (PEG) and polysaccharides (dextran and Ficoll) are often used as in vitro crowding agents to mimic in vivo crowding conditions, with a common observation that high concentrations of these polymers stabilize the more compact biomolecule conformation. However, such an analysis can be distorted by differences in polymer interactions with the folded vs unfolded conformers, requiring temperature-dependent analysis of the thermodynamics to reliably assess competing enthalpic vs entropic contributions and thus the explicit role of excluded volume. In this work, temperature-controlled single-molecule fluorescence resonance energy transfer (smFRET) is used to characterize the thermodynamic interaction between nucleic acids and common polymer crowders PEG, dextran, and Ficoll. The results reveal that PEG promotes secondary and tertiary nucleic acid folding by simultaneously increasing the folding rate while decreasing the unfolding rate, with temperature-dependent studies confirming that the source of PEG stabilization is predominantly entropic and consistent with a true excluded volume crowding mechanism. By way of contrast, neither dextran nor Ficoll induces any significant concentration-dependent change in nucleic acid folding stability at room temperature, but instead, stabilization effects gradually appear with a temperature increase. Such a thermal response indicates that both folding enthalpies and entropies are impacted by dextran and Ficoll. A detailed thermodynamic analysis of the kinetics suggests that, instead of true entropic molecular crowding, dextran and Ficoll associate preferentially with the unfolded vs folded nucleic acid conformer as a result of larger solvent accessible surface area, thereby skewing the free energy landscapes through both significant entropic/enthalpic contributions that compete and fortuitously cancel near room temperature.

In vitro selected GUAA tetraloop-binding receptors with structural plasticity and evolvability towards natural RNA structural modules

GNRA tetraloop-binding receptor interactions are key components in the macromolecular assembly of a variety of functional RNAs. In nature, there is an apparent bias for GAAA/11nt receptor and GYRA/helix interactions, with the former interaction being thermodynamically more stable than the latter. While past in vitro selections allowed isolation of novel GGAA and GUGA receptors, we report herein an in vitro selection that revealed several novel classes of specific GUAA receptors with binding affinities comparable to those from natural GAAA/11nt interactions. These GUAA receptors have structural homology with double-locked bulge RNA modules naturally occurring in ribosomal RNAs. They display mutational robustness that enables exploration of the sequence/phenotypic space associated to GNRA/receptor interactions through epistasis. Their thermodynamic self-assembly fitness landscape is characterized by a rugged neutral network with possible evolutionary trajectories toward natural GNRA/receptor interactions. High throughput sequencing analysis revealed synergetic mutations located away from the tertiary interactions that positively contribute to assembly fitness. Our study suggests that the repertoire of GNRA/receptor interactions is much larger than initially thought from the analysis of natural stable RNA molecules and also provides clues for their evolution towards natural GNRA/receptors.

Effect of GCAA stabilizing loops on three- and four-way intramolecular junctions

Tetraloops are a common way of changing the melting behavior of a DNA or RNA structure without changing the sequence of the stem. Because of the ubiquitous nature of tetraloops, our goal is to understand the effect a GCAA tetraloop, which belongs to the GNRA family of tetraloops, has on the unfolding thermodynamics of intramolecular junctions. Specifically, we have described the melting behavior of intramolecular three-way and four-way junctions where a T₅ loop has been replaced with a GCAA tetraloops in different positions. Their thermodynamic profiles, including Δn_Na⁺ and Δn_W, were analyzed based on the position of the tetraloop. We obtained between -16.7 and -27.5 kcal mol^-1 for all junctions studied. The experimental data indicates the influence of the GCAA tetraloop is primarily dictated by the native unfolding of the junction; if the tetraloop is placed on a stem that unfolds as a single domain when the tetraloop is not present, it will unfold as a single domain when the tetraloop is present but with a higher thermal stability. Conversely, if the tetraloop is placed on a stem which unfolds cooperatively with other stems when the tetraloop is not present, the tetraloop will increase the thermal stability of all the stems in the melting domain. The oligonucleotide structure and not the tetraloop itself affects ion uptake; three-way junctions do not gain an increase in ion uptake, but four-way junctions do. This is not the case for water immobilization, where the position of the tetraloop dictates the amount of water immobilized.

Thermodynamics and kinetics of RNA tertiary structure formation in the junctionless hairpin ribozyme

The hairpin ribozyme consists of two RNA internal loops that interact to form the catalytically active structure. This docking transition is a rare example of intermolecular formation of RNA tertiary structure without coupling to helix annealing. We have used temperature-dependent surface plasmon resonance (SPR) to characterize the thermodynamics and kinetics of RNA tertiary structure formation for the junctionless form of the ribozyme, in which loops A and B reside on separate molecules. We find docking to be strongly enthalpy-driven and to be accompanied by substantial activation barriers for association and dissociation, consistent with the structural reorganization of both internal loops upon complex formation. Comparisons with the parallel analysis of a ribozyme variant carrying a 2'-O-methyl modification at the self-cleavage site and with published data in other systems reveal a surprising diversity of thermodynamic signatures, emphasizing the delicate balance of contributions to the free energy of formation of RNA tertiary structure.

Amino Acid Specific Effects on RNA Tertiary Interactions: Single-Molecule Kinetic and Thermodynamic Studies

In light of the current models for an early RNA-based universe, the potential influence of simple amino acids on tertiary folding of ribozymal RNA into biochemically competent structures is speculated to be of significant evolutionary importance. In the present work, the folding-unfolding kinetics of a ubiquitous tertiary interaction motif, the GAAA tetraloop-tetraloop receptor (TL-TLR), is investigated by single-molecule fluorescence resonance energy transfer spectroscopy in the presence of natural amino acids both with (e.g., lysine, arginine) and without (e.g., glycine) protonated side chain residues. By way of control, we also investigate the effects of a special amino acid (e.g., proline) and amino acid mimetic (e.g., betaine) that contain secondary or quaternary amine groups rather than a primary amine group. This combination permits systematic study of amino acid induced (or amino acid like) RNA folding dynamics as a function of side chain complexity, pK_a, charge state, and amine group content. Most importantly, each of the naturally occurring amino acids is found to destabilize the TL-TLR tertiary folding equilibrium, the kinetic origin of which is dominated by a decrease in the folding rate constant (k_dock), also affected by a strongly amino acid selective increase in the unfolding rate constant (k_undock). To further elucidate the underlying thermodynamics, single-molecule equilibrium constants (K_eq) for TL-TLR folding have been probed as a function of temperature, which reveal an amino acid dependent decrease in both overall exothermicity (ΔΔH° > 0) and entropic cost (-TΔΔS° < 0) for the overall folding process. Temperature-dependent studies on the folding/unfolding kinetic rate constants reveal analogous amino acid specific changes in both enthalpy (ΔΔH^⧧) and entropy (ΔΔS^⧧) for accessing the transition state barrier. The maximum destabilization of the TL-TLR tertiary interaction is observed for arginine, which is consistent with early studies of arginine and guanidine ion-inhibited self-splicing kinetics for the full Tetrahymena ribozyme [ Yarus , M. ; Christian , E. L. Nature 1989 , 342 , 349 - 350 ; Yarus , M. Science 1988 , 240 , 1751 - 1758 ].

Entropic stabilization of folded RNA in crowded solutions measured by SAXS

Non-coding RNAs must fold into specific structures that are stabilized by metal ions and other co-solutes in the cell's interior. Large crowder molecules such as PEG stabilize a bacterial group I ribozyme so that the RNA folds in low Mg²⁺ concentrations typical of the cell's interior. To understand the thermodynamic origins of stabilization by crowder molecules, small angle X-ray scattering was used to measure the folding and helix assembly of a bacterial group I ribozyme at different temperatures and in different MgCl₂ and polyethylene glycol (PEG) concentrations. The resulting phase diagrams show that perturbations to folding by each variable do not overlap. A favorable enthalpy change drives the formation of compact, native-like structures, but requires Mg²⁺ ions at all temperatures studied (5–55°C). PEG reduces the entropic cost of helix assembly and increases correlations between RNA segments at all temperatures. The phase diagrams also revealed a semi-compact intermediate between the unfolded and folded ensemble that is locally more flexible than the unfolded state, as judged by SHAPE modification. These results suggest that environmental variables such as temperature and solute density will favor different types of RNA structures.

Isothermal Titration Calorimetry: Assisted Crystallization of RNA-Ligand Complexes

The success rate of nucleic acids/ligands co-crystallization can be significantly improved by performing preliminary biophysical analyses. Among suitable biophysical approaches, isothermal titration calorimetry (ITC) is certainly a method of choice. ITC can be used in a wide range of experimental conditions to monitor in real time the formation of the RNA- or DNA-ligand complex, with the advantage of providing in addition the complete binding profile of the interaction. Following the ITC experiment, the complex is ready to be concentrated for crystallization trials. This chapter describes a detailed experimental protocol for using ITC as a tool for monitoring RNA/small molecule binding, followed by co-crystallization.

Measuring the Kinetics of Molecular Association by Isothermal Titration Calorimetry

The real-time power response inherent in an isothermal titration calorimetry (ITC) experiment provides an opportunity to directly analyze association kinetics, which, together with the conventional measurement of thermodynamic quantities, can provide an incredibly rich description of molecular binding in a single experiment. Here, we detail our application of this method, in which interactions occurring with relaxation times ranging from slightly below the instrument response time constant (12.5 s in this case) to as large as 600 s can be fully detailed in terms of both the thermodynamics and kinetics. In a binding titration scenario, in the most general case an injection can reveal an association rate constant (kon). Under more restrictive conditions, the instrument time constant-corrected power decay following each injection is simply an exponential decay described by a composite rate constant (kobs), from which both kon and the dissociation rate constant (koff) can be extracted. The data also support the viability of this exponential approach, for kon only, for a slightly larger set of conditions. Using a bimolecular RNA folding model and a protein-ligand interaction, we demonstrate and have internally validated this approach to experiment design, data processing, and error analysis. An updated guide to thermodynamic and kinetic regimes accessible by ITC is provided.

Dynamics of riboswitches: Molecular simulations

Riboswitch RNAs play key roles in bacterial metabolism and represent a promising new class of antibiotic targets for treatment of infectious disease. While many studies of riboswitches have been performed, the exact mechanism of riboswitch operation is still not fully understood at the atomistic level of detail. Molecular dynamics simulations are useful for interpreting existing experimental data and producing predictions for new experiments. Here, a wide range of computational studies on riboswitches is reviewed. By elucidating the key principles of riboswitch operation, computation may aid in the effort to design more specific antibiotics with affinities greater than those of the native ligand. Such a detailed understanding may be required to improve efficacy and reduce side effects. These studies are laying the groundwork for understanding the action mechanism of new compounds that inhibit riboswitch activity. Future directions such as magnesium effects, large-scale conformational changes, expression platforms and co-transcriptional folding are also discussed. This article is part of a Special Issue entitled: Riboswitches.

Exploring the electrostatic energy landscape for tetraloop-receptor docking

It has long been appreciated that Mg(2+) is essential for the stabilization of RNA tertiary structure. However, the problem of quantitative prediction for the ion effect in tertiary structure folding remains. By using the virtual bond RNA folding model (Vfold) to generate RNA conformations and the newly improved tightly bound ion model (TBI) to treat ion-RNA interactions, we investigate Mg(2+)-facilitated tetraloop-receptor docking. For the specific construct of the tetraloop-receptor system, the theoretical analysis shows that the Mg(2+)-induced stabilizing force for the docked state is predominantly entropic and the major contribution comes from the entropy of the diffusive ions. Furthermore, our results show that Mg(2+) ions promote tetraloop-receptor docking mainly through the entropy of the diffusive ions. The theoretical prediction agrees with experimental analysis. The method developed in this paper, which combines the theory for the (Mg(2+)) ion effects in RNA folding and RNA conformational sampling, may provide a useful framework for studying the ion effect in the folding of more complex RNA structures.

An RNA folding motif: GNRA tetraloop-receptor interactions

Nearly two decades after Westhof and Michel first proposed that RNA tetraloops may interact with distal helices, tetraloop–receptor interactions have been recognized as ubiquitous elements of RNA tertiary structure. The unique architecture of GNRA tetraloops (N=any nucleotide, R=purine) enables interaction with a variety of receptors, e.g., helical minor grooves and asymmetric internal loops. The most common example of the latter is the GAAA tetraloop–11 nt tetraloop receptor motif. Biophysical characterization of this motif provided evidence for the modularity of RNA structure, with applications spanning improved crystallization methods to RNA tectonics. In this review, we identify and compare types of GNRA tetraloop–receptor interactions. Then we explore the abundance of structural, kinetic, and thermodynamic information on the frequently occurring and most widely studied GAAA tetraloop–11 nt receptor motif. Studies of this interaction have revealed powerful paradigms for structural assembly of RNA, as well as providing new insights into the roles of cations, transition states and protein chaperones in RNA folding pathways. However, further research will clearly be necessary to characterize other tetraloop–receptor and long-range tertiary binding interactions in detail – an important milestone in the quantitative prediction of free energy landscapes for RNA folding.

Intermolecular domain docking in the hairpin ribozyme: metal dependence, binding kinetics and catalysis

The hairpin ribozyme is a prototype small, self-cleaving RNA motif. It exists naturally as a four-way RNA junction containing two internal loops on adjoining arms. These two loops interact in a cation-driven docking step prior to chemical catalysis to form a tightly integrated structure, with dramatic changes occurring in the conformation of each loop upon docking. We investigate the thermodynamics and kinetics of the docking process using constructs in which loop A and loop B reside on separate molecules. Using a novel CD difference assay to isolate the effects of metal ions linked to domain docking, we find the intermolecular docking process to be driven by sub-millimolar concentrations of the exchange-inert Co(NH₃)₆³⁺. RNA self-cleavage requires binding of lower-affinity ions with greater apparent cooperativity than the docking process itself, implying that, even in the absence of direct coordination to RNA, metal ions play a catalytic role in hairpin ribozyme function beyond simply driving loop-loop docking. Surface plasmon resonance assays reveal remarkably slow molecular association, given the relatively tight loop-loop interaction. This observation is consistent with a “double conformational capture” model in which only collisions between loop A and loop B molecules that are simultaneously in minor, docking-competent conformations are productive for binding.

Thermodynamic origins of monovalent facilitated RNA folding

Cations have long been associated with formation of native RNA structure and are commonly thought to stabilize the formation of tertiary contacts by favorably interacting with the electrostatic potential of the RNA, giving rise to an "ion atmosphere". A significant amount of information regarding the thermodynamics of structural transitions in the presence of an ion atmosphere has accumulated and suggests stabilization is dominated by entropic terms. This work provides an analysis of how RNA-cation interactions affect the entropy and enthalpy associated with an RNA tertiary transition. Specifically, temperature-dependent single-molecule fluorescence resonance energy transfer studies have been exploited to determine the free energy (ΔG°), enthalpy (ΔH°), and entropy (ΔS°) of folding for an isolated tetraloop-receptor tertiary interaction as a function of Na(+) concentration. Somewhat unexpectedly, increasing the Na(+) concentration changes the folding enthalpy from a strongly exothermic process [e.g., ΔH° = -26(2) kcal/mol at 180 mM] to a weakly exothermic process [e.g., ΔH° = -4(1) kcal/mol at 630 mM]. As a direct corollary, it is the strong increase in folding entropy [Δ(ΔS°) > 0] that compensates for this loss of exothermicity for the achievement of more favorable folding [Δ(ΔG°) < 0] at higher Na(+) concentrations. In conjunction with corresponding measurements of the thermodynamics of the transition state barrier, these data provide a detailed description of the folding pathway associated with the GAAA tetraloop-receptor interaction as a function of Na(+) concentration. The results support a potentially universal mechanism for monovalent facilitated RNA folding, whereby an increasing monovalent concentration stabilizes tertiary structure by reducing the entropic penalty for folding.

Thermodynamic and kinetic analysis of an RNA kissing interaction and its resolution into an extended duplex

Kissing hairpin interactions form when the loop residues of two hairpins have Watson-Crick complementarity. In a unimolecular context, kissing interactions are important for tertiary folding and pseudoknot formation, whereas in a bimolecular context, they provide a basis for molecular recognition. In some cases, kissing complexes can be a prelude to strand displacement reactions where the two hairpins resolve to form a stable extended intermolecular duplex. The kinetics and thermodynamics of kissing-complex formation and their subsequent strand-displacement reactions are poorly understood. Here, biophysical techniques including isothermal titration calorimetry, surface plasmon resonance, and single-molecule fluorescence have been employed to probe the factors that govern the stability of kissing complexes and their subsequent structural rearrangements. We show that the general understanding of RNA duplex formation can be extended to kissing complexes but that kissing complexes display an unusual level of stability relative to simple duplexes of the same sequence. These interactions form and break many times at room temperature before becoming committed to a slow, irreversible forward transition to the strand-displaced form. Furthermore, using smFRET we show that the primary difference between stable and labile kissing complexes is based almost completely on their off rates. Both stable and labile complexes form at the same rate within error, but less stable species dissociate rapidly, allowing us to understand how these complexes can help generate specificity along a folding pathway or during a gene regulation event.

Characterization of the kinetic and thermodynamic landscape of RNA folding using a novel application of isothermal titration calorimetry

A novel isothermal titration calorimetry (ITC) method was applied to investigate RNA helical packing driven by the GAAA tetraloop-receptor interaction in magnesium and potassium solutions. Both the kinetics and thermodynamics were obtained in individual ITC experiments, and analysis of the kinetic data over a range of temperatures provided Arrhenius activation energies (ΔH(‡)) and Eyring transition state entropies (ΔS(‡)). The resulting rich dataset reveals strongly contrasting kinetic and thermodynamic profiles for this RNA folding system when stabilized by potassium versus magnesium. In potassium, association is highly exothermic (ΔH(25°C) = -41.6 ± 1.2 kcal/mol in 150 mM KCl) and the transition state is enthalpically barrierless (ΔH(‡) = -0.6 ± 0.5). These parameters are significantly positively shifted in magnesium (ΔH(25°C) = -20.5 ± 2.1 kcal/mol, ΔH(‡) = 7.3 ± 2.2 kcal/mol in 0.5 mM MgCl(2)). Mixed salt solutions approximating physiological conditions exhibit an intermediate thermodynamic character. The cation-dependent thermodynamic landscape may reflect either a salt-dependent unbound receptor conformation, or alternatively and more generally, it may reflect a small per-cation enthalpic penalty associated with folding-coupled magnesium uptake.

Entropic origin of Mg2+-facilitated RNA folding

Mg(2+) is essential for the proper folding and function of RNA, though the effect of Mg(2+) concentration on the free energy, enthalpy, and entropy landscapes of RNA folding is unknown. This work exploits temperature-controlled single-molecule FRET methods to address the thermodynamics of RNA folding pathways by probing the intramolecular docking/undocking kinetics of the ubiquitous GAAA tetraloop-receptor tertiary interaction as a function of [Mg(2+)]. These measurements yield the barrier and standard state enthalpies, entropies, and free energies for an RNA tertiary transition, in particular, revealing the thermodynamic origin of [Mg(2+)]-facilitated folding. Surprisingly, these studies reveal that increasing [Mg(2+)] promotes tetraloop-receptor interaction by reducing the entropic barrier (-TΔS(++)(dock)) and the overall entropic penalty (-TΔS(+) (dock)) for docking, with essentially negligible effects on both the activation enthalpy (ΔH(++)(dock)) and overall exothermicity (ΔH(+)(dock)). These observations contrast with the conventional notion that increasing [Mg(2+)] facilitates folding by minimizing electrostatic repulsion of opposing RNA helices, which would incorrectly predict a decrease in ΔH(++)(dock)) and ΔH(+)(dock)) with [Mg(2+)]. Instead we propose that higher [Mg(2+)] can aid RNA folding by decreasing the entropic penalty of counterion uptake and by reducing disorder of the unfolded conformational ensemble.

The molecular interactions that stabilize RNA tertiary structure: RNA motifs, patterns, and networks

RNA molecules adopt specific three-dimensional structures critical to their function. Many essential metabolic processes, including protein synthesis and RNA splicing, are carried out by RNA molecules with elaborate tertiary structures (e.g. 3QIQ, right). Indeed, the ribosome and self-splicing introns are complex RNA machines. But even the coding regions in messenger RNAs and viral RNAs are flanked by highly structured untranslated regions, which provide regulatory information necessary for gene expression. RNA tertiary structure is defined as the three-dimensional arrangement of RNA building blocks, which include helical duplexes, triple-stranded structures, and other components that are held together through connections collectively termed RNA tertiary interactions. The structural diversity of these interactions is now a subject of intense investigation, involving the techniques of NMR, X-ray crystallography, chemical genetics, and phylogenetic analysis. At the same time, many investigators are using biophysical techniques to elucidate the driving forces for tertiary structure formation and the mechanisms for its stabilization. RNA tertiary folding is promoted by maximization of base stacking, much like the hydrophobic effect that drives protein folding. RNA folding also requires electrostatic stabilization, both through charge screening and site binding of metals, and it is enhanced by desolvation of the phosphate backbone. In this Account, we provide an overview of the features that specify and stabilize RNA tertiary structure. A major determinant for overall tertiary RNA architecture is local conformation in secondary-structure junctions, which are regions from which two or more duplexes project. At junctions and other structures, such as pseudoknots and kissing loops, adjacent helices stack on one another, and these coaxial stacks play a major role in dictating the overall architectural form of an RNA molecule. In addition to RNA junction topology, a second determinant for RNA tertiary structure is the formation of sequence-specific interactions. Networks of triple helices, tetraloop-receptor interactions, and other sequence-specific contacts establish the framework for the overall tertiary fold. The third determinant of tertiary structure is the formation of stabilizing stacking and backbone interactions, and many are not sequence specific. For example, ribose zippers allow 2'-hydroxyl groups on different RNA strands to form networks of interdigitated hydrogen bonds, serving to seal strands together and thereby stabilize adjacent substructures. These motifs often require monovalent and divalent cations, which can interact diffusely or through chelation to specific RNA functional groups. As we learn more about the components of RNA tertiary structure, we will be able to predict the structures of RNA molecules from their sequences, thereby obtaining key information about biological function. Understanding and predicting RNA structure is particularly important given the recent discovery that although most of our genome is transcribed into RNA molecules, few of them have a known function. The prevalence of RNA viruses and pathogens with RNA genomes makes RNA drug discovery an active area of research. Finally, knowledge of RNA structure will facilitate the engineering of supramolecular RNA structures, which can be used as nanomechanical components for new materials. But all of this promise depends on a better understanding of the RNA parts list, and how the pieces fit together.

Conformation effects of base modification on the anticodon stem-loop of Bacillus subtilis tRNA(Tyr)

tRNA molecules contain 93 chemically unique nucleotide base modifications that expand the chemical and biophysical diversity of RNA and contribute to the overall fitness of the cell. Nucleotide modifications of tRNA confer fidelity and efficiency to translation and are important in tRNA-dependent RNA-mediated regulatory processes. The three-dimensional structure of the anticodon is crucial to tRNA-mRNA specificity, and the diverse modifications of nucleotide bases in the anticodon region modulate this specificity. We have determined the solution structures and thermodynamic properties of Bacillus subtilis tRNA(Tyr) anticodon arms containing the natural base modifications N(6)-dimethylallyl adenine (i(6)A(37)) and pseudouridine (ψ(39)). UV melting and differential scanning calorimetry indicate that the modifications stabilize the stem and may enhance base stacking in the loop. The i(6)A(37) modification disrupts the hydrogen bond network of the unmodified anticodon loop including a C(32)-A(38)(+) base pair and an A(37)-U(33) base-base interaction. Although the i(6)A(37) modification increases the dynamic nature of the loop nucleotides, metal ion coordination reestablishes conformational homogeneity. Interestingly, the i(6)A(37) modification and Mg(2+) are sufficient to promote the U-turn fold of the anticodon loop of Escherichia coli tRNA(Phe), but these elements do not result in this signature feature of the anticodon loop in tRNA(Tyr).

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.

Enthalpy-driven RNA folding: single-molecule thermodynamics of tetraloop-receptor tertiary interaction

RNA folding thermodynamics are crucial for structure prediction, which requires characterization of both enthalpic and entropic contributions of tertiary motifs to conformational stability. We explore the temperature dependence of RNA folding due to the ubiquitous GAAA tetraloop-receptor docking interaction, exploiting immobilized and freely diffusing single-molecule fluorescence resonance energy transfer (smFRET) methods. The equilibrium constant for intramolecular docking is obtained as a function of temperature (T = 21-47 degrees C), from which a van't Hoff analysis yields the enthalpy (DeltaH degrees) and entropy (DeltaS degrees) of docking. Tetraloop-receptor docking is significantly exothermic and entropically unfavorable in 1 mM MgCl(2) and 100 mM NaCl, with excellent agreement between immobilized (DeltaH degrees = -17.4 +/- 1.6 kcal/mol, and DeltaS degrees = -56.2 +/- 5.4 cal mol(-1) K(-1)) and freely diffusing (DeltaH degrees = -17.2 +/- 1.6 kcal/mol, and DeltaS degrees = -55.9 +/- 5.2 cal mol(-1) K(-1)) species. Kinetic heterogeneity in the tetraloop-receptor construct is unaffected over the temperature range investigated, indicating a large energy barrier for interconversion between the actively docking and nondocking subpopulations. Formation of the tetraloop-receptor interaction can account for approximately 60% of the DeltaH degrees and DeltaS degrees of P4-P6 domain folding in the Tetrahymena ribozyme, suggesting that it may act as a thermodynamic clamp for the domain. Comparison of the isolated tetraloop-receptor and other tertiary folding thermodynamics supports a theme that enthalpy- versus entropy-driven folding is determined by the number of hydrogen bonding and base stacking interactions.

Fluorescence competition assay measurements of free energy changes for RNA pseudoknots

RNA pseudoknots have important functions, and thermodynamic stability is a key to predicting pseudoknots in RNA sequences and to understanding their functions. Traditional methods, such as UV melting and differential scanning calorimetry, for measuring RNA thermodynamics are restricted to temperature ranges around the melting temperature for a pseudoknot. Here, we report RNA pseudoknot free energy changes at 37 °C measured by fluorescence competition assays. Sequence-dependent studies for the loop 1−stem 2 region reveal (1) the individual nearest-neighbor hydrogen bonding (INN-HB) model provides a reasonable estimate for the free energy change when a Watson−Crick base pair in stem 2 is changed, (2) the loop entropy can be estimated by a statistical polymer model, although some penalty for certain loop sequences is necessary, and (3) tertiary interactions can significantly stabilize pseudoknots and extending the length of stem 2 may alter tertiary interactions such that the INN-HB model does not predict the net effect of adding a base pair. The results can inform writing of algorithms for predicting and/or designing RNA secondary structures.

Studying RNA-RNA and RNA-protein interactions by isothermal titration calorimetry

Isothermal Titration Calorimetry (ITC) provides a sensitive and accurate means by which to study the thermodynamics of RNA folding, RNA binding to small molecules, and RNA-protein interactions. The advent of extremely sensitive instrumentation and the increasing availability of ITC in shared facilities have made it increasingly valuable as a tool for RNA biochemistry. As an isothermal measurement, it allows analysis at a defined temperature, distinguishing it from thermal melting approaches (UV melting and differential scanning calorimetry, for instance) that provide thermodynamic information specific to the melting temperature. Residual structures at low temperature in the unfolded state and heat capacity changes lead to potential differences between thermodynamic values measured by ITC and those derived from melting studies. This article describes how ITC can be put to use in the study of RNA biochemistry.