ScanFold: an approach for genome-wide discovery of local RNA structural elements-applications to Zika virus and HIV

A survey of RNA secondary structural propensity encoded within human herpesvirus genomes: global comparisons and local motifs

There are nine herpesviruses known to infect humans, of which Epstein-Barr virus (EBV) is the most widely distributed (>90% of adults infected). This ubiquitous virus is implicated in a variety of cancers and autoimmune diseases. Previous analyses of the EBV genome revealed numerous regions with evidence of generating unusually stable and conserved RNA secondary structures and led to the discovery of a novel class of EBV non-coding (nc)RNAs: the stable intronic sequence (sis)RNAs. To gain a better understanding of the roles of RNA structure in EBV biology and pathogenicity, we revisit EBV using recently developed tools for genome-wide motif discovery and RNA structural characterization. This corroborated previous results and revealed novel motifs with potential functionality; one of which has been experimentally validated. Additionally, since many herpesviruses increasingly rival the seroprevalence of EBV (VZV, HHV-6 and HHV-7 being the most notable), analyses were expanded to include all sequenced human Herpesvirus RefSeq genomes, allowing for genomic comparisons. In total 10 genomes were analyzed, for EBV (types 1 and 2), HCMV, HHV-6A, HHV-6B, HHV-7, HSV-1, HSV-2, KSHV, and VZV. All resulting data were archived in the RNAStructuromeDB (https://structurome.bb.iastate.edu/herpesvirus) to make them available to a wide array of researchers.

RNA genome conservation and secondary structure in SARS-CoV-2 and SARS-related viruses: a first look

As the COVID-19 outbreak spreads, there is a growing need for a compilation of conserved RNA genome regions in the SARS-CoV-2 virus along with their structural propensities to guide development of antivirals and diagnostics. Here we present a first look at RNA sequence conservation and structural propensities in the SARS-CoV-2 genome. Using sequence alignments spanning a range of betacoronaviruses, we rank genomic regions by RNA sequence conservation, identifying 79 regions of length at least 15 nt as exactly conserved over SARS-related complete genome sequences available near the beginning of the COVID-19 outbreak. We then confirm the conservation of the majority of these genome regions across 739 SARS-CoV-2 sequences subsequently reported from the COVID-19 outbreak, and we present a curated list of 30 "SARS-related-conserved" regions. We find that known RNA structured elements curated as Rfam families and in prior literature are enriched in these conserved genome regions, and we predict additional conserved, stable secondary structures across the viral genome. We provide 106 "SARS-CoV-2-conserved-structured" regions as potential targets for antivirals that bind to structured RNA. We further provide detailed secondary structure models for the extended 5' UTR, frameshifting stimulation element, and 3' UTR. Lastly, we predict regions of the SARS-CoV-2 viral genome that have low propensity for RNA secondary structure and are conserved within SARS-CoV-2 strains. These 59 "SARS-CoV-2-conserved-unstructured" genomic regions may be most easily accessible by hybridization in primer-based diagnostic strategies.

Plasmodium falciparum translational machinery condones polyadenosine repeats

Plasmodium falciparum is a causative agent of human malaria. Sixty percent of mRNAs from its extremely AT-rich (81%) genome harbor long polyadenosine (polyA) runs within their ORFs, distinguishing the parasite from its hosts and other sequenced organisms. Recent studies indicate polyA runs cause ribosome stalling and frameshifting, triggering mRNA surveillance pathways and attenuating protein synthesis. Here, we show that P. falciparum is an exception to this rule. We demonstrate that both endogenous genes and reporter sequences containing long polyA runs are efficiently and accurately translated in P. falciparum cells. We show that polyA runs do not elicit any response from No Go Decay (NGD) or result in the production of frameshifted proteins. This is in stark contrast to what we observe in human cells or T. thermophila, an organism with similar AT-content. Finally, using stalling reporters we show that Plasmodium cells evolved not to have a fully functional NGD pathway.

An in silico map of the SARS-CoV-2 RNA Structurome

SARS-CoV-2 is a positive-sense single-stranded RNA virus that has exploded throughout the global human population. This pandemic coronavirus strain has taken scientists and public health researchers by surprise and knowledge of its basic biology (e.g. structure/function relationships in its genomic, messenger and template RNAs) and modes for therapeutic intervention lag behind that of other human pathogens. In this report we used a recently-developed bioinformatics approach, ScanFold, to deduce the RNA structural landscape of the SARS-CoV-2 transcriptome. We recapitulate known elements of RNA structure and provide a model for the folding of an essential frameshift signal. Our results find that the SARS-CoV-2 is greatly enriched in unusually stable and likely evolutionarily ordered RNA structure, which provides a huge reservoir of potential drug targets for RNA-binding small molecules. Our results also predict regions that are accessible for intermolecular interactions, which can aid in the design of antisense therapeutics. All results are made available via a public database (the RNAStructuromeDB) where they may hopefully drive drug discovery efforts to inhibit SARS-CoV-2 pathogenesis.

Translation of the intrinsically disordered protein α-synuclein is inhibited by a small molecule targeting its structured mRNA

Many proteins are refractory to targeting because they lack small-molecule binding pockets. An alternative to drugging these proteins directly is to target the messenger (m)RNA that encodes them, thereby reducing protein levels. We describe such an approach for the difficult-to-target protein α-synuclein encoded by the SNCA gene. Multiplication of the SNCA gene locus causes dominantly inherited Parkinson's disease (PD), and α-synuclein protein aggregates in Lewy bodies and Lewy neurites in sporadic PD. Thus, reducing the expression of α-synuclein protein is expected to have therapeutic value. Fortuitously, the SNCA mRNA has a structured iron-responsive element (IRE) in its 5' untranslated region (5' UTR) that controls its translation. Using sequence-based design, we discovered small molecules that target the IRE structure and inhibit SNCA translation in cells, the most potent of which is named Synucleozid. Both in vitro and cellular profiling studies showed Synucleozid directly targets the α-synuclein mRNA 5' UTR at the designed site. Mechanistic studies revealed that Synucleozid reduces α-synuclein protein levels by decreasing the amount of SNCA mRNA loaded into polysomes, mechanistically providing a cytoprotective effect in cells. Proteome- and transcriptome-wide studies showed that the compound's selectivity makes Synucleozid suitable for further development. Importantly, transcriptome-wide analysis of mRNAs that encode intrinsically disordered proteins revealed that each has structured regions that could be targeted with small molecules. These findings demonstrate the potential for targeting undruggable proteins at the level of their coding mRNAs. This approach, as applied to SNCA, is a promising disease-modifying therapeutic strategy for PD and other α-synucleinopathies.

Mapping the RNA structural landscape of viral genomes

Functional RNA structures are prevalent in viral genomes, and have been shown to play roles in almost every aspect of their biology. However, the majority of viral RNA remains structurally uncharacterized. This is likely to remain true as the cost of sequencing decreases much faster than the cost of structural characterizations. Because of this, there is a need for rapid, inexpensive methods to highlight regions of viral RNA which are ideal candidates for structure-function analyses. The ScanFold method was developed as a single sequence alternative to traditional RNA structural motif discovery pipelines, which rely heavily on well curated sequence alignments to identify conserved RNA structures. ScanFold focuses on identifying (based on their more stable than expected folding energies) the most likely functional structures encoded within a single large RNA sequence, while allowing predicted motifs to be tested for evidence of structural conservation later. Decoupling these processes can be a benefit to researchers studying viruses lacking the ideal phylogenetic depth to yield evidence of structural conservation. Here, we demonstrate how the most significant ScanFold predicted structures correspond to higher base pairing probabilities, SHAPE reactivities, and predict known functional structures within the ZIKV and HIV-1 genomes with accuracy. Best practices and examples are also shown to aid users in utilizing ScanFold for their own systems of interest. ScanFold is available as a Webserver (https://mosslabtools.bb.iastate.edu/scanfold) or can be downloaded (https://github.com/moss-lab/ScanFold) and run locally.

Probing RNA structure in vivo

RNA structure underpins many essential functions in biology. New chemical reagents and techniques for probing RNA structure in living cells have emerged in recent years. High-throughput, genome-wide techniques such as Structure-seq2 and DMS-MaPseq exploit nucleobase modification by dimethylsulfate (DMS) to obtain complete structuromes, and are applicable to multiple domains of life and conditions. New reagents such as 1-ethyl-3-(3-dimethylaminopropyl)carbodiimide (EDC), glyoxal, and nicotinoyl azide (NAz) greatly expand the capabilities of nucleobase probing in cells. Additionally, ribose-targeting reagents in selective 2'-hydroxyl acylation and primer extension (SHAPE) detect RNA flexibility in vivo. These techniques, coupled with crosslinking nucleobases in psoralen analysis of RNA interactions and structures (PARIS), provide new and diverse ways to elucidate RNA secondary and tertiary structure in vivo and genome-wide.

RNA structural analysis of the MYC mRNA reveals conserved motifs that affect gene expression

The MYC gene encodes a human transcription factor and proto-oncogene that is dysregulated in over half of all known cancers. To better understand potential post-transcriptional regulatory features affecting MYC expression, we analyzed secondary structures in the MYC mRNA using a program that is optimized for finding small locally-folded motifs with a high propensity for function. This was accomplished by calculating folding metrics across the MYC sequence using a sliding analysis window and generating unique consensus base pairing models weighted by their lower-than-random predicted folding energy. A series of 30 motifs were identified, primarily in the 5' and 3' untranslated regions, which show evidence of structural conservation and compensating mutations across vertebrate MYC homologs. This analysis was able to recapitulate known elements found within an internal ribosomal entry site, as well as discover a novel element in the 3' UTR that is unusually stable and conserved. This novel motif was shown to affect MYC expression, potentially via the modulation of miRNA target accessibility or other trans-regulatory factors. In addition to providing basic insights into mechanisms that regulate MYC expression, this study provides numerous, potentially druggable RNA targets for the MYC gene, which is considered “undruggable” at the protein level.

Computational approaches for the discovery of splicing regulatory RNA structures

Global RNA structure and local functional motifs mediate interactions important in determining the rates and patterns of mRNA splicing. In this review, we overview approaches for the computational prediction of RNA secondary structure with a special emphasis on the discovery of motifs important to RNA splicing. The process of identifying and modeling potential splicing regulatory structures is illustrated using a recently-developed approach for RNA structural motif discovery, the ScanFold pipeline, which is applied to the identification of a known splicing regulatory structure in influenza virus.