Thermodynamic origins of monovalent facilitated RNA folding

High-throughput determination of RNA tertiary contact thermodynamics by quantitative DMS chemical mapping

Structured RNAs often contain long-range tertiary contacts that are critical to their function. Despite the importance of tertiary contacts, methods to measure their thermodynamics are low throughput or require specialized instruments. Here, we introduce a new quantitative chemical mapping method (qMaPseq) to measure Mg2+-induced formation of tertiary contact thermodynamics in a high-throughput manner using standard biochemistry equipment. With qMaPseq, we measured the ΔG of 98 unique tetraloop/tetraloop receptor (TL/TLR) variants in a one-pot reaction. These results agree well with measurements from specialized instruments (R2= 0.64). Furthermore, the DMS reactivity of the TL directly correlates to the stability of the contact (R2= 0.68), the first direct evidence that a single DMS reactivity measurement reports on thermodynamics. Combined with structure prediction, DMS reactivity allowed the development of experimentally accurate 3D models of TLR mutants. These results demonstrate that qMaPseq is broadly accessible, high-throughput and directly links DMS reactivity to thermodynamics.

Localized Plasmonic Heating for Single-Molecule DNA Rupture Measurements in Optical Tweezers

To date, studies on the thermodynamic and kinetic processes that underlie biological function and nanomachine actuation in biological- and biology-inspired molecular constructs have primarily focused on photothermal heating of ensemble systems, highlighting the need for probes that are localized within the molecular construct and capable of resolving single-molecule response. Here we present an experimental demonstration of wavelength-selective, localized heating at the single-molecule level using the surface plasmon resonance of a 15 nm gold nanoparticle (AuNP). Our approach is compatible with force-spectroscopy measurements and can be applied to studies of the single-molecule thermodynamic properties of DNA origami nanomachines as well as biomolecular complexes. We further demonstrate wavelength selectivity and establish the temperature dependence of the reaction coordinate for base-pair disruption in the shear-rupture geometry, demonstrating the utility and flexibility of this approach for both fundamental studies of local (nanometer-scale) temperature gradients and rapid and multiplexed nanomachine actuation.

Factors Affecting Stability of RNA - Temperature, Length, Concentration, pH, and Buffering Species

RNA is prone to both chemical degradation and/or physical instability. Some of the factors affecting stability of RNA in solution are its length, 3' poly A tail and 5' cap integrity, excipients, buffering species, pH of the solution, nucleases, and divalent cations. In this work, we showed the effect of temperature, messenger RNA (mRNA) length, buffering species, pH of the solution, and the concentration of mRNA on its chemical and physical stability. Our thermodynamic analysis of a 4000 nucleotide-long mRNA measured an activation energy of 31.5 kcal/mol normalized per phosphodiester backbone. We found mRNA length to be negatively correlated to its stability. Buffering species and pH of the solution affected mRNA integrity along with affecting the onset temperature of melting obtained by Differential Scanning Calorimetry (DSC) thermograms. It was also found that increasing the concentration of mRNA in solution increased its stability.

Ionic Cooperativity between Lysine and Potassium in the Lysine Riboswitch: Single-Molecule Kinetic and Thermodynamic Studies

