Cryo-EM structures of full-length Tetrahymena ribozyme at 3.1 Å resolution

Automated detection and de novo structure modeling of nucleic acids from cryo-EM maps

Cryo-electron microscopy (cryo-EM) is one of the most powerful experimental methods for macromolecular structure determination. However, accurate DNA/RNA structure modeling from cryo-EM maps is still challenging especially for protein-DNA/RNA or multi-chain DNA/RNA complexes. Here we propose a deep learning-based method for accurate de novo structure determination of DNA/RNA from cryo-EM maps at  <5 Å resolutions, which is referred to as EM2NA. EM2NA is extensively evaluated on a diverse test set of 50 experimental maps at 2.0-5.0 Å resolutions, and compared with state-of-the-art methods including CryoREAD, ModelAngelo, and phenix.map_to_model. On average, EM2NA achieves a residue coverage of 83.15%, C4' RMSD of 1.06 Å, and sequence recall of 46.86%, which outperforms the existing methods. Moreover, EM2NA is applied to build the DNA/RNA structures with 10 to 5347 nt from an EMDB-wide data set of 263 unmodeled raw maps, demonstrating its ability in the blind model building of DNA/RNA from cryo-EM maps. EM2NA is fast and can normally build a DNA/RNA structure of  <500 nt within 10 minutes.

Dynamic conformation: Marching toward circular RNA function and application

Circular RNA is a group of covalently closed, single-stranded transcripts with unique biogenesis, stability, and conformation that play distinct roles in modulating cellular functions and also possess a great potential for developing circular RNA-based therapies. Importantly, due to its circular conformation, circular RNA generates distinct intramolecular base pairing that is different from the linear transcript. In this perspective, we review how circular RNA conformation can affect its turnover and modes of action, as well as what factors can modulate circular RNA conformation. We also discuss how understanding circular RNA conformation can facilitate learning about their functions as well as the remaining technological issues to further address their conformation. These efforts will ultimately inform the design of circular RNA-based platforms for biomedical applications.

Cryo-EM: A window into the dynamic world of RNA molecules

RNAs are critical for complex cellular functions, characterized by their structural versatility and ability to undergo conformational transitions in response to cellular cues. The elusive structures of RNAs are being unraveled with unprecedented clarity, thanks to the technological advancements in structural biology, including nuclear magnetic resonance (NMR), X-ray crystallography, cryo-electron microscopy (cryo-EM) etc. This review focuses on examining the revolutionary impact of cryo-EM on our comprehension of RNA structural dynamics, underscoring the technique's contributions to structural biology and envisioning the future trajectory of this rapidly evolving field.

Challenges, advances, and opportunities in RNA structural biology by Cryo-EM

RNAs are remarkably versatile molecules that can fold into intricate three-dimensional (3D) structures to perform diverse cellular and viral functions. Despite their biological importance, relatively few RNA 3D structures have been solved, and our understanding of RNA structure-function relationships remains in its infancy. This limitation partly arises from challenges posed by RNA's complex conformational landscape, characterized by structural flexibility, formation of multiple states, and a propensity to misfold. Recently, cryo-electron microscopy (cryo-EM) has emerged as a powerful tool for the visualization of conformationally dynamic RNA-only 3D structures. However, RNA's characteristics continue to pose challenges. We discuss experimental methods developed to overcome these hurdles, including the engineering of modular modifications that facilitate the visualization of small RNAs, improve particle alignment, and validate structural models.

Structural and biophysical dissection of RNA conformational ensembles

RNA's ability to form and interconvert between multiple secondary and tertiary structures is critical to its functional versatility and the traditional view of RNA structures as static entities has shifted towards understanding them as dynamic conformational ensembles. In this review we discuss RNA structural ensembles and their dynamics, highlighting the concept of conformational energy landscapes as a unifying framework for understanding RNA processes such as folding, misfolding, conformational changes, and complex formation. Ongoing advancements in cryo-electron microscopy and chemical probing techniques are significantly enhancing our ability to investigate multiple structures adopted by conformationally dynamic RNAs, while traditional methods such as nuclear magnetic resonance spectroscopy continue to play a crucial role in providing high-resolution, quantitative spatial and temporal information. We discuss how these methods, when used synergistically, can provide a comprehensive understanding of RNA conformational ensembles, offering new insights into their regulatory functions.

Building molecular model series from heterogeneous CryoEM structures using Gaussian mixture models and deep neural networks

Cryogenic electron microscopy (CryoEM) produces structures of macromolecules at near-atomic resolution. However, building molecular models with good stereochemical geometry from those structures can be challenging and time-consuming, especially when many structures are obtained from datasets with conformational heterogeneity. Here we present a model refinement protocol that automatically generates series of molecular models from CryoEM datasets, which describe the dynamics of the macromolecular system and have near-perfect geometry scores.

Efficient circularization of protein-encoding RNAs via a novel cis-splicing system

Circular RNAs (circRNAs) have emerged as a promising alternative to linear mRNA, owing to their unique properties and potential therapeutic applications, driving the development of novel approaches for their production. This study introduces a cis-splicing system that efficiently produces circRNAs by incorporating a ribozyme core at one end of the precursor, thereby eliminating the need for additional spacer elements between the ribozyme and the gene of interest (GOI). In this cis-splicing system, sequences resembling homologous arms at both ends of the precursor are crucial for forming the P9.0 duplex, which in turn facilitates effective self-splicing and circularization. We demonstrate that the precise recognition of the second transesterification site depends more on the structural characteristics of P9.0 adjacent to the ωG position than on the nucleotide composition of the P9.0-ωG itself. Further optimization of structural elements, like P10 and P1-ex, significantly improves circularization efficiency. The circRNAs generated through the cis-splicing system exhibit prolonged protein expression and minimal activation of the innate immune response. This study provides a comprehensive exploration of circRNA generation via a novel strategy and offers valuable insights into the structural engineering of RNA, paving the way for future advancements in circRNA-based applications.

