Single-molecule studies highlight conformational heterogeneity in the early folding steps of a large ribozyme

The conformational space of RNase P RNA in solution

RNA conformational diversity has fundamental biological roles1-5, but direct visualization of its full conformational space in solution has not been possible using traditional biophysical techniques. Using solution atomic force microscopy, a deep neural network and statistical analyses, we show that the ribonuclease P (RNase P) RNA adopts heterogeneous conformations consisting of a conformationally invariant core and highly flexible peripheral structural elements that sample a broad conformational space, with amplitudes as large as 20-60 Å in a multitude of directions, with very low net energy cost. Increasing Mg2+ drives compaction and enhances enzymatic activity, probably by narrowing the conformational space. Moreover, analyses of the correlations and anticorrelations between spatial flexibility and sequence conservation suggest that the functional roles of both the structure and dynamics of key regions are embedded in the primary sequence. These findings reveal the structure-dynamics basis for the embodiment of both enzymatic precision and substrate promiscuity in the RNA component of the RNase P. Mapping the conformational space of the RNase P RNA demonstrates a new general approach to studying RNA structure and dynamics.

Determining structures of RNA conformers using AFM and deep neural networks

Much of the human genome is transcribed into RNAs1, many of which contain structural elements that are important for their function. Such RNA molecules-including those that are structured and well-folded2-are conformationally heterogeneous and flexible, which is a prerequisite for function3,4, but this limits the applicability of methods such as NMR, crystallography and cryo-electron microscopy for structure elucidation. Moreover, owing to the lack of a large RNA structure database, and no clear correlation between sequence and structure, approaches such as AlphaFold5 for protein structure prediction do not apply to RNA. Therefore, determining the structures of heterogeneous RNAs remains an unmet challenge. Here we report holistic RNA structure determination method using atomic force microscopy, unsupervised machine learning and deep neural networks (HORNET), a novel method for determining three-dimensional topological structures of RNA using atomic force microscopy images of individual molecules in solution. Owing to the high signal-to-noise ratio of atomic force microscopy, this method is ideal for capturing structures of large RNA molecules in distinct conformations. In addition to six benchmark cases, we demonstrate the utility of HORNET by determining multiple heterogeneous structures of RNase P RNA and the HIV-1 Rev response element (RRE) RNA. Thus, our method addresses one of the major challenges in determining heterogeneous structures of large and flexible RNA molecules, and contributes to the fundamental understanding of RNA structural biology.

Universal cold RNA phase transitions

RNA's diversity of structures and functions impacts all life forms since primordia. We use calorimetric force spectroscopy to investigate RNA folding landscapes in previously unexplored low-temperature conditions. We find that Watson-Crick RNA hairpins, the most basic secondary structure elements, undergo a glass-like transition below [Formula: see text]C where the heat capacity abruptly changes and the RNA folds into a diversity of misfolded structures. We hypothesize that an altered RNA biochemistry, determined by sequence-independent ribose-water interactions, outweighs sequence-dependent base pairing. The ubiquitous ribose-water interactions lead to universal RNA phase transitions below TG, such as maximum stability at [Formula: see text]C where water density is maximum, and cold denaturation at [Formula: see text]C. RNA cold biochemistry may have a profound impact on RNA function and evolution.

Watching ion-driven kinetics of ribozyme folding and misfolding caused by energetic and topological frustration one molecule at a time

Folding of ribozymes into well-defined tertiary structures usually requires divalent cations. How Mg2+ ions direct the folding kinetics has been a long-standing unsolved problem because experiments cannot detect the positions and dynamics of ions. To address this problem, we used molecular simulations to dissect the folding kinetics of the Azoarcus ribozyme by monitoring the path each molecule takes to reach the folded state. We quantitatively establish that Mg2+ binding to specific sites, coupled with counter-ion release of monovalent cations, stimulate the formation of secondary and tertiary structures, leading to diverse pathways that include direct rapid folding and trapping in misfolded structures. In some molecules, key tertiary structural elements form when Mg2+ ions bind to specific RNA sites at the earliest stages of the folding, leading to specific collapse and rapid folding. In others, the formation of non-native base pairs, whose rearrangement is needed to reach the folded state, is the rate-limiting step. Escape from energetic traps, driven by thermal fluctuations, occurs readily. In contrast, the transition to the native state from long-lived topologically trapped native-like metastable states is extremely slow. Specific collapse and formation of energetically or topologically frustrated states occur early in the assembly process.

Determining structures of individual RNA conformers using atomic force microscopy images and deep neural networks