Functionality in many biological systems, including proteins and nucleic acid structures, including protein and nucleic acid riboswitch structures, can depend on cooperative kinetic behavior between multiple small molecule ligands. In this work, single-molecule FRET data on the Bacillus subtilis lysine riboswitch reveals that affinity for the cognate lysine ligand increases significantly with K⁺, providing evidence for synergism between lysine/K⁺ binding to the aptamer and successful folding of the riboswitch. To describe/interpret this more complex kinetic scenario, we explore the conventional 4-state ("square") model for aptamer binding as a function of K⁺. Extension into this additional dimension generates a novel "cube" model for riboswitch folding dynamics with respect to lysine/K⁺ binding, revealing that riboswitch folding (k_fold) and unfolding (k_unfold) rate constants increase and decrease dramatically with K⁺, respectively. Furthermore, temperature-dependent single-molecule kinetic studies indicate that the presence of K⁺ entropically enhances the transition state barrier to folding but partially compensates for this by increasing the overall exothermicity for lysine binding. We rationalize this behavior as evidence that K⁺ facilitates hydrogen bonding between the negatively charged carboxyl group of lysine and the RNA, increasing structural rigidity and lowering entropy in the binding pocket. Finally, we explore the effects of cation size with Na⁺ and Cs⁺ studies to demonstrate that K⁺ is optimally suited for bridging interactions between lysine and the riboswitch aptamer domain. Regulation of lysine production and transport, dictated by the riboswitch's ability to recognize and bind lysine, is therefore intimately tied to the presence of K⁺ in the binding pocket and is strongly modulated by local cation conditions. The results suggest an increase in lysine riboswitch functionality by sensitivity to additional species in the cellular riboswitch environment.

Entropy Driving the Mg2+-Induced Folding of TPP Riboswitch RNA

Mg²⁺ is well known to facilitate the structural folding of RNA. However, the thermodynamic and dynamic roles of Mg²⁺ in RNA folding remain elusive. Here, we exploit single-molecule fluorescence resonance energy transfer (smFRET) and isothermal titration calorimetry (ITC) to study the mechanism of Mg²⁺ in facilitating the folding of thiamine pyrophosphate (TPP) riboswitch RNA. The results of smFRET identify that the presence of Mg²⁺ compacts the RNA and enlarges the conformational dispersity among individual RNA molecules, resulting in a large gain of entropy. The compact yet flexible conformations triggered by Mg²⁺ may help the riboswitch recognize its specific ligand and further fold. This is supported by the ITC experiments, in which the Mg²⁺-induced RNA folding is driven by entropy (ΔS) instead of enthalpy (ΔH). Our results complement the understanding of the Mg²⁺-induced RNA folding. The strategy developed in this work can be used to model other RNAs' folding under different conditions.

Sensitive and Low-Bias Transcriptome Sequencing Using Agarose PCR

Transcriptome sequencing has emerged as an important research tool for exploring the mysteries of life at the single-cell level. However, its wide application is limited by the bias associated with the amplification reactions which is essential for library building of trace RNA. In this study, low-melting-point agarose was added to the amplification reactions to take advantage of its molecular crowding effect and polymer cross-linked structure to improve the sensitivity of the reactions and reduce bias. To further evaluate the performance of the method, it was applied to transcriptome sequencing of microregion samples from brain tissue sections of mice with Parkinson's disease at the single cell level. The results showed that agarose PCR had better performance than in-tube PCR. Further application of agarose PCR to transcriptome library sequencing could obtain data closer to that of unamplified. With the addition of low melting point agarose, the sensitivity of the amplification reaction was significantly increased, while homogeneity was increased by approximately 2-fold. Not only that, but this work also provides 11% sensitivity improvement for spatial transcriptomic study on Parkinson's disease-associated gene detection. The agarose PCR provides a new tool for efficient and homogeneous amplification of trace samples and can be widely used for spatial transcriptome library sequencing and studies.

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.

Apolar chemical environments compact unfolded RNAs and can promote folding

It is well documented that the structure, and thus function, of nucleic acids depends on the chemical environment surrounding them, which often includes potential proteinaceous binding partners. The nonpolar amino acid side chains of these proteins will invariably alter the polarity of the local chemical environment around the nucleic acid. However, we are only beginning to understand how environmental polarity generally influences the structural and energetic properties of RNA folding. Here, we use a series of aqueous-organic cosolvent mixtures to systematically modulate the solvent polarity around two different RNA folding constructs that can form either secondary or tertiary structural elements. Using single-molecule Förster resonance energy transfer spectroscopy to simultaneously monitor the structural and energetic properties of these RNAs, we show that the unfolded conformations of both model RNAs become more compact in apolar environments characterized by dielectric constants less than that of pure water. In the case of tertiary structure formation, this compaction also gives rise to more energetically favorable folding. We propose that these physical changes arise from an enhanced accumulation of counterions in the low dielectric environment surrounding the unfolded RNA.