Outcomes of the EMDataResource cryo-EM Ligand Modeling Challenge

The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein-nucleic acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: Escherichia coli beta-galactosidase with inhibitor, SARS-CoV-2 virus RNA-dependent RNA polymerase with covalently bound nucleotide analog and SARS-CoV-2 virus ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. The quality of submitted ligand models and surrounding atoms were analyzed by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics and contact scores. A composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.

Scaffold-enabled high-resolution cryo-EM structure determination of RNA

Cryo-EM structure determination of protein-free RNAs has remained difficult with most attempts yielding low to moderate resolution and lacking nucleotide-level detail. These difficulties are compounded for small RNAs as cryo-EM is inherently more difficult for lower molecular weight macromolecules. Here we present a strategy for fusing small RNAs to a group II intron that yields high resolution structures of the appended RNA, which we demonstrate with the 86-nucleotide thiamine pyrophosphate (TPP) riboswitch, and visualizing the riboswitch ligand binding pocket at 2.5 Å resolution. We also determined the structure of the ligand-free apo state and observe that the aptamer domain of the riboswitch undergoes a large-scale conformational change upon ligand binding, illustrating how small molecule binding to an RNA can induce large effects on gene expression. This study both sets a new standard for cryo-EM riboswitch visualization and offers a versatile strategy applicable to a broad range of small to moderate-sized RNAs, which were previously intractable for high-resolution cryo-EM studies.

All-atom RNA structure determination from cryo-EM maps

Many methods exist for determining protein structures from cryogenic electron microscopy maps, but this remains challenging for RNA structures. Here we developed EMRNA, a method for accurate, automated determination of full-length all-atom RNA structures from cryogenic electron microscopy maps. EMRNA integrates deep learning-based detection of nucleotides, three-dimensional backbone tracing and scoring with consideration of sequence and secondary structure information, and full-atom construction of the RNA structure. We validated EMRNA on 140 diverse RNA maps ranging from 37 to 423 nt at 2.0-6.0 Å resolutions, and compared EMRNA with auto-DRRAFTER, phenix.map_to_model and CryoREAD on a set of 71 cases. EMRNA achieves a median accuracy of 2.36 Å root mean square deviation and 0.86 TM-score for full-length RNA structures, compared with 6.66 Å and 0.58 for auto-DRRAFTER. EMRNA also obtains a high residue coverage and sequence match of 93.30% and 95.30% in the built models, compared with 58.20% and 42.20% for phenix.map_to_model and 56.45% and 52.3% for CryoREAD. EMRNA is fast and can build an RNA structure of 100 nt within 3 min.

Outcomes of the EMDataResource Cryo-EM Ligand Modeling Challenge

The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein/nucleic-acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: E. coli beta-galactosidase with inhibitor, SARS-CoV-2 RNA-dependent RNA polymerase with covalently bound nucleotide analog, and SARS-CoV-2 ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. We found that (1) the quality of submitted ligand models and surrounding atoms varied, as judged by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics, and contact scores, and (2) a composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.

Tracking Native Tetrahymena Ribozyme Folding with Fluorescence

Folding of the Tetrahymena group I intron ribozyme and other structured RNAs has been measured using a catalytic activity assay to monitor the native state formation by cleavage of a radiolabeled oligonucleotide substrate. While highly effective, the assay has inherent limitations present in any radioactivity- and gel-based assay. Administrative and safety considerations arise from the radioisotope, and data collection is laborious due to the use of polyacrylamide gels. Here we describe a fluorescence-based, solution assay that allows for more efficient data acquisition. The substrate is labeled with 6-carboxyfluorescein (6FAM) fluorophore and black hole quencher (BHQ1) at the 5' and 3' ends, respectively. Substrate cleavage results in release of the quencher, increasing the fluorescence signal by an average of 30-fold. A side-by-side comparison with the radioactivity-based assay shows good agreement in monitoring Tetrahymena ribozyme folding from a misfolded conformation to the native state, albeit with increased uncertainty. The lower precision of the fluorescence assay is compensated for by the relative ease and efficiency of the workflow. In addition, this assay will allow institutions that do not use radioactive materials to monitor native folding of the Tetrahymena ribozyme, and the same strategy should be amenable to native folding of other ribozymes.

A generalizable scaffold-based approach for structure determination of RNAs by cryo-EM

Single-particle cryo-electron microscopy (cryo-EM) can reveal the structures of large and often dynamic molecules, but smaller biomolecules (≤50 kDa) remain challenging targets due to their intrinsic low signal to noise ratio. Methods to help resolve small proteins have been applied but development of similar approaches to aid in structural determination of small, structured RNA elements have lagged. Here, we present a scaffold-based approach that we used to recover maps of sub-25 kDa RNA domains to 4.5-5.0 Å. While lacking the detail of true high-resolution maps, these maps are suitable for model building and preliminary structure determination. We demonstrate this method helped faithfully recover the structure of several RNA elements of known structure, and that it promises to be generalized to other RNAs without disturbing their native fold. This approach may streamline the sample preparation process and reduce the optimization required for data collection. This first-generation scaffold approach provides a robust system to aid in RNA structure determination by cryo-EM and lays the groundwork for further scaffold optimization to achieve higher resolution.

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.

Impact of Ion-Mixing Entropy on Orientational Preferences of DNA Helices: FRET Measurements and Computer Simulations