The vast percentage of the human genome is transcribed into RNA, many of which contain various structural elements and are important for functions. RNA molecules are conformationally heterogeneous and functionally dyanmics1, even when they are structured and well-folded2, which limit the applicability of methods such as NMR, crystallography, or cryo-EM. Moreover, because of the lack of a large structure RNA database, and no clear correlation between sequence and structure, approaches like AlphaFold3 for protein structure prediction, do not apply to RNA. Therefore determining the structures of heterogeneous RNA is an unmet challenge. Here we report a novel method of determining RNA three-dimensional topological structures using deep neural networks and atomic force microscopy (AFM) images of individual RNA molecules in solution. Owing to the high signal-to-noise ratio of AFM, our method is ideal for capturing structures of individual conformationally heterogeneous RNA. We show that our method can determine 3D topological structures of any large folded RNA conformers, from ~ 200 to ~ 420 residues, the size range that most functional RNA structures or structural elements fall into. Thus our method addresses one of the major challenges in frontier RNA structural biology and may impact our fundamental understanding of RNA structure.

Advances in optical counting and imaging of micro/nano single-entity reactors for biomolecular analysis

Ultrasensitive detection of biomarkers is of paramount importance in various fields. Superior to the conventional ensemble measurement-based assays, single-entity assays, especially single-entity detection-based digital assays, not only can reach ultrahigh sensitivity, but also possess the potential to examine the heterogeneities among the individual target molecules within a population. In this review, we summarized the current biomolecular analysis methods that based on optical counting and imaging of the micro/nano-sized single entities that act as the individual reactors (e.g., micro-/nanoparticles, microemulsions, and microwells). We categorize the corresponding techniques as analog and digital single-entity assays and provide detailed information such as the design principles, the analytical performance, and their implementation in biomarker analysis in this work. We have also set critical comments on each technique from these aspects. At last, we reflect on the advantages and limitations of the optical single-entity counting and imaging methods for biomolecular assay and highlight future opportunities in this field.

Role of ergodicity, aging, and Gaussianity in resolving the origins of biomolecule subdiffusion

The internal motions of biomolecules are essential to their function. Although biological macromolecules conventionally show subdiffusive dynamics, only recently has subdiffusion been associated with non-ergodicity. These findings have stimulated new questions in biophysics and statistical mechanics. Is non-ergodic subdiffusion a general strategy shared by biomolecules? What underlying mechanisms are responsible for it? Here, we performed extensive molecular dynamics (MD) simulations to characterize the internal dynamics of six different biomolecules, ranging from single or double-stranded DNA, a single domain protein (KRAS), two globular proteins (PGK and SHP2), to an intrinsically disordered protein (SNAP-25). We found that the subdiffusive behavior of these biomolecules falls into two classes. The internal motion of the first three cases is ergodic subdiffusion and can be interpreted by fractional Brownian motion (FBM), while the latter three cases involve non-ergodic subdiffusion and can be modeled by mixed origins of continuous-time random walk (CTRW) and FBM.

Non-ergodicity of a globular protein extending beyond its functional timescale

Internal motions of folded proteins have been assumed to be ergodic, i.e., that the dynamics of a single protein molecule averaged over a very long time resembles that of an ensemble. Here, by performing single-molecule fluorescence resonance energy transfer (smFRET) experiments and molecular dynamics (MD) simulations of a multi-domain globular protein, cytoplasmic protein-tyrosine phosphatase (SHP2), we demonstrate that the functional inter-domain motion is observationally non-ergodic over the time spans 10⁻¹² to 10⁻⁷ s and 10⁻¹ to 10² s. The difference between observational non-ergodicity and simple non-convergence is discussed. In comparison, a single-strand DNA of similar size behaves ergodically with an energy landscape resembling a one-dimensional linear chain. The observed non-ergodicity results from the hierarchical connectivity of the high-dimensional energy landscape of the protein molecule. As the characteristic time for the protein to conduct its dephosphorylation function is ∼10 s, our findings suggest that, due to the non-ergodicity, individual, seemingly identical protein molecules can be dynamically and functionally different.

Internal motions of folded proteins have been assumed to be ergodic, i.e., that the dynamics of a single protein molecule averaged over a very long time resembles that of an ensemble.

Snapshotting the transient conformations and tracing the multiple pathways of single peptide folding using a solid-state nanopore

