Single-molecule studies of group II intron ribozymes

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.

Single-Molecule Study of DNAzyme Reveals Its Intrinsic Conformational Dynamics

DNAzyme is a class of DNA molecules that can perform catalytic functions with high selectivity towards specific metal ions. Due to its potential applications for biosensors and medical therapeutics, DNAzyme has been extensively studied to characterize the relationships between its biochemical properties and functions. Similar to protein enzymes and ribozymes, DNAzymes have been found to undergo conformational changes in a metal-ion-dependent manner for catalysis. Despite the important role the conformation plays in the catalysis process, such structural and dynamic information might not be revealed by conventional approaches. Here, by using the single-molecule fluorescence resonance energy transfer (smFRET) technique, we were able to investigate the detailed conformational dynamics of a uranyl-specific DNAzyme 39E. We observed conformation switches of 39E to a folded state with the addition of Mg²⁺ and to an extended state with the addition of UO₂²⁺. Furthermore, 39E can switch to a more compact configuration with or without divalent metal ions. Our findings reveal that 39E can undergo conformational changes spontaneously between different configurations.

Combining Coarse-Grained Simulations and Single Molecule Analysis Reveals a Three-State Folding Model of the Guanidine-II Riboswitch

Riboswitch RNAs regulate gene expression by conformational changes induced by environmental conditions and specific ligand binding. The guanidine-II riboswitch is proposed to bind the small molecule guanidinium and to subsequently form a kissing loop interaction between the P1 and P2 hairpins. While an interaction was shown for isolated hairpins in crystallization and electron paramagnetic resonance experiments, an intrastrand kissing loop formation has not been demonstrated. Here, we report the first evidence of this interaction in cis in a ligand and Mg²⁺ dependent manner. Using single-molecule FRET spectroscopy and detailed structural information from coarse-grained simulations, we observe and characterize three interconvertible states representing an open and kissing loop conformation as well as a novel Mg²⁺ dependent state for the guanidine-II riboswitch from E. coli. The results further substantiate the proposed switching mechanism and provide detailed insight into the regulation mechanism for the guanidine-II riboswitch class. Combining single molecule experiments and coarse-grained simulations therefore provides a promising perspective in resolving the conformational changes induced by environmental conditions and to yield molecular insights into RNA regulation.

Single-molecule study of the effects of temperature, pH, and RNA base on the stepwise enzyme kinetics of 10–23 deoxyribozyme

We investigated how the stepwise enzyme kinetics of 10–23 deoxyribozyme was affected by temperature, pH, and RNA residue of the substrate at the single-molecule level. A deoxyribozyme-substrate system was employed to temporally categorize a single-turnover reaction into four distinct steps: binding, cleavage, dissociation of one of the cleaved fragments, and dissociation of the other fragment. The dwell time of each step was measured as the temperature was varied from 26 to 34 °C, to which the transition state theory was applied to obtain the enthalpy and entropy of activation for individual steps. In addition, we found that only the cleavage step was significantly affected by pH, indicating that it involves deprotonation of a single proton. We also found that different RNA residues specifically affect the cleavage step and cause the dwell time to change by as much as 5 times.

We investigated how the stepwise enzyme kinetics of 10–23 deoxyribozyme was affected by temperature, pH, and RNA residue of the substrate at the single-molecule level.

Pseudoknot length modulates the folding, conformational dynamics, and robustness of Xrn1 resistance of flaviviral xrRNAs

To understand how RNA dynamics is regulated and connected to its function, we investigate the folding, conformational dynamics and robustness of Xrn1 resistance of a set of flaviviral xrRNAs using SAXS, smFRET and in vitro enzymatic assays. Flaviviral xrRNAs form discrete ring-like 3D structures, in which the length of a conserved long-range pseudoknot (PK2) ranges from 2 bp to 7 bp. We find that xrRNAs' folding, conformational dynamics and Xrn1 resistance are strongly correlated and highly Mg²⁺-dependent, furthermore, the Mg²⁺-dependence is modulated by PK2 length variations. xrRNAs with long PK2 require less Mg²⁺ to stabilize their folding, exhibit reduced conformational dynamics and strong Xrn1 resistance even at low Mg²⁺, and tolerate mutations at key tertiary motifs at high Mg²⁺, which generally are destructive to xrRNAs with short PK2. These results demonstrate an unusual regulatory mechanism of RNA dynamics providing insights into the functions and future biomedical applications of xrRNAs.

Site-Specific Dual-Color Labeling of Long RNAs

Labeling of large RNAs with reporting entities, e.g., fluorophores, has significant impact on RNA studies in vitro and in vivo. Here, we describe a minimally invasive RNA labeling method featuring nucleotide and position selectivity, which solves the long-standing challenge of how to achieve accurate site-specific labeling of large RNAs with a least possible influence on folding and/or function. We use a custom-designed reactive DNA strand to hybridize to the RNA and transfer the alkyne group onto the targeted adenine or cytosine. Simultaneously, the 3'-terminus of RNA is converted to a dialdehyde moiety under the experimental condition applied. The incorporated functionalities at the internal and the 3'-terminal sites can then be conjugated with reporting entities via bioorthogonal chemistry. This method is particularly valuable for, but not limited to, single-molecule fluorescence applications. We demonstrate the method on an RNA construct of 275 nucleotides, the btuB riboswitch of Escherichia coli.