Biological processes require DNA and RNA helices to pack together in specific interhelical orientations. While electrostatic repulsion between backbone charges is expected to be maximized when helices are in parallel alignment, such orientations are commonplace in nature. To better understand how the repulsion is overcome, we used experimental and computational approaches to investigate how the orientational preferences of DNA helices depend on the concentration and valence of mobile cations. We used Förster resonance energy transfer (FRET) to probe the relative orientations of two 24-bp helices held together via a freely rotating PEG linker. At low cation concentrations, the helices preferred more "cross"-like orientations over those closer to parallel, and this preference was reduced with increasing salt concentrations. The results were in good quantitative agreement with Poisson-Boltzmann (PB) calculations for monovalent salt (Na+). However, PB underestimated the ability of mixtures of monovalent and divalent ions (Mg2+) to reduce the conformational preference. As a complementary approach, we performed all-atom molecular dynamics (MD) simulations and found better agreement with the experimental results. While MD and PB predict similar electrostatic forces, MD predicts a greater accumulation of Mg2+ in the ion atmosphere surrounding the DNA. Mg2+ occupancy is predicted to be greater in conformations close to the parallel orientation than in conformations close to the crossed orientation, enabling a greater release of Na+ ions and providing an entropic gain (one bound ion for two released). MD predicts an entropy gain larger than that of PB because of the increased Mg2+ occupancy. The entropy changes have a negligible effect at low Mg2+ concentrations because the free energies are dominated by electrostatic repulsion. However, as the Mg2+ concentration increases, charge screening is more effective and the mixing entropy produces readily detectable changes in packing preferences. Our results underline the importance of mixing entropy of counterions in nucleic acid interactions and provide a new understanding on the impact of a mixed ion atmosphere on the packing of DNA helices.

Capturing heterogeneous conformers of cobalamin riboswitch by cryo-EM

RNA conformational heterogeneity often hampers its high-resolution structure determination, especially for large and flexible RNAs devoid of stabilizing proteins or ligands. The adenosylcobalamin riboswitch exhibits heterogeneous conformations under 1 mM Mg2+ concentration and ligand binding reduces conformational flexibility. Among all conformers, we determined one apo (5.3 Å) and four holo cryo-electron microscopy structures (overall 3.0-3.5 Å, binding pocket 2.9-3.2 Å). The holo dimers exhibit global motions of helical twisting and bending around the dimer interface. A backbone comparison of the apo and holo states reveals a large structural difference in the P6 extension position. The central strand of the binding pocket, junction 6/3, changes from an 'S'- to a 'U'-shaped conformation to accommodate ligand. Furthermore, the binding pocket can partially form under 1 mM Mg2+ and fully form under 10 mM Mg2+ within the bound-like structure in the absence of ligand. Our results not only demonstrate the stabilizing ligand-induced conformational changes in and around the binding pocket but may also provide further insight into the role of the P6 extension in ligand binding and selectivity.

In Vitro Selection of M2+-Independent, Fast-Responding Acidic Deoxyribozymes for Bacterial Detection

We report on the first efforts to isolate acidic RNA-cleaving DNAzymes (aRCDs) from a random-sequence DNA pool by in vitro selection that are activated by a microbe Escherichia coli (E. coli), at pH 5.3. Importantly, these E. coli-responsive aRCDs only require monovalent metal ions as cofactors for cleaving a fluorogenic chimeric DNA/RNA substrate. Such characteristics can be used to efficiently protect RCDs from both intrinsic chemical instability and external enzymatic degradation. One remarkable DNAzyme, aRCD-EC1, is specific for E. coli, and its target is likely a protein. Furthermore, truncated aRCD-EC1 had significantly improved catalytic activity with an observed rate constant (kobs) of 1.18 min-1, making it the fastest bacteria-responding RCD reported to date. Clinical evaluation of this aRCD-based fluorescent assay using 40 patient urine samples demonstrated a diagnostic sensitivity of 100% and a specificity of 100% at a total analysis time of 50 min without a bacterial culture. This work can expand the repertoire of DNAzymes that are active under nonphysiological conditions, thus facilitating the development of diverse DNAzyme-based biosensors in clinical diagnosis.

RNA structure determination: From 2D to 3D

RNA molecules serve a wide range of functions that are closely linked to their structures. The basic structural units of RNA consist of single- and double-stranded regions. In order to carry out advanced functions such as catalysis and ligand binding, certain types of RNAs can adopt higher-order structures. The analysis of RNA structures has progressed alongside advancements in structural biology techniques, but it comes with its own set of challenges and corresponding solutions. In this review, we will discuss recent advances in RNA structure analysis techniques, including structural probing methods, X-ray crystallography, nuclear magnetic resonance, cryo-electron microscopy, and small-angle X-ray scattering. Often, a combination of multiple techniques is employed for the integrated analysis of RNA structures. We also survey important RNA structures that have been recently determined using various techniques.

Advances in chaperone-assisted RNA crystallography using synthetic antibodies

RNA molecules play essential roles in many biological functions, from gene expression regulation, cellular growth, and metabolism to catalysis. They frequently fold into three-dimensional structures to perform their functions. Therefore, determining RNA structure represents a key step for understanding the structure-function relationships and developing RNA-targeted therapeutics. X-ray crystallography remains a method of choice for determining high-resolution RNA structures, but it has been challenging due to difficulties associated with RNA crystallization and phasing. Several natural and synthetic RNA binding proteins have been used to facilitate RNA crystallography. Having unique properties to help crystal packing and phasing, synthetic antibody fragments, specifically the Fabs, have emerged as promising RNA crystallization chaperones, and so far, over a dozen of RNA structures have been solved using this strategy. Nevertheless, multiple steps in this approach need to be improved, including the recombinant expression of these anti-RNA Fabs, to warrant the full potential of these synthetic Fabs as RNA crystallization chaperones. This review highlights the nuts and bolts and recent advances in the chaperone-assisted RNA crystallography approach, specifically emphasizing the Fab antibody fragments as RNA crystallization chaperones.

A precise and efficient circular RNA synthesis system based on a ribozyme derived from Tetrahymena thermophila