A fundamental question relating to protein folding/unfolding is the time evolution of the folding of a protein into its precisely defined native structure. The proper identification of transition conformations is essential for accurately describing the dynamic protein folding/unfolding pathways. Owing to the rapid transitions and sub-nm conformation differences involved, the acquisition of the transient conformations and dynamics of proteins is difficult due to limited instrumental resolution. Using the electrochemical confinement effect of a solid-state nanopore, we were able to snapshot the transient conformations and trace the multiple transition pathways of a single peptide inside a nanopore. By combining the results with a Markov chain model, this new single-molecule technique is applied to clarify the transition pathways of the β-hairpin peptide, which shows nonequilibrium fluctuations among several blockage current stages. This method enables the high-throughput investigation of transition pathways experimentally to access previously obscure peptide dynamics, which is significant for understanding the folding/unfolding mechanisms and misfolding of peptides or proteins.

A solid-state nanopore based method is described for resolving protein-folding-related problems via snapshotting the folding intermediates and characterizing the kinetics of a single peptide.

High-throughput determination of RNA structures

RNA performs and regulates a diverse range of cellular processes, with new functional roles being uncovered at a rapid pace. Interest is growing in how these functions are linked to RNA structures that form in the complex cellular environment. A growing suite of technologies that use advances in RNA structural probes, high-throughput sequencing and new computational approaches to interrogate RNA structure at unprecedented throughput are beginning to provide insights into RNA structures at new spatial, temporal and cellular scales.

Life under the Microscope: Single-Molecule Fluorescence Highlights the RNA World

The emergence of single-molecule (SM) fluorescence techniques has opened up a vast new toolbox for exploring the molecular basis of life. The ability to monitor individual biomolecules in real time enables complex, dynamic folding pathways to be interrogated without the averaging effect of ensemble measurements. In parallel, modern biology has been revolutionized by our emerging understanding of the many functions of RNA. In this comprehensive review, we survey SM fluorescence approaches and discuss how the application of these tools to RNA and RNA-containing macromolecular complexes in vitro has yielded significant insights into the underlying biology. Topics covered include the three-dimensional folding landscapes of a plethora of isolated RNA molecules, their assembly and interactions in RNA-protein complexes, and the relation of these properties to their biological functions. In all of these examples, the use of SM fluorescence methods has revealed critical information beyond the reach of ensemble averages.

Target Confinement in Small Reaction Volumes Using Microfluidic Technologies: A Smart Approach for Single-Entity Detection and Analysis

Over the last decades, the study of cells, nucleic acid molecules, and proteins has evolved from ensemble measurements to so-called single-entity studies. The latter offers huge benefits, not only as biological research tools to examine heterogeneities among individual entities within a population, but also as biosensing tools for medical diagnostics, which can reach the ultimate sensitivity by detecting single targets. Whereas various techniques for single-entity detection have been reported, this review focuses on microfluidic systems that physically confine single targets in small reaction volumes. We categorize these techniques as droplet-, microchamber-, and nanostructure-based and provide an overview of their implementation for studying single cells, nucleic acids, and proteins. We furthermore reflect on the advantages and limitations of these techniques and highlight future opportunities in the field.

Quantifying kinetics from time series of single-molecule Förster resonance energy transfer efficiency histograms

Single-molecule fluorescence spectroscopy is a powerful approach for probing biomolecular structure and dynamics, including protein folding. For the investigation of nonequilibrium kinetics, Förster resonance energy transfer combined with confocal multiparameter detection has proven particularly versatile, owing to the large number of observables and the broad range of accessible timescales, especially in combination with rapid microfluidic mixing. However, a comprehensive kinetic analysis of the resulting time series of transfer efficiency histograms and complementary observables can be challenging owing to the complexity of the data. Here we present and compare three different methods for the analysis of such kinetic data: singular value decomposition, multivariate curve resolution with alternating least square fitting, and model-based peak fitting, where an explicit model of both the transfer efficiency histogram of each species and the kinetic mechanism of the process is employed. While each of these methods has its merits for specific applications, we conclude that model-based peak fitting is most suitable for a quantitative analysis and comparison of kinetic mechanisms.

Single-molecule imaging of telomerase activity via linear plasmon rulers

A strategy for real-time monitoring of the extension of the telomerase primer based on plasmon rulers was demonstrated for the first time.

Triple helix conformation-specific blinking of Cy3 in DNA

We report that Cy3 undergoes triple helix conformation-specific blinking in DNA. Blinking patterns were affected by the stabilization of the Hoogsteen base-pair, suggesting that not only the presence but also the fluctuating behaviour of the triple helix can be monitored by the changes in the Cy3 blinking patterns.

The assembly dynamics of the cytolytic pore toxin ClyA

