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

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

Structural basis of rotavirus RNA chaperone displacement and RNA annealing.

Rotavirus genomes are distributed between 11 distinct RNA molecules, all of which must be selectively copackaged during virus assembly. This likely occurs through sequence-specific RNA interactions facilitated by the RNA chaperone NSP2. Here, we report that NSP2 autoregulates its chaperone activity through its C-terminal region (CTR) that promotes RNA-RNA interactions by limiting its helix-unwinding activity. Unexpectedly, structural proteomics data revealed that the CTR does not directly interact with RNA, while accelerating RNA release from NSP2. Cryo-electron microscopy reconstructions of an NSP2-RNA complex reveal a highly conserved acidic patch on the CTR, which is poised toward the bound RNA. Virus replication was abrogated by charge-disrupting mutations within the acidic patch but completely restored by charge-preserving mutations. Mechanistic similarities between NSP2 and the unrelated bacterial RNA chaperone Hfq suggest that accelerating RNA dissociation while promoting intermolecular RNA interactions may be a widespread strategy of RNA chaperone recycling.

Cis regulation within a cluster of viral microRNAs.

MicroRNAs (miRNAs) are small regulatory RNAs involved in virtually all biological processes. Although many of them are co-expressed from clusters, little is known regarding the impact of this organization on the regulation of their accumulation. In this study, we set to decipher a regulatory mechanism controlling the expression of the ten clustered pre-miRNAs from Kaposi's sarcoma associated herpesvirus (KSHV). We measured in vitro the efficiency of cleavage of each individual pre-miRNA by the Microprocessor and found that pre-miR-K1 and -K3 were the most efficiently cleaved pre-miRNAs. A mutational analysis showed that, in addition to producing mature miRNAs, they are also important for the optimal expression of the whole set of miRNAs. We showed that this feature depends on the presence of a canonical pre-miRNA at this location since we could functionally replace pre-miR-K1 by a heterologous pre-miRNA. Further in vitro processing analysis suggests that the two stem-loops act in cis and that the cluster is cleaved in a sequential manner. Finally, we exploited this characteristic of the cluster to inhibit the expression of the whole set of miRNAs by targeting the pre-miR-K1 with LNA-based antisense oligonucleotides in cells either expressing a synthetic construct or latently infected with KSHV.

Complex Conformational Dynamics of the Heart Failure-Associated Pre-miRNA-377 Hairpin Revealed by Single-Molecule Optical Tweezers.

Pre-miRNA-377 is a hairpin-shaped regulatory RNA associated with heart failure. Here, we use single-molecule optical tweezers to unzip pre-miRNA-377 and study its stability and dynamics. We show that magnesium ions have a strong stabilizing effect, and that sodium ions stabilize the hairpin more than potassium ions. The hairpin unfolds in a single step, regardless of buffer composition. Interestingly, hairpin folding occurs either in a single step (type 1) or through the formation of intermediates, in multiple steps (type 2) or gradually (type 3). Type 3 occurs only in the presence of both sodium and magnesium, while type 1 and 2 take place in all buffers, with type 1 being the most prevalent. By reducing the size of the native hairpin loop from fourteen to four nucleotides, we demonstrate that the folding heterogeneity originates from the large size of the hairpin loop. Further, while efficient pre-miRNA-377 binders are lacking, we demonstrate that the recently developed C2 ligand displays bimodal activity: it enhances the mechanical stability of the pre-miRNA-377 hairpin and perturbs its folding. The knowledge regarding pre-miRNA stability and dynamics that we provide is important in understanding its regulatory function and how it can be modulated to achieve a therapeutic effect, e.g., in heart failure treatment.

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.

Computer-Aided Design of Active Pseudoknotted Hammerhead Ribozymes.

Pseudoknots are important motifs for stabilizing the structure of functional RNAs. As an example, pseudoknotted hammerhead ribozymes are highly active compared to minimal ribozymes. The design of new RNA sequences that retain the function of a model RNA structure includes taking in account pseudoknots presence in the structure, which is usually a challenge for bioinformatics tools. Our method includes using "Enzymer," a software for designing RNA sequences with desired secondary structures that may include pseudoknots. Enzymer implements an efficient stochastic search and optimization algorithm to sample RNA sequences from low ensemble defect mutational landscape of an initial design template to generate an RNA sequence that is predicted to fold into the desired target structure.

Inverse RNA Folding Workflow to Design and Test Ribozymes that Include Pseudoknots.