Salt dependent mesoscopic model for RNA at multiple strand concentrations

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

Chirality-Dependent Amino Acid Modulation of RNA Folding

The preponderance of a specific d- or l-chirality in fats, sugars, amino acids, nucleic acids, and so on is ubiquitous in nature, yet the biological origin of such chiral dominance (i.e., with one enantiomer overwhelmingly present) remains an open question. One plausible proposal for the predominance of l-chirality in amino acids could be through evolutionary templating of chiral RNA-folding via chaperone activity. To help evaluate this possibility, single molecule fluorescence experiments have been performed that measure the chiral dependence of chaperone folding dynamics for the simple tetraloop-tetraloop receptor (TL-TLR) tertiary binding motif in the presence of a series of chiral amino acids. Specifically, d- vs l-arginine is found to accelerate the unfolding of this RNA motif in a chirally selective fashion, with temperature-dependent studies of the kinetics performed to extract free energy, enthalpy, and entropy landscapes for the underlying thermodynamics. Furthermore, all-atom molecular dynamics (MD) simulations are pursued to provide additional physical insight into this chiral sensitivity, which reveal enantiomer-specific sampling of nucleic acid surfaces by d- vs l-arginine and support a putative mechanism for chirally specific denaturation of RNA tertiary structure by arginine but not other amino acids.

Simple and Analytical Model of RNA Collapse

An analytical model for the free energy change during collapse of an RNA molecule from an extended to a compact conformation is proposed. It considers explicit binding of water and ion molecules to the RNA and the exchange of these molecules with the aqueous solution. Microscopic states of the system are captured on a two-dimensional square lattice and evaluated using contact energies. We compute the free energy as a function of a collapse variable and the number of ions bound to the RNA. The major driving force to the collapse of the RNA chain is the gain in water entropy once expelled from the surface of the RNA molecule illustrated by decomposing the free energy into species contributions and their energy and entropy components. The sensitivity of the conclusions of the model to variations in parameters is computed and appears to be weak.

Correction: High pressure single-molecule FRET studies of the lysine riboswitch: cationic and osmolytic effects on pressure induced denaturation

Correction for ‘High pressure single-molecule FRET studies of the lysine riboswitch: cationic and osmolytic effects on pressure induced denaturation’ by Hsuan-Lei Sung et al., Phys. Chem. Chem. Phys., 2020, DOI: 10.1039/d0cp01921f.

Multiple prebiotic metals mediate translation

Today, Mg2+ is an essential cofactor with diverse structural and functional roles in life's oldest macromolecular machine, the translation system. We tested whether ancient Earth conditions (low O2, high Fe2+, and high Mn2+) can revert the ribosome to a functional ancestral state. First, SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension) was used to compare the effect of Mg2+, Fe2+, and Mn2+ on the tertiary structure of rRNA. Then, we used in vitro translation reactions to test whether Fe2+ or Mn2+ could mediate protein production, and quantified ribosomal metal content. We found that (i) Mg2+, Fe2+, and Mn2+ had strikingly similar effects on rRNA folding; (ii) Fe2+ and Mn2+ can replace Mg2+ as the dominant divalent cation during translation of mRNA to functional protein; and (iii) Fe and Mn associate extensively with the ribosome. Given that the translation system originated and matured when Fe2+ and Mn2+ were abundant, these findings suggest that Fe2+ and Mn2+ played a role in early ribosomal evolution.

Amino Acid Stabilization of Nucleic Acid Secondary Structure: Kinetic Insights from Single-Molecule Studies

