Discovery of a large-scale, cell-state-responsive allosteric switch in the 7SK RNA using DANCE-MaP.

7SK is a conserved noncoding RNA that regulates transcription by sequestering the transcription factor P-TEFb. 7SK function entails complex changes in RNA structure, but characterizing RNA dynamics in cells remains an unsolved challenge. We developed a single-molecule chemical probing strategy, DANCE-MaP (deconvolution and annotation of ribonucleic conformational ensembles), that defines per-nucleotide reactivity, direct base pairing interactions, tertiary interactions, and thermodynamic populations for each state in RNA structural ensembles from a single experiment. DANCE-MaP reveals that 7SK RNA encodes a large-scale structural switch that couples dissolution of the P-TEFb binding site to structural remodeling at distal release factor binding sites. The 7SK structural equilibrium shifts in response to cell growth and stress and can be targeted to modulate expression of P-TEFbresponsive genes. Our study reveals that RNA structural dynamics underlie 7SK function as an integrator of diverse cellular signals to control transcription and establishes the power of DANCE-MaP to define RNA dynamics in cells.

Genomic RNA Elements Drive Phase Separation of the SARS-CoV-2 Nucleocapsid.

We report that the SARS-CoV-2 nucleocapsid protein (N-protein) undergoes liquid-liquid phase separation (LLPS) with viral RNA. N-protein condenses with specific RNA genomic elements under physiological buffer conditions and condensation is enhanced at human body temperatures (33°C and 37°C) and reduced at room temperature (22°C). RNA sequence and structure in specific genomic regions regulate N-protein condensation while other genomic regions promote condensate dissolution, potentially preventing aggregation of the large genome. At low concentrations, N-protein preferentially crosslinks to specific regions characterized by single-stranded RNA flanked by structured elements and these features specify the location, number, and strength of N-protein binding sites (valency). Liquid-like N-protein condensates form in mammalian cells in a concentration-dependent manner and can be altered by small molecules. Condensation of N-protein is RNA sequence and structure specific, sensitive to human body temperature, and manipulatable with small molecules, and therefore presents a screenable process for identifying antiviral compounds effective against SARS-CoV-2.

Template switching enables chemical probing of native RNA structures.

RNAs are often studied in non-native sequence contexts to facilitate structural studies. However, seemingly innocuous changes to an RNA sequence may perturb the native structure and generate inaccurate or ambiguous structural models. To facilitate the investigation of native RNA secondary structure by selective 2' hydroxyl acylation analyzed by primer extension (SHAPE), we engineered an approach that couples minimal enzymatic steps to RNA chemical probing and mutational profiling (MaP) reverse transcription (RT) methods - a process we call template switching and mutational profiling (Switch-MaP). In Switch-MaP, RT templates and additional library sequences are added post-probing through ligation and template switching, capturing reactivities for every nucleotide. For a candidate SAM-I riboswitch, we compared RNA structure models generated by the Switch-MaP approach to those of traditional primer-based MaP, including RNAs with or without appended structure cassettes. Primer-based MaP masked reactivity data in the 5' and 3' ends of the RNA, producing ambiguous ensembles inconsistent with the conserved SAM-I riboswitch secondary structure. Structure cassettes enabled unambiguous modeling of an aptamer-only construct but introduced non-native interactions in the full length riboswitch. In contrast, Switch-MaP provided reactivity data for all nucleotides in each RNA and enabled unambiguous modeling of secondary structure, consistent with the conserved SAM-I fold. Switch-MaP is a straightforward alternative approach to primer-based and cassette-based chemical probing methods that precludes primer masking and the formation of alternative secondary structures due to non-native sequence elements.

The TRIM-NHL RNA-binding protein Brain Tumor coordinately regulates expression of the glycolytic pathway and vacuolar ATPase complex.