Ribozymes are RNAs that catalyze reactions. They occur in nature, and can also be evolved in vitro to catalyze novel reactions. This chapter provides detailed protocols for using inverse folding software to design a ribozyme sequence that will fold to a known ribozyme secondary structure and for testing the catalytic activity of the sequence experimentally. This protocol is able to design sequences that include pseudoknots, which is important as all naturally occurring full-length ribozymes have pseudoknots. The starting point is the known pseudoknot-containing secondary structure of the ribozyme and knowledge of any nucleotides whose identity is required for function. The output of the protocol is a set of sequences that have been tested for function. Using this protocol, we were previously successful at designing highly active double-pseudoknotted HDV ribozymes.

Using an L7Ae-Tethered, Hydroxyl Radical-Mediated Footprinting Strategy to Identify and Validate Kink-Turns in RNAs.

Kink-turns are important RNA structural modules that facilitate long-range tertiary interactions and form binding sites for members of the L7Ae family of proteins. Present in a wide variety of functional RNAs, kink-turns play key organizational roles in many RNA-based cellular processes, including translation, modification, and tRNA biogenesis. It is important to determine the contribution of kink-turns to the overall architecture of resident RNAs, as these modules dictate ribonucleoprotein (RNP) assembly and function. This chapter describes a site-directed, hydroxyl radical-mediated footprinting strategy that utilizes L7Ae-tethered chemical nucleases to experimentally validate computationally identified kink-turns in any RNA and under a wide variety of conditions. The work plan described here uses the catalytic RNase P RNA as an example to provide a blueprint for using this footprinting method to map RNA-protein interactions in other RNP complexes.

RNA structure-wide discovery of functional interactions with multiplexed RNA motif library.

Biochemical assays and computational analyses have discovered RNA structures throughout various transcripts. However, the roles of these structures are mostly unknown. Here we develop folded RNA element profiling with structure library (FOREST), a multiplexed affinity assay system to identify functional interactions from transcriptome-wide RNA structure datasets. We generate an RNA structure library by extracting validated or predicted RNA motifs from gene-annotated RNA regions. The RNA structure library with an affinity enrichment assay allows for the comprehensive identification of target-binding RNA sequences and structures in a high-throughput manner. As a proof-of-concept, FOREST discovers multiple RNA-protein interaction networks with quantitative scores, including translational regulatory elements that function in living cells. Moreover, FOREST reveals different binding landscapes of RNA G-quadruplex (rG4) structures-binding proteins and discovers rG4 structures in the terminal loops of precursor microRNAs. Overall, FOREST serves as a versatile platform to investigate RNA structure-function relationships on a large scale.

Structure-Based Mechanisms of a Molecular RNA Polymerase/Chaperone Machine Required for Ribosome Biosynthesis.

Bacterial ribosomal RNAs are synthesized by a dedicated, conserved transcription-elongation complex that transcribes at high rates, shields RNA polymerase from premature termination, and supports co-transcriptional RNA folding, modification, processing, and ribosomal subunit assembly by presently unknown mechanisms. We have determined cryo-electron microscopy structures of complete Escherichia coli ribosomal RNA transcription elongation complexes, comprising RNA polymerase; DNA; RNA bearing an N-utilization-site-like anti-termination element; Nus factors A, B, E, and G; inositol mono-phosphatase SuhB; and ribosomal protein S4. Our structures and structure-informed functional analyses show that fast transcription and anti-termination involve suppression of NusA-stabilized pausing, enhancement of NusG-mediated anti-backtracking, sequestration of the NusG C-terminal domain from termination factor ρ, and the ρ blockade. Strikingly, the factors form a composite RNA chaperone around the RNA polymerase RNA-exit tunnel, which supports co-transcriptional RNA folding and annealing of distal RNA regions. Our work reveals a polymerase/chaperone machine required for biosynthesis of functional ribosomes.

Real-time monitoring of single ZTP riboswitches reveals a complex and kinetically controlled decision landscape.

RNAs begin to fold and function during transcription. Riboswitches undergo cotranscriptional switching in the context of transcription elongation, RNA folding, and ligand binding. To investigate how these processes jointly modulate the function of the folate stress-sensing Fusobacterium ulcerans ZTP riboswitch, we apply a single-molecule vectorial folding (VF) assay in which an engineered superhelicase Rep-X sequentially releases fluorescently labeled riboswitch RNA from a heteroduplex in a 5'-to-3' direction, at ~60 nt s-1 [comparable to the speed of bacterial RNA polymerase (RNAP)]. We demonstrate that the ZTP riboswitch is kinetically controlled and that its activation is favored by slower unwinding, strategic pausing between but not before key folding elements, or a weakened transcription terminator. Real-time single-molecule monitoring captures folding riboswitches in multiple states, including an intermediate responsible for delayed terminator formation. These results show how individual nascent RNAs occupy distinct channels within the folding landscape that controls the fate of the riboswitch.

