Secondary Structure Prediction of Single Sequences Using RNAstructure

How Parameters Influence SHAPE-Directed Predictions

The structure of an RNA sequence encodes information about its biological function. Dynamic programming algorithms are often used to predict the conformation of an RNA molecule from its sequence alone, and adding experimental data as auxiliary information improves prediction accuracy. This auxiliary data is typically incorporated into the nearest neighbor thermodynamic model22 by converting the data into pseudoenergies. Here, we look at how much of the space of possible structures auxiliary data allows prediction methods to explore. We find that for a large class of RNA sequences, auxiliary data shifts the predictions significantly. Additionally, we find that predictions are highly sensitive to the parameters which define the auxiliary data pseudoenergies. In fact, the parameter space can typically be partitioned into regions where different structural predictions predominate.

Rapid and Quantitative Detection of Lung Cancer Biomarker ENOX2 Using a Novel Aptamer in an Electrochemical DNA-Based (E-DNA) Biosensor

To overcome early cancer detection challenges, diagnostic tools enabling more sensitive, rapid, and noninvasive detection are necessary. An attractive cancer target for diagnostic blood tests is human Ecto-NOX disulfide-thiol exchanger 2 (ENOX2), expressed in most human cancer types and regularly shed into blood sera. Here, we developed an electrochemical DNA-based (E-DNA) biosensor that rapidly detects physiologically relevant levels of ENOX2. To identify ENOX2-binding aptamers that could potentially be used in a biosensor, recombinantly expressed ENOX2 was used as a binding target in an oligonucleotide library pull-down that generated a highly enriched ENOX2-binding aptamer. This candidate aptamer sensitively bound ENOX2 via gel mobility shift assays. To enable this aptamer to function in an ENOX2 E-DNA biosensor, the aptamer sequence was modified to adopt two conformations, one capable of ENOX2 binding, and one with disrupted ENOX2 binding. Upon ENOX2 introduction, a conformational shift to the ENOX2 binding state resulted in changed dynamics of a redox reporter molecule, which generated a rapid, significant, and target-specific electrical current readout change. ENOX2 biosensor sensitivity was at or below the diagnostic range. The ENOX2 E-DNA biosensor design presented here may enable the development of more sensitive, rapid, diagnostic tools for early cancer detection.

The Cynosure of CtBP: Evolution of a Bilaterian Transcriptional Corepressor

Evolution of sequence-specific transcription factors clearly drives lineage-specific innovations, but less is known about how changes in the central transcriptional machinery may contribute to evolutionary transformations. In particular, transcriptional regulators are rich in intrinsically disordered regions that appear to be magnets for evolutionary innovation. The C-terminal Binding Protein (CtBP) is a transcriptional corepressor derived from an ancestral lineage of alpha hydroxyacid dehydrogenases; it is found in mammals and invertebrates, and features a core NAD-binding domain as well as an unstructured C-terminus (CTD) of unknown function. CtBP can act on promoters and enhancers to repress transcription through chromatin-linked mechanisms. Our comparative phylogenetic study shows that CtBP is a bilaterian innovation whose CTD of about 100 residues is present in almost all orthologs. CtBP CTDs contain conserved blocks of residues and retain a predicted disordered property, despite having variations in the primary sequence. Interestingly, the structure of the C-terminus has undergone radical transformation independently in certain lineages including flatworms and nematodes. Also contributing to CTD diversity is the production of myriad alternative RNA splicing products, including the production of "short" tailless forms of CtBP in Drosophila. Additional diversity stems from multiple gene duplications in vertebrates, where up to five CtBP orthologs have been observed. Vertebrate lineages show fewer major modifications in the unstructured CTD, possibly because gene regulatory constraints of the vertebrate body plan place specific constraints on this domain. Our study highlights the rich regulatory potential of this previously unstudied domain of a central transcriptional regulator.

Sense and antisense RNA products of the uxuR gene can affect motility and chemotaxis acting independent of the UxuR protein

