Group II Intron Self-Splicing

MatK impacts differential chloroplast translation by limiting spliced tRNA-K(UUU) abundance

The protein levels of chloroplast photosynthetic genes and genes related to the chloroplast genetic apparatus vary to adapt to different conditions. However, the underlying mechanisms governing these variations remain unclear. The chloroplast intron Maturase K is encoded within the trnK intron and has been suggested to be required for splicing several group IIA introns, including the trnK intron. In this study, we used RNA immunoprecipitation followed by high-throughput sequencing (RIP-Seq) to identify MatK's preference for binding to group IIA intron domains I and VI within target transcripts. Importantly, these domains are crucial for splice site selection, and we discovered alternative 5'-splice sites in three MatK target introns. The resulting alternative trnK lariat structure showed increased accumulation during heat acclimation. The cognate codon of tRNA-K(UUU) is highly enriched in mRNAs encoding ribosomal proteins and a trnK-matK over-expressor exhibited elevated levels of the spliced tRNA-K(UUU). Ribosome profiling analysis of the overexpressor revealed a significant up-shift in the translation of ribosomal proteins compared to photosynthetic genes. Our findings suggest the existence of a novel regulatory mechanism linked to the abundance of tRNA-K(UUU), enabling the differential expression of functional chloroplast gene groups.

The discovery of a catalytic RNA within RNase P and its legacy

Sidney Altman's discovery of the processing of one RNA by another RNA that acts like an enzyme was revolutionary in biology and the basis for his sharing the 1989 Nobel Prize in Chemistry with Thomas Cech. These breakthrough findings support the key role of RNA in molecular evolution, where replicating RNAs (and similar chemical derivatives) either with or without peptides functioned in protocells during the early stages of life on Earth, an era referred to as the RNA world. Here, we cover the historical background highlighting the work of Altman and his colleagues and the subsequent efforts of other researchers to understand the biological function of RNase P and its catalytic RNA subunit and to employ it as a tool to downregulate gene expression. We primarily discuss bacterial RNase P-related studies but acknowledge that many groups have significantly contributed to our understanding of archaeal and eukaryotic RNase P, as reviewed in this special issue and elsewhere.

Plastid DNA is a major source of nuclear genome complexity and of RNA genes in the orphan crop moringa

BACKGROUND

Unlike Transposable Elements (TEs) and gene/genome duplication, the role of the so-called nuclear plastid DNA sequences (NUPTs) in shaping the evolution of genome architecture and function remains poorly studied. We investigate here the functional and evolutionary fate of NUPTs in the orphan crop Moringa oleifera (moringa), featured by the highest fraction of plastid DNA found so far in any plant genome, focusing on (i) any potential biases in their distribution in relation to specific nuclear genomic features, (ii) their contribution to the emergence of new genes and gene regions, and (iii) their impact on the expression of target nuclear genes.

RESULTS

In agreement with their potential mutagenic effect, NUPTs are underrepresented among structural genes, although their overall transcription levels and broadness were only lower when involved exonic regions; the occurrence of plastid DNA generally did not result in a broader expression, except among those affected in introns by older NUPTs. In contrast, we found a strong enrichment of NUPTs among specific superfamilies of retrotransposons and several classes of RNA genes, including those participating in the protein biosynthetic machinery (i.e., rRNA and tRNA genes) and a specific class of regulatory RNAs. A significant fraction of NUPT RNA genes was found to be functionally expressed, thus potentially contributing to the nuclear pool.

CONCLUSIONS

Our results complete our view of the molecular factors driving the evolution of nuclear genome architecture and function, and support plastid DNA in moringa as a major source of (i) genome complexity and (ii) the nuclear pool of RNA genes.

Highly Reactive Group I Introns Ubiquitous in Pathogenic Fungi

Systemic fungal infections are a growing public health threat, and yet viable antifungal drug targets are limited as fungi share a similar proteome with humans. However, features of RNA metabolism and the noncoding transcriptomes in fungi are distinctive. For example, fungi harbor highly structured RNA elements that humans lack, such as self-splicing introns within key housekeeping genes in the mitochondria. However, the location and function of these mitochondrial riboregulatory elements has largely eluded characterization. Here we used an RNA-structure-based bioinformatics pipeline to identify the group I introns interrupting key mitochondrial genes in medically relevant fungi, revealing their fixation within a handful of genetic hotspots and their ubiquitous presence across divergent phylogenies of fungi, including all highest priority pathogens such as Candida albicans, Candida auris, Aspergillus fumigatus and Cryptococcus neoformans. We then biochemically characterized two representative introns from C. albicans and C. auris, demonstrating their exceptionally efficient splicing catalysis relative to previously-characterized group I introns. Indeed, the C. albicans mitochondrial intron displays extremely rapid catalytic turnover, even at ambient temperatures and physiological magnesium ion concentrations. Our results unmask a significant new set of players in the RNA metabolism of pathogenic fungi, suggesting a promising new type of antifungal drug target.

Long-read direct RNA sequencing of the mitochondrial transcriptome of Saccharomyces cerevisiae reveals condition-dependent intron abundance

Mitochondria fulfil many essential roles and have their own genome, which is expressed as polycistronic transcripts that undergo co- or posttranscriptional processing and splicing. Due to the inherent complexity and limited technical accessibility of the mitochondrial transcriptome, fundamental questions regarding mitochondrial gene expression and splicing remain unresolved, even in the model eukaryote Saccharomyces cerevisiae. Long-read sequencing could address these fundamental questions. Therefore, a method for the enrichment of mitochondrial RNA and sequencing using Nanopore technology was developed, enabling the resolution of splicing of polycistronic genes and the quantification of spliced RNA. This method successfully captured the full mitochondrial transcriptome and resolved RNA splicing patterns with single-base resolution and was applied to explore the transcriptome of S. cerevisiae grown with glucose or ethanol as the sole carbon source, revealing the impact of growth conditions on mitochondrial RNA expression and splicing. This study uncovered a remarkable difference in the turnover of Group II introns between yeast grown in either mostly fermentative or fully respiratory conditions. Whether this accumulation of introns in glucose medium has an impact on mitochondrial functions remains to be explored. Combined with the high tractability of the model yeast S. cerevisiae, the developed method enables to monitor mitochondrial transcriptome responses in a broad range of relevant contexts, including oxidative stress, apoptosis and mitochondrial diseases.

Hydrolytic endonucleolytic ribozyme (HYER) is programmable for sequence-specific DNA cleavage