FRET-based dynamic structural biology: Challenges, perspectives and an appeal for open-science practices

Single-molecule FRET (smFRET) has become a mainstream technique for studying biomolecular structural dynamics. The rapid and wide adoption of smFRET experiments by an ever-increasing number of groups has generated significant progress in sample preparation, measurement procedures, data analysis, algorithms and documentation. Several labs that employ smFRET approaches have joined forces to inform the smFRET community about streamlining how to perform experiments and analyze results for obtaining quantitative information on biomolecular structure and dynamics. The recent efforts include blind tests to assess the accuracy and the precision of smFRET experiments among different labs using various procedures. These multi-lab studies have led to the development of smFRET procedures and documentation, which are important when submitting entries into the archiving system for integrative structure models, PDB-Dev. This position paper describes the current 'state of the art' from different perspectives, points to unresolved methodological issues for quantitative structural studies, provides a set of 'soft recommendations' about which an emerging consensus exists, and lists openly available resources for newcomers and seasoned practitioners. To make further progress, we strongly encourage 'open science' practices.

Role of the central junction in folding topology of the protein-free human U2-U6 snRNA complex

U2 and U6 small nuclear (sn)RNAs are the only snRNAs directly implicated in catalyzing the splicing of pre-mRNA, but assembly and rearrangement steps prior to catalysis require numerous proteins. Previous studies have shown that the protein-free U2-U6 snRNA complex adopts two conformations in equilibrium, characterized by four and three helices surrounding a central junction. The four-helix conformer is strongly favored in the in vitro protein-free state, but the three-helix conformer predominates in spliceosomes. To analyze the role of the central junction in positioning elements forming the active site, we derived three-dimensional models of the two conformations from distances measured between fluorophores at selected locations in constructs representing the protein-free human U2-U6 snRNA complex by time-resolved fluorescence resonance energy transfer. Data describing four angles in the four-helix conformer suggest tetrahedral geometry; addition of Mg²⁺ results in shortening of the distances between neighboring helices, indicating compaction of the complex around the junction. In contrast, the three-helix conformer shows a closer approach between helices bearing critical elements, but the addition of Mg²⁺ widens the distance between them; thus in neither conformer are the critical helices positioned to favor the proposed triplex interaction. The presence of Mg²⁺ also enhances the fraction of the three-helix conformer, as does incubation with the Prp19-related protein RBM22, which has been implicated in the remodeling of the U2-U6 snRNA complex to render it catalytically active. These data suggest that although the central junction assumes a significant role in orienting helices, spliceosomal proteins and Mg²⁺ facilitate formation of the catalytically active conformer.

Encapsulation of Fluorescently Labeled RNAs into Surface-Tethered Vesicles for Single-Molecule FRET Studies in TIRF Microscopy

Imaging fluorescently labeled biomolecules on a single-molecule level is a well-established technique to follow intra- and intermolecular processes in time, usually hidden in the ensemble average. The classical approach comprises surface immobilization of the molecule of interest, which increases the risk of restricting the natural behavior due to surface interactions. Encapsulation of such biomolecules into surface-tethered phospholipid vesicles enables to follow one molecule at a time, freely diffusing and without disturbing surface interactions. Further, the encapsulation allows to keep reaction partners (reactants and products) in close proximity and enables higher temperatures otherwise leading to desorption of the direct immobilized biomolecules.Here, we describe a detailed protocol for the encapsulation of a catalytically active RNA starting from surface passivation over RNA encapsulation to data evaluation of single-molecule FRET experiments in TIRF microscopy. We present an optimized procedure that preserves RNA functionality and applies to investigations of, e.g., large ribozymes and RNAs, where direct immobilization is structurally not possible.

Site-specific dual-color labeling of long RNAs for single-molecule spectroscopy

Labeling of long RNA molecules in a site-specific yet generally applicable manner is integral to many spectroscopic applications. Here we present a novel covalent labeling approach that is site-specific and scalable to long intricately folded RNAs. In this approach, a custom-designed DNA strand that hybridizes to the RNA guides a reactive group to target a preselected adenine residue. The functionalized nucleotide along with the concomitantly oxidized 3'-terminus can subsequently be conjugated to two different fluorophores via bio-orthogonal chemistry. We validate this modular labeling platform using a regulatory RNA of 275 nucleotides, the btuB riboswitch of Escherichia coli, demonstrate its general applicability by modifying a base within a duplex, and show its site-selectivity in targeting a pair of adjacent adenines. Native folding and function of the RNA is confirmed on the single-molecule level by using FRET as a sensor to visualize and characterize the conformational equilibrium of the riboswitch upon binding of its cofactor adenosylcobalamin. The presented labeling strategy overcomes size and site constraints that have hampered routine production of labeled RNA that are beyond 200 nt in length.

Optimal molecular crowding accelerates group II intron folding and maximizes catalysis