Pore-forming toxins are protein assemblies used by many organisms to disrupt the membranes of target cells. They are expressed as soluble monomers that assemble spontaneously into multimeric pores. However, owing to their complexity, the assembly processes have not been resolved in detail for any pore-forming toxin. To determine the assembly mechanism for the ring-shaped, homododecameric pore of the bacterial cytolytic toxin ClyA, we collected a diverse set of kinetic data using single-molecule spectroscopy and complementary techniques on timescales from milliseconds to hours, and from picomolar to micromolar ClyA concentrations. The entire range of experimental results can be explained quantitatively by a surprisingly simple mechanism. First, addition of the detergent n-dodecyl-β-D-maltopyranoside to the soluble monomers triggers the formation of assembly-competent toxin subunits, accompanied by the transient formation of a molten-globule-like intermediate. Then, all sterically compatible oligomers contribute to assembly, which greatly enhances the efficiency of pore formation compared with simple monomer addition.

Dark-field-based observation of single-nanoparticle dynamics on a supported lipid bilayer for in situ analysis of interacting molecules and nanoparticles

Observation of single plasmonic nanoparticles in reconstituted biological systems allows us to obtain snapshots of dynamic processes between molecules and nanoparticles with unprecedented spatiotemporal resolution and single-molecule/single-particle-level data acquisition. This Concept is intended to introduce nanoparticle-tethered supported lipid bilayer platforms that allow for the dynamic confinement of nanoparticles on a two-dimensional fluidic surface. The dark-field-based long-term, stable, real-time observation of freely diffusing plasmonic nanoparticles on a lipid bilayer enables one to extract a broad range of information about interparticle and molecular interactions throughout the entire reaction period. Herein, we highlight important developments in this context to provide ideas on how molecular interactions can be interpreted by monitoring dynamic behaviors and optical signals of laterally mobile nanoparticles.

Molecular plasmonics for nanoscale spectroscopy

Surface- and tip-enhanced Raman and LSPR spectroscopies have developed over the past 15 years as unique tools for uncovering the properties of single particles and single molecules that are unobservable in ensemble measurements. Measurements of individual events provide insight into the distribution of molecular properties that are averaged over in ensemble experiments. Raman and LSPR spectroscopy can provide detailed information on the identity of molecular species and changes in the local environment, respectively. In this review a detailed discussion is presented on single-molecule and single-particle Raman and LSPR spectroscopy focusing on the major developments in the fields and applications of the techniques.

Out-of-equilibrium biomolecular interactions monitored by picosecond fluorescence in microfluidic droplets

We developed a new experimental approach combining Time-Resolved Fluorescence (TRF) spectroscopy and Droplet Microfluidics (DμF) to investigate the relaxation dynamics of structurally heterogeneous biomolecular systems. Here DμF was used to produce with minimal material consumption an out-of-equilibrium, fluorescently labeled biomolecular complex by rapid mixing within the droplets. TRF detection was implemented with a streak camera to monitor the time evolution of the structural heterogeneity of the complex along its relaxation towards equilibrium while it propagates inside the microfluidic channel. The approach was validated by investigating the fluorescence decay kinetics of a model interacting system of bovine serum albumin and Patent Blue V. Fluorescence decay kinetics are acquired with very good signal-to-noise ratio and allow for global, multicomponent fluorescence decay analysis, evidencing heterogeneous structural relaxation over several 100 ms.

Protein-guided RNA dynamics during early ribosome assembly

The assembly of 30S ribosomes requires the precise addition of 20 proteins to the 16S ribosomal RNA. How early binding proteins change the rRNA structure so that later proteins may join the complex is poorly understood. Here we use single molecule fluorescence resonance energy transfer (smFRET) to observe real-time encounters between ribosomal protein S4 and the 16S 5′ domain RNA at an early stage of 30S assembly. Dynamic initial S4-RNA complexes pass through a stable non-native intermediate before converting to the native complex, showing that non-native structures can offer a low free energy path to protein-RNA recognition. Three-color FRET and molecular dynamics (MD) simulations reveal how S4 changes the frequency and direction of RNA helix motions, guiding a conformational switch that enforces the hierarchy of protein addition. This protein-guided dynamics offers an alternative explanation for induced fit in RNA-protein complexes.

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.

Hidden complexity in the isomerization dynamics of Holliday junctions

A plausible consequence of the rugged folding energy landscapes inherent to biomolecules is that there may be more than one functionally competent folded state. Indeed, molecule-to-molecule variations in the folding dynamics of enzymes and ribozymes have recently been identified in single-molecule experiments, but without systematic quantification or an understanding of their structural origin. Here, using concepts from glass physics and complementary clustering analysis, we provide a quantitative method to analyse single-molecule fluorescence resonance energy transfer (smFRET) data, thereby probing the isomerization dynamics of Holliday junctions, which display such heterogeneous dynamics over a long observation time (T(obs) ≈ 40 s). We show that the ergodicity of Holliday junction dynamics is effectively broken and that their conformational space is partitioned into a folding network of kinetically disconnected clusters. Theory suggests that the persistent heterogeneity of Holliday junction dynamics is a consequence of internal multiloops with varying sizes and flexibilities frozen by Mg(2+) ions. An annealing experiment using Mg(2+) pulses lends support to this idea by explicitly showing that interconversions between trajectories with different patterns can be induced.