Ribozymes are catalytic RNAs with diverse functions including self-splicing and polymerization. This work aims to discover natural ribozymes that behave as hydrolytic and sequence-specific DNA endonucleases, which could be repurposed as DNA manipulation tools. Focused on bacterial group II-C introns, we found that many systems without intron-encoded protein propagate multiple copies in their resident genomes. These introns, named HYdrolytic Endonucleolytic Ribozymes (HYERs), cleaved RNA, single-stranded DNA, bubbled double-stranded DNA (dsDNA), and plasmids in vitro. HYER1 generated dsDNA breaks in the mammalian genome. Cryo-electron microscopy analysis revealed a homodimer structure for HYER1, where each monomer contains a Mg2+-dependent hydrolysis pocket and captures DNA complementary to the target recognition site (TRS). Rational designs including TRS extension, recruiting sequence insertion, and heterodimerization yielded engineered HYERs showing improved specificity and flexibility for DNA manipulation.

Structures, functions and adaptations of the human LINE-1 ORF2 protein

The LINE-1 (L1) retrotransposon is an ancient genetic parasite that has written around one-third of the human genome through a 'copy and paste' mechanism catalysed by its multifunctional enzyme, open reading frame 2 protein (ORF2p)1. ORF2p reverse transcriptase (RT) and endonuclease activities have been implicated in the pathophysiology of cancer2,3, autoimmunity4,5 and ageing6,7, making ORF2p a potential therapeutic target. However, a lack of structural and mechanistic knowledge has hampered efforts to rationally exploit it. We report structures of the human ORF2p 'core' (residues 238-1061, including the RT domain) by X-ray crystallography and cryo-electron microscopy in several conformational states. Our analyses identified two previously undescribed folded domains, extensive contacts to RNA templates and associated adaptations that contribute to unique aspects of the L1 replication cycle. Computed integrative structural models of full-length ORF2p show a dynamic closed-ring conformation that appears to open during retrotransposition. We characterize ORF2p RT inhibition and reveal its underlying structural basis. Imaging and biochemistry show that non-canonical cytosolic ORF2p RT activity can produce RNA:DNA hybrids, activating innate immune signalling through cGAS/STING and resulting in interferon production6-8. In contrast to retroviral RTs, L1 RT is efficiently primed by short RNAs and hairpins, which probably explains cytosolic priming. Other biochemical activities including processivity, DNA-directed polymerization, non-templated base addition and template switching together allow us to propose a revised L1 insertion model. Finally, our evolutionary analysis demonstrates structural conservation between ORF2p and other RNA- and DNA-dependent polymerases. We therefore provide key mechanistic insights into L1 polymerization and insertion, shed light on the evolutionary history of L1 and enable rational drug development targeting L1.

Structural insights into intron catalysis and dynamics during splicing

The group II intron ribonucleoprotein is an archetypal splicing system with numerous mechanistic parallels to the spliceosome, including excision of lariat introns1,2. Despite the importance of branching in RNA metabolism, structural understanding of this process has remained elusive. Here we present a comprehensive analysis of three single-particle cryogenic electron microscopy structures captured along the splicing pathway. They reveal the network of molecular interactions that specifies the branchpoint adenosine and positions key functional groups to catalyse lariat formation and coordinate exon ligation. The structures also reveal conformational rearrangements of the branch helix and the mechanism of splice site exchange that facilitate the transition from branching to ligation. These findings shed light on the evolution of splicing and highlight the conservation of structural components, catalytic mechanism and dynamical strategies retained through time in premessenger RNA splicing machines.

The DEAD-box RNA helicase ZmRH48 is required for the splicing of multiple mitochondrial introns, mitochondrial complex biosynthesis, and seed development in maize

RNA helicases participate in nearly all aspects of RNA metabolism by rearranging RNAs or RNA-protein complexes in an adenosine triphosphate-dependent manner. Due to the large RNA helicase families in plants, the precise roles of many RNA helicases in plant physiology and development remain to be clarified. Here, we show that mutations in maize (Zea mays) DEAD-box RNA helicase 48 (ZmRH48) impair the splicing of mitochondrial introns, mitochondrial complex biosynthesis, and seed development. Loss of ZmRH48 function severely arrested embryogenesis and endosperm development, leading to defective kernel formation. ZmRH48 is targeted to mitochondria, where its deficiency dramatically reduced the splicing efficiency of five cis-introns (nad5 intron 1; nad7 introns 1, 2, and 3; and ccmFc intron 1) and one trans-intron (nad2 intron 2), leading to lower levels of mitochondrial complexes I and III. ZmRH48 interacts with two unique pentatricopeptide repeat (PPR) proteins, PPR-SMR1 and SPR2, which are required for the splicing of over half of all mitochondrial introns. PPR-SMR1 interacts with SPR2, and both proteins interact with P-type PPR proteins and Zm-mCSF1 to facilitate intron splicing. These results suggest that ZmRH48 is likely a component of a splicing complex and is critical for mitochondrial complex biosynthesis and seed development.

In Silico Discovery of Group II Intron RNA Splicing Inhibitors

Here, we describe the discovery of compounds that inhibit self-splicing in group II introns. Using docking calculations, we targeted the catalytic active site within the Oceanobacillus iheyensis group IIC intron and virtually screened a library of lead-like compounds. From this initial virtual screen, we identified three unique scaffolds that inhibit splicing in vitro. Additional tests revealed that an analog of the lead scaffold inhibits splicing in an intron-dependent manner. Furthermore, this analog exhibited activity against the group II intron from a different class: the yeast ai5γ IIB intron. The splicing inhibitors we identified could serve as chemical tools for developing group II intron-targeted antifungals, and, more broadly, our results highlight the potential of in silico techniques for identifying bioactive hits against structured and functionally complex RNAs.

Arena: Rapid and Accurate Reconstruction of Full Atomic RNA Structures From Coarse-grained Models

RNA tertiary structures from experiments or computational predictions often contain missing atoms, which prevent analyses requiring full atomic structures. Current programs for RNA reconstruction can be slow, inaccurate, and/or require specific atoms to be present in the input. We present Arena (Atomic Reconstruction of RNA), which reconstructs a full atomic RNA structure from residues that can have as few as one atom. Arena first fills in missing atoms and then iteratively refines their placement to reduce nonideal geometries. We benchmarked Arena on a dataset of 361 RNA structures, where Arena achieves high accuracy and speed compared to other structure reconstruction programs. For example, Arena was used to reconstruct full atomic structures from a single phosphorus atom per nucleotide to, on average, within 3.63 Å RMSD of the experimental structure, while virtually removing all clashes and running in <3 s, which is 353× and 46× faster than state-of-the-art programs PDBFixer and C2A, respectively. The Arena source code is available at https://github.com/pylelab/Arena and the webserver at https://zhanggroup.org/Arena/.