The essential Drosophila RNA-binding protein Brain Tumor (Brat) represses specific genes to control embryogenesis and differentiation of stem cells. In the brain, Brat functions as a tumor suppressor that diminishes neural stem cell proliferation while promoting differentiation. Though important Brat-regulated target mRNAs have been identified in these contexts, the full impact of Brat on gene expression remains to be discovered. Here, we identify the network of Brat-regulated mRNAs by performing RNA sequencing (RNA-seq) following depletion of Brat from cultured cells. We identify 158 mRNAs, with high confidence, that are repressed by Brat. De novo motif analysis identified a functionally enriched RNA motif in the 3' untranslated regions (UTRs) of Brat-repressed mRNAs that matches the biochemically defined Brat binding site. Integrative data analysis revealed a high-confidence list of Brat-repressed and Brat-bound mRNAs containing 3'UTR Brat binding motifs. Our RNA-seq and reporter assays show that multiple 3'UTR motifs promote the strength of Brat repression, whereas motifs in the 5'UTR are not functional. Strikingly, we find that Brat regulates expression of glycolytic enzymes and the vacuolar ATPase complex, providing new insight into its role as a tumor suppressor and the coordination of metabolism and intracellular pH.

Structural analysis of MALAT1 long noncoding RNA in cells and in evolution.

Although not canonically polyadenylated, the long noncoding RNA MALAT1 (metastasis-associated lung adenocarcinoma transcript 1) is stabilized by a highly conserved 76-nt triple helix structure on its 3' end. The entire MALAT1 transcript is over 8000 nt long in humans. The strongest structural conservation signal in MALAT1 (as measured by covariation of base pairs) is in the triple helix structure. Primary sequence analysis of covariation alone does not reveal the degree of structural conservation of the entire full-length transcript, however. Furthermore, RNA structure is often context dependent; RNA binding proteins that are differentially expressed in different cell types may alter structure. We investigate here the in-cell and cell-free structures of the full-length human and green monkey (Chlorocebus sabaeus) MALAT1 transcripts in multiple tissue-derived cell lines using SHAPE chemical probing. Our data reveal levels of uniform structural conservation in different cell lines, in cells and cell-free, and even between species, despite significant differences in primary sequence. The uniformity of the structural conservation across the entire transcript suggests that, despite seeing covariation signals only in the triple helix junction of the lncRNA, the rest of the transcript's structure is remarkably conserved, at least in primates and across multiple cell types and conditions.

Functional Validation of SAM Riboswitch Element A from Listeria monocytogenes.

SreA is one of seven candidate S-adenosyl methionine (SAM) class I riboswitches identified in Listeria monocytogenes, a saprophyte and opportunistic foodborne pathogen. SreA precedes genes encoding a methionine ATP-binding cassette (ABC) transporter, which imports methionine and is presumed to regulate transcription of its downstream genes in a SAM-dependent manner. The proposed role of SreA in controlling the transcription of genes encoding an ABC transporter complex may have important implications for how the bacteria senses and responds to the availability of the metabolite SAM in the diverse environments in which L. monocytogenes persists. Here we validate SreA as a functional SAM-I riboswitch through ligand binding studies, structure characterization, and transcription termination assays. We determined that SreA has both a structure and SAM binding properties similar to those of other well-characterized SAM-I riboswitches. Despite the apparent structural similarities to previously described SAM-I riboswitches, SreA induces transcription termination in response to comparatively lower (nanomolar) ligand concentrations. Furthermore, SreA is a leaky riboswitch that permits some transcription of the downstream gene even in the presence of millimolar SAM, suggesting that L. monocytogenes may "dampen" the expression of genes for methionine import but likely does not turn them "OFF".

Single-Molecule Correlated Chemical Probing: A Revolution in RNA Structure Analysis.