Classic strategies for circular RNA (circRNA) preparation always introduce large numbers of linear transcripts or extra nucleotides to the circularized product. In this study, we aimed to develop an efficient system for circRNA preparation based on a self-splicing ribozyme derived from an optimized Tetrahymena thermophila group Ⅰ intron. The target RNA sequence was inserted downstream of the ribozyme and a complementary antisense region was added upstream of the ribozyme to assist cyclization. Then, we compared the circularization efficiency of ribozyme or flanking intronic complementary sequence (ICS)-mediated methods through the DNMT1, CDR1as, FOXO3, and HIPK3 genes and found that the efficiency of our system was remarkably higher than that of flanking ICS-mediated method. Consequently, the circularized products mediated by ribozyme are not introduced with additional nucleotides. Meanwhile, the overexpressed circFOXO3 maintained its biological functions in regulating cell proliferation, migration, and apoptosis. Finally, a ribozyme-based circular mRNA expression system was demonstrated with a split green fluorescent protein (GFP) using an optimized Coxsackievirus B3 (CVB3) internal ribosome entry site (IRES) sequence, and this system achieved successful translation of circularized mRNA. Therefore, this novel, convenient, and rapid engineering RNA circularization system can be applied for the functional study and large-scale preparation of circular RNA in the future.

Making RNA: Using T7 RNA polymerase to produce high yields of RNA from DNA templates

RNA is playing an ever-growing role in molecular biology and biomedicine due to the many ways it influences gene expression and its increasing use in modern therapeutics. Hence, production of RNA molecules in large quantity and high purity has become essential for advancing basic scientific research and for developing next-generation therapeutics. T7 RNA polymerase (RNAP) is a DNA-dependent RNA polymerase of bacteriophage origin and it is the most widely-utilized tool enzyme for producing RNA. Here we describe a set of robust methods for in vitro transcribing RNA molecules from DNA templates using T7 RNAP, along with a set of subsequent RNA purification schemes. In the first part of this chapter, we provide the general method for T7 RNAP-based in vitro transcription and technical notes for troubleshooting failed or inefficient transcription. We also provide modified protocols for preparing specialized RNA transcripts. In the second part, we provide two purification methods using either gel-based denaturing purification or size exclusion column-based non-denaturing purification for isolating high-purity RNA products from transcription reaction mixtures and preparing them for downstream applications. This chapter is designed to provide researchers with versatile ways to efficiently generate RNA molecules of interest and a troubleshooting guide should they encounter problems while working with in vitro transcription using T7 RNAP.

A Generalizable Scaffold-Based Approach for Structure Determination of RNAs by Cryo-EM

Single-particle cryo-electron microscopy (cryo-EM) can reveal the structures of large and often dynamic molecules, but smaller biomolecules remain challenging targets due to their intrinsic low signal to noise ratio. Methods to resolve small proteins have been applied but development of similar approaches for small structured RNA elements have lagged. Here, we present a scaffold-based approach that we used to recover maps of sub-25 kDa RNA domains to 4.5 - 5.0 Å. While lacking the detail of true high-resolution maps, these are suitable for model building and preliminary structure determination. We demonstrate this method faithfully recovers the structure of several RNA elements of known structure, and it promises to be generalized to other RNAs without disturbing their native fold. This approach may streamline the sample preparation process and reduce the optimization required for data collection. This first-generation scaffold approach provides a system for RNA structure determination by cryo-EM and lays the groundwork for further scaffold optimization to achieve higher resolution.

Structure, folding and flexibility of co-transcriptional RNA origami

RNA origami is a method for designing RNA nanostructures that can self-assemble through co-transcriptional folding with applications in nanomedicine and synthetic biology. However, to advance the method further, an improved understanding of RNA structural properties and folding principles is required. Here we use cryogenic electron microscopy to study RNA origami sheets and bundles at sub-nanometre resolution revealing structural parameters of kissing-loop and crossover motifs, which are used to improve designs. In RNA bundle designs, we discover a kinetic folding trap that forms during folding and is only released after 10 h. Exploration of the conformational landscape of several RNA designs reveal the flexibility of helices and structural motifs. Finally, sheets and bundles are combined to construct a multidomain satellite shape, which is characterized by individual-particle cryo-electron tomography to reveal the domain flexibility. Together, the study provides a structural basis for future improvements to the design cycle of genetically encoded RNA nanodevices.

RNA ligase ribozymes with a small catalytic core

Catalytic RNAs, or ribozymes, catalyze diverse chemical reactions that could have sustained primordial life in the hypothetical RNA world. Many natural ribozymes and laboratory evolved ribozymes exhibit efficient catalysis mediated by elaborate catalytic cores within complex tertiary structures. However, such complex RNA structures and sequences are unlikely to have emerged by chance during the earliest phase of chemical evolution. Here, we explored simple and small ribozyme motifs capable of ligating two RNA fragments in a template-directed fashion (ligase ribozymes). One-round selection of small ligase ribozymes followed by deep sequencing revealed a ligase ribozyme motif comprising a three-nucleotide loop opposite to the ligation junction. The observed ligation was magnesium(II) dependent and appears to form a 2'-5' phosphodiester linkage. The fact that such a small RNA motif can function as a catalyst supports a scenario in which RNA or other primordial nucleic acids played a central role in chemical evolution of life.

Pairwise Engineering of Tandemly Aligned Self-Splicing Group I Introns for Analysis and Control of Their Alternative Splicing

Alternative splicing is an important mechanism in the process of eukaryotic nuclear mRNA precursors producing multiple protein products from a single gene. Although group I self-splicing introns usually perform regular splicing, limited examples of alternative splicing have also been reported. The exon-skipping type of splicing has been observed in genes containing two group I introns. To characterize splicing patterns (exon-skipping/exon-inclusion) of tandemly aligned group I introns, we constructed a reporter gene containing two Tetrahymena introns flanking a short exon. To control splicing patterns, we engineered the two introns in a pairwise manner to design pairs of introns that selectively perform either exon-skipping or exon-inclusion splicing. Through pairwise engineering and biochemical characterization, the structural elements important for the induction of exon-skipping splicing were elucidated.

Ribozyme-mediated RNA synthesis and replication in a model Hadean microenvironment