Amino acid and nucleic acid interactions are central in biology and may have played a role in the evolutionary development of protein-based life from an early "RNA Universe." To explore the possible role of single amino acids in promoting nucleic acid folding, single-molecule Förster resonance energy transfer experiments have been implemented with a DNA hairpin construct (7 nucleotide double strand with a 40A loop) as a simple model for secondary structure formation. Exposure to positively charged amino acids (arginine and lysine) is found to clearly stabilize the secondary structure. Kinetically, each amino acid promotes folding by generating a large increase in the folding rate with little change in the unfolding rate. From analysis as a function of temperature, arginine and lysine are found to significantly increase the overall exothermicity of folding while imposing only a small entropic penalty on the folding process. Detailed investigations into the kinetics and thermodynamics of this amino acid-induced folding stability reveal arginine and lysine to interact with nucleic acids in a manner reminiscent of monovalent cations. Specifically, these observations are interpreted in the context of an ion atmosphere surrounding the nucleic acid, in which amino acid salts stabilize folding qualitatively like small monovalent cations but also exhibit differences because of the composition of their side chains.

Single-Molecule Fluorescence Reveals Commonalities and Distinctions among Natural and in Vitro-Selected RNA Tertiary Motifs in a Multistep Folding Pathway

Decades of study of the RNA folding problem have revealed that diverse and complex structured RNAs are built from a common set of recurring structural motifs, leading to the perspective that a generalizable model of RNA folding may be developed from understanding of the folding properties of individual structural motifs. We used single-molecule fluorescence to dissect the kinetic and thermodynamic properties of a set of variants of a common tertiary structural motif, the tetraloop/tetraloop-receptor (TL/TLR). Our results revealed a multistep TL/TLR folding pathway in which preorganization of the ubiquitous AA-platform submotif precedes the formation of the docking transition state and tertiary A-minor hydrogen bond interactions form after the docking transition state. Differences in ion dependences between TL/TLR variants indicated the occurrence of sequence-dependent conformational rearrangements prior to and after the formation of the docking transition state. Nevertheless, varying the junction connecting the TL/TLR produced a common kinetic and ionic effect for all variants, suggesting that the global conformational search and compaction electrostatics are energetically independent from the formation of the tertiary motif contacts. We also found that in vitro-selected variants, despite their similar stability at high Mg²⁺ concentrations, are considerably less stable than natural variants under near-physiological ionic conditions, and the occurrence of the TL/TLR sequence variants in Nature correlates with their thermodynamic stability in isolation. Overall, our findings are consistent with modular but complex energetic properties of RNA structural motifs and will aid in the eventual quantitative description of RNA folding from its secondary and tertiary structural elements.

Biophysical Insights from Temperature-Dependent Single-Molecule Förster Resonance Energy Transfer

Single-molecule fluorescence microscopy techniques can be used in combination with micrometer length-scale temperature control and Förster resonance energy transfer (FRET) in order to gain detailed information about fundamental biophysical phenomena. In particular, this combination of techniques has helped foster the development of remarkable quantitative tools for studying both time- and temperature-dependent structural kinetics of biopolymers. Over the past decade, multiple research efforts have successfully incorporated precise spatial and temporal control of temperature into single-molecule FRET (smFRET)-based experiments, which have uncovered critical thermodynamic information on a wide range of biological systems such as conformational dynamics of nucleic acids. This review provides an overview of various temperature-dependent smFRET approaches from our laboratory and others, highlighting efforts in which such methods have been successfully applied to studies of single-molecule nucleic acid folding.

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 ].

Probing the kinetic and thermodynamic consequences of the tetraloop/tetraloop receptor monovalent ion-binding site in P4-P6 RNA by smFRET

Structured RNA molecules play roles in central biological processes and understanding the basic forces and features that govern RNA folding kinetics and thermodynamics can help elucidate principles that underlie biological function. Here we investigate one such feature, the specific interaction of monovalent cations with a structured RNA, the P4-P6 domain of the Tetrahymena ribozyme. We employ single molecule FRET (smFRET) approaches as these allow determination of folding equilibrium and rate constants over a wide range of stabilities and thus allow direct comparisons without the need for extrapolation. These experiments provide additional evidence for specific binding of monovalent cations, Na+ and K+, to the RNA tetraloop-tetraloop receptor (TL-TLR) tertiary motif. These ions facilitate both folding and unfolding, consistent with an ability to help order the TLR for binding and further stabilize the tertiary contact subsequent to attainment of the folding transition state.