Unlike in vivo conditions, group II intron ribozymes are known to require high magnesium(II) concentrations ([Mg2+]) and high temperatures (42 °C) for folding and catalysis in vitro. A possible explanation for this difference is the highly crowded cellular environment, which can be mimicked in vitro by macromolecular crowding agents. Here, we combined bulk activity assays and single-molecule Förster Resonance Energy Transfer (smFRET) to study the influence of polyethylene glycol (PEG) on catalysis and folding of the ribozyme. Our activity studies reveal that PEG reduces the [Mg2+] required, and we found an "optimum" [PEG] that yields maximum activity. smFRET experiments show that the most compact state population, the putative active state, increases with increasing [PEG]. Dynamic transitions between folded states also increase. Therefore, this study shows that optimal molecular crowding concentrations help the ribozyme not only to reach the native fold but also to increase its in vitro activity to approach that in physiological conditions.

Reliable State Identification and State Transition Detection in Fluorescence Intensity-Based Single-Molecule Förster Resonance Energy-Transfer Data

Single-molecule Förster resonance energy transfer (smFRET) is a powerful technique to probe biomolecular structure and dynamics. A popular implementation of smFRET consists of recording fluorescence intensity time traces of surface-immobilized, chromophore-tagged molecules. This approach generates large and complex data sets, the analysis of which is to date not standardized. Here, we address a key challenge in smFRET data analysis: the generation of thermodynamic and kinetic models that describe with statistical rigor the behavior of FRET trajectories recorded from surface-tethered biomolecules in terms of the number of FRET states, the corresponding mean FRET values, and the kinetic rates at which they interconvert. For this purpose, we first perform Monte Carlo simulations to generate smFRET trajectories, in which a relevant space of experimental parameters is explored. Then, we provide an account on current strategies to achieve such model selection, as well as a quantitative assessment of their performances. Specifically, we evaluate the performance of each algorithm (change-point analysis, STaSI, HaMMy, vbFRET, and ebFRET) with respect to accuracy, reproducibility, and computing time, which yields a range of algorithm-specific referential benchmarks for various data qualities. Data simulation and analysis were performed with our MATLAB-based multifunctional analysis software for handling smFRET data (MASH-FRET).

Bulky Lesion Bypass Requires Dpo4 Binding in Distinct Conformations

Translesion DNA synthesis is an essential process that helps resume DNA replication at forks stalled near bulky adducts on the DNA. Benzo[a]pyrene (B[a]P) is a polycyclic aromatic hydrocarbon (PAH) that can be metabolically activated to benzo[a]pyrene diol epoxide (BPDE), which then can react with DNA to form carcinogenic DNA adducts. Here, we have used single-molecule florescence resonance energy transfer (smFRET) experiments, classical molecular dynamics simulations, and nucleotide incorporation assays to investigate the mechanism by which the model Y-family polymerase, Dpo4, bypasses a (+)-cis-B[a]P-N²-dG adduct in DNA. Our data show that when (+)-cis-B[a]P-N²-dG is the templating base, the B[a]P moiety is in a non-solvent exposed conformation stacked within the DNA helix, where it effectively blocks nucleotide incorporation across the adduct by Dpo4. However, when the media contains a small amount of dimethyl sulfoxide (DMSO), the adduct is able to move to a solvent-exposed conformation, which enables error-prone DNA replication past the adduct. When the primer terminates across from the adduct position, the addition of DMSO leads to the formation of an insertion complex capable of accurate nucleotide incorporation.

Group II Intron Self-Splicing

Group II introns are large, autocatalytic ribozymes that catalyze RNA splicing and retrotransposition. Splicing by group II introns plays a major role in the metabolism of plants, fungi, and yeast and contributes to genetic variation in many bacteria. Group II introns have played a major role in genome evolution, as they are likely progenitors of spliceosomal introns, retroelements, and other machinery that controls genetic variation and stability. The structure and catalytic mechanism of group II introns have recently been elucidated through a combination of genetics, chemical biology, solution biochemistry, and crystallography. These studies reveal a dynamic machine that cycles progressively through multiple conformations as it stimulates the various stages of splicing. A central active site, containing a reactive metal ion cluster, catalyzes both steps of self-splicing. These studies provide insights into RNA structure, folding, and catalysis, as they raise new questions about the behavior of RNA machines.

Pseudoknot Formation Seeds the Twister Ribozyme Cleavage Reaction Coordinate

The twister RNA is a recently discovered nucleolytic ribozyme that is present in both bacteria and eukarya. While its biological role remains unclear, crystal structure analyses and biochemical approaches have revealed critical features of its catalytic mechanism. Here, we set out to explore dynamic aspects of twister RNA folding along the cleavage reaction coordinate. To do so, we have employed both bulk and single-molecule fluorescence resonance energy transfer (FRET) methods to investigate a set of twister RNAs with labels strategically positioned at communicating segments. The data reveal that folding of the central pseudoknot (T1), the most crucial structural determinant to promote cleavage, exhibits reversible opening and closing dynamics at physiological Mg2+ concentration. Uncoupled folding, in which T1 formation precedes structuring for closing of stem P1, was confirmed using pre-steady-state three-color smFRET experiments initiated by Mg2+ injection. This finding suggests that the folding path of twister RNA requires proper orientation of the substrate prior to T1 closure such that the U5-A6 cleavage site becomes embraced to achieve its cleavage competent conformation. We also find that the cleaved 3'-fragment retains its compacted pseudoknot fold, despite the absence of the phylogenetically conserved stem P1, rationalizing the poor turnover efficiency of the twister ribozyme.

Sequence-Specific Post-Synthetic Oligonucleotide Labeling for Single-Molecule Fluorescence Applications