Rock, scissors, paper: How RNA structure informs function

RNA can fold back on itself to adopt a wide range of structures. These range from relatively simple hairpins to intricate 3D folds and can be accompanied by regulatory interactions with both metabolites and macromolecules. The last 50 yr have witnessed elucidation of an astonishing array of RNA structures including transfer RNAs, ribozymes, riboswitches, the ribosome, the spliceosome, and most recently entire RNA structuromes. These advances in RNA structural biology have deepened insight into fundamental biological processes including gene editing, transcription, translation, and structure-based detection and response to temperature and other environmental signals. These discoveries reveal that RNA can be relatively static, like a rock; that it can have catalytic functions of cutting bonds, like scissors; and that it can adopt myriad functional shapes, like paper. We relate these extraordinary discoveries in the biology of RNA structure to the plant way of life. We trace plant-specific discovery of ribozymes and riboswitches, alternative splicing, organellar ribosomes, thermometers, whole-transcriptome structuromes and pan-structuromes, and conclude that plants have a special set of RNA structures that confer unique types of gene regulation. We finish with a consideration of future directions for the RNA structure-function field.

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.

Categorizing 161 plant (streptophyte) mitochondrial group II introns into 29 families of related paralogues finds only limited links between intron mobility and intron-borne maturases

Group II introns are common in the two endosymbiotic organelle genomes of the plant lineage. Chloroplasts harbor 22 positionally conserved group II introns whereas their occurrence in land plant (embryophyte) mitogenomes is highly variable and specific for the seven major clades: liverworts, mosses, hornworts, lycophytes, ferns, gymnosperms and flowering plants. Each plant group features "signature selections" of ca. 20-30 paralogues from a superset of altogether 105 group II introns meantime identified in embryophyte mtDNAs, suggesting massive intron gains and losses along the backbone of plant phylogeny. We report on systematically categorizing plant mitochondrial group II introns into "families", comprising evidently related paralogues at different insertion sites, which may even be more similar than their respective orthologues in phylogenetically distant taxa. Including streptophyte (charophyte) algae extends our sampling to 161 and we sort 104 streptophyte mitochondrial group II introns into 25 core families of related paralogues evidently arising from retrotransposition events. Adding to discoveries of only recently created intron paralogues, hypermobile introns and twintrons, our survey led to further discoveries including previously overlooked "fossil" introns in spacer regions or e.g., in the rps8 pseudogene of lycophytes. Initially excluding intron-borne maturase sequences for family categorization, we added an independent analysis of maturase phylogenies and find a surprising incongruence between intron mobility and the presence of intron-borne maturases. Intriguingly, however, we find that several examples of nuclear splicing factors meantime characterized simultaneously facilitate splicing of independent paralogues now placed into the same intron families. Altogether this suggests that plant group II intron mobility, in contrast to their bacterial counterparts, is not intimately linked to intron-encoded maturases.

Importance of mobile genetic element immunity in numerically abundant Trichodesmium clades

The colony-forming cyanobacteria Trichodesmium spp. are considered one of the most important nitrogen-fixing genera in the warm, low nutrient ocean. Despite this central biogeochemical role, many questions about their evolution, physiology, and trophic interactions remain unanswered. To address these questions, we describe Trichodesmium pangenomic potential via significantly improved genomic assemblies from two isolates and 15 new >50% complete Trichodesmium metagenome-assembled genomes from hand-picked, Trichodesmium colonies spanning the Atlantic Ocean. Phylogenomics identified ~four N2 fixing clades of Trichodesmium across the transect, with T. thiebautii dominating the colony-specific reads. Pangenomic analyses showed that all T. thiebautii MAGs are enriched in COG defense mechanisms and encode a vertically inherited Type III-B Clustered Regularly Interspaced Short Palindromic Repeats and associated protein-based immunity system (CRISPR-Cas). Surprisingly, this CRISPR-Cas system was absent in all T. erythraeum genomes, vertically inherited by T. thiebautii, and correlated with increased signatures of horizontal gene transfer. Additionally, the system was expressed in metaproteomic and transcriptomic datasets and CRISPR spacer sequences with 100% identical hits to field-assembled, putative phage genome fragments were identified. While the currently CO2-limited T. erythraeum is expected to be a 'winner' of anthropogenic climate change, their genomic dearth of known phage resistance mechanisms, compared to T. thiebautii, could put this outcome in question. Thus, the clear demarcation of T. thiebautii maintaining CRISPR-Cas systems, while T. erythraeum does not, identifies Trichodesmium as an ecologically important CRISPR-Cas model system, and highlights the need for more research on phage-Trichodesmium interactions.

A New EBS2b-IBS2b Base Paring (A-8/T-8) Improved the Gene-Targeting Efficiency of Thermotargetron in Escherichia coli

Thermophilic group II intron is one type of retrotransposon composed of intron RNA and intron-encoded protein (IEP), which can be utilized in gene targeting by harnessing their novel ribozyme-based DNA integration mechanism termed "retrohoming." It is mediated by a ribonucleoprotein (RNP) complex that contains the excised intron lariat RNA and an IEP with reverse transcriptase (RT) activity. The RNP recognizes targeting sites by exon-binding sequences 2 (EBS2)/intron-binding sequences 2 (IBS2), EBS1/IBS1, and EBS3/IBS3 bases pairing. Previously, we developed the TeI3c/4c intron as a thermophilic gene targeting system-Thermotargetron (TMT). However, we found that the targeting efficiency of TMT varies significantly at different targeting sites, which leads to a relatively low success rate. To further improve the success rate and gene-targeting efficiency of TMT, we constructed a Random Gene-targeting Plasmids Pool (RGPP) to analyze the sequence recognition preference of TMT. A new base pairing, located at the -8 site between EBS2/IBS2 and EBS1/IBS1 (named EBS2b-IBS2b), increased the success rate (2.45- to 5.07-fold) and significantly improved gene-targeting efficiency of TMT. A computer algorithm (TMT 1.0), based on the newly discovered sequence recognition roles, was also developed to facilitate the design of TMT gene-targeting primers. The present work could essentially expand the practicalities of TMT in the genome engineering of heat-tolerance mesophilic and thermophilic bacteria. IMPORTANCE The randomized base pairing in the interval of IBS2 and IBS1 of Tel3c/4c intron (-8 and -7 sites) in Thermotargetron (TMT) results in a low success rate and gene-targeting efficiency in bacteria. In the present work, we constructed a randomized gene-targeting plasmids pool (RGPP) to study whether there is a base preference in target sequences. Among all the successful "retrohoming" targets, we found that a new EBS2b-IBS2b base paring (A-8/T-8) significantly increased TMT's gene-targeting efficiency, and the concept is also applicable to other gene targets in redesigned gene-targeting plasmids pool in E. coli. The improved TMT is a promising tool for the genetic engineering of bacteria and could promote metabolic engineering and synthetic biology research in valuable microbes that recalcitrance for genetic manipulation.