Small non-coding and antisense RNAs are widespread in all kingdoms of life, however, the diversity of their functions in bacteria is largely unknown. Here, we study RNAs synthesised from divergent promoters located in the 3'-end of the uxuR gene, encoding transcription factor regulating hexuronate metabolism in Escherichia coli. These overlapping promoters were predicted in silico with rather high scores, effectively bound RNA polymerase in vitro and in vivo and were capable of initiating transcription in sense and antisense directions. The genome-wide correlation between in silico promoter scores and RNA polymerase binding in vitro and in vivo was higher for promoters located on the antisense strands of the genes, however, sense promoters within the uxuR gene were more active. Both regulatory RNAs synthesised from the divergent promoters inhibited expression of genes associated with the E. coli motility and chemotaxis independent of a carbon source on which bacteria had been grown. Direct effects of these RNAs were confirmed for the fliA gene encoding σ²⁸ subunit of RNA polymerase. In addition to intracellular sRNAs, promoters located within the uxuR gene could initiate synthesis of transcripts found in the fraction of RNAs secreted in the extracellular medium. Their profile was also carbon-independent suggesting that intragenic uxuR transcripts have a specific regulatory role not directly related to the function of the protein in which gene they are encoded.

Sensitive and reproducible cell-free methylome quantification with synthetic spike-in controls

Cell-free methylated DNA immunoprecipitation sequencing (cfMeDIP-seq) identifies genomic regions with DNA methylation, using a protocol adapted to work with low-input DNA samples and with cell-free DNA (cfDNA). We developed a set of synthetic spike-in DNA controls for cfMeDIP-seq to provide a simple and inexpensive reference for quantitative normalization. We designed 54 DNA fragments with combinations of methylation status (methylated and unmethylated), fragment length (80 bp, 160 bp, 320 bp), G + C content (35%, 50%, 65%), and fraction of CpG dinucleotides within the fragment (1/80 bp, 1/40 bp, 1/20 bp). Using 0.01 ng of spike-in controls enables training a generalized linear model that absolutely quantifies methylated cfDNA in MeDIP-seq experiments. It mitigates batch effects and corrects for biases in enrichment due to known biophysical properties of DNA fragments and other technical biases.

A map of the SARS-CoV-2 RNA structurome