Enzyme-catalyzed replication of nucleic acid sequences is a prerequisite for the survival and evolution of biological entities. Before the advent of protein synthesis, genetic information was most likely stored in and replicated by RNA. However, experimental systems for sustained RNA-dependent RNA-replication are difficult to realise, in part due to the high thermodynamic stability of duplex products and the low chemical stability of catalytic RNAs. Using a derivative of a group I intron as a model for an RNA replicase, we show that heated air-water interfaces that are exposed to a plausible CO₂-rich atmosphere enable sense and antisense RNA replication as well as template-dependent synthesis and catalysis of a functional ribozyme in a one-pot reaction. Both reactions are driven by autonomous oscillations in salt concentrations and pH, resulting from precipitation of acidified dew droplets, which transiently destabilise RNA duplexes. Our results suggest that an abundant Hadean microenvironment may have promoted both replication and synthesis of functional RNAs.

Snapshots of the second-step self-splicing of Tetrahymena ribozyme revealed by cryo-EM

Group I introns are catalytic RNAs that coordinate two consecutive transesterification reactions for self-splicing. To understand how the group I intron promotes catalysis and coordinates self-splicing reactions, we determine the structures of L-16 Tetrahymena ribozyme in complex with a 5'-splice site analog product and a 3'-splice site analog substrate using cryo-EM. We solve six conformations from a single specimen, corresponding to different splicing intermediates after the first ester-transfer reaction. The structures reveal dynamics during self-splicing, including large conformational changes of the internal guide sequence and the J5/4 junction as well as subtle rearrangements of active-site metals and the hydrogen bond formed between the 2'-OH group of A261 and the N2 group of guanosine substrate. These results help complete a detailed structural and mechanistic view of this paradigmatic group I intron undergoing the second step of self-splicing.

Cryo-EM-guided engineering of T-box-tRNA modules with enhanced selectivity and sensitivity in translational regulation

Riboswitches are non-coding RNA elements that play vital roles in regulating gene expression. Their specific ligand-dependent structural reorganization facilitates their use as templates for design of engineered RNA switches for therapeutics, nanotechnology and synthetic biology. T-box riboswitches bind tRNAs to sense aminoacylation and control gene expression via transcription attenuation or translation inhibition. Here we determine the cryo-EM structure of the wild-type Mycobacterium smegmatis ileS T-box in complex with its cognate tRNA ^Ile . This structure shows a very flexible antisequestrator region that tolerates both 3'-OH and 2',3'-cyclic phosphate modification at the 3' end of tRNA ^Ile . Elongation of one helical turn (11-base pair) in both the tRNA acceptor arm and T-box Stem III maintains T-box-tRNA complex formation and increases the selectivity for tRNA 3' end modification. Moreover, elongation of Stem III results in ∼6-fold tighter binding to tRNA, which leads to increased sensitivity of downstream translational regulation indicated by precedent translation. Our results demonstrate that cryo-EM can guide RNA engineering to design improved riboswitch modules for translational regulation, and potentially a variety of additional functions.

A split ribozyme that links detection of a native RNA to orthogonal protein outputs

Individual RNA remains a challenging signal to synthetically transduce into different types of cellular information. Here, we describe Ribozyme-ENabled Detection of RNA (RENDR), a plug-and-play strategy that uses cellular transcripts to template the assembly of split ribozymes, triggering splicing reactions that generate orthogonal protein outputs. To identify split ribozymes that require templating for splicing, we use laboratory evolution to evaluate the activities of different split variants of the Tetrahymena thermophila ribozyme. The best design delivers a 93-fold dynamic range of splicing with RENDR controlling fluorescent protein production in response to an RNA input. We further resolve a thermodynamic model to guide RENDR design, show how input signals can be transduced into diverse outputs, demonstrate portability across different bacteria, and use RENDR to detect antibiotic-resistant bacteria. This work shows how transcriptional signals can be monitored in situ and converted into different types of biochemical information using RNA synthetic biology.

Snapshots of the first-step self-splicing of Tetrahymena ribozyme revealed by cryo-EM

Tetrahymena ribozyme is a group I intron, whose self-splicing is the result of two sequential ester-transfer reactions. To understand how it facilitates catalysis in the first self-splicing reaction, we used cryogenic electron microscopy (cryo-EM) to resolve the structures of L-16 Tetrahymena ribozyme complexed with a 11-nucleotide 5'-splice site analog substrate. Four conformations were achieved to 4.14, 3.18, 3.09 and 2.98 Å resolutions, respectively, corresponding to different splicing intermediates during the first enzymatic reaction. Comparison of these structures reveals structural alterations, including large conformational changes in IGS/IGSext (P1-P1ext duplex) and J5/4, as well as subtle local rearrangements in the G-binding site. These structural changes are required for the enzymatic activity of the Tetrahymena ribozyme. Our study demonstrates the ability of cryo-EM to capture dynamic RNA structural changes, ushering in a new era in the analysis of RNA structure-function by cryo-EM.

Auto-DRRAFTER: Automated RNA Modeling Based on Cryo-EM Density

RNA three-dimensional structures provide rich and vital information for understanding their functions. Recent advances in cryogenic electron microscopy (cryo-EM) allow structure determination of RNAs and ribonucleoprotein (RNP) complexes. However, limited global and local resolutions of RNA cryo-EM maps pose great challenges in tracing RNA coordinates. The Rosetta-based "auto-DRRAFTER" method builds RNA models into moderate-resolution RNA cryo-EM density as part of the Ribosolve pipeline. Here, we describe a step-by-step protocol for auto-DRRAFTER using a glycine riboswitch from Fusobacterium nucleatum as an example. Successful implementation of this protocol allows automated RNA modeling into RNA cryo-EM density, accelerating our understanding of RNA structure-function relationships. Input and output files are being made available at https://github.com/auto-DRRAFTER/springer-chapter .

Near-Atomic Resolution Cryo-EM Image Reconstruction of RNA