Sodium and Potassium Interactions with Nucleic Acids

Metal ions are essential cofactors for the structure and functions of nucleic acids. Yet, the early discovery in the 70s of the crucial role of Mg(2+) in stabilizing tRNA structures has occulted for a long time the importance of monovalent cations. Renewed interest in these ions was brought in the late 90s by the discovery of specific potassium metal ions in the core of a group I intron. Their importance in nucleic acid folding and catalytic activity is now well established. However, detection of K(+) and Na(+) ions is notoriously problematic and the question about their specificity is recurrent. Here we review the different methods that can be used to detect K(+) and Na(+) ions in nucleic acid structures such as X-ray crystallography, nuclear magnetic resonance or molecular dynamics simulations. We also discuss specific versus non-specific binding to different structures through various examples.

Kinetic and thermodynamic origins of osmolyte-influenced nucleic acid folding

The influential role of monovalent and divalent metal cations in facilitating conformational transitions in both RNA and DNA has been a target of intense biophysical research efforts. However, organic neutrally charged cosolutes can also significantly alter nucleic acid conformational transitions. For example, highly soluble small molecules such as trimethylamine N-oxide (TMAO) and urea are occasionally utilized by organisms to regulate cellular osmotic pressure. Ensemble studies have revealed that these so-called osmolytes can substantially influence the thermodynamics of nucleic acid conformational transitions. In the present work, we exploit single-molecule FRET (smFRET) techniques to measure, for first time, the kinetic origins of these osmolyte-induced changes to the folding free energy. In particular, we focus on smFRET RNA and DNA constructs designed as model systems for secondary and tertiary structure formation. These findings reveal that TMAO preferentially stabilizes both secondary and tertiary interactions by increasing kfold and decreasing kunfold, whereas urea destabilizes both conformational transitions, resulting in the exact opposite shift in kinetic rate constants (i.e., decreasing kfold and increasing kunfold). Complementary temperature-dependent smFRET experiments highlight a thermodynamic distinction between the two different mechanisms responsible for TMAO-facilitated conformational transitions, while only a single mechanism is seen for the destabilizing osmolyte urea. Finally, these results are interpreted in the context of preferential interactions between osmolytes, and the solvent accessible surface area (SASA) associated with the (i) nucleobase, (ii) sugar, and (iii) phosphate groups of nucleic acids in order to map out structural changes that occur during the conformational transitions.

Pulsed IR heating studies of single-molecule DNA duplex dissociation kinetics and thermodynamics

Single-molecule fluorescence spectroscopy is a powerful technique that makes it possible to observe the conformational dynamics associated with biomolecular processes. The addition of precise temperature control to these experiments can yield valuable thermodynamic information about equilibrium and kinetic rate constants. To accomplish this, we have developed a microscopy technique based on infrared laser overtone/combination band absorption to heat small (≈10(-11) liter) volumes of water. Detailed experimental characterization of this technique reveals three major advantages over conventional stage heating methods: 1), a larger range of steady-state temperatures (20-100°C); 2), substantially superior spatial (≤20 μm) control; and 3), substantially superior temporal (≈1 ms) control. The flexibility and breadth of this spatial and temporally resolved laser-heating approach is demonstrated in single-molecule fluorescence assays designed to probe the dissociation of a 21 bp DNA duplex. These studies are used to support a kinetic model based on nucleic acid end fraying that describes dissociation for both short (<10 bp) and long (>10 bp) DNA duplexes. These measurements have been extended to explore temperature-dependent kinetics for the 21 bp construct, which permit determination of single-molecule activation enthalpies and entropies for DNA duplex dissociation.

Molecular-crowding effects on single-molecule RNA folding/unfolding thermodynamics and kinetics