RNA molecules convey biological information both in their linear sequence and in their base-paired secondary and tertiary structures. Chemical probing experiments, which involve treating an RNA with a reagent that modifies conformationally dynamic nucleotides, have broadly enabled examination of short- and long-range RNA structure in diverse contexts, including in living cells. For decades, chemical probing experiments have been interpreted in a per-nucleotide way, such that the reactivity measured at each nucleotide reports the average structure at a position over all RNA molecules within a sample. However, there are numerous important cases where per-nucleotide chemical probing falls short, including for RNAs that are bound by proteins, RNAs that form complex higher order structures, and RNAs that sample multiple conformations.Recent experimental and computational innovations have started a revolution in RNA structure analysis by transforming chemical probing into a massively parallel, single-molecule experiment. Enabled by a specialized reverse transcription strategy called mutational profiling (MaP), multiple chemical modification events can be measured within individual RNA molecules. Nucleotides that communicate structurally through direct base pairing or large-scale folding-unfolding transitions will react with chemical probes in a correlated manner, thereby revealing structural complexity hidden to conventional approaches. These single-molecule correlated chemical probing (smCCP) experiments can be interpreted to directly identify nucleotides that base pair (the PAIR-MaP strategy) and to reveal long-range, through-space structural communication (RING-MaP). Correlated probing can also define the thermodynamic populations of complex RNA ensembles (DANCE-MaP). Complex RNA-protein networks can be interrogated by cross-linking proteins to RNA and measuring correlations between cross-linked positions (RNP-MaP).smCCP thus visualizes RNA secondary and higher-order structure with unprecedented accuracy, defining novel structures, RNA-protein interaction networks, time-resolved dynamics, and allosteric structural switches. These strategies are not mutually exclusive; in favorable cases, multiple levels of RNA structure ─ base pairing, through-space structural communication, and equilibrium ensembles ─ can be resolved concurrently. The physical experimentation required for smCCP is profoundly simple, and experiments are readily performed in cells on RNAs of any size, including large noncoding RNAs and mRNAs. Single-molecule correlated chemical probing is paving the way for a new generation of biophysical studies on RNA in living systems.

Global 5'-UTR RNA structure regulates translation of a SERPINA1 mRNA.

SERPINA1 mRNAs encode the protease inhibitor α-1-antitrypsin and are regulated through post-transcriptional mechanisms. α-1-antitrypsin deficiency leads to chronic obstructive pulmonary disease (COPD) and liver cirrhosis, and specific variants in the 5'-untranslated region (5'-UTR) are associated with COPD. The NM_000295.4 transcript is well expressed and translated in lung and blood and features an extended 5'-UTR that does not contain a competing upstream open reading frame (uORF). We show that the 5'-UTR of NM_000295.4 folds into a well-defined multi-helix structural domain. We systematically destabilized mRNA structure across the NM_000295.4 5'-UTR, and measured changes in (SHAPE quantified) RNA structure and cap-dependent translation relative to a native-sequence reporter. Surprisingly, despite destabilizing local RNA structure, most mutations either had no effect on or decreased translation. Most structure-destabilizing mutations retained native, global 5'-UTR structure. However, those mutations that disrupted the helix that anchors the 5'-UTR domain yielded three groups of non-native structures. Two of these non-native structure groups refolded to create a stable helix near the translation initiation site that decreases translation. Thus, in contrast to the conventional model that RNA structure in 5'-UTRs primarily inhibits translation, complex folding of the NM_000295.4 5'-UTR creates a translation-optimized message by promoting accessibility at the translation initiation site.

Direct Duplex Detection: An Emerging Tool in the RNA Structure Analysis Toolbox.

While a variety of powerful tools exists for analyzing RNA structure, identifying long-range and intermolecular base-pairing interactions has remained challenging. Recently, three groups introduced a high-throughput strategy that uses psoralen-mediated crosslinking to directly identify RNA-RNA duplexes in cells. Initial application of these methods highlights the preponderance of long-range structures within and between RNA molecules and their widespread structural dynamics.

Guidelines for SHAPE Reagent Choice and Detection Strategy for RNA Structure Probing Studies.