The rapid development of cryogenic electron microscopy (cryo-EM) enables the structure determination of macromolecules without the need for crystallization. Protein, protein-lipid, and protein-nucleic acid complexes can now be routinely resolved by cryo-EM single-particle analysis (SPA) to near-atomic or atomic resolution. Here we describe the structure determination of pure RNAs by SPA, from cryo-specimen preparation to data collection and 3D reconstruction. This protocol is useful to yield many cryo-EM structures of RNA, here exemplified by the Tetrahymena L-21 ScaI ribozyme at 3.1-Å resolution.

Brief considerations on targeting RNA with small molecules

For more than three decades, RNA has been known to be a relevant and attractive macromolecule to target but figuring out which RNA should be targeted and how remains challenging. Recent years have seen the confluence of approaches for screening, drug optimization, and target validation that have led to the approval of a few RNA-targeting therapeutics for clinical applications. This focused perspective aims to highlight - but not exhaustively review - key factors accounting for these successes while pointing at crucial aspects worth considering for further breakthroughs.

Moderate activity of RNA chaperone maximizes the yield of self-spliced pre-RNA in vivo

CYT-19 is a DEAD-box protein whose adenosine-triphosphate (ATP)-dependent helicase activity facilitates the folding of group I introns in precursor RNA (pre-RNA) of Neurospora crassa (N. crassa). In the process, they consume a substantial amount of ATP. While much of the mechanistic insight into CYT-19 activity has been gained through the studies on the folding of Tetrahymena group I intron ribozyme, the more biologically relevant issue, namely the effect of CYT-19 on the self-splicing of pre-RNA, remains largely unexplored. Here, we employ a kinetic network model, based on the generalized iterative annealing mechanism (IAM), to investigate the relation between CYT-19 activity, rate of ribozyme folding, and the kinetics of the self-splicing reaction. The network rate parameters are extracted by analyzing the recent biochemical data for CYT-19-facilitated folding of Tetrahymena ribozyme. We then build extended models to explore the metabolism of pre-RNA. We show that the timescales of chaperone-mediated folding of group I ribozyme and self-splicing reaction compete with each other. As a consequence, in order to maximize the self-splicing yield of group I introns in pre-RNA, the chaperone activity must be sufficiently large to unfold the misfolded structures, but not too large to unfold the native structures prior to the self-splicing event. We discover that despite the promiscuous action on structured RNAs, the helicase activity of CYT-19 on group I ribozyme gives rise to self-splicing yields that are close to the maximum.

Emerging technologies for genetic control systems in cellular therapies

Progress in synthetic biology has enabled the construction of designer cells that sense biological inputs, and, in response, activate user-defined biomolecular programs. Such engineered cells provide unique opportunities for treating a wide variety of diseases. Current strategies mostly rely on cell-surface receptor systems engineered to convert binding interactions into activation of a transcriptional program. Genetic control systems are emerging as an appealing alternative to receptor-based sensors as they overcome the need for receptor engineering and result in cellular behaviors that operate over therapeutically relevant timescales. Genetic control systems include synthetic gene networks, RNA-based sensors, and post-translational tools. These technologies present fundamental challenges, including the requirement for precise integration with innate pathways, the need for parts orthogonal to existing circuitries, and the metabolic burden induced by such complex cell engineering endeavors. This review discusses the challenges in the design of genetic control systems for cellular therapies and their translational applications.

Mobile group I introns at nuclear rDNA position L2066 harbor sense and antisense homing endonuclease genes intervened by spliceosomal introns

Background

Mobile group I introns encode homing endonucleases that confer intron mobility initiated by a double-strand break in the intron-lacking allele at the site of insertion. Nuclear ribosomal DNA of some fungi and protists contain mobile group I introns harboring His-Cys homing endonuclease genes (HEGs). An intriguing question is how protein-coding genes embedded in nuclear ribosomal DNA become expressed. To address this gap of knowledge we analyzed nuclear L2066 group I introns from myxomycetes and ascomycetes.

Results

A total of 34 introns were investigated, including two identified mobile-type introns in myxomycetes with HEGs oriented in sense or antisense directions. Intriguingly, both HEGs are interrupted by spliceosomal introns. The intron in Didymium squamulosum, which harbors an antisense oriented HEG, was investigated in more detail. The group I intron RNA self-splices in vitro, thus generating ligated exons and full-length intron circles. The intron HEG is expressed in vivo in Didymium cells, which involves removal of a 47-nt spliceosomal intron (I-47) and 3′ polyadenylation of the mRNA. The D. squamulosum HEG (lacking the I-47 intron) was over-expressed in E. coli, and the corresponding protein was purified and shown to confer endonuclease activity. The homing endonuclease was shown to cleave an intron-lacking DNA and to produce a pentanucleotide 3′ overhang at the intron insertion site.

Conclusions

The L2066 family of nuclear group I introns all belong to the group IE subclass. The D. squamulosum L2066 intron contains major hallmarks of a true mobile group I intron by encoding a His-Cys homing endonuclease that generates a double-strand break at the DNA insertion site. We propose a potential model to explain how an antisense HEG becomes expressed from a nuclear ribosomal DNA locus.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-022-00280-4.

Structure of the OMEGA nickase IsrB in complex with ωRNA and target DNA

RNA-guided systems, such as CRISPR-Cas, combine programmable substrate recognition with enzymatic function, a combination that has been used advantageously to develop powerful molecular technologies1,2. Structural studies of these systems have illuminated how the RNA and protein jointly recognize and cleave their substrates, guiding rational engineering for further technology development3. Recent work identified a new class of RNA-guided systems, termed OMEGA, which include IscB, the likely ancestor of Cas9, and the nickase IsrB, a homologue of IscB lacking the HNH nuclease domain4. IsrB consists of only around 350 amino acids, but its small size is counterbalanced by a relatively large RNA guide (roughly 300-nt ωRNA). Here, we report the cryogenic-electron microscopy structure of Desulfovirgula thermocuniculi IsrB (DtIsrB) in complex with its cognate ωRNA and a target DNA. We find the overall structure of the IsrB protein shares a common scaffold with Cas9. In contrast to Cas9, however, which uses a recognition (REC) lobe to facilitate target selection, IsrB relies on its ωRNA, part of which forms an intricate ternary structure positioned analogously to REC. Structural analyses of IsrB and its ωRNA as well as comparisons to other RNA-guided systems highlight the functional interplay between protein and RNA, advancing our understanding of the biology and evolution of these diverse systems.