Ultraviolet shadowing of RNA can cause significant chemical damage in seconds

Chemical purity of RNA samples is important for high-precision studies of RNA folding and catalytic behavior, but photodamage accrued during ultraviolet (UV) shadowing steps of sample preparation can reduce this purity. Here, we report the quantitation of UV-induced damage by using reverse transcription and single-nucleotide-resolution capillary electrophoresis. We found photolesions in a dozen natural and artificial RNAs; across multiple sequence contexts, dominantly at but not limited to pyrimidine doublets; and from multiple lamps recommended for UV shadowing. Irradiation time-courses revealed detectable damage within a few seconds of exposure for 254 nm lamps held at a distance of 5 to 10 cm from 0.5-mm thickness gels. Under these conditions, 200-nucleotide RNAs subjected to 20 seconds of UV shadowing incurred damage to 16-27% of molecules; and, due to a 'skin effect', the molecule-by-molecule distribution of lesions gave 4-fold higher variance than a Poisson distribution. Thicker gels, longer wavelength lamps, and shorter exposure times reduced but did not eliminate damage. These results suggest that RNA biophysical studies should report precautions taken to avoid artifactual heterogeneity from UV shadowing.

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.

Single-molecule spectroscopic study of dynamic nanoscale DNA bending behavior of HIV-1 nucleocapsid protein

We have studied the conformational dynamics associated with the nanoscale DNA bending induced by human immunodeficiency virus type 1 (HIV-1) nucleocapsid (NC) protein using single-molecule Förster resonance energy transfer (SM-FRET). To gain molecular-level insights into how the HIV-1 NC locally distorts the structures of duplexed DNA segments, the dynamics, reversibility, and sequence specificity of the DNA bending behavior of NC have been systematically studied. We have performed SM-FRET measurements on a series of duplexed DNA segments with varying sequences, lengths, and local structures in the presence of the wide-type HIV-1 NC and NC mutants lacking either the basic N-terminal domain or the zinc fingers. On the basis of the SM-FRET results, we have proposed a possible mechanism for the NC-induced DNA bending in which both NC's zinc fingers and N-terminal domain are found to play crucial roles. The SM-FRET results reported here add new mechanistic insights into the biological behaviors and functions of HIV-1 NC as a retroviral DNA-architectural protein which may play critical roles in the compaction, nuclear import, and integration of the proviral DNA during the retroviral life cycle.

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.

Single Molecule Analysis Research Tool (SMART): an integrated approach for analyzing single molecule data

Single molecule studies have expanded rapidly over the past decade and have the ability to provide an unprecedented level of understanding of biological systems. A common challenge upon introduction of novel, data-rich approaches is the management, processing, and analysis of the complex data sets that are generated. We provide a standardized approach for analyzing these data in the freely available software package SMART: Single Molecule Analysis Research Tool. SMART provides a format for organizing and easily accessing single molecule data, a general hidden Markov modeling algorithm for fitting an array of possible models specified by the user, a standardized data structure and graphical user interfaces to streamline the analysis and visualization of data. This approach guides experimental design, facilitating acquisition of the maximal information from single molecule experiments. SMART also provides a standardized format to allow dissemination of single molecule data and transparency in the analysis of reported data.

EXPLORING THE FOLDING FREE ENERGY LANDSCAPE OF SMALL RNA MOLECULES BY SINGLE-PAIR FÖRSTER RESONANCE ENERGY TRANSFER

Proteins and RNA are biological macromolecules built from linear polymers. The process by which they fold into compact, well-defined, three-dimensional architectures to perform their functional tasks is still not well understood. It can be visualized by Brownian motion of an ensemble of molecules through a rugged energy landscape in search of an energy minimum corresponding to the native state. To explore the conformational energy landscape of small RNAs, single pair Förster resonance energy transfer (spFRET) experiments on solutions as well as on surface-immobilized samples have provided new insights. In this review, we focus on our recent work on two FRET-labeled small RNAs, the Diels-Alderase (DAse) ribozyme and the human mitochondrial tRNA Lys . For both RNAs, three different conformational states can be distinguished, and the associated mean FRET efficiencies provide clues about their structural properties. The systematic variation of their free energies with the concentration of Mg 2+ counterions was analyzed quantitatively by using a thermodynamic model that separates conformational changes from Mg 2+ binding. Furthermore, time-resolved spFRET studies on immobilized DAse reveal slow interconversions between intermediate and folded states on the time scale of ~ 100 ms. The quantitative data obtained from spFRET experiments may likely assist in the further development of theories and models addressing the folding dynamics and (counterion-dependent) energetics of RNA molecules.