Chemical probing is an important tool for characterizing the complex folded structures of RNA molecules, many of which play key cellular roles. Electrophilic SHAPE reagents create adducts at the 2'-hydroxyl position on the RNA backbone of flexible ribonucleotides with relatively little dependence on nucleotide identity. Strategies for adduct detection such as mutational profiling (MaP) allow accurate, automated calculation of relative adduct frequencies for each nucleotide in a given RNA or group of RNAs. A number of alternative reagents and adduct detection strategies have been proposed, especially for use in living cells. Here we evaluate five SHAPE reagents: three previously well-validated reagents 1M7 (1-methyl-7-nitroisatoic anhydride), 1M6 (1-methyl-6-nitroisatoic anhydride), and NMIA ( N-methylisatoic anhydride), one more recently proposed NAI (2-methylnicotinic acid imidazolide), and one novel reagent 5NIA (5-nitroisatoic anhydride). We clarify the importance of carefully designed software in reading out SHAPE experiments using massively parallel sequencing approaches. We examine SHAPE modification in living cells in diverse cell lines, compare MaP and reverse transcription-truncation as SHAPE adduct detection strategies, make recommendations for SHAPE reagent choice, and outline areas for future development.

SHAPE Probing Reveals Human rRNAs Are Largely Unfolded in Solution.

Chemical probing experiments, coupled with empirically determined free energy change relationships, can enable accurate modeling of the secondary structures of diverse and complex RNAs. A current frontier lies in modeling large and structurally heterogeneous transcripts, including complex eukaryotic RNAs. To validate and improve on experimentally driven approaches for modeling large transcripts, we obtained high-quality SHAPE data for the protein-free human 18S and 28S ribosomal RNAs (rRNAs). To our surprise, SHAPE-directed structure models for the human rRNAs poorly matched accepted structures. Analysis of predicted rRNA structures based on low-SHAPE and low-entropy (lowSS) metrics revealed that, whereas ∼75% of Escherichia coli rRNA sequences form well-determined lowSS secondary structure, only ∼40% of the human rRNAs do. Critically, regions of the human rRNAs that specifically fold into well-determined lowSS structures were modeled to high accuracy using SHAPE data. This work reveals that eukaryotic rRNAs are more unfolded than are those of prokaryotic rRNAs and indeed are largely unfolded overall, likely reflecting increased protein dependence for eukaryotic ribosome structure. In addition, those regions and substructures that are well-determined can be identified de novo and successfully modeled by SHAPE-directed folding.

Integrated analysis of RNA-binding protein complexes using in vitro selection and high-throughput sequencing and sequence specificity landscapes (SEQRS).

RNA-binding proteins (RBPs) collaborate to control virtually every aspect of RNA function. Tremendous progress has been made in the area of global assessment of RBP specificity using next-generation sequencing approaches both in vivo and in vitro. Understanding how protein-protein interactions enable precise combinatorial regulation of RNA remains a significant problem. Addressing this challenge requires tools that can quantitatively determine the specificities of both individual proteins and multimeric complexes in an unbiased and comprehensive way. One approach utilizes in vitro selection, high-throughput sequencing, and sequence-specificity landscapes (SEQRS). We outline a SEQRS experiment focused on obtaining the specificity of a multi-protein complex between Drosophila RBPs Pumilio (Pum) and Nanos (Nos). We discuss the necessary controls in this type of experiment and examine how the resulting data can be complemented with structural and cell-based reporter assays. Additionally, SEQRS data can be integrated with functional genomics data to uncover biological function. Finally, we propose extensions of the technique that will enhance our understanding of multi-protein regulatory complexes assembled onto RNA.

Combinatorial control of messenger RNAs by Pumilio, Nanos and Brain Tumor Proteins.