Mitochondrial reactivity following acute exposure to experimental pain testing in people with HIV and chronic pain

Background: Physical stressors can cause a physiological response that can contribute to an increase in mitochondrial dysfunction and Mitochondrial DNA damage (mtDNA damage). People living with HIV (PWH) are more likely to suffer from chronic pain and may be more susceptible to mitochondrial dysfunction following exposure to a stressor. We used Quantitative Sensory Testing (QST) as an acute painful stressor in order to investigate whether PWH with/without chronic pain show differential mitochondrial physiological responses. Methods: The current study included PWH with (n = 26), and without (n = 29), chronic pain. Participants completed a single session that lasted approximately 180 min, including QST. Blood was taken prior to and following the QST battery for assays measuring mtDNA damage, mtDNA copy number, and mtDNA damage-associated molecular pattern (DAMP) levels (i.e., ND1 and ND6). Results: We examined differences between those with and without pain on various indicators of mitochondrial reactivity following exposure to QST. However, only ND6 and mtDNA damage were shown to be statistically significant between pain groups. Conclusion: PWH with chronic pain showed greater mitochondrial reactivity to laboratory stressors. Consequently, PWH and chronic pain may be more susceptible to conditions in which mitochondrial damage/dysfunction play a central role, such as cognitive decline.

Structures of a mobile intron retroelement poised to attack its structured DNA target

Group II introns are ribozymes that catalyze their self-excision and function as retroelements that invade DNA. As retrotransposons, group II introns form ribonucleoprotein (RNP) complexes that roam the genome, integrating by reversal of forward splicing. Here we show that retrotransposition is achieved by a tertiary complex between a structurally elaborate ribozyme, its protein mobility factor, and a structured DNA substrate. We solved cryo-electron microscopy structures of an intact group IIC intron-maturase retroelement that was poised for integration into a DNA stem-loop motif. By visualizing the RNP before and after DNA targeting, we show that it is primed for attack and fits perfectly with its DNA target. This study reveals design principles of a prototypical retroelement and reinforces the hypothesis that group II introns are ancient elements of genetic diversification.

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.

Unique and Repeated Stwintrons (Spliceosomal Twin Introns) in the Hypoxylaceae

Introns are usually non-coding sequences interrupting open reading frames in pre-mRNAs [D1,2]. Stwintrons are nested spliceosomal introns, where an internal intron splits a second donor sequence into two consecutive splicing reactions leading to mature mRNA. In Hypoxylon sp. CO27-5, 36 highly sequence-similar [D1,2] stwintrons are extant (sister stwintrons). An additional 81 [D1,2] sequence-unrelated stwintrons are described here. Most of them are located at conserved gene positions rooted deep in the Hypoxylaceae. Absence of exonic sequence bias at the exon–stwintron junctions and a very similar phase distribution were noted for both groups. The presence of an underlying sequence symmetry in all 117 stwintrons was striking. This symmetry, more pronounced near the termini of most of the full-length sister stwintrons, may lead to a secondary structure that brings into close proximity the most distal splice sites, the donor of the internal and the acceptor of the external intron. The Hypoxylon stwintrons were overwhelmingly excised by consecutive splicing reactions precisely removing the whole intervening sequence, whereas one excision involving the distal splice sites led to a frameshift. Alternative (mis)splicing took place for both sister and uniquely occurring stwintrons. The extraordinary symmetry of the sister stwintrons thus seems dispensable for the infrequent, direct utilisation of the distal splice sites.

Distinct Expansion of Group II Introns During Evolution of Prokaryotes and Possible Factors Involved in Its Regulation

Group II introns (G2Is) are ribozymes that have retroelement characteristics in prokaryotes. Although G2Is are suggested to have been an important evolutionary factor in the prokaryote-to-eukaryote transition, comprehensive analyses of these introns among the tens of thousands of prokaryotic genomes currently available are still limited. Here, we developed a bioinformatic pipeline that systematically collects G2Is and applied it to prokaryotic genomes. We found that in bacteria, 25% (447 of 1,790) of the total representative genomes had an average of 5.3 G2Is, and in archaea, 9% (28 of 296) of the total representative genomes had an average of 3.0 G2Is. The greatest number of G2Is per genome was 101 in Arthrospira platensis (phylum Cyanobacteriota). A comprehensive sequence analysis of the intron-encoded protein (IEP) in each G2I sequence was conducted and resulted in the addition of three new IEP classes (U1-U3) to the previous classification. This analysis suggested that about 30% of all IEPs are non-canonical IEPs. The number of G2Is per genome was defined almost at the phylum level, and at least in the following two phyla, Firmicutes, and Cyanobacteriota, the type of IEP was largely associated as a factor in the G2I increase, i.e., there was an explosive increase in G2Is with bacterial C-type IEPs, mainly in the phylum Firmicutes, and in G2Is with CL-type IEPs, mainly in the phylum Cyanobacteriota. We also systematically analyzed the relationship between genomic signatures and the mechanism of these increases in G2Is. This is the first study to systematically characterize G2Is in the prokaryotic phylogenies.

Chemical Dual End-Labeling of Large Ribozymes

Fast and efficient site-specific labeling of long RNAs is one of the main bottlenecks limiting distance measurements by means of Förster resonance energy transfer (FRET) or electron paramagnetic resonance (EPR) spectroscopy. Here, we present an optimized protocol for dual end-labeling with different fluorophores at the same time meeting the restrictions of highly labile and degradation-sensitive RNAs. We describe in detail the dual-labeling of a catalytically active wild-type group II intron as a typical representative of long functional RNAs. The modular procedure chemically activates the 5'-phosphate and the 3'-ribose for bioconjugation with a pair of fluorophores, as shown herein, or with spin labels. The mild reaction conditions preserve the structural and functional integrity of the biomacromolecule and results in covalent, dual-labeled RNA in its pre-catalytic state in yields suitable for both ensemble and single-molecule FRET experiments.