Selection of RNA aptamers against mouse embryonic stem cells

Embryonic stem cells (ESCs) are capable of unlimited self-renewal and differentiation into multiple cell types. Recent large-scale analyses have identified various cell surface molecules on ESCs. Some of them are considered to be beneficial markers for characterization of cellular phenotypes and/or play an essential role for regulating the differentiation state. Thus, it is desired to efficiently produce affinity reagents specific to these molecules. In this study, to develop such reagents for mouse ESCs (mESCs), we selected RNA aptamers against intact, live mESCs using several selection strategies. The initial selection provided us with several anti-mESC aptamers of distinct sequences, which unexpectedly react with the same molecule on mESCs. Then, to isolate aptamers against different surface markers on mESCs, one of the selected aptamers was used as a competitor in the subsequent selections. In addition, one of the selections further employed negative selection against differentiated mouse cells. Consequently, we successfully isolated three classes of anti-mESC aptamers that do not compete with one another. The isolated aptamers were shown to distinguish mESCs from differentiated mouse cell lines and trace the differentiation process of mESCs. These aptamers could prove useful for developing molecular probes and manipulation tools for mESCs.

Metal-ion-dependent folding of a uranyl-specific DNAzyme: insight into function from fluorescence resonance energy transfer studies

Fluorescence resonance energy transfer (FRET) has been used to study the global folding of an uranyl (UO(2)(2+))-specific 39E DNAzyme in the presence of Mg(2+), Zn(2+), Pb(2+), or UO(2)(2+). At pH 5.5 and physiological ionic strength (100 mM Na(+)), two of the three stems in this DNAzyme folded into a compact structure in the presence of Mg(2+) or Zn(2+). However, no folding occurred in the presence of Pb(2+) or UO(2)(2+); this is analogous to the "lock-and-key" catalysis mode first observed in the Pb(2+)-specific 8-17 DNAzyme. However, Mg(2+) and Zn(2+) exert different effects on the 8-17 and 39E DNAzymes. Whereas Mg(2+) or Zn(2+)-dependent folding promoted 8-17 DNAzyme activity, the 39E DNAzyme folding induced by Mg(2+) or Zn(2+) inhibited UO(2)(2+)-specific activity. Group IIA series of metal ions (Mg(2+), Ca(2+), Sr(2+)) also caused global folding of the 39E DNAzyme, for which the apparent binding affinity between these metal ions and the DNAzyme decreases as the ionic radius of the metal ions increases. Because the ionic radius of Sr(2+) (1.12 Å) is comparable to that of Pb(2+) (1.20 Å), but contrary to Pb(2+), Sr(2+) induces the DNAzyme to fold under identical conditions, ionic size alone cannot account for the unique folding behaviors induced by Pb(2+) and UO(2)(2+). Under low ionic strength (30 mM Na(+)), all four metal ions (Mg(2+), Zn(2+), Pb(2+), and UO(2)(2+)), caused 39E DNAzyme folding, suggesting that metal ions can neutralize the negative charge of DNA-backbone phosphates in addition to playing specific catalytic roles. Mg(2+) at low (<2 mM) concentration promoted UO(2)(2+)-specific activity, whereas Mg(2+) at high (>2 mM) concentration inhibited the UO(2)(2+)-specific activity. Therefore, the lock-and-key mode of DNAzymes depends on ionic strength, and the 39E DNAzyme is in the lock-and-key mode only at ionic strengths of 100 mM or greater.

Removal of covalent heterogeneity reveals simple folding behavior for P4-P6 RNA