The sequence-specific fluorescence labeling of nucleic acids is a prerequisite for various methods including single-molecule Förster resonance energy transfer (smFRET) for the detailed study of nucleic acid folding and function. Such nucleic acid derivatives are commonly obtained by solid-phase methods; however, yields decrease rapidly with increasing length and restrict the practicability of this approach for long strands. Here, we report a new labeling strategy for the postsynthetic incorporation of a bioorthogonal group into single stranded regions of both DNA and RNA of unrestricted length. A 12-alkyne-etheno-adenine modification is sequence-selectively formed using DNA-templated synthesis, followed by conjugation of the fluorophore Cy3 via a copper-catalyzed azide-alkyne cycloaddition (CuAAC). Evaluation of the labeled strands in smFRET measurements shows that the strategy developed here has the potential to be used for the study of long functional nucleic acids by (single-molecule) fluorescence or other methods. To prove the universal use of the method, its application was successfully extended to the labeling of a short RNA single strand. As a proof-of-concept, also the labeling of a large RNA molecule in form of a 633 nucleotide long construct derived from the Saccharomyces cerevisiae group II intron Sc.ai5γ was performed, and covalent attachment of the Cy3 fluorophore was shown with gel electrophoresis.

Native purification and labeling of RNA for single molecule fluorescence studies

The recent discovery that non-coding RNAs are considerably more abundant and serve a much wider range of critical cellular functions than recognized over previous decades of research into molecular biology has sparked a renewed interest in the study of structure-function relationships of RNA. To perform their functions in the cell, RNAs must dominantly adopt their native conformations, avoiding deep, non-productive kinetic traps that may exist along a frustrated (rugged) folding free energy landscape. Intracellularly, RNAs are synthesized by RNA polymerase and fold co-transcriptionally starting from the 5' end, sometimes with the aid of protein chaperones. By contrast, in the laboratory RNAs are commonly generated by in vitro transcription or chemical synthesis, followed by purification in a manner that includes the use of high concentrations of urea, heat and UV light (for detection), resulting in the denaturation and subsequent refolding of the entire RNA. Recent studies into the nature of heterogeneous RNA populations resulting from this process have underscored the need for non-denaturing (native) purification methods that maintain the co-transcriptional fold of an RNA. Here, we present protocols for the native purification of an RNA after its in vitro transcription and for fluorophore and biotin labeling methods designed to preserve its native conformation for use in single molecule fluorescence resonance energy transfer (smFRET) inquiries into its structure and function. Finally, we present methods for taking smFRET data and for analyzing them, as well as a description of plausible overall preparation schemes for the plethora of non-coding RNAs.

Strategy for Internal Labeling of Large RNAs with Minimal Perturbation by Using Fluorescent PNA

Fluorescence techniques for the investigation of biomolecules and their folding pathways require an efficient labeling strategy. A common method to internally label large RNAs involves the introduction of long loops for hybridization of fluorophore-carrying DNA strands. Such loops often disturb the structure, and thus the functionality, of the RNA. Here we show, in a proof of concept study with a >600 nucleotide group II intron ribozyme, that the usage of the nucleic acid analogue peptide nucleic acid (PNA) is more efficient in several aspects, minimizing the required structural modifications of the RNA. We demonstrate by various methods, including smFRET, that much smaller concentrations and shorter PNAs can be applied, compared to DNA, for rapid and specific internal RNA labeling. The folding pathway and catalytic activity of this large ribozyme is only minimally affected by the PNA, but the background signal is significantly reduced.

Cation-induced kinetic heterogeneity of the intron-exon recognition in single group II introns

RNA is commonly believed to undergo a number of sequential folding steps before reaching its functional fold, i.e., the global minimum in the free energy landscape. However, there is accumulating evidence that several functional conformations are often in coexistence, corresponding to multiple (local) minima in the folding landscape. Here we use the 5'-exon-intron recognition duplex of a self-splicing ribozyme as a model system to study the influence of Mg(2+) and Ca(2+) on RNA tertiary structure formation. Bulk and single-molecule spectroscopy reveal that near-physiological M(2+) concentrations strongly promote interstrand association. Moreover, the presence of M(2+) leads to pronounced kinetic heterogeneity, suggesting the coexistence of multiple docked and undocked RNA conformations. Heterogeneity is found to decrease at saturating M(2+) concentrations. Using NMR, we locate specific Mg(2+) binding pockets and quantify their affinity toward Mg(2+). Mg(2+) pulse experiments show that M(2+) exchange occurs on the timescale of seconds. This unprecedented combination of NMR and single-molecule Förster resonance energy transfer demonstrates for the first time to our knowledge that a rugged free energy landscape coincides with incomplete occupation of specific M(2+) binding sites at near-physiological M(2+) concentrations. Unconventional kinetics in nucleic acid folding frequently encountered in single-molecule experiments are therefore likely to originate from a spectrum of conformations that differ in the occupation of M(2+) binding sites.

The role of Mg(II) in DNA cleavage site recognition in group II intron ribozymes: solution structure and metal ion binding sites of the RNA-DNA complex