Natural circularly permuted group II introns in bacteria produce RNA circles

Group II self-splicing introns are large structured RNAs that remove themselves from transcripts while simultaneously sealing the resulting gaps. Some representatives can subsequently reverse splice into DNA, accounting for their pervasive distribution in bacteria. The catalytically active tertiary structure of each group II intron is assembled from six domains that are arranged in a conserved order. Here, we report structural isomers of group II introns, called CP group II ribozymes, wherein the characteristic order of domains has been altered. Domains five and six, which normally reside at the 3' end of group II introns, instead occupy the 5' end to form circularly permuted variants. These unusual group II intron derivatives are catalytically active and generate large linear branched and small circular RNAs, reaction products that are markedly different from those generated by canonical group II introns. The biological role of CP group II ribozymes is currently unknown.

Development of an In Vitro Chloroplast Splicing System: Sequences Required for Correct pre-mRNA Splicing

Chloroplast genomes in land plants include approximately 20 intron-containing genes. Most of the introns are similar to the group II introns found in fungi, algae and some bacteria, but no self-splicing has been reported. To analyze splicing reactions in chloroplasts, we developed a tobacco (Nicotiana tabacum) chloroplast-based in vitro system. We optimized the splicing reaction using atpF precursor messenger RNA (pre-mRNA). Our system requires a high ATP concentration, whereas ATP is not necessary for self-splicing group II introns. Self-splicing group II introns possess two exon-binding sites (EBS1 and 2) complementary to two intron-binding sites (IBS1 and 2) in the 3' end of 5' exons, which are involved in 5' splice-site selection. Using our in vitro system and atpF pre-mRNA, we analyzed short sequences corresponding to the above EBSs and IBSs. Mutation analyses revealed that EBS1-IBS1 pairing is essential, while EBS2-IBS2 pairing is important but not crucial for splicing. The first 3' exon nucleotide determines the 3' splice sites of self-splicing introns. However, mutations to this nucleotide in atpF pre-mRNA did not affect splicing. This result suggests that the mechanism underlying chloroplast pre-mRNA splicing differs partly from that mediating the self-splicing of group II introns.

The pentatricopeptide repeat protein EMP603 is required for the splicing of mitochondrial Nad1 intron 2 and seed development in maize

Intron splicing is an essential event in post-transcriptional RNA processing in plant mitochondria, which requires the participation of diverse nuclear-encoded splicing factors. However, it is presently unclear how these proteins cooperatively take part in the splicing of specific introns. In this study, we characterized a nuclear-encoded mitochondrial P-type pentatricopeptide repeat (PPR) protein named EMP603. This protein is essential for splicing of intron 2 in the Nad1 gene and interacts with the mitochondria-localized DEAD-box RNA helicase PMH2-5140, the RAD52-like proteins ODB1-0814 and ODB1-5061, and the CRM domain-containing protein Zm-mCSF1. Further study revealed that the N-terminal region of EMP603 interacts with the DEAD-box of PMH2-5140, the CRM domain of Zm-mCSF1, and OBD1-5061, but not with OBD1-0814, whereas the PPR domain of EMP603 can interact with ODB1-0814, ODB1-5061, and PMH2-5140, but not with Zm-mCSF1. Defects in EMP603 severely disrupt the assembly and activity of mitochondrial complex I, leading to impaired mitochondrial function, and delayed seed development. The interactions revealed between EMP603 and PMH2-5140, ODB1-0814, ODB1-5061, and Zm-mCSF1 indicate a possible involvement of a dynamic 'spliceosome-like' complex in intron splicing, and may accelerate the elucidation of the intron splicing mechanism in plant mitochondria.

Getting to the bottom of lncRNA mechanism: structure-function relationships

While long non-coding RNAs are known to play key roles in disease and development, relatively few structural studies have been performed for this important class of RNAs. Here, we review functional studies of long non-coding RNAs and expose the need for high-resolution 3-D structural studies, discussing the roles of long non-coding RNAs in the cell and how structure–function relationships might be used to elucidate further understanding. We then describe structural studies of other classes of RNAs using chemical probing, nuclear magnetic resonance, small-angle X-ray scattering, X-ray crystallography, and cryogenic electron microscopy (cryo-EM). Next, we review early structural studies of long non-coding RNAs to date and describe the way forward for the structural biology of long non-coding RNAs in terms of cryo-EM.

Mechanisms of catalytic RNA molecules

Ribozymes are folded catalytic RNA molecules that perform important biological functions. Since the discovery of the first RNA with catalytic activity in 1982, a large number of ribozymes have been reported. While most catalytic RNA molecules act alone, some RNA-based catalysts, such as RNase P, the ribosome, and the spliceosome, need protein components to perform their functions in the cell. In the last decades, the structure and mechanism of several ribozymes have been studied in detail. Aside from the ribosome, which catalyzes peptide bond formation during protein synthesis, the majority of known ribozymes carry out mostly phosphoryl transfer reactions, notably trans-esterification or hydrolysis reactions. In this review, we describe the main features of the mechanisms of various types of ribozymes that can function with or without the help of proteins to perform their biological functions.

The emerging role of circular RNAs in cardiovascular diseases

Cardiovascular disease (CVD) is one of the vital causes of morbidity and mortality, and the number of deaths from CVD has increased worldwide. Circular RNAs (circRNAs) is a novel type of endogenous noncoding RNA, which can form covalent closed continuous rings and are highly expressed in the eukaryotic transcriptome. In recent years, research on circRNAs have been increasing and the researchers have also become cumulatively aware of the association between circRNAs and CVD. This review highlights the biogenesis and functions of circRNAs and the role in cardiovascular diseases.

Social Networking of Quasi-Species Consortia drive Virolution via Persistence

The emergence of cooperative quasi-species consortia (QS-C) thinking from the more accepted quasispecies equations of Manfred Eigen, provides a conceptual foundation from which concerted action of RNA agents can now be understood. As group membership becomes a basic criteria for the emergence of living systems, we also start to understand why the history and context of social RNA networks become crucial for survival and function. History and context of social RNA networks also lead to the emergence of a natural genetic code. Indeed, this QS-C thinking can also provide us with a transition point between the chemical world of RNA replicators and the living world of RNA agents that actively differentiate self from non-self and generate group identity with membership roles. Importantly the social force of a consortia to solve complex, multilevel problems also depend on using opposing and minority functions. The consortial action of social networks of RNA stem-loops subsequently lead to the evolution of cellular organisms representing a tree of life.

RNA-Targeting Splicing Modifiers: Drug Development and Screening Assays