The effects of "molecular crowding" on elementary biochemical processes due to high solute concentrations are poorly understood and yet clearly essential to the folding of nucleic acids and proteins into correct, native structures. The present work presents, to our knowledge, first results on the single-molecule kinetics of solute molecular crowding, specifically focusing on GAAA tetraloop-receptor folding to isolate a single RNA tertiary interaction using time-correlated single-photon counting and confocal single-molecule FRET microscopy. The impact of crowding by high-molecular-weight polyethylene glycol on the RNA folding thermodynamics is dramatic, with up to ΔΔG° ∼ -2.5 kcal/mol changes in free energy and thus >60-fold increase in the folding equilibrium constant (Keq) for excluded volume fractions of 15%. Most importantly, time-correlated single-molecule methods permit crowding effects on the kinetics of RNA folding/unfolding to be explored for the first time (to our knowledge), which reveal that this large jump in Keq is dominated by a 35-fold increase in tetraloop-receptor folding rate, with only a modest decrease in the corresponding unfolding rate. This is further explored with temperature-dependent single-molecule RNA folding measurements, which identify that crowding effects are dominated by entropic rather than enthalpic contributions to the overall free energy change. Finally, a simple "hard-sphere" treatment of the solute excluded volume is invoked to model the observed kinetic trends, and which predict ΔΔG° ∼ -5 kcal/mol free-energy stabilization at excluded volume fractions of 30%.

Single-molecule fluorescence resonance energy transfer studies of the human telomerase RNA pseudoknot: temperature-/urea-dependent folding kinetics and thermodynamics

The ribonucleoprotein telomerase is an RNA-dependent DNA polymerase that catalyzes the repetitive addition of a short, species-specific, DNA sequence to the ends of linear eukaryotic chromosomes. The single RNA component of telomerase contains both the template sequence for DNA synthesis and a functionally critical pseudoknot motif, which can also exist as a less stable hairpin. Here we use a minimal version of the human telomerase RNA pseudoknot to study this hairpin-pseudoknot structural equilibrium using temperature-controlled single-molecule fluorescence resonance energy transfer (smFRET) experiments. The urea dependence of these experiments aids in determination of the folding kinetics and thermodynamics. The wild-type pseudoknot behavior is compared and contrasted to a mutant pseudoknot sequence implicated in a genetic disorder-dyskeratosis congenita. These findings clearly identify that this 2nt noncomplementary mutation destabilizes the folding of the wild-type pseudoknot by substantially reducing the folding rate constant (≈ 400-fold) while only nominally increasing the unfolding rate constant (≈ 5-fold). Furthermore, the urea dependence of the equilibrium and rate constants is used to develop a free energy landscape for this unimolecular equilibrium and propose details about the structure of the transition state. Finally, the urea-dependent folding experiments provide valuable physical insights into the mechanism for destabilization of RNA pseudoknots by such chemical denaturants.

Single-molecule kinetics reveal cation-promoted DNA duplex formation through ordering of single-stranded helices

In this work, the kinetics of short, fully complementary oligonucleotides are investigated at the single-molecule level. Constructs 6-9 bp in length exhibit single exponential kinetics over 2 orders of magnitude time for both forward (kon, association) and reverse (koff, dissociation) processes. Bimolecular rate constants for association are weakly sensitive to the number of basepairs in the duplex, with a 2.5-fold increase between 9 bp (k'on = 2.1(1) × 10(6) M(-1) s(-1)) and 6 bp (k'on = 5.0(1) × 10(6) M(-1) s(-1)) sequences. In sharp contrast, however, dissociation rate constants prove to be exponentially sensitive to sequence length, varying by nearly 600-fold over the same 9 bp (koff = 0.024 s(-1)) to 6 bp (koff = 14 s(-1)) range. The 8 bp sequence is explored in more detail, and the NaCl dependence of kon and koff is measured. Interestingly, kon increases by >40-fold (kon = 0.10(1) s(-1) to 4.0(4) s(-1) between [NaCl] = 25 mM and 1 M), whereas in contrast, koff decreases by fourfold (0.72(3) s(-1) to 0.17(7) s(-1)) over the same range of conditions. Thus, the equilibrium constant (Keq) increases by ≈160, largely due to changes in the association rate, kon. Finally, temperature-dependent measurements reveal that increased [NaCl] reduces the overall exothermicity (ΔΔH° > 0) of duplex formation, albeit by an amount smaller than the reduction in entropic penalty (-TΔΔS° < 0). This reduced entropic cost is attributed to a cation-facilitated preordering of the two single-stranded species, which lowers the association free-energy barrier and in turn accelerates the rate of duplex formation.

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.