The Sac10b homolog from Sulfolobus islandicus is an RNA chaperone.

Nucleic acid-binding proteins of the Sac10b family, also known as Alba, are widely distributed in Archaea. However, the physiological roles of these proteins have yet to be clarified. Here, we show that Sis10b, a member of the Sac10b family from the hyperthermophilic archaeon Sulfolobus islandicus, was active in RNA strand exchange, duplex RNA unwinding in vitro and RNA unfolding in a heterologous host cell. This protein exhibited temperature-dependent binding preference for ssRNA over dsRNA and was more efficient in RNA unwinding and RNA unfolding at elevated temperatures. Notably, alanine substitution of a highly conserved basic residue (K) at position 17 in Sis10b drastically reduced the ability of this protein to catalyse RNA strand exchange and RNA unwinding. Additionally, the preferential binding of Sis10b to ssRNA also depended on the presence of K17 or R17. Furthermore, normal growth was restored to a slow-growing Sis10b knockdown mutant by overproducing wild-type Sis10b but not by overproducing K17A in this mutant strain. Our results indicate that Sis10b is an RNA chaperone that likely functions most efficiently at temperatures optimal for the growth of S. islandicus, and K17 is essential for the chaperone activity of the protein.

Same fold, different properties: polarizable molecular dynamics simulations of telomeric and TERRA G-quadruplexes.

DNA and RNA sequences rich in guanine can fold into noncanonical structures called G-quadruplexes (GQs), which exhibit a common stem structure of Hoogsteen hydrogen-bonded guanine tetrads and diverse loop structures. GQ sequence motifs are overrepresented in promoters, origins of replication, telomeres, and untranslated regions in mRNA, suggesting roles in modulating gene expression and preserving genomic integrity. Given these roles and unique aspects of different structures, GQs are attractive targets for drug design, but greater insight into GQ folding pathways and the interactions stabilizing them is required. Here, we performed molecular dynamics simulations to study two bimolecular GQs, a telomeric DNA GQ and the analogous telomeric repeat-containing RNA (TERRA) GQ. We applied the Drude polarizable force field, which we show outperforms the additive CHARMM36 force field in both ion retention and maintenance of the GQ folds. The polarizable simulations reveal that the GQs bind bulk K+ ions differently, and that the TERRA GQ accumulates more K+ ions, suggesting different ion interactions stabilize these structures. Nucleobase dipole moments vary as a function of position and also contribute to ion binding. Finally, we show that the TERRA GQ is more sensitive than the telomeric DNA GQ to water-mediated modulation of ion-induced dipole-dipole interactions.

The Biogenesis of SRP RNA Is Modulated by an RNA Folding Intermediate Attained during Transcription.

The signal recognition particle (SRP), responsible for co-translational protein targeting and delivery to cellular membranes, depends on the native long-hairpin fold of its RNA to confer functionality. Since RNA initiates folding during its synthesis, we used high-resolution optical tweezers to follow in real time the co-transcriptional folding of SRP RNA. Surprisingly, SRP RNA folding is robust to transcription rate changes and the presence or absence of its 5'-precursor sequence. The folding pathway also reveals the obligatory attainment of a non-native hairpin intermediate (H1) that eventually rearranges into the native fold. Furthermore, H1 provides a structural platform alternative to the native fold for RNase P to bind and mature SRP RNA co-transcriptionally. Delays in attaining the final native fold are detrimental to the cell, altogether showing that a co-transcriptional folding pathway underpins the proper biogenesis of function-essential SRP RNA.

METTL15 introduces N4-methylcytidine into human mitochondrial 12S rRNA and is required for mitoribosome biogenesis.

Post-transcriptional RNA modifications, the epitranscriptome, play important roles in modulating the functions of RNA species. Modifications of rRNA are key for ribosome production and function. Identification and characterization of enzymes involved in epitranscriptome shaping is instrumental for the elucidation of the functional roles of specific RNA modifications. Ten modified sites have been thus far identified in the mammalian mitochondrial rRNA. Enzymes responsible for two of these modifications have not been characterized. Here, we identify METTL15, show that it is the main N4-methylcytidine (m4C) methyltransferase in human cells and demonstrate that it is responsible for the methylation of position C839 in mitochondrial 12S rRNA. We show that the lack of METTL15 results in a reduction of the mitochondrial de novo protein synthesis and decreased steady-state levels of protein components of the oxidative phosphorylation system. Without functional METTL15, the assembly of the mitochondrial ribosome is decreased, with the late assembly components being unable to be incorporated efficiently into the small subunit. We speculate that m4C839 is involved in the stabilization of 12S rRNA folding, therefore facilitating the assembly of the mitochondrial small ribosomal subunits. Taken together our data show that METTL15 is a novel protein necessary for efficient translation in human mitochondria.