Eukaryotes possess a vast array of RNA-binding proteins (RBPs) that affect mRNAs in diverse ways to control protein expression. Combinatorial regulation of mRNAs by RBPs is emerging as the rule. No example illustrates this as vividly as the partnership of 3 Drosophila RBPs, Pumilio, Nanos and Brain Tumor, which have overlapping functions in development, stem cell maintenance and differentiation, fertility and neurologic processes. Here we synthesize 30 y of research with new insights into their molecular functions and mechanisms of action. First, we provide an overview of the key properties of each RBP. Next, we present a detailed analysis of their collaborative regulatory mechanism using a classic example of the developmental morphogen, hunchback, which is spatially and temporally regulated by the trio during embryogenesis. New biochemical, structural and functional analyses provide insights into RNA recognition, cooperativity, and regulatory mechanisms. We integrate these data into a model of combinatorial RNA binding and regulation of translation and mRNA decay. We then use this information, transcriptome wide analyses and bioinformatics predictions to assess the global impact of Pumilio, Nanos and Brain Tumor on gene regulation. Together, the results support pervasive, dynamic post-transcriptional control.

The RNA binding domain of Pumilio antagonizes poly-adenosine binding protein and accelerates deadenylation.

PUF proteins are potent repressors that serve important roles in stem cell maintenance, neurological processes, and embryonic development. These functions are driven by PUF protein recognition of specific binding sites within the 3' untranslated regions of target mRNAs. In this study, we investigated mechanisms of repression by the founding PUF, Drosophila Pumilio, and its human orthologs. Here, we evaluated a previously proposed model wherein the Pumilio RNA binding domain (RBD) binds Argonaute, which in turn blocks the translational activity of the eukaryotic elongation factor 1A. Surprisingly, we found that Argonautes are not necessary for repression elicited by Drosophila and human PUFs in vivo. A second model proposed that the RBD of Pumilio represses by recruiting deadenylases to shorten the mRNA's polyadenosine tail. Indeed, the RBD binds to the Pop2 deadenylase and accelerates deadenylation; however, this activity is not crucial for regulation. Rather, we determined that the poly(A) is necessary for repression by the RBD. Our results reveal that poly(A)-dependent repression by the RBD requires the poly(A) binding protein, pAbp. Furthermore, we show that repression by the human PUM2 RBD requires the pAbp ortholog, PABPC1. Pumilio associates with pAbp but does not disrupt binding of pAbp to the mRNA. Taken together, our data support a model wherein the Pumilio RBD antagonizes the ability of pAbp to promote translation. Thus, the conserved function of the PUF RBD is to bind specific mRNAs, antagonize pAbp function, and promote deadenylation.

Measuring Proximity-Mediated Function of mRNA Regulatory Proteins by Engineered Tethering.

A powerful approach for studying the functional consequences of site-specific RNA-protein interactions is to artificially tether a protein to a messenger (or noncoding) RNA through a selective, high-affinity interaction. We share a strategy for evaluating the contribution of protein positioning within an mRNA on gene expression. We introduced an RNA hairpin recognition site for the MS2 coat protein into the untranslated regions or coding sequence of mRNAs expressing a luminescent reporter protein, NanoLuc. Effector proteins fused to the MS2 coat protein could thus be targeted to distinct regions across the mRNA. We illustrate this approach using ZFP36L2, which recruits the CCR4-NOT complex for poly(A) tail deadenylation. Tethering ZFP36L2 to the 3'-UTR decreased NanoLuc expression, as expected, given the known interaction of this adapter protein with adenine uridine-rich elements (AREs). Intriguingly, ZFP36L2 also decreased NanoLuc expression when bound within the coding sequence, revealing that ZFP36L2-and potentially many other mRNA regulatory proteins-can function when targeted to diverse locations within an mRNA. This multi-target tethering strategy enables exploration of the interplay between mRNA-protein proximity and gene expression.

Human Pumilio proteins recruit multiple deadenylases to efficiently repress messenger RNAs.