Group II intron ribozymes catalyze the cleavage of (and their reinsertion into) DNA and RNA targets using a Mg2(+)-dependent reaction. The target is cleaved 3' to the last nucleotide of intron binding site 1 (IBS1), one of three regions that form base pairs with the intron's exon binding sites (EBS1 to -3).We solved the NMR solution structure of the d3' hairpin of the Sc.ai5γ intron containing EBS1 in its 11-nucleotide loop in complex with the dIBS1 DNA 7-mer and compare it with the analogous RNA-RNA contact. The EBS1-dIBS1 helix is slightly flexible and non-symmetric. NMR data reveal two major groove binding sites for divalent metal ions at the EBS1-dIBS1 helix, and surface plasmon resonance experiments show that low concentrations of Mg2(+) considerably enhance the affinity of dIBS1 for EBS1. Our results indicate that identification of both RNA and DNA IBS1 targets, presentation of the scissile bond, and stabilization of the structure by metal ions are governed by the overall structure of EBS1-dIBS1 and the surrounding loop nucleotides but are irrespective of different EBS1-(d)IBS1 geometries and interstrand affinities.

Direct observation of the growth and shrinkage of microtubules by single molecule Förster resonance energy transfer

Single molecule Förster resonance energy transfer (FRET) has been applied to monitor the growth and the shrinkage of the dynamic microtubules.

Simple Method for Constructing RNA Triangle, Square, Pentagon by Tuning Interior RNA 3WJ Angle from 60° to 90° or 108°

Precise shape control of architectures at the nanometer scale is an intriguing but extremely challenging facet. RNA has recently emerged as a unique material and thermostable building block for use in nanoparticle construction. Here, we describe a simple method from design to synthesis of RNA triangle, square, and pentagon by stretching RNA 3WJ native angle from 60° to 90° and 108°, using the three-way junction (3WJ) of the pRNA from bacteriophage phi29 dsDNA packaging motor. These methods for the construction of elegant polygons can be applied to other RNA building blocks including the utilization and application of RNA 4-way, 5-way, and other multi-way junctions.

RNA folding dynamics using laser-assisted single-molecule refolding

RNA folding pathways can be complex and even include kinetic traps or misfolded intermediates that can be slow to resolve. Characterizing these pathways is critical to understanding how RNA molecules acquire their biological function. We have previously developed a novel approach to help characterize such misfolded intermediates. Laser-assisted single-molecule refolding (LASR) is a powerful technique that combines temperature-jump (T-jump) kinetics with single-molecule detection. In a typical LASR experiment, the temperature is rapidly increased and conformational dynamics are characterized, in real-time, at the single-molecule level using single-molecule fluorescence resonance energy transfer (smFRET). Here, we provide detailed protocols for performing LASR experiments including sample preparation, temperature calibration, and data analysis.

Single-molecule fluorescence of nucleic acids

Single-molecule fluorescence studies of nucleic acids are revolutionizing our understanding of fundamental cellular processes related to DNA and RNA processing mechanisms. Detailed molecular insights into DNA repair, replication, transcription, and RNA folding and function are continuously being uncovered by using the full repertoire of single-molecule fluorescence techniques. The fundamental reason behind the stunning growth in the application of single-molecule techniques to study nucleic acid structure and dynamics is the unmatched ability of single-molecule fluorescence, and mostly single-molecule FRET, to resolve heterogeneous static and dynamic populations and identify transient and low-populated states without the need for sample synchronization. New advances in DNA and RNA synthesis, post-synthetic dye-labeling methods, immobilization and passivation strategies, improved dye photophysics, and standardized analysis methods have enabled the implementation of single-molecule techniques beyond specialized laboratories. In this chapter, we introduce the practical aspects of applying single-molecule techniques to investigate nucleic acid structure, dynamics, and function.

A chemogenetic approach to study the structural basis of protein-facilitated RNA folding

Large RNA molecules play important roles in all aspects of cellular metabolism ranging from mRNA splicing and protein biosynthesis to regulation of gene expression. In order to correctly perform its function in the cell, an RNA molecule must fold into a complex tertiary structure. Folding of many large RNAs is slow either due to formation of stable misfolded intermediates or due to high contact order or instability of obligate folding intermediates. Therefore many RNAs use protein cofactors to facilitate their folding in vivo. Folding of the yeast mitochondrial group II intron ai5γ to the native state under physiological conditions is facilitated by the protein cofactor Mss116. This chapter describes the use of Nucleotide Analog Interference Mapping (NAIM) to identify specific substructures within the intron molecule that are directly affected by the protein.

Post-transcriptional modifications modulate conformational dynamics in human U2-U6 snRNA complex

The spliceosome catalyzes precursor-mRNA splicing in all eukaryotes. It consists of over 100 proteins and five small nuclear RNAs (snRNAs), including U2 and U6 snRNAs, which are essential for catalysis. Human and yeast snRNAs share structural similarities despite the fact that human snRNAs contain numerous post-transcriptional modifications. Although functions for these modifications have been proposed, their exact roles are still not well understood. To help elucidate these roles in pre-mRNA splicing, we have used single-molecule fluorescence to characterize the effect of several post-transcriptional modifications in U2 snRNA on the conformation and dynamics of the U2-U6 complex in vitro. Consistent with yeast, the human U2-U6 complex reveals the presence of a magnesium-dependent dynamic equilibrium among three conformations. Interestingly, our data show that modifications in human U2 stem I modulate the dynamic equilibrium of the U2-U6 complex by stabilizing the four-helix structure. However, the small magnitude of this effect suggests that post-transcriptional modifications in human snRNAs may have a primary role in the mediation of specific RNA-protein interactions in vivo.