Effects of salt, polyethylene glycol, and locked nucleic acids on the thermodynamic stabilities of consecutive terminal adenosine mismatches in RNA duplexes

Consecutive terminal mismatches add thermodynamic stability to RNA duplexes and occur frequently in microRNA-mRNA interactions. Accurate thermodynamic stabilities of consecutive terminal mismatches contribute to the development of specific, high-affinity siRNA therapeutics. Consecutive terminal adenosine mismatches (TAMS) are studied at different salt concentrations, with polyethylene glycol cosolutes, and with locked nucleic acid (LNA) substitutions. These measurements provide benchmarks for the application of thermodynamic predictions to different physiological or therapeutic conditions. The salt dependence for RNA duplex stability is similar for TAMS, internal loops, and Watson-Crick duplexes. A unified model for predicting the free energy of an RNA duplex with or without loops and mismatches at lower sodium concentrations is presented. The destabilizing effects of PEG 200 are larger for TAMS than internal loops or Watson-Crick duplexes, which may result from different base stacking conformations, dynamics, and water hydration. In contrast, LNA substitutions stabilize internal loops much more than TAMS. Surprisingly, the average per adenosine increase in stability for LNA substitutions in internal loops is -1.82 kcal/mol and only -0.20 kcal/mol for TAMS. The stabilities of TAMS and internal loops with LNA substitutions have similar favorable free energies. Thus, the unfavorable free energy of adenosine internal loops is largely an entropic effect. The favorable stabilities of TAMS result mainly from base stacking. The ability of RNA duplexes to form extended terminal mismatches in the absence of proteins such as argonaute and identifying the enthalpic contributions to terminal mismatch stabilities provide insight into the physical basis of microRNA-mRNA molecular recognition and specificity.

The role of counterion valence and size in GAAA tetraloop-receptor docking/undocking kinetics

For RNA to fold into compact, ordered structures, it must overcome electrostatic repulsion between negatively charged phosphate groups by counterion recruitment. A physical understanding of the counterion-assisted folding process requires addressing how cations kinetically and thermodynamically control the folding equilibrium for each tertiary interaction in a full-length RNA. In this work, single-molecule FRET (fluorescence resonance energy transfer) techniques are exploited to isolate and explore the cation-concentration-dependent kinetics for formation of a ubiquitous RNA tertiary interaction, that is, the docking/undocking of a GAAA tetraloop with its 11-nt receptor. Rate constants for docking (k(dock)) and undocking (k(undock)) are obtained as a function of cation concentration, size, and valence, specifically for the series Na(+), K(+), Mg(2+), Ca(2+), Co(NH(3))(6)(3+), and spermidine(3+). Increasing cation concentration acceleratesk(dock)dramatically but achieves only a slight decrease in k(undock). These results can be kinetically modeled using parallel cation-dependent and cation-independent docking pathways, which allows for isolation of the folding kinetics from the interaction energetics of the cations with the undocked and docked states, respectively. This analysis reveals a preferential interaction of the cations with the transition state and docked state as compared to the undocked RNA, with the ion-RNA interaction strength growing with cation valence. However, the corresponding number of cations that are taken up by the RNA upon folding decreases with charge density of the cation. The only exception to these behaviors is spermidine(3+), whose weaker influence on the docking equilibria with respect to Co(NH(3))(6)(3+) can be ascribed to steric effects preventing complete neutralization of the RNA phosphate groups.