In silico structure analysis of alphaviral RNA genomes shows diversity in the evasion of IFIT1-mediated innate immunity.

The IFIT (interferon-induced proteins with tetratricopeptide repeats) family constitutes a major arm of the antiviral function of type I interferon (IFN). Human IFIT1, the earliest discovered member of this family, inhibits several viruses of positivestrand RNA genome. IFIT1 specifically recognizes single-stranded RNAwith canonical 7-methylguanylate cap at the 50 end (Cap0), and inhibits their translation by competing with eIF4E (eukaryotic initiation factor 4E), an essential factor for 50Cap recognition. Recently, a novel viral mechanism of IFIT1 suppression was reported, in which an RNA hairpin in the 50 untranslated region (50UTR) of the viral genome prevented recognition by IFIT1 and enhanced virus growth. Here, I have analyzed the in silico predicted structures in the 50UTR of the genomes of the Alphaviruses, a large group of enveloped RNA virus with positive-sense single-stranded genome. The results uncovered a large ensemble of RNA secondary structures of diverse size and shape in the different viruses, which showed little correspondence to the phylogeny of the viruses. Unexpectedly, the 50UTR of several viral genomes in this family did not fold into any structure, suggesting either their extreme sensitivity to IFIT1 or the existence of alternative viral mechanisms of subverting IFIT1 function.

Conformation and mechanical property of rpoS mRNA inhibitory stem studied by optical tweezers and X-ray scattering.

3' downstream inhibitory stem plays a crucial role in locking rpoS mRNA 5' untranslated region in a self-inhibitory state. Here, we used optical tweezers to study the unfolding/refolding of rpoS inhibitory stem in the absence and presence of Mg2+. We found adding Mg2+ decreased the free energy of the RNA junction without re-arranging its secondary structure, through confirming that this RNA formed a canonical RNA three-way junction. We suspected increased free energy might change the relative orientation of different stems of rpoS and confirmed this by small angle X-ray scattering. Such changed conformation may improve Hfq-bridged annealing between sRNA and rpoS RNA inhibitory stem. We established a convenient route to analyze the changes of RNA conformation and folding dynamics by combining optical tweezers with X-ray scattering methods. This route can be easily applied in the studies of other RNA structure and ligand-RNA.

Efficient inhibition of RNA self-primed extension by addition of competing 3'-capture DNA-improved RNA synthesis by T7 RNA polymerase.

In vitro synthesized RNA is used widely in studies of RNA biology, biotechnology and RNA therapeutics. However, in vitro synthesized RNA often contains impurities, such as RNAs with lengths shorter and longer than the expected runoff RNA. We have recently confirmed that longer RNA products are formed predominantly via cis self-primed extension, in which released runoff RNA folds back on itself to prime its own RNA-templated extension. In the current work, we demonstrate that addition of a DNA oligonucleotide (capture DNA) that is complementary to the 3' end of the expected runoff RNA effectively prevents self-primed extension, even under conditions commonly used for high RNA yields. Moreover, the presence of this competing capture DNA during 'high yield' transcription, leads to an increase in the yield of expected runoff RNA by suppressing the formation of undesired longer RNA byproducts.

A pan-cancer analysis of synonymous mutations.

Synonymous mutations have been viewed as silent mutations, since they only affect the DNA and mRNA, but not the amino acid sequence of the resulting protein. Nonetheless, recent studies suggest their significant impact on splicing, RNA stability, RNA folding, translation or co-translational protein folding. Hence, we compile 659194 synonymous mutations found in human cancer and characterize their properties. We provide the user-friendly, comprehensive resource for synonymous mutations in cancer, SynMICdb ( http://SynMICdb.dkfz.de ), which also contains orthogonal information about gene annotation, recurrence, mutation loads, cancer association, conservation, alternative events, impact on mRNA structure and a SynMICdb score. Notably, synonymous and missense mutations are depleted at the 5'-end of the coding sequence as well as at the ends of internal exons independent of mutational signatures. For patient-derived synonymous mutations in the oncogene KRAS, we indicate that single point mutations can have a relevant impact on expression as well as on mRNA secondary structure.