SARS-CoV-2 has exploded throughout the human population. To facilitate efforts to gain insights into SARS-CoV-2 biology and to target the virus therapeutically, it is essential to have a roadmap of likely functional regions embedded in its RNA genome. In this report, we used a bioinformatics approach, ScanFold, to deduce the local RNA structural landscape of the SARS-CoV-2 genome with the highest likelihood of being functional. We recapitulate previously-known elements of RNA structure and provide a model for the folding of an essential frameshift signal. Our results find that SARS-CoV-2 is greatly enriched in unusually stable and likely evolutionarily ordered RNA structure, which provides a large reservoir of potential drug targets for RNA-binding small molecules. Results are enhanced via the re-analyses of publicly-available genome-wide biochemical structure probing datasets that are broadly in agreement with our models. Additionally, ScanFold was updated to incorporate experimental data as constraints in the analysis to facilitate comparisons between ScanFold and other RNA modelling approaches. Ultimately, ScanFold was able to identify eight highly structured/conserved motifs in SARS-CoV-2 that agree with experimental data, without explicitly using these data. All results are made available via a public database (the RNAStructuromeDB: https://structurome.bb.iastate.edu/sars-cov-2) and model comparisons are readily viewable at https://structurome.bb.iastate.edu/sars-cov-2-global-model-comparisons.

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.

Expansion segments in bacterial and archaeal 5S ribosomal RNAs

The large ribosomal RNAs of eukaryotes frequently contain expansion sequences that add to the size of the rRNAs but do not affect their overall structural layout and are compatible with major ribosomal function as an mRNA translation machine. The expansion of prokaryotic ribosomal RNAs is much less explored. In order to obtain more insight into the structural variability of these conserved molecules, we herein report the results of a comprehensive search for the expansion sequences in prokaryotic 5S rRNAs. Overall, 89 expanded 5S rRNAs of 15 structural types were identified in 15 archaeal and 36 bacterial genomes. Expansion segments ranging in length from 13 to 109 residues were found to be distributed among 17 insertion sites. The strains harboring the expanded 5S rRNAs belong to the bacterial orders Clostridiales, Halanaerobiales, Thermoanaerobacterales, and Alteromonadales as well as the archael order Halobacterales When several copies of a 5S rRNA gene are present in a genome, the expanded versions may coexist with normal 5S rRNA genes. The insertion sequences are typically capable of forming extended helices, which do not seemingly interfere with folding of the conserved core. The expanded 5S rRNAs have largely been overlooked in 5S rRNA databases.

Selection of High-Affinity RNA Aptamers That Distinguish between Doxycycline and Tetracycline

Oligonucleotide aptamers are found in prokaryotes and eukaryotes, and they can be selected from large synthetic libraries to bind protein or small-molecule ligands with high affinities and specificities. Aptamers can function as biosensors, as protein recognition elements, and as components of riboswitches allowing ligand-dependent control of gene expression. One of the best studied laboratory-selected aptamers binds the antibiotic tetracycline, but it binds with a much lower affinity to the closely related but more bioavailable antibiotic doxycycline. Here we report enrichment of doxycycline binding aptamers from a selectively randomized library of tetracycline aptamer variants over four selection rounds. Selected aptamers distinguish between doxycycline, which they bind with dissociation constants of approximately 7 nM, and tetracycline, which they bind undetectably. They thus function as orthogonal complements to the original tetracycline aptamer. Unexpectedly, doxycycline aptamers adopt a conformation distinct from that of the tetracycline aptamer and depend on constant regions originally installed as primer binding sites. We show that the fluorescence emission intensity of doxycycline increases upon aptamer binding, permitting their use as biosensors. This new class of aptamers can be used in multiple contexts where doxycycline detection, or doxycycline-mediated regulation, is necessary.

Excessive excision of correct nucleotides during DNA synthesis explained by replication hurdles

The proofreading exonuclease activity of replicative DNA polymerase excises misincorporated nucleotides during DNA synthesis, but these events are rare. Therefore, we were surprised to find that T7 replisome excised nearly 7% of correctly incorporated nucleotides during leading and lagging strand syntheses. Similar observations with two other DNA polymerases establish its generality. We show that excessive excision of correctly incorporated nucleotides is not due to events such as processive degradation of nascent DNA or spontaneous partitioning of primer‐end to the exonuclease site as a “cost of proofreading”. Instead, we show that replication hurdles, including secondary structures in template, slowed helicase, or uncoupled helicase–polymerase, increase DNA reannealing and polymerase backtracking, and generate frayed primer‐ends that are shuttled to the exonuclease site and excised efficiently. Our studies indicate that active‐site shuttling occurs at a high frequency, and we propose that it serves as a proofreading mechanism to protect primer‐ends from mutagenic extensions.

Interplay of primary sequence, position and secondary RNA structure determines alternative splicing of LMNA in a pre-mature aging syndrome

Aberrant splicing in exon 11 of the LMNA gene causes the premature aging disorder Hutchinson-Gilford Progeria Syndrome. A de novo C1824T mutation activates an internal alternative 5' splice site, resulting in formation of the disease-causing progerin protein. The underlying mechanism for this 5' splice site selection is unknown. Here, we have applied a combination of targeted mutational analysis in a cell-based system and structural mapping by SHAPE-MaP to comprehensively probe the contributions of primary sequence, secondary RNA structure and linear splice site position in determining in vivo mechanisms of splice site choice in LMNA. While splice site choice is in part defined by sequence complementarity to U1 snRNA, we identify RNA secondary structural elements near the alternative 5' splice sites and show that splice site choice is significantly influenced by the structural context of the available splice sites. Furthermore, relative positioning of the competing sites within the primary sequence of the pre-mRNA is a predictor of 5' splice site usage, with the distal position favored over the proximal, regardless of sequence composition. Together, these results demonstrate that 5' splice site selection in LMNA is determined by an intricate interplay among RNA sequence, secondary structure and splice site position.

MicroRNA and Nonsense Transcripts as Putative Viral Evasion Mechanisms

Viral proteins encode numerous antiviral activities to modify the host immunity. In this article, we hypothesize that viral genomes and gene transcripts interfere with host gene expression using passive mechanisms to deregulate host microRNA (miRNA) activity. We postulate that various RNA viruses mimic or block binding between a host miRNA and its target transcript, a phenomenon mediated by the miRNA seed site at the 5′ end of miRNA. Virus-encoded miRNA seed sponges (vSSs) can potentially bind to host miRNA seed sites and prevent interaction with their native targets thereby relieving native miRNA suppression. In contrast, virus-encoded miRNA seed mimics (vSMs) may mediate considerable downregulation of host miRNA activity. We analyzed genomes from diverse RNA viruses for vSS and vSM signatures and found an abundance of these motifs indicating that this may be a mechanism of deceiving host immunity. Employing respiratory syncytial virus and measles virus as models, we reveal that regions surrounding vSS or vSM motifs have features characteristics of pre-miRNA templates and show that RSV viral transcripts are processed into small RNAs that may behave as vSS or vSM effectors. These data suggest that complex molecular interactions likely occur at the host-virus interface. Identifying the mechanisms in the network of interactions between the host and viral transcripts can help uncover ways to improve vaccine efficacy, therapeutics, and potentially mitigate the adverse events that may be associated with some vaccines.

A cohabiting bacterium alters the spectrum of short RNAs secreted by Escherichia coli

Recently, it has been found that bacteria secrete short RNAs able to affect gene expression in eukaryotic cells, while certain mammalian microRNAs shape the gut microbiome altering bacterial transcriptome. The involvement of bacterial RNAs in communication with other bacteria is also expected, but has not been documented yet. Here, we compared the fractions of extremely short (12-22 nucleotides) RNAs secreted by Escherichia coli grown in a pure culture and jointly with bacteria of the Paenibacillus genus. Besides fragments of rRNAs and tRNAs, abundant in all samples, secreted oligonucleotides (exoRNAs) predominantly contained GC-rich fragments of messenger and antisense RNAs processed from regions with stable secondary structures. They differed in composition from oligonucleotides of intracellular fraction, where fragments of small regulatory RNAs were prevalent. Both fractions contained RNAs capable of forming complementary duplexes, while for exoRNA samples a higher percentage of 3΄-end modified RNAs and different endonuclease cleavage were detected. The presence of a cohabiting bacterium altered the spectrum of E. coli exoRNAs, indicating a population-dependent control over their composition. Possible mechanisms of this effect are discussed.

Two ribosome recruitment sites direct multiple translation events within HIV1 Gag open reading frame

In the late phase of the HIV virus cycle, the unspliced genomic RNA is exported to the cytoplasm for the necessary translation of the Gag and Gag-pol polyproteins. Three distinct translation initiation mechanisms ensuring Gag production have been described with little rationale for their multiplicity. The Gag-IRES has the singularity to be located within Gag ORF and to directly interact with ribosomal 40S. Aiming at elucidating the specificity and the relevance of this interaction, we probed HIV-1 Gag-IRES structure and developed an innovative integrative modelling strategy to take into account all the gathered information. We propose a novel Gag-IRES secondary structure strongly supported by all experimental data. We further demonstrate the presence of two regions within Gag-IRES that independently and directly interact with the ribosome. Importantly, these binding sites are functionally relevant to Gag translation both in vitro and ex vivo. This work provides insight into the Gag-IRES molecular mechanism and gives compelling evidence for its physiological importance. It allows us to propose original hypotheses about the IRES physiological role and conservation among primate lentiviruses.