RNA folding landscapes have been described alternately as simple and as complex. The limited diversity of RNA residues and the ability of RNA to form stable secondary structures prior to adoption of a tertiary structure would appear to simplify folding relative to proteins. Nevertheless, there is considerable evidence for long-lived misfolded RNA states, and these observations have suggested rugged energy landscapes. Recently, single molecule fluorescence resonance energy transfer (smFRET) studies have exposed heterogeneity in many RNAs, consistent with deeply furrowed rugged landscapes. We turned to an RNA of intermediate complexity, the P4-P6 domain from the Tetrahymena group I intron, to address basic questions in RNA folding. P4-P6 exhibited long-lived heterogeneity in smFRET experiments, but the inability to observe exchange in the behavior of individual molecules led us to probe whether there was a non-conformational origin to this heterogeneity. We determined that routine protocols in RNA preparation and purification, including UV shadowing and heat annealing, cause covalent modifications that alter folding behavior. By taking measures to avoid these treatments and by purifying away damaged P4-P6 molecules, we obtained a population of P4-P6 that gave near-uniform behavior in single molecule studies. Thus, the folding landscape of P4-P6 lacks multiple deep furrows that would trap different P4-P6 molecules in different conformations and contrasts with the molecular heterogeneity that has been seen in many smFRET studies of structured RNAs. The simplicity of P4-P6 allowed us to reliably determine the thermodynamic and kinetic effects of metal ions on folding and to now begin to build more detailed models for RNA folding behavior.

Geometrical Analysis of Cytochrome c Unfolding

We have developed a geometrical model to study the unfolding of iso-1 cytochrome c. The model draws on the crystallographic data reported for this protein. These data were used to calculate the distance between specific residues in the folded state, and in a sequence of extended states defined by n= 3, 5, 7, 9, 11, 13, and 15 residue units. Exact calculations carried out for each of the 103 residues in the polypeptide chain demonstrate that different regions of the chain have different unfolding histories. Regions where there is a persistence of compact structures can be identified, and this geometrical characterization is fully consistent with analyses of time-resolved fluorescence energy-transfer (TrFET) data using dansyl-derivatized cysteine side-chain probes at positions 39, 50, 66, 85, and 99. Our calculations were carried out assuming that different regions of the polypeptide chain unfold synchronously. To test this assumption, we performed lattice Monte Carlo simulations to study systematically the possible importance of asynchronicity. Our calculations show that small departures from synchronous dynamics can arise if displacements of residues in the main body of the chain are much more sluggish than near-terminal residues.

A rearrangement of the guanosine-binding site establishes an extended network of functional interactions in the Tetrahymena group I ribozyme active site

Protein enzymes appear to use extensive packing and hydrogen bonding interactions to precisely position catalytic groups within active sites. Because of their inherent backbone flexibility and limited side chain repertoire, RNA enzymes face additional challenges relative to proteins in precisely positioning substrates and catalytic groups. Here, we use the group I ribozyme to probe the existence, establishment, and functional consequences of an extended network of interactions in an RNA active site. The group I ribozyme catalyzes a site-specific attack of guanosine on an oligonucleotide substrate. We previously determined that the hydrogen bond between the exocyclic amino group of guanosine and the 2'-hydroxyl group at position A261 of the Tetrahymena group I ribozyme contributes to overall catalysis. We now use functional data, aided by double mutant cycles, to probe this hydrogen bond in the individual reaction steps of the catalytic cycle. Our results indicate that this hydrogen bond is not formed upon guanosine binding to the ribozyme but instead forms at a later stage of the catalytic cycle. Formation of this hydrogen bond is correlated with other structural rearrangements in the ribozyme's active site that are promoted by docking of the oligonucleotide substrate into the ribozyme's active site, and disruption of this interaction has deleterious consequences for the chemical transformation within the ternary complex. These results, combined with earlier results, provide insight into the nature of the multiple conformational steps used by the Tetrahymena group I ribozyme to achieve its active structure and reveal an intricate, extended network of interactions that is used to establish catalytic interactions within this RNA's active site.

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.

Multivector fluorescence analysis of the xpt guanine riboswitch aptamer domain and the conformational role of guanine

Purine riboswitches are RNA regulatory elements that control purine metabolism in response to intracellular concentrations of the purine ligands. Conformational changes of the guanine riboswitch aptamer domain induced by guanine binding lead to transcriptional regulation of genes involved in guanine biosynthesis. The guanine riboswitch aptamer domain has three RNA helices designated P1, P2, and P3. An overall model for the Mg²⁺- and guanine-dependent relative orientations and dynamics of P1, P2, and P3 has not been reported, and the conformational role of guanine under physiologically relevant conditions has not been fully elucidated. In this study, an ensemble and single-molecule fluorescence resonance energy transfer (FRET) study was performed on three orthogonally labeled variants of the xpt guanine riboswitch aptamer domain. The combined FRET data support a model in which the unfolded state of the aptamer domain has a highly dynamic P2 helix that switches rapidly between two orientations relative to nondynamic P1 and P3. At ≪1 mM Mg²⁺ (in the presence of a saturating level of guanine) or ≥1 mM Mg²⁺ (in the absence of guanine), the riboswitch starts to adopt a folded conformation in which loop−loop interactions lock P2 and P3 into place. At >5 mM Mg²⁺, further compaction occurs in which P1 more closely approaches P3. Our data help to explain the biological role of guanine as stabilizing the globally folded aptamer domain conformation at physiologically relevant Mg²⁺ concentrations (≤1 mM), whereas in the absence of guanine, much higher Mg²⁺ concentrations are required to induce this folding event.