RNA splicing is an essential step in producing mature messenger RNA (mRNA) and other RNA species. Harnessing RNA splicing modifiers as a new pharmacological modality is promising for the treatment of diseases caused by aberrant splicing. This drug modality can be used for infectious diseases by disrupting the splicing of essential pathogenic genes. Several antisense oligonucleotide splicing modifiers were approved by the U.S. Food and Drug Administration (FDA) for the treatment of spinal muscular atrophy (SMA) and Duchenne muscular dystrophy (DMD). Recently, a small-molecule splicing modifier, risdiplam, was also approved for the treatment of SMA, highlighting small molecules as important warheads in the arsenal for regulating RNA splicing. The cellular targets of these approved drugs are all mRNA precursors (pre-mRNAs) in human cells. The development of novel RNA-targeting splicing modifiers can not only expand the scope of drug targets to include many previously considered “undruggable” genes but also enrich the chemical-genetic toolbox for basic biomedical research. In this review, we summarized known splicing modifiers, screening methods for novel splicing modifiers, and the chemical space occupied by the small-molecule splicing modifiers.

SHAPE Profiling to Probe Group II Intron Conformational Dynamics During Splicing

Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) is a widely used technique for studying the structure and function of RNA molecules. It characterizes the flexibility of single nucleotides in the context of the local RNA structure. Here we describe the application of SHAPE-MaP (mutational profiling) to study different conformational states of the group II intron during the self-splicing reaction.

Group II Introns Generate Functional Chimeric Relaxase Enzymes with Modified Specificities through Exon Shuffling at Both the RNA and DNA Level

Group II introns are large self-splicing RNA enzymes with a broad but somewhat irregular phylogenetic distribution. These ancient retromobile elements are the proposed ancestors of approximately half the human genome, including the abundant spliceosomal introns and non-long terminal repeat retrotransposons. In contrast to their eukaryotic derivatives, bacterial group II introns have largely been considered as harmful selfish mobile retroelements that parasitize the genome of their host. As a challenge to this view, we recently uncovered a new intergenic trans-splicing pathway that generates an assortment of mRNA chimeras. The ability of group II introns to combine disparate mRNA fragments was proposed to increase the genetic diversity of the bacterial host by shuffling coding sequences. Here, we show that the Ll.LtrB and Ef.PcfG group II introns from Lactococcus lactis and Enterococcus faecalis respectively can both use the intergenic trans-splicing pathway to catalyze the formation of chimeric relaxase mRNAs and functional proteins. We demonstrated that some of these compound relaxase enzymes yield gain-of-function phenotypes, being significantly more efficient than their precursor wild-type enzymes at supporting bacterial conjugation. We also found that relaxase enzymes with shuffled functional domains are produced in biologically relevant settings under natural expression levels. Finally, we uncovered examples of lactococcal chimeric relaxase genes with junctions exactly at the intron insertion site. Overall, our work demonstrates that the genetic diversity generated by group II introns, at the RNA level by intergenic trans-splicing and at the DNA level by recombination, can yield new functional enzymes with shuffled exons, which can lead to gain-of-function phenotypes.

Alternative RNA structures formed during transcription depend on elongation rate and modify RNA processing

Most RNA processing occurs co-transcriptionally. We interrogated nascent pol II transcripts by chemical and enzymatic probing and determined how the "nascent RNA structureome" relates to splicing, A-I editing and transcription speed. RNA folding within introns and steep structural transitions at splice sites are associated with efficient co-transcriptional splicing. A slow pol II mutant elicits extensive remodeling into more folded conformations with increased A-I editing. Introns that become more structured at their 3' splice sites get co-transcriptionally excised more efficiently. Slow pol II altered folding of intronic Alu elements where cryptic splicing and intron retention are stimulated, an outcome mimicked by UV, which decelerates transcription. Slow transcription also remodeled RNA folding around alternative exons in distinct ways that predict whether skipping or inclusion is favored, even though it occurs post-transcriptionally. Hence, co-transcriptional RNA folding modulates post-transcriptional alternative splicing. In summary, the plasticity of nascent transcripts has widespread effects on RNA processing.

Engineering highly efficient backsplicing and translation of synthetic circRNAs

Circular RNAs (circRNAs) are highly stable RNA molecules that are attractive templates for expression of therapeutic proteins and non-coding RNAs. In eukaryotes, circRNAs are primarily generated by the spliceosome through backsplicing. Here, we interrogate different molecular elements including intron type and length, Alu repeats, internal ribosome entry sites (IRESs), and exon length essential for circRNA formation and exploit this information to engineer robust backsplicing and circRNA expression. Specifically, we leverage the finding that the downstream intron can tolerate large inserts without affecting splicing to achieve tandem expression of backspliced circRNAs and tRNA intronic circRNAs from the same template. Further, truncation of selected intronic regions markedly increased circRNA formation in different cell types in vitro as well as AAV-mediated circRNA expression in cardiac and skeletal muscle tissue in vivo. We also observed that different IRES elements and exon length influenced circRNA expression and translation, revealing an exonic contribution to splicing, as evidenced by different RNA species produced. Taken together, these data provide new insight into improving the design and expression of synthetic circRNAs. When combined with AAV capsid and promoter technologies, the backsplicing introns and IRES elements constituting this modular platform significantly expand the gene expression toolkit.

Conjugation as a Highly Sensitive Assay to Study Group II Intron Splicing In Vivo

Group II introns are noncoding sequences that interrupt genes, and that must be removed or spliced-out at the RNA level during gene expression. Following the transcription of interrupted genes, group II introns self-splice while concurrently ligating their flanking exons to generate mature mRNAs ready for translation. Ll.LtrB, the model group II intron from the gram-positive bacterium Lactococcus lactis, interrupts the gene coding for a relaxase enzyme that initiates the transfer of mobile elements by conjugation. This functional link between group II intron splicing and conjugative transfer enabled us to engineer highly sensitive splicing assays using the native biological context of Ll.LtrB. The splicing efficiency/conjugation assay was developed to determine the splicing competence of various Ll.LtrB mutants, whereas the splicing selection/conjugation assay was established to isolate splicing-proficient variants from a randomly generated bank of mutated introns.

Site-specific RNA methylation by a methyltransferase ribozyme