The promise of cryo-EM to explore RNA structural dynamics

Conformational dynamics are essential to macromolecular function. This is certainly true of RNA, whose ability to undergo programmed conformational dynamics is essential to create and regulate complex biological processes. However, methods to easily and simultaneously interrogate both the structure and conformational dynamics of fully functional RNAs in isolation and in complex with proteins have not historically been available. Due to its ability to image and classify single particles, cryogenic electron microscopy (cryo-EM) has the potential to address this gap and may be particularly amenable to exploring structural dynamics within the three-dimensional folds of biologically active RNAs. We discuss the possibilities and current limitations of applying cryo-EM to simultaneously study RNA structure and conformational dynamics, and present one example that illustrates this (as of yet) not fully realized potential.

Topological crossing in the misfolded Tetrahymena ribozyme resolved by cryo-EM

Taking advantage of single-particle cryogenic electron microscopy (cryo-EM) to analyze highly heterogeneous or flexible samples, we obtained long-awaited three-dimensional (3D) structures of the misfolded Tetrahymena ribozyme. These structures provide clear evidence for a previously proposed topological isomer model, in which the stereochemically impossible crossing of two core RNA strands prevents rapid rearrangement of the misfolded state to the native state. Topological isomers may be widespread in misfolding of complex RNA, and these cryo-EM structures set a foundation for dissecting their detailed kinetic mechanisms and functional consequences in a paradigmatic model system.

The Tetrahymena group I intron has been a key system in the understanding of RNA folding and misfolding. The molecule folds into a long-lived misfolded intermediate (M) in vitro, which has been known to form extensive native-like secondary and tertiary structures but is separated by an unknown kinetic barrier from the native state (N). Here, we used cryogenic electron microscopy (cryo-EM) to resolve misfolded structures of the Tetrahymena L-21 ScaI ribozyme. Maps of three M substates (M1, M2, M3) and one N state were achieved from a single specimen with overall resolutions of 3.5 Å, 3.8 Å, 4.0 Å, and 3.0 Å, respectively. Comparisons of the structures reveal that all the M substates are highly similar to N, except for rotation of a core helix P7 that harbors the ribozyme’s guanosine binding site and the crossing of the strands J7/3 and J8/7 that connect P7 to the other elements in the ribozyme core. This topological difference between the M substates and N state explains the failure of 5′-splice site substrate docking in M, supports a topological isomer model for the slow refolding of M to N due to a trapped strand crossing, and suggests pathways for M-to-N refolding.

Cryo-electron tomography related radiation-damage parameters for individual-molecule 3D structure determination

To understand the dynamic structure-function relationship of soft- and biomolecules, the determination of the three-dimensional (3D) structure of each individual molecule (nonaveraged structure) in its native state is sought-after. Cryo-electron tomography (cryo-ET) is a unique tool for imaging an individual object from a series of tilted views. However, due to radiation damage from the incident electron beam, the tolerable electron dose limits image contrast and the signal-to-noise ratio (SNR) of the data, preventing the 3D structure determination of individual molecules, especially at high-resolution. Although recently developed technologies and techniques, such as the direct electron detector, phase plate, and computational algorithms, can partially improve image contrast/SNR at the same electron dose, the high-resolution structure, such as tertiary structure of individual molecules, has not yet been resolved. Here, we review the cryo-electron microscopy (cryo-EM) and cryo-ET experimental parameters to discuss how these parameters affect the extent of radiation damage. This discussion can guide us in optimizing the experimental strategy to increase the imaging dose or improve image SNR without increasing the radiation damage. With a higher dose, a higher image contrast/SNR can be achieved, which is crucial for individual-molecule 3D structure. With 3D structures determined from an ensemble of individual molecules in different conformations, the molecular mechanism through their biochemical reactions, such as self-folding or synthesis, can be elucidated in a straightforward manner.

Cryo-EM reveals an entangled kinetic trap in the folding of a catalytic RNA

Functional RNAs fold through complex pathways that can contain misfolded "kinetic traps." A complete model of RNA folding requires understanding the formation of these misfolded states, but they are difficult to characterize because of their transient and potentially conformationally dynamic nature. We used cryo-electron microscopy (cryo-EM) to visualize a long-lived misfolded state in the folding pathway of the Tetrahymena thermophila group I intron, a paradigmatic RNA structure-function model system. The structure revealed how this state forms native-like secondary structure and tertiary contacts but contains two incorrectly crossed strands, consistent with a previous model. This incorrect topology mispositions a critical catalytic domain and cannot be resolved locally as extensive refolding is required. This work provides a structural framework for interpreting decades of biochemical and functional studies and demonstrates the power of cryo-EM for the exploration of RNA folding pathways.

Box-shaped ribozyme octamer formed by face-to-face dimerization of a pair of square-shaped ribozyme tetramers

Naturally occurring ribozymes with defined three-dimensional (3D) structures serve as promising platforms for the design and construction of artificial RNA nanostructures. We constructed a hexameric ribozyme nanostructure by face-to-face dimerization of a pair of triangular ribozyme trimers, unit RNAs of which were derived from the Tetrahymena group I ribozyme. In this study, we have expanded the dimerization strategy to a square-shaped ribozyme tetramer by introducing four pillar units. The resulting box-shaped nanostructures, which contained eight ribozyme units, can be assembled from either four or two components of their unit RNAs.