BOBA FRET: bootstrap-based analysis of single-molecule FRET data

Time-binned single-molecule Förster resonance energy transfer (smFRET) experiments with surface-tethered nucleic acids or proteins permit to follow folding and catalysis of single molecules in real-time. Due to the intrinsically low signal-to-noise ratio (SNR) in smFRET time traces, research over the past years has focused on the development of new methods to extract discrete states (conformations) from noisy data. However, limited observation time typically leads to pronounced cross-sample variability, i.e., single molecules display differences in the relative population of states and the corresponding conversion rates. Quantification of cross-sample variability is necessary to perform statistical testing in order to assess whether changes observed in response to an experimental parameter (metal ion concentration, the presence of a ligand, etc.) are significant. However, such hypothesis testing has been disregarded to date, precluding robust biological interpretation. Here, we address this problem by a bootstrap-based approach to estimate the experimental variability. Simulated time traces are presented to assess the robustness of the algorithm in conjunction with approaches commonly used in thermodynamic and kinetic analysis of time-binned smFRET data. Furthermore, a pair of functionally important sequences derived from the self-cleaving group II intron Sc.ai5γ (d3'EBS1*/IBS1*) is used as a model system. Through statistical hypothesis testing, divalent metal ions are shown to have a statistically significant effect on both thermodynamic and kinetic aspects of their interaction. The Matlab source code used for analysis (bootstrap-based analysis of smFRET data, BOBA FRET), as well as a graphical user interface, is available via http://www.aci.uzh.ch/rna/.

Crystal structure of 3WJ core revealing divalent ion-promoted thermostability and assembly of the Phi29 hexameric motor pRNA

The bacteriophage phi29 DNA packaging motor, one of the strongest biological motors characterized to date, is geared by a packaging RNA (pRNA) ring. When assembled from three RNA fragments, its three-way junction (3WJ) motif is highly thermostable, is resistant to 8 M urea, and remains associated at extremely low concentrations in vitro and in vivo. To elucidate the structural basis for its unusual stability, we solved the crystal structure of this pRNA 3WJ motif at 3.05 Å. The structure revealed two divalent metal ions that coordinate 4 nt of the RNA fragments. Single-molecule fluorescence resonance energy transfer (smFRET) analysis confirmed a structural change of 3WJ upon addition of Mg²⁺. The reported pRNA 3WJ conformation is different from a previously published construct that lacks the metal coordination sites. The phi29 DNA packaging motor contains a dodecameric connector at the vertex of the procapsid, with a central pore for DNA translocation. This portal connector serves as the foothold for pRNA binding to procapsid. Subsequent modeling of a connector/pRNA complex suggests that the pRNA of the phi29 DNA packaging motor exists as a hexameric complex serving as a sheath over the connector. The model of hexameric pRNA on the connector agrees with AFM images of the phi29 pRNA hexamer acquired in air and matches all distance parameters obtained from cross-linking, complementary modification, and chemical modification interference.

Dissecting single-molecule signal transduction in carbon nanotube circuits with protein engineering

Single-molecule experimental methods have provided new insights into biomolecular function, dynamic disorder, and transient states that are all invisible to conventional measurements. A novel, nonfluorescent single-molecule technique involves attaching single molecules to single-walled carbon nanotube field-effective transistors (SWNT FETs). These ultrasensitive electronic devices provide long-duration, label-free monitoring of biomolecules and their dynamic motions. However, generalization of the SWNT FET technique first requires design rules that can predict the success and applicability of these devices. Here, we report on the transduction mechanism linking enzymatic processivity to electrical signal generation by a SWNT FET. The interaction between SWNT FETs and the enzyme lysozyme was systematically dissected using eight different lysozyme variants synthesized by protein engineering. The data prove that effective signal generation can be accomplished using a single charged amino acid, when appropriately located, providing a foundation to widely apply SWNT FET sensitivity to other biomolecular systems.

Complex formation of cadmium with sugar residues, nucleobases, phosphates, nucleotides, and nucleic acids

Cadmium(II), commonly classified as a relatively soft metal ion, prefers indeed aromatic-nitrogen sites (e.g., N7 of purines) over oxygen sites (like sugar-hydroxyl groups). However, matters are not that simple, though it is true that the affinity of Cd(2+) towards ribose-hydroxyl groups is very small; yet, a correct orientation brought about by a suitable primary binding site and a reduced solvent polarity, as it is expected to occur in a folded nucleic acid, may facilitate metal ion-hydroxyl group binding very effectively. Cd(2+) prefers the guanine(N7) over the adenine(N7), mainly because of the steric hindrance of the (C6)NH(2) group in the adenine residue. This Cd(2+)-(N7) interaction in a guanine moiety leads to a significant acidification of the (N1)H meaning that the deprotonation reaction occurs now in the physiological pH range. N3 of the cytosine residue, together with the neighboring (C2)O, is also a remarkable Cd(2+) binding site, though replacement of (C2)O by (C2)S enhances the affinity towards Cd(2+) dramatically, giving in addition rise to the deprotonation of the (C4)NH(2) group. The phosphodiester bridge is only a weak binding site but the affinity increases further from the mono- to the di- and the triphosphate. The same also holds for the corresponding nucleotides. Complex stability of the pyrimidine-nucleotides is solely determined by the coordination tendency of the phosphate group(s), whereas in the case of purine-nucleotides macrochelate formation takes place by the interaction of the phosphate-coordinated Cd(2+) with N7. The extents of the formation degrees of these chelates are summarized and the effect of a non-bridging sulfur atom in a thiophosphate group (versus a normal phosphate group) is considered. Mixed ligand complexes containing a nucleotide and a further mono- or bidentate ligand are covered and it is concluded that in these species N7 is released from the coordination sphere of Cd(2+). In the case that the other ligand contains an aromatic residue (e.g., 2,2'-bipyridine or the indole ring of tryptophanate) intramolecular stack formation takes place. With buffers like Tris or Bistris mixed ligand complexes are formed. Cd(2+) coordination to dinucleotides and to dinucleoside monophosphates provides some insights regarding the interaction between Cd(2+) and nucleic acids. Cd(2+) binding to oligonucleotides follows the principles of coordination to its units. The available crystal studies reveal that N7 of purines is the prominent binding site followed by phosphate oxygens and other heteroatoms in nucleic acids. Due to its high thiophilicity, Cd(2+) is regularly used in so-called thiorescue experiments, which lead to the identification of a direct involvement of divalent metal ions in ribozyme catalysis.

Helicase-mediated changes in RNA structure at the single-molecule level

RNA helicases are a diverse group of RNA-dependent ATPases known to play a large number of biological roles inside the cell, such as RNA unwinding, remodeling, export and degradation. Understanding how helicases mediate changes in RNA structure is therefore of fundamental interest. The advent of single-molecule spectroscopic techniques has unveiled with unprecedented detail the interplay of RNA helicases with their substrates. In this review, we describe the characterization of helicase-RNA interactions by single-molecule approaches. State-of-the-art techniques are presented, followed by a discussion of recent advancements in this exciting field.

The structural stabilization of the κ three-way junction by Mg(II) represents the first step in the folding of a group II intron

Folding of group II introns is characterized by a first slow compaction of domain 1 (D1) followed by the rapid docking of other domains to this scaffold. D1 compaction initiates in a small subregion encompassing the κ and ζ elements. These two tertiary elements are also the major interaction sites with domain 5 to form the catalytic core. Here, we provide the first characterization of the structure adopted at an early folding step and show that the folding control element can be narrowed down to the three-way junction with the κ motif. In our nuclear magnetic resonance studies of this substructure derived from the yeast mitochondrial group II intron Sc.ai5γ, we show that a high affinity Mg(II) ion stabilizes the κ element and enables coaxial stacking between helices d′ and d′′, favoring a rigid duplex across the three-way junction. The κ-element folds into a stable GAAA-tetraloop motif and engages in A-minor interactions with helix d′. The addition of cobalt(III)hexammine reveals three distinct binding sites. The Mg(II)-promoted structural rearrangement and rigidification of the D1 core can be identified as the first micro-step of D1 folding.

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.

The brace for a growing scaffold: Mss116 protein promotes RNA folding by stabilizing an early assembly intermediate

The ai5γ group II intron requires a protein cofactor to facilitate native folding in the cell. Yeast protein Mss116 greatly accelerates intron folding under near-physiological conditions both in vivo and in vitro. Although the effect of Mss116 on the kinetics of ai5γ ribozyme folding and catalysis has been extensively studied, the precise structural role and interaction sites of Mss116 have been elusive. Using Nucleotide Analog Interference Mapping to study the folding of splicing precursor constructs, we have identified specific intron functional groups that participate in Mss116-facilitated folding and we have determined their role in the folding mechanism. The data indicate that Mss116 stabilizes an early, obligate folding intermediate within intron domain 1, thereby laying the foundation for productive folding to the native state. In addition, the data reveal an important role for the IBS2 exon sequence and for the terminus of domain 6, during the folding of self-splicing group IIB intron constructs.

Revealing -1 programmed ribosomal frameshifting mechanisms by single-molecule techniques and computational methods

Programmed ribosomal frameshifting (PRF) serves as an intrinsic translational regulation mechanism employed by some viruses to control the ratio between structural and enzymatic proteins. Most viral mRNAs which use PRF adapt an H-type pseudoknot to stimulate −1 PRF. The relationship between the thermodynamic stability and the frameshifting efficiency of pseudoknots has not been fully understood. Recently, single-molecule force spectroscopy has revealed that the frequency of −1 PRF correlates with the unwinding forces required for disrupting pseudoknots, and that some of the unwinding work dissipates irreversibly due to the torsional restraint of pseudoknots. Complementary to single-molecule techniques, computational modeling provides insights into global motions of the ribosome, whose structural transitions during frameshifting have not yet been elucidated in atomic detail. Taken together, recent advances in biophysical tools may help to develop antiviral therapies that target the ubiquitous −1 PRF mechanism among viruses.

RNA catalysis as a probe for chaperone activity of DEAD-box helicases

DEAD-box proteins are vitally important to cellular processes and make up the largest class of helicases. Many DEAD-box proteins function as RNA chaperones by accelerating structural transitions of RNA, which can result in the resolution of misfolded conformers or conversion between functional structures. While the biological importance of chaperone proteins is clear, their mechanisms are incompletely understood. Here, we illustrate how the catalytic activity of certain RNAs can be used to measure RNA chaperone activity. By measuring the amount of substrate converted to product, the fraction of catalytically active molecules is measured over time, providing a quantitative measure of the formation or loss of native RNA. The assays are described with references to group I and group II introns and their ribozyme derivatives, and examples are included that illustrate potential complications and indicate how catalytic activity measurements can be combined with physical approaches to gain insights into the mechanisms of DEAD-box proteins as RNA chaperones.

Single molecule FRET characterization of large ribozyme folding

A procedure to investigate the folding of group II intron by single molecule Fluorescence Resonance Energy Transfer (smFRET) using total internal reflection fluorescence microscopy (TIRFM) is described in this chapter. Using our previous studies on the folding and dynamics of a large ribozyme in the presence of metal ions (i.e., Mg(2+) and Ca(2+)) and/or the DEAD-box protein Mss116 as an example, we here describe step-by-step procedures to perform experiments. smFRET allows the investigation of individual molecules, thus, providing kinetic and mechanistic information hidden in ensemble averaged experiments.

Kinetic characterization of group II intron folding and splicing

Group II introns are large self-splicing ribozymes found in bacterial genomes, in organelles of plants and fungi, and even in some animal organisms. Many organellar group II introns interrupt important housekeeping genes; therefore, their splicing is critical for the survival of the host organism. Group II introns are versatile catalytic RNAs: they facilitate their own excision from a pre-mRNA, they promote ligation of exons to form a translation-competent mature mRNA; they can act like mobile genomic elements and insert themselves into RNA and DNA targets with remarkable precision, which makes them attractive tools for genetic engineering. The first step in characterization of any group II intron is the evaluation of its catalytic activity and its ability to properly fold into the native functionally active structure. This chapter describes kinetic assays used to characterize folding and catalytic properties of group II intron-derived ribozymes.

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.

Direct structural analysis of modified RNA by fluorescent in-line probing

Chemical probing is a common method for the structural characterization of RNA. Typically, RNA is radioactively end-labelled, subjected to probing conditions, and the cleavage fragment pattern is analysed by gel electrophoresis. In recent years, many chemical modifications, like fluorophores, were introduced into RNA, but methods are lacking that detect the influence of the modification on the RNA structure with single-nucleotide resolution. Here, we first demonstrate that a 5′-terminal ³²P label can be replaced by a dye label for in-line probing of riboswitch RNAs. Next, we show that small, highly structured FRET-labelled Diels–Alderase ribozymes can be directly probed, using the internal or terminal FRET dyes as reporters. The probing patterns indeed reveal whether or not the attachment of the dyes influences the structure. The existence of two dye labels in typical FRET constructs is found to be beneficial, as ‘duplexing’ allows observation of the complete RNA on a single gel. Structural information can be derived from the probing gels by deconvolution of the superimposed band patterns. Finally, we use fluorescent in-line probing to experimentally validate the structural consequences of photocaging, unambiguously demonstrating the intentional destruction of selected elements of secondary or tertiary structure.

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.

Stability and structure of mixed-ligand metal ion complexes that contain Ni2+, Cu2+, or Zn2+, and Histamine, as well as adenosine 5'-triphosphate (ATP4-) or uridine 5'-triphosphate (UTP(4-): an intricate network of equilibria

With a view on protein-nucleic acid interactions in the presence of metal ions we studied the "simple" mixed-ligand model systems containing histamine (Ha), the metal ions Ni(2+), Cu(2+), or Zn(2+) (M(2+)), and the nucleotides adenosine 5'-triphosphate (ATP(4-)) or uridine 5'-triphosphate (UTP(4-)), which will both be referred to as nucleoside 5'-triphosphate (NTP(4-)). The stability constants of the ternary M(NTP)(Ha)(2-) complexes were determined in aqueous solution by potentiometric pH titrations. We show for both ternary-complex types, M(ATP)(Ha)(2-) and M(UTP)(Ha)(2-), that intramolecular stacking between the nucleobase and the imidazole residue occurs and that the stacking intensity is approximately the same for a given M(2+) in both types of complexes: The formation degree of the intramolecular stacks is estimated to be 20 to 50%. Consequently, in protein-nucleic acid interactions imidazole-nucleobase stacks may well be of relevance. Furthermore, the well-known formation of macrochelates in binary M(2+) complexes of purine nucleotides, that is, the phosphate-coordinated M(2+) interacts with N7, is confirmed for the M(ATP)(2-) complexes. It is concluded that upon formation of the mixed-ligand complexes the M(2+)-N7 bond is broken and the energy needed for this process corresponds to the stability differences determined for the M(UTP)(Ha)(2-) and M(ATP)(Ha)(2-) complexes. It is, therefore, possible to calculate from these stability differences of the ternary complexes the formation degrees of the binary macrochelates: The closed forms amount to (65±10)%, (75±8)%, and (31±14) % for Ni(ATP)(2-), Cu(ATP)(2-), and Zn(ATP)(2-), respectively, and these percentages agree excellently with previous results obtained by different methods, confirming thus the internal validity of the data and the arguments used in the evaluation processes. Based on the overall results it is suggested that M(ATP)(2-) species, when bound to an enzyme, may exist in a closed macrochelated form only, if no enzyme groups coordinate directly to the metal ion.