Nearly all classes of coding and non-coding RNA undergo post-transcriptional modification, including RNA methylation. Methylated nucleotides are among the evolutionarily most-conserved features of transfer (t)RNA and ribosomal (r)RNA1,2. Many contemporary methyltransferases use the universal cofactor S-adenosylmethionine (SAM) as a methyl-group donor. SAM and other nucleotide-derived cofactors are considered to be evolutionary leftovers from an RNA world, in which ribozymes may have catalysed essential metabolic reactions beyond self-replication3. Chemically diverse ribozymes seem to have been lost in nature, but may be reconstructed in the laboratory by in vitro selection. Here we report a methyltransferase ribozyme that catalyses the site-specific installation of 1-methyladenosine in a substrate RNA, using O6-methylguanine as a small-molecule cofactor. The ribozyme shows a broad RNA-sequence scope, as exemplified by site-specific adenosine methylation in various RNAs. This finding provides fundamental insights into the catalytic abilities of RNA, serves a synthetic tool to install 1-methyladenosine in RNA and may pave the way to in vitro evolution of other methyltransferase and demethylase ribozymes.

Molecular crowding and RNA catalysis

RNA enzymes or ribozymes catalyze some of the most important reactions in biology and are thought to have played a central role in the origin and evolution of life on earth. Catalytic function in RNA has evolved in crowded cellular environments that are different from dilute solutions in which most in vitro assays are performed. The presence of molecules such as amino acids, polypeptides, alcohols, and sugars in the cell introduces forces that modify the kinetics and thermodynamics of ribozyme-catalyzed reactions. Synthetic molecules are routinely used in in vitro studies to better approximate the properties of biomolecules under in vivo conditions. This review discusses the various forces that operate within simulated crowded solutions in the context of RNA structure, folding, and catalysis. It also explores ideas about how crowding could have been beneficial to the evolution of functional RNAs and the development of primitive cellular systems in a prebiotic milieu.

The molecular structure of long non-coding RNAs: emerging patterns and functional implications

Long non-coding RNAs (lncRNAs) are recently-discovered transcripts that regulate vital cellular processes and are crucially connected to diseases. Despite their unprecedented molecular complexity, it is emerging that lncRNAs possess distinct structural motifs. Remarkably, the 3D shape and topology of full-length, native lncRNAs have been visualized for the first time in the last year. These studies reveal that lncRNA structures dictate lncRNA functions. Here, we review experimentally determined lncRNA structures and emphasize that lncRNA structural characterization requires synergistic integration of computational, biochemical and biophysical approaches. Based on these emerging paradigms, we discuss how to overcome the challenges posed by the complex molecular architecture of lncRNAs, with the goal of obtaining a detailed understanding of lncRNA functions and molecular mechanisms in the future.

DEK46 performs C-to-U editing of a specific site in mitochondrial nad7 introns that is critical for intron splicing and seed development in maize

The self-splicing of group II introns during RNA processing depends on their catalytic structure and is influenced by numerous factors that promote the formation of that structure through direct binding. Here we report that C-to-U editing at a specific position in two nad7 introns is essential to splicing, which also implies that the catalytic activity of non-functional group II introns could be restored by editing. We characterized a maize (Zea mays) mutant, dek46, with a defective kernel phenotype; Dek46 encodes a pentatricopeptide repeat DYW protein exclusively localized in mitochondria. Analyses of the coding regions of mitochondrial transcripts did not uncover differences in RNA editing between dek46 mutant and wild-type maize, but showed that splicing of nad7 introns 3 and 4 is severely reduced in the mutant. Furthermore, editing at nucleotide 22 of domain 5 (D5-C22) of both introns is abolished in dek46. We constructed chimeric introns by swapping D5 of P.li.LSUI2 with D5 of nad7 intron 3. In vitro splicing assays indicated that the chimeric intron containing D5-U22 can be self-spliced, but the one containing D5-C22 cannot. These results indicate that DEK46 functions in the C-to-U editing of D5-C22 of both introns, and the U base at this position is critical to intron splicing.

How Is Precursor Messenger RNA Spliced by the Spliceosome?

Splicing of the precursor messenger RNA, involving intron removal and exon ligation, is mediated by the spliceosome. Together with biochemical and genetic investigations of the past four decades, structural studies of the intact spliceosome at atomic resolution since 2015 have led to mechanistic delineation of RNA splicing with remarkable insights. The spliceosome is proven to be a protein-orchestrated metalloribozyme. Conserved elements of small nuclear RNA (snRNA) constitute the splicing active site with two catalytic metal ions and recognize three conserved intron elements through duplex formation, which are delivered into the splicing active site for branching and exon ligation. The protein components of the spliceosome stabilize the conformation of the snRNA, drive spliceosome remodeling, orchestrate the movement of the RNA elements, and facilitate the splicing reaction. The overall organization of the spliceosome and the configuration of the splicing active site are strictly conserved between human and yeast.

Mg2+ vs Ca2+ bound active site of group II intron- A MD study

Group II introns are enzymes which undergo self-splicing and remove itself from pre-messenger RNA. X-ray structures of group II intron of Oceanobacillus iheyensis at various stages of the self-splicing pathway (pre-hydrolytic, post-hydrolytic, and ligand-free state) revealed intricate atomic interaction network in the active site of the intron. It has been confirmed that a heteronuclear metal ion cluster consisting of four metal ions (K1, K2 sites with K⁺ and M1, M2 sites with Mg²⁺) are crucial for function. Substitution of Mg²⁺ by Ca²⁺ results in loss of enzymatic activity. The X-ray structures not only opens up the possibility of modelling Mg²⁺ and Ca²⁺ bound active site of group II intron and quantitatively estimate the energetics of Mg²⁺ vs Ca²⁺ preference but also explore the relative structural and dynamical differences in response to divalent metal ion substitution. Thus, using X-ray structures as a template we performed molecular dynamics simulations to compare structural and dynamical differences between Mg²⁺ and Ca²⁺ bound active site of group II intron at various stages of the splicing pathway (i.e., pre-hydrolytic, post-hydrolytic, and ligand-free state). Quantitative estimation of Mg²⁺ vs Ca²⁺ selectivity at the M1, M2 sites confirmed Mg²⁺ preference at intron active sites relative to Ca²⁺. Ca²⁺ is relatively more hydrated in the intron active site relative to Mg²⁺. The local environment (bound nucleophilic water, interaction with scissile phosphate) around M2 is strikingly different between Mg²⁺ and Ca²⁺ bound pre-hydrolytic state. In the post-hydrolytic state, the exon part of the hydrolysis product is involved in direct interaction with the M1, whereas the intron part is highly flexible in our MD trajectories. Solvent exposure of M1, M2 sites are least in the pre-hydrolytic state, highest in the ligand-free state, and intermediate in the post-hydrolytic state.