Single-molecule Förster resonance energy transfer studies of RNA structure, dynamics and function

Single-molecule fluorescence microscopy experiments on RNA molecules brought to light the highly complex dynamics of key biological processes, including RNA folding, catalysis of ribozymes, ligand sensing of riboswitches and aptamers, and protein synthesis in the ribosome. By using highly advanced biophysical spectroscopy techniques in combination with sophisticated biochemical synthesis approaches, molecular dynamics of individual RNA molecules can be observed in real time and under physiological conditions in unprecedented detail that cannot be achieved with bulk experiments. Here, we review recent advances in RNA folding and functional studies of RNA and RNA-protein complexes addressed by using single-molecule Förster (fluorescence) resonance energy transfer (smFRET) technique.

Methods of site-specific labeling of RNA with fluorescent dyes

Single molecule fluorescence techniques offer unique insights into mechanisms of conformational changes of RNA. Knowing how to make fluorescently labeled RNA molecules and understanding potential limitations of different labeling strategies is essential for successful implementation of single molecule fluorescence techniques. This chapter offers a step by step overview of the process of obtaining RNA constructs ready for single molecule measurements. Several alternative methods are described for each step, and ways of troubleshooting the most common problems, in particular, splinted RNA ligation, are suggested.

Continuous imaging of plasmon rulers in live cells reveals early-stage caspase-3 activation at the single-molecule level

The use of plasmon coupling in metal nanoparticles has shown great potential for the optical characterization of many biological processes. Recently, we have demonstrated the use of "plasmon rulers" to observe conformational changes of single biomolecules in vitro. Plasmon rulers provide robust signals without photobleaching or blinking. Here, we show the first application of plasmon rulers to in vivo studies to observe very long trajectories of single biomolecules in live cells. We present a unique type of plasmon ruler comprised of peptide-linked gold nanoparticle satellites around a core particle, which was used as a probe to optically follow cell-signaling pathways in vivo at the single-molecule level. These "crown nanoparticle plasmon rulers" allowed us to continuously monitor trajectories of caspase-3 activity in live cells for over 2 h, providing sufficient time to observe early-stage caspase-3 activation, which was not possible by conventional ensemble analyses.

Global stabilization of rRNA structure by ribosomal proteins S4, S17, and S20

Ribosomal proteins stabilize the folded structure of the ribosomal RNA and enable the recruitment of further proteins to the complex. Quantitative hydroxyl radical footprinting was used to measure the extent to which three different primary assembly proteins, S4, S17, and S20, stabilize the three-dimensional structure of the Escherichia coli 16S 5' domain. The stability of the complexes was perturbed by varying the concentration of MgCl(2). Each protein influences the stability of the ribosomal RNA tertiary interactions beyond its immediate binding site. S4 and S17 stabilize the entire 5' domain, while S20 has a more local effect. Multistage folding of individual helices within the 5' domain shows that each protein stabilizes a different ensemble of structural intermediates that include nonnative interactions at low Mg(2+) concentration. We propose that the combined interactions of S4, S17, and S20 with different helical junctions bias the free-energy landscape toward a few RNA conformations that are competent to add the secondary assembly protein S16 in the next step of assembly.

Single-Molecule Fluorescence Studies of RNA: A Decade's Progress

Over the past decade, single-molecule fluorescence studies have elucidated the structure-function relationship of RNA molecules. The real-time observation of individual RNAs by single-molecule fluorescence has unveiled the dynamic behavior of complex RNA systems in unprecedented detail, revealing the presence of transient intermediate states and their kinetic pathways. This review provides an overview of how single-molecule fluorescence has been used to explore the dynamics of RNA folding and catalysis.

Single-molecule observations of ribosome function

Single-molecule investigations promise to greatly advance our understanding of basic and regulated ribosome functions during the process of translation. Here, recent progress towards directly imaging the elemental translation elongation steps using fluorescence resonance energy transfer (FRET)-based imaging methods is discussed, which provide striking evidence of the highly dynamic nature of the ribosome. In this view, global rates and fidelities of protein synthesis reactions may be regulated by interactions of the ribosome with mRNA, tRNA, translation factors and potentially many other cellular ligands that modify intrinsic conformational equilibria in the translating particle. Future investigations probing this model must aim to visualize translation processes from multiple structural and kinetic perspectives simultaneously, to provide direct correlations between factor binding and conformational events.