DNA-catalysed alternative RNA splicing

We report DNA-catalysed alternative RNA splicing in vitro. Using modular DNA catalysts with RNA endonuclease and RNA ligase activities, we show that DNA can modulate RNA structure and activity. Furthermore, we illustrate that such DNA-catalysed reactions can yield, from a common precursor, different splicing isoforms with distinct functions.

How does RNA fold dynamically?

Recent advances in interrogating RNA folding dynamics have shown the classical model of RNA folding to be incomplete. Here, we pose three prominent questions for the field that are at the forefront of our understanding of the importance of RNA folding dynamics for RNA function. The first centers on the most appropriate biophysical framework to describe changes to the RNA folding energy landscape that a growing RNA chain encounters during transcriptional elongation. The second focuses on the potential ubiquity of strand displacement - a process by which RNA can rapidly change conformations - and how this process may be generally present in broad classes of seemingly different RNAs. The third raises questions about the potential importance and roles of cellular protein factors in RNA conformational switching. Answers to these questions will greatly improve our fundamental knowledge of RNA folding and function, drive biotechnological advances that utilize engineered RNAs, and potentially point to new areas of biology yet to be discovered.

Structural Organization of S516 Group I Introns in Myxomycetes

Group I introns are mobile genetic elements encoding self-splicing ribozymes. Group I introns in nuclear genes are restricted to ribosomal DNA of eukaryotic microorganisms. For example, the myxomycetes, which represent a distinct protist phylum with a unique life strategy, are rich in nucleolar group I introns. We analyzed and compared 75 group I introns at position 516 in the small subunit ribosomal DNA from diverse and distantly related myxomycete taxa. A consensus secondary structure revealed a conserved group IC1 ribozyme core, but with a surprising RNA sequence complexity in the peripheral regions. Five S516 group I introns possess a twintron organization, where a His-Cys homing endonuclease gene insertion was interrupted by a small spliceosomal intron. Eleven S516 introns contained direct repeat arrays with varying lengths of the repeated motif, a varying copy number, and different structural organizations. Phylogenetic analyses of S516 introns and the corresponding host genes revealed a complex inheritance pattern, with both vertical and horizontal transfers. Finally, we reconstructed the evolutionary history of S516 nucleolar group I introns from insertion of mobile-type introns at unoccupied cognate sites, through homing endonuclease gene degradation and loss, and finally to the complete loss of introns. We conclude that myxomycete S516 introns represent a family of genetic elements with surprisingly dynamic structures despite a common function in RNA self-splicing.

Isotope Labels Combined with Solution NMR Spectroscopy Make Visible the Invisible Conformations of Small-to-Large RNAs

RNA is central to the proper function of cellular processes important for life on earth and implicated in various medical dysfunctions. Yet, RNA structural biology lags significantly behind that of proteins, limiting mechanistic understanding of RNA chemical biology. Fortunately, solution NMR spectroscopy can probe the structural dynamics of RNA in solution at atomic resolution, opening the door to their functional understanding. However, NMR analysis of RNA, with only four unique ribonucleotide building blocks, suffers from spectral crowding and broad linewidths, especially as RNAs grow in size. One effective strategy to overcome these challenges is to introduce NMR-active stable isotopes into RNA. However, traditional uniform labeling methods introduce scalar and dipolar couplings that complicate the implementation and analysis of NMR measurements. This challenge can be circumvented with selective isotope labeling. In this review, we outline the development of labeling technologies and their application to study biologically relevant RNAs and their complexes ranging in size from 5 to 300 kDa by NMR spectroscopy.

Sub-3-Å cryo-EM structure of RNA enabled by engineered homomeric self-assembly

High-resolution structural studies are essential for understanding the folding and function of diverse RNAs. Herein, we present a nanoarchitectural engineering strategy for efficient structural determination of RNA-only structures using single-particle cryogenic electron microscopy (cryo-EM). This strategy-ROCK (RNA oligomerization-enabled cryo-EM via installing kissing loops)-involves installing kissing-loop sequences onto the functionally nonessential stems of RNAs for homomeric self-assembly into closed rings with multiplied molecular weights and mitigated structural flexibility. ROCK enables cryo-EM reconstruction of the Tetrahymena group I intron at 2.98-Å resolution overall (2.85 Å for the core), allowing de novo model building of the complete RNA, including the previously unknown peripheral domains. ROCK is further applied to two smaller RNAs-the Azoarcus group I intron and the FMN riboswitch, revealing the conformational change of the former and the bound ligand in the latter. ROCK holds promise to greatly facilitate the use of cryo-EM in RNA structural studies.

NMR-derived secondary structure of the full-length Ox40 mRNA 3'UTR and its multivalent binding to the immunoregulatory RBP Roquin

Control of posttranscriptional mRNA decay is a crucial determinant of cell homeostasis and differentiation. mRNA lifetime is governed by cis-regulatory elements in their 3' untranslated regions (UTR). Despite ongoing progress in the identification of cis elements we have little knowledge about the functional and structural integration of multiple elements in 3'UTR regulatory hubs and their recognition by mRNA-binding proteins (RBPs). Structural analyses are complicated by inconsistent mapping and prediction of RNA fold, by dynamics, and size. We here, for the first time, provide the secondary structure of a complete mRNA 3'UTR. We use NMR spectroscopy in a divide-and-conquer strategy complemented with SAXS, In-line probing and SHAPE-seq applied to the 3'UTR of Ox40 mRNA, which encodes a T-cell co-receptor repressed by the protein Roquin. We provide contributions of RNA elements to Roquin-binding. The protein uses its extended bi-modal ROQ domain to sequentially engage in a 2:1 stoichiometry with a 3'UTR core motif. We observe differential binding of Roquin to decay elements depending on their structural embedment. Our data underpins the importance of studying RNA regulation in a full sequence and structural context. This study serves as a paradigm for an approach in analysing structured RNA-regulatory hubs and their binding by RBPs.