PPR20 Is Required for the cis-Splicing of Mitochondrial nad2 Intron 3 and Seed Development in Maize

Pentatricopeptide repeat (PPR) proteins are helical repeat RNA-binding proteins that function in RNA processing by conferring sequence-specific RNA-binding activity. Owing to the lethality of PPR mutants, functions of many PPR proteins remain obscure. In this study, we report the function of PPR20 in intron splicing in mitochondria and its role in maize seed development. PPR20 is a P-type PPR protein targeted to mitochondria. The ppr20 mutants display slow embryo and endosperm development. Null mutation of PPR20 severely reduces the cis-splicing of mitochondrial nad2 intron 3, resulting in reduction in the assembly and activity of mitochondrial complex I. The ppr20-35 allele with a Mu insertion in the N-terminal region shows a much weaker phenotype. Molecular analyses revealed that the mutant produces a truncated transcript, coding for PPR20ΔN120 lacking the N-terminal 120 amino acids. Subcellular localization revealed that PPR20ΔN120:GFP is able to target to mitochondria as well, suggesting the sequence diversity of the mitochondrial targeting peptides. Another mutant zm_mterf15 was also found to be impaired in the splicing of mitochondrial nad2 intron 3. Further analyses are required to identify the exact function of PPR20 and Zm_mTERF15 in the splicing of nad2 intron 3.

Metal ions and sugar puckering balance single-molecule kinetic heterogeneity in RNA and DNA tertiary contacts

The fidelity of group II intron self-splicing and retrohoming relies on long-range tertiary interactions between the intron and its flanking exons. By single-molecule FRET, we explore the binding kinetics of the most important, structurally conserved contact, the exon and intron binding site 1 (EBS1/IBS1). A comparison of RNA-RNA and RNA-DNA hybrid contacts identifies transient metal ion binding as a major source of kinetic heterogeneity which typically appears in the form of degenerate FRET states. Molecular dynamics simulations suggest a structural link between heterogeneity and the sugar conformation at the exon-intron binding interface. While Mg2+ ions lock the exon in place and give rise to long dwell times in the exon bound FRET state, sugar puckering alleviates this structural rigidity and likely promotes exon release. The interplay of sugar puckering and metal ion coordination may be an important mechanism to balance binding affinities of RNA and DNA interactions in general.

Two Pentatricopeptide Repeat Proteins Are Required for the Splicing of nad5 Introns in Maize

Mitochondrial genes in flowering plants contain predominantly group II introns that require precise splicing before translation into functional proteins. Splicing of these introns is facilitated by various nucleus-encoded splicing factors. Due to lethality of mutants, functions of many splicing factors have not been revealed. Here, we report the function of two P-type PPR proteins PPR101 and PPR231, and their role in maize seed development. PPR101 and PPR231 are targeted to mitochondria. Null mutation of PPR101 and PPR231 arrests embryo and endosperm development, generating empty pericarp and small kernel phenotype, respectively, in maize. Loss-of-function in PPR101 abolishes the splicing of nad5 intron 2, and reduces the splicing of nad5 intron 1. Loss-of-function in PPR231 reduces the splicing of nad5 introns 1, 2, 3 and nad2 intron 3. The absence of Nad5 protein eliminates assembly of complex I, and activates the expression of alternative oxidase AOX2. These results indicate that both PPR101 and PPR231 are required for mitochondrial nad5 introns 1 and 2 splicing, while PPR231 is also required for nad5 intron 3 and nad2 intron 3. Both genes are essential to complex I assembly, mitochondrial function, and maize seed development. This work reveals that the splicing of a single intron involves multiple PPRs.

Visualizing Disordered Single-Stranded RNA: Connecting Sequence, Structure, and Electrostatics

Disordered homopolymeric regions of single-stranded RNA, such as U or A tracts, are found within functional RNAs where they play distinct roles in defining molecular structure and facilitating recognition by partners. Despite this prominence, details of conformational and biophysical properties of these regions have not yet been resolved. We apply a number of experimental techniques to investigate the conformations of these biologically important motifs and provide quantitative measurements of their ion atmospheres. Single strands of RNA display pronounced sequence-dependent conformations that relate to the unique ion atmospheres each attracts. Chains of rU bases are relatively unstructured under all conditions, while chains of rA bases display distinct ordering through stacking or clustering motifs, depending on the composition of the surrounding solution. These dramatic structural differences are consistent with the measured disparity in ion composition and atmospheres around each homopolymer, revealing a complex interplay of base, ion, and single-strand ordering. The unique structural and ionic signatures of homopolymer ssRNAs explains their role(s) in folding structured RNAs and may explain their distinct recognition by protein partners.

Base-pair ambiguity and the kinetics of RNA folding

Background

A pairings of nucleotide sequences. Given this forbidding free-energy landscape, mechanisms have evolved that contribute to a directed and efficient folding process, including catalytic proteins and error-detecting chaperones. Among structural RNA molecules we make a distinction between “bound” molecules, which are active as part of ribonucleoprotein (RNP) complexes, and “unbound,” with physiological functions performed without necessarily being bound in RNP complexes. We hypothesized that unbound molecules, lacking the partnering structure of a protein, would be more vulnerable than bound molecules to kinetic traps that compete with native stem structures. We defined an “ambiguity index”—a normalized function of the primary and secondary structure of an individual molecule that measures the number of kinetic traps available to nucleotide sequences that are paired in the native structure, presuming that unbound molecules would have lower indexes. The ambiguity index depends on the purported secondary structure, and was computed under both the comparative (“gold standard”) and an equilibrium-based prediction which approximates the minimum free energy (MFE) structure. Arguing that kinetically accessible metastable structures might be more biologically relevant than thermodynamic equilibrium structures, we also hypothesized that MFE-derived ambiguities would be less effective in separating bound and unbound molecules.

Results

We have introduced an intuitive and easily computed function of primary and secondary structures that measures the availability of complementary sequences that could disrupt the formation of native stems on a given molecule—an ambiguity index. Using comparative secondary structures, the ambiguity index is systematically smaller among unbound than bound molecules, as expected. Furthermore, the effect is lost when the presumably more accurate comparative structure is replaced instead by the MFE structure.

Conclusions

A statistical analysis of the relationship between the primary and secondary structures of non-coding RNA molecules suggests that stem-disrupting kinetic traps are substantially less prevalent in molecules not participating in RNP complexes. In that this distinction is apparent under the comparative but not the MFE secondary structure, the results highlight a possible deficiency in structure predictions when based upon assumptions of thermodynamic equilibrium.