PUF proteins are a conserved family of eukaryotic RNA-binding proteins that regulate specific mRNAs: they control many processes including stem cell proliferation, fertility, and memory formation. PUFs repress protein expression from their target mRNAs but the mechanism by which they do so remains unclear, especially for humans. Humans possess two PUF proteins, PUM1 and PUM2, which exhibit similar RNA binding specificities. Here we report new insights into their regulatory activities and mechanisms of action. We developed functional assays to measure sequence-specific repression by PUM1 and PUM2. Both robustly inhibit translation and promote mRNA degradation. Purified PUM complexes were found to contain subunits of the CCR4-NOT (CNOT) complex, which contains multiple enzymes that catalyze mRNA deadenylation. PUMs interact with the CNOT deadenylase subunits in vitro. We used three approaches to determine the importance of deadenylases for PUM repression. First, dominant-negative mutants of CNOT7 and CNOT8 reduced PUM repression. Second, RNA interference depletion of the deadenylases alleviated PUM repression. Third, the poly(A) tail was necessary for maximal PUM repression. These findings demonstrate a conserved mechanism of PUF-mediated repression via direct recruitment of the CCR4-POP2-NOT deadenylase leading to translational inhibition and mRNA degradation. A second, deadenylation independent mechanism was revealed by the finding that PUMs repress an mRNA that lacks a poly(A) tail. Thus, human PUMs are repressors capable of deadenylation-dependent and -independent modes of repression.

Distinct MUNC lncRNA structural domains regulate transcription of different promyogenic factors.

Many lncRNAs have been discovered using transcriptomic data; however, it is unclear what fraction of lncRNAs is functional and what structural properties affect their phenotype. MUNC lncRNA (also known as DRReRNA) acts as an enhancer RNA for the Myod1 gene in cis and stimulates the expression of other promyogenic genes in trans by recruiting the cohesin complex. Here, experimental probing of the RNA structure revealed that MUNC contains multiple structural domains not detected by prediction algorithms in the absence of experimental information. We show that these specific and structurally distinct domains are required for induction of promyogenic genes, for binding genomic sites and gene expression regulation, and for binding the cohesin complex. Myod1 induction and cohesin interaction comprise only a subset of MUNC phenotype. Our study reveals unexpectedly complex, structure-driven functions for the MUNC lncRNA and emphasizes the importance of experimentally determined structures for understanding structure-function relationships in lncRNAs.

Analysis of RNA-protein networks with RNP-MaP defines functional hubs on RNA.

RNA-protein interaction networks govern many biological processes but are difficult to examine comprehensively. We devised ribonucleoprotein networks analyzed by mutational profiling (RNP-MaP), a live-cell chemical probing strategy that maps cooperative interactions among multiple proteins bound to single RNA molecules at nucleotide resolution. RNP-MaP uses a hetero-bifunctional crosslinker to freeze interacting proteins in place on RNA and then maps multiple bound proteins on single RNA strands by read-through reverse transcription and DNA sequencing. RNP-MaP revealed that RNase P and RMRP, two sequence-divergent but structurally related non-coding RNAs, share RNP networks and that network hubs define functional sites in these RNAs. RNP-MaP also identified protein interaction networks conserved between mouse and human XIST long non-coding RNAs and defined protein communities whose binding sites colocalize and form networks in functional regions of XIST. RNP-MaP enables discovery and efficient validation of functional protein interaction networks on long RNAs in living cells.

Targeting the Oncogenic Long Non-coding RNA SLNCR1 by Blocking Its Sequence-Specific Binding to the Androgen Receptor.

Long non-coding RNAs (lncRNAs) are critical regulators of numerous physiological processes and diseases, especially cancers. However, development of lncRNA-based therapies is limited because the mechanisms of many lncRNAs are obscure, and interactions with functional partners, including proteins, remain uncharacterized. The lncRNA SLNCR1 binds to and regulates the androgen receptor (AR) to mediate melanoma invasion and proliferation in an androgen-independent manner. Here, we use biochemical analyses coupled with selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) RNA structure probing to show that the N-terminal domain of AR binds a pyrimidine-rich motif in an unstructured region of SLNCR1. This motif is predictive of AR binding, as we identify an AR-binding motif in lncRNA HOXA11-AS-203. Oligonucleotides that bind either the AR N-terminal domain or the AR RNA motif block the SLNCR1-AR interaction and reduce SLNCR1-mediated melanoma invasion. Delivery of oligos that block SLNCR1-AR interaction thus represent a plausible therapeutic strategy.