Structural characterization of naturally occurring RNA single mismatches

RNA tertiary structure and conformational dynamics revealed by BASH MaP

The functional effects of an RNA can arise from complex three-dimensional folds known as tertiary structures. However, predicting the tertiary structure of an RNA and whether an RNA adopts distinct tertiary conformations remains challenging. To address this, we developed BASH MaP, a single-molecule dimethyl sulfate (DMS) footprinting method and DAGGER, a computational pipeline, to identify alternative tertiary structures adopted by different molecules of RNA. BASH MaP utilizes potassium borohydride to reveal the chemical accessibility of the N7 position of guanosine, a key mediator of tertiary structures. We used BASH MaP to identify diverse conformational states and dynamics of RNA G-quadruplexes, an important RNA tertiary motif, in vitro and in cells. BASH MaP and DAGGER analysis of the fluorogenic aptamer Spinach reveals that it adopts alternative tertiary conformations which determine its fluorescence states. BASH MaP thus provides an approach for structural analysis of RNA by revealing previously undetectable tertiary structures.

An RNA excited conformational state at atomic resolution

Sparse and short-lived excited RNA conformational states are essential players in cell physiology, disease, and therapeutic development, yet determining their 3D structures remains challenging. Combining mutagenesis, NMR spectroscopy, and computational modeling, we determined the 3D structural ensemble formed by a short-lived (lifetime ~2.1 ms) lowly-populated (~0.4%) conformational state in HIV-1 TAR RNA. Through a strand register shift, the excited conformational state completely remodels the 3D structure of the ground state (RMSD from the ground state = 7.2 ± 0.9 Å), forming a surprisingly more ordered conformational ensemble rich in non-canonical mismatches. The structure impedes the formation of the motifs recognized by Tat and the super elongation complex, explaining why this alternative TAR conformation cannot activate HIV-1 transcription. The ability to determine the 3D structures of fleeting RNA states using the presented methodology holds great promise for our understanding of RNA biology, disease mechanisms, and the development of RNA-targeting therapeutics.

Mass Spectrometry of Nucleic Acid Noncovalent Complexes

Nucleic acids have been among the first targets for antitumor drugs and antibiotics. With the unveiling of new biological roles in regulation of gene expression, specific DNA and RNA structures have become very attractive targets, especially when the corresponding proteins are undruggable. Biophysical assays to assess target structure as well as ligand binding stoichiometry, affinity, specificity, and binding modes are part of the drug development process. Mass spectrometry offers unique advantages as a biophysical method owing to its ability to distinguish each stoichiometry present in a mixture. In addition, advanced mass spectrometry approaches (reactive probing, fragmentation techniques, ion mobility spectrometry, ion spectroscopy) provide more detailed information on the complexes. Here, we review the fundamentals of mass spectrometry and all its particularities when studying noncovalent nucleic acid structures, and then review what has been learned thanks to mass spectrometry on nucleic acid structures, self-assemblies (e.g., duplexes or G-quadruplexes), and their complexes with ligands.

Electronic Circular Dichroism of the Cas9 Protein and gRNA:Cas9 Ribonucleoprotein Complex

The Streptococcus pyogenes Cas9 protein (SpCas9), a component of CRISPR-based immune system in microbes, has become commonly utilized for genome editing. This nuclease forms a ribonucleoprotein (RNP) complex with guide RNA (gRNA) which induces Cas9 structural changes and triggers its cleavage activity. Here, electronic circular dichroism (ECD) spectroscopy was used to confirm the RNP formation and to determine its individual components. The ECD spectra had characteristic features differentiating Cas9 and gRNA, the former showed a negative/positive profile with maxima located at 221, 209 and 196 nm, while the latter revealed positive/negative/positive/negative pattern with bands observed at 266, 242, 222 and 209 nm, respectively. For the first time, the experimental ECD spectrum of the gRNA:Cas9 RNP complex is presented. It exhibits a bisignate positive/negative ECD couplet with maxima at 273 and 235 nm, and it differs significantly from individual spectrum of each RNP components. Additionally, the Cas9 protein and RNP complex retained biological activity after ECD measurements and they were able to bind and cleave DNA in vitro. Hence, we conclude that ECD spectroscopy can be considered as a quick and non-destructive method of monitoring conformational changes of the Cas9 protein as a result of Cas9 and gRNA interaction, and identification of the gRNA:Cas9 RNP complex.

RNA CoSSMos 2.0: an improved searchable database of secondary structure motifs in RNA three-dimensional structures

The RNA Characterization of Secondary Structure Motifs, RNA CoSSMos, database is a freely accessible online database that allows users to identify secondary structure motifs among RNA 3D structures and explore their structural features. RNA CoSSMos 2.0 now requires two closing base pairs for all RNA loop motifs to create a less redundant database of secondary structures. Furthermore, RNA CoSSMos 2.0 represents an upgraded database with new features that summarize search findings and aid in the search for 3D structural patterns among RNA secondary structure motifs. Previously, users were limited to viewing search results individually, with no built-in tools to compare search results. RNA CoSSMos 2.0 provides two new features, allowing users to summarize, analyze and compare their search result findings. A function has been added to the website that calculates the average and representative structures of the search results. Additionally, users can now view a summary page of their search results that reports percentages of each structural feature found, including sugar pucker, glycosidic linkage, hydrogen bonding patterns and stacking interactions. Other upgrades include a newly embedded NGL structural viewer, the option to download the clipped structure coordinates in *.pdb format and improved NMR structure results. RNA CoSSMos 2.0 is no longer simply a search engine for a structure database; it now has the capability of analyzing, comparing and summarizing search results.

Database URL: http://rnacossmos.com

Identification of key sequence features required for microRNA biogenesis in plants

MicroRNAs (miRNAs) are endogenous small RNAs of ∼21 nt that regulate multiple biological pathways in multicellular organisms. They derive from longer transcripts that harbor an imperfect stem-loop structure. In plants, the ribonuclease type III DICER-LIKE1 assisted by accessory proteins cleaves the precursor to release the mature miRNA. Numerous studies highlight the role of the precursor secondary structure during plant miRNA biogenesis; however, little is known about the relevance of the precursor sequence. Here, we analyzed the sequence composition of plant miRNA primary transcripts and found specifically located sequence biases. We show that changes in the identity of specific nucleotides can increase or abolish miRNA biogenesis. Most conspicuously, our analysis revealed that the identity of the nucleotides at unpaired positions of the precursor plays a crucial role during miRNA biogenesis in Arabidopsis.

The secondary structure of miRNA precursor sequences is known to affect processing by DICER-like proteins. Here Rojas et al. show that additional sequence features also play a regulatory role in plants with nucleotide identity at unpaired positions substantially impacting processing efficiency.

Identification and Characterization of New RNA Tetraloop Sequence Families

There is an abundance of RNA sequence information available due to the efforts of sequencing projects. However, current techniques implemented to solve the tertiary structures of RNA, such as NMR and X-ray crystallography, are difficult and time-consuming. Therefore, biophysical techniques are not able to keep pace with the abundance of sequence information available. Because of this, there is a need to develop quick and efficient ways to predict RNA tertiary structure from sequence. One promising approach is to identify structural patterns within previously solved 3D structures and apply these patterns to new sequences. RNA tetraloops are one of the most common naturally occurring secondary structure motifs. Here, we use RNA Characterization of Secondary Structure Motifs (CoSSMos), Dissecting the Spatial Structure of RNA (DSSR), and a bioinformatic approach to search for and characterize tertiary structure patterns among tetraloops. Not surprising, we identified the well-known GNRA and UNCG tetraloops, as well as the previously identified RNYA tetraloop. However, some previously identified characteristics of these families were not observed in this data set, and some new characteristics were identified. In addition, we also identified and characterized three new tetraloop sequence families: YGAR, UGGU, and RMSA. This new structural information sheds light on the tertiary structure of tetraloops and contributes to the efforts of RNA tertiary structure prediction from sequence.

Evidence of Intra-individual SSU Polymorphisms in Dark-spored Myxomycetes (Amoebozoa)

The nuclear small subunit rRNA gene (SSU or 18S) is a marker frequently used in phylogenetic and barcoding studies in Amoebozoa, including Myxomycetes. Despite its common usage and the confirmed existence of divergent copies of ribosomal genes in other protists, the potential presence of intra-individual SSU variability in Myxomycetes has never been studied before. Here we investigated the pattern of nucleotide polymorphism in the 5' end fragment of SSU by cloning and sequencing a total of 238 variants from eight specimens, each representing a species of the dark-spored orders Stemonitidales and Physarales. After excluding singletons, a relatively low SSU intra-individual variability was found but our data indicate that this might be a widely distributed phenomenon in Myxomycetes as all samples analyzed possessed various ribotypes. To determine if the occurrence of multiple SSU variants within a single specimen has a negative effect on the circumscription of species boundaries, we conducted phylogenetic analyses that revealed that clone variation may be detrimental for inferring phylogenetic relationships among some of the specimens analyzed. Despite that intra-individual variability should be assessed in additional taxa, our results indicate that special care should be taken for species identification when working with closely related species.

Single-Molecule Mechanical Folding and Unfolding of RNA Hairpins: Effects of Single A-U to A·C Pair Substitutions and Single Proton Binding and Implications for mRNA Structure-Induced -1 Ribosomal Frameshifting

A wobble A·C pair can be protonated at near physiological pH to form a more stable wobble A⁺·C pair. Here, we constructed an RNA hairpin (rHP) and three mutants with one A-U base pair substituted with an A·C mismatch on the top (near the loop, U22C), middle (U25C), and bottom (U29C) positions of the stem, respectively. Our results on single-molecule mechanical (un)folding using optical tweezers reveal the destabilization effect of A-U to A·C pair substitution and protonation-dependent enhancement of mechanical stability facilitated through an increased folding rate, or decreased unfolding rate, or both. Our data show that protonation may occur rapidly upon the formation of an apparent mechanical folding transition state. Furthermore, we measured the bulk -1 ribosomal frameshifting efficiencies of the hairpins by a cell-free translation assay. For the mRNA hairpins studied, -1 frameshifting efficiency correlates with mechanical unfolding force at equilibrium and folding rate at around 15 pN. U29C has a frameshifting efficiency similar to that of rHP (∼2%). Accordingly, the bottom 2-4 base pairs of U29C may not form under a stretching force at pH 7.3, which is consistent with the fact that the bottom base pairs of the hairpins may be disrupted by ribosome at the slippery site. U22C and U25C have a similar frameshifting efficiency (∼1%), indicating that both unfolding and folding rates of an mRNA hairpin in a crowded environment may affect frameshifting. Our data indicate that mechanical (un)folding of RNA hairpins may mimic how mRNAs unfold and fold in the presence of translating ribosomes.

Modulation of the Thermodynamic Signatures of an RNA Thermometer by Osmolytes and Salts

Folding of ribonucleic acids (RNAs) is driven by several factors, such as base pairing and stacking, chain entropy, and ion-mediated electrostatics, which have been studied in great detail. However, the power of background molecules in the cellular milieu is often neglected. Herein, we study the effect of common osmolytes on the folding equilibrium of a hairpin-structured RNA and, using pressure perturbation, provide novel thermodynamic and volumetric insights into the modulation mechanism. The presence of TMAO causes an increased thermal stability and a more positive volume change for the helix-to-coil transition, whereas urea destabilizes the hairpin and leads to an increased expansibility of the unfolded state. Further, we find a strong interplay between water, salt, and osmolyte in driving the thermodynamics and defining the temperature and pressure stability limit of the RNA. Our results support a universal working mechanism of TMAO and urea to (de)stabilize proteins and the RNA.

Enhancing potency of siRNA targeting fusion genes by optimization outside of target sequence

Canonical siRNA design algorithms have become remarkably effective at predicting favorable binding regions within a target mRNA, but in some cases (e.g., a fusion junction site) region choice is restricted. In these instances, alternative approaches are necessary to obtain a highly potent silencing molecule. Here we focus on strategies for rational optimization of two siRNAs that target the junction sites of fusion oncogenes BCR-ABL and TMPRSS2-ERG. We demonstrate that modifying the termini of these siRNAs with a terminal G-U wobble pair or a carefully selected pair of terminal asymmetry-enhancing mismatches can result in an increase in potency at low doses. Importantly, we observed that improvements in silencing at the mRNA level do not necessarily translate to reductions in protein level and/or cell death. Decline in protein level is also heavily influenced by targeted protein half-life, and delivery vehicle toxicity can confound measures of cell death due to silencing. Therefore, for BCR-ABL, which has a long protein half-life that is difficult to overcome using siRNA, we also developed a nontoxic transfection vector: poly(lactic-coglycolic acid) nanoparticles that release siRNA over many days. We show that this system can achieve effective killing of leukemic cells. These findings provide insights into the implications of siRNA sequence for potency and suggest strategies for the design of more effective therapeutic siRNA molecules. Furthermore, this work points to the importance of integrating studies of siRNA design and delivery, while heeding and addressing potential limitations such as restricted targetable mRNA regions, long protein half-lives, and nonspecific toxicities.

A Tandem Oligonucleotide Approach for SNP-Selective RNA Degradation Using Modified Antisense Oligonucleotides

Antisense oligonucleotides have been studied for many years as a tool for gene silencing. One of the most difficult cases of selective RNA silencing involves the alleles of single nucleotide polymorphisms, in which the allele sequence is differentiated by a single nucleotide. A new approach to improve the performance of allele selectivity for antisense oligonucleotides is proposed. It is based on the simultaneous application of two oligonucleotides. One is complementary to the mutated form of the targeted RNA and is able to activate RNase H to cleave the RNA. The other oligonucleotide, which is complementary to the wild type allele of the targeted RNA, is able to inhibit RNase H cleavage. Five types of SNPs, C/G, G/C, G/A, A/G, and C/U, were analyzed within the sequence context of genes associated with neurodegenerative disorders such as Alzheimer’s disease, Parkinson’s disease, ALS (Amyotrophic Lateral Sclerosis), and Machado-Joseph disease. For most analyzed cases, the application of the tandem approach increased allele-selective RNA degradation 1.5–15 fold relative to the use of a single antisense oligonucleotide. The presented study proves that differentiation between single substitution is highly dependent on the nature of the SNP and surrounding nucleotides. These variables are crucial for determining the proper length of the inhibitor antisense oligonucleotide. In the tandem approach, the comparison of thermodynamic stability of the favorable duplexes WT RNA-inhibitor and Mut RNA-gapmer with the other possible duplexes allows for the evaluation of chances for the allele-selective degradation of RNA. A larger difference in thermodynamic stability between favorable duplexes and those that could possibly form, usually results in the better allele selectivity of RNA degradation.

Hybridization Properties of RNA Containing 8-Methoxyguanosine and 8-Benzyloxyguanosine

Modified nucleobase analogues can serve as powerful tools for changing physicochemical and biological properties of DNA or RNA. Guanosine derivatives containing bulky substituents at 8 position are known to adopt syn conformation of N-glycoside bond. On the contrary, in RNA the anti conformation is predominant in Watson-Crick base pairing. In this paper two 8-substituted guanosine derivatives, 8-methoxyguanosine and 8-benzyloxyguanosine, were synthesized and incorporated into oligoribonucleotides to investigate their influence on the thermodynamic stability of RNA duplexes. The methoxy and benzyloxy substituents are electron-donating groups, decreasing the rate of depurination in the monomers, as confirmed by N-glycoside bond stability assessments. Thermodynamic stability studies indicated that substitution of guanosine by 8-methoxy- or 8-benzyloxyguanosine significantly decreased the thermodynamic stability of RNA duplexes. Moreover, the presence of 8-substituted guanosine derivatives decreased mismatch discrimination. Circular dichroism spectra of modified RNA duplexes exhibited patterns typical for A-RNA geometry.

Phosphoribosylphosphate and phosphoribosylnicotinate pairing with phosphoribosylamine at the origin of the RNA world

The gap between prebiotic chemistry and the RNA origin of life lies in the nature of the pre-ribonucleotides that initiated replication-like activity. The present paper illustrates how the constraints for self-catalytic synthesis of the building blocks point to ionic pairing of the first pre-ribonucleotides and point to a central role for nicotinate.

Hierarchy of RNA functional dynamics

RNA dynamics play a fundamental role in many cellular functions. However, there is no general framework to describe these complex processes, which typically consist of many structural maneuvers that occur over timescales ranging from picoseconds to seconds. Here, we classify RNA dynamics into distinct modes representing transitions between basins on a hierarchical free-energy landscape. These transitions include large-scale secondary-structural transitions at >0.1-s timescales, base-pair/tertiary dynamics at microsecond-to-millisecond timescales, stacking dynamics at timescales ranging from nanoseconds to microseconds, and other "jittering" motions at timescales ranging from picoseconds to nanoseconds. We review various modes within these three different tiers, the different mechanisms by which they are used to regulate function, and how they can be coupled together to achieve greater functional complexity.

Effect of sodium ions on RNA duplex stability

The standard sodium concentration for RNA optical melting experiments is 1.021 M. Algorithms that predict Tm, ΔG°37, and secondary structure from sequence generally rely on parameters derived from optical melting experiments performed in 1.021 M sodium. Physiological monovalent cation concentrations are much lower than 1.021 M. In fact, many molecular biology techniques require buffers containing monovalent cation concentrations other than 1.021 M. Predictions based on the 1.021 M Na(+) parameters may not be accurate when the monovalent cation concentration is not 1.021 M. Here, we report thermodynamic data from optical melting experiments for a set of 18 RNA duplexes, each melted over a wide range of sodium ion concentrations (71, 121, 221, and 621 mM). Using these data and previously published data for the same sequences melted in 1.021 M Na(+), we report Tm and ΔG°37 correction factors to scale the standard 1.021 M Na(+) RNA parameters to other sodium ion concentrations. The recommended Tm correction factor predicts the melting temperature within 0.7 °C, and the recommended ΔG°37 correction factor predicts the free energy within 0.14 kcal/mol. These correction factors can be incorporated into prediction algorithms that predict RNA secondary structure from sequence and provide Tm and ΔG°37 values for RNA duplexes.

A dynamic structural model of expanded RNA CAG repeats: a refined X-ray structure and computational investigations using molecular dynamics and umbrella sampling simulations

One class of functionally important RNA is repeating transcripts that cause disease through various mechanisms. For example, expanded CAG repeats can cause Huntington's and other disease through translation of toxic proteins. Herein, a crystal structure of r[5'UUGGGC(CAG)3GUCC]2, a model of CAG expanded transcripts, refined to 1.65 Å resolution is disclosed that shows both anti-anti and syn-anti orientations for 1 × 1 nucleotide AA internal loops. Molecular dynamics (MD) simulations using AMBER force field in explicit solvent were run for over 500 ns on the model systems r(5'GCGCAGCGC)2 (MS1) and r(5'CCGCAGCGG)2 (MS2). In these MD simulations, both anti-anti and syn-anti AA base pairs appear to be stable. While anti-anti AA base pairs were dynamic and sampled multiple anti-anti conformations, no syn-anti ↔ anti-anti transformations were observed. Umbrella sampling simulations were run on MS2, and a 2D free energy surface was created to extract transformation pathways. In addition, an explicit solvent MD simulation over 800 ns was run on r[5'GGGC(CAG)3GUCC]2, which closely represents the refined crystal structure. One of the terminal AA base pairs (syn-anti conformation), transformed to anti-anti conformation. The pathway followed in this transformation was the one predicted by umbrella sampling simulations. Further analysis showed a binding pocket near AA base pairs in syn-anti conformations. Computational results combined with the refined crystal structure show that global minimum conformation of 1 × 1 nucleotide AA internal loops in r(CAG) repeats is anti-anti but can adopt syn-anti depending on the environment. These results are important to understand RNA dynamic-function relationships and to develop small molecules that target RNA dynamic ensembles.

Ensemble analysis of primary microRNA structure reveals an extensive capacity to deform near the Drosha cleavage site

Most noncoding RNAs function properly only when folded into complex three-dimensional (3D) structures, but the experimental determination of these structures remains challenging. Understanding of primary microRNA (miRNA) maturation is currently limited by a lack of determined structures for nonprocessed forms of the RNA. SHAPE chemistry efficiently determines RNA secondary structural information with single-nucleotide resolution, providing constraints suitable for input into MC-Pipeline for refinement of 3D structure models. Here we combine these approaches to analyze three structurally diverse primary microRNAs, revealing deviations from canonical double-stranded RNA structure in the stem adjacent to the Drosha cut site for all three. The necessity of these deformable sites for efficient processing is demonstrated through Drosha processing assays. The structure models generated herein support the hypothesis that deformable sequences spaced roughly once per turn of A-form helix, created by noncanonical structure elements, combine with the necessary single-stranded RNA-double-stranded RNA junction to define the correct Drosha cleavage site.

RNA CoSSMos: Characterization of Secondary Structure Motifs--a searchable database of secondary structure motifs in RNA three-dimensional structures

RNA secondary structure is important for designing therapeutics, understanding protein–RNA binding and predicting tertiary structure of RNA. Several databases and downloadable programs exist that specialize in the three-dimensional (3D) structure of RNA, but none focus specifically on secondary structural motifs such as internal, bulge and hairpin loops. The RNA Characterization of Secondary Structure Motifs (RNA CoSSMos) database is a freely accessible and searchable online database and website of 3D characteristics of secondary structure motifs. To create the RNA CoSSMos database, 2156 Protein Data Bank (PDB) files were searched for internal, bulge and hairpin loops, and each loop's structural information, including sugar pucker, glycosidic linkage, hydrogen bonding patterns and stacking interactions, was included in the database. False positives were defined, identified and reclassified or omitted from the database to ensure the most accurate results possible. Users can search via general PDB information, experimental parameters, sequence and specific motif and by specific structural parameters in the subquery page after the initial search. Returned results for each search can be viewed individually or a complete set can be downloaded into a spreadsheet to allow for easy comparison. The RNA CoSSMos database is automatically updated weekly and is available at http://cossmos.slu.edu.

A mutate-and-map strategy accurately infers the base pairs of a 35-nucleotide model RNA

We present a rapid experimental strategy for inferring base pairs in structured RNAs via an information-rich extension of classic chemical mapping approaches. The mutate-and-map method, previously applied to a DNA/RNA helix, systematically searches for single mutations that enhance the chemical accessibility of base-pairing partners distant in sequence. To test this strategy for structured RNAs, we have carried out mutate-and-map measurements for a 35-nt hairpin, called the MedLoop RNA, embedded within an 80-nt sequence. We demonstrate the synthesis of all 105 single mutants of the MedLoop RNA sequence and present high-throughput DMS, CMCT, and SHAPE modification measurements for this library at single-nucleotide resolution. The resulting two-dimensional data reveal visually clear, punctate features corresponding to RNA base pair interactions as well as more complex features; these signals can be qualitatively rationalized by comparison to secondary structure predictions. Finally, we present an automated, sequence-blind analysis that permits the confident identification of nine of the 10 MedLoop RNA base pairs at single-nucleotide resolution, while discriminating against all 1460 false-positive base pairs. These results establish the accuracy and information content of the mutate-and-map strategy and support its feasibility for rapidly characterizing the base-pairing patterns of larger and more complex RNA systems.

Molecular recognition of 6'-N-5-hexynoate kanamycin A and RNA 1x1 internal loops containing CA mismatches

In our previous study to identify the RNA internal loops that bind an aminoglycoside derivative, we determined that 6'-N-5-hexynoate kanamycin A prefers to bind 1x1 nucleotide internal loops containing C·A mismatches. In this present study, the molecular recognition between a variety of RNAs that are mutated around the C·A loop and the ligand was investigated. Studies show that both loop nucleotides and loop closing pairs affect binding affinity. Most interestingly, it was shown that there is a correlation between the thermodynamic stability of the C·A internal loops and ligand affinity. Specifically, C·A loops that had relatively high or low stability bound the ligand most weakly whereas loops with intermediate stability bound the ligand most tightly. In contrast, there is no correlation between the likelihood that a loop forms a C-A(+) pair at lower pH and ligand affinity. It was also found that a 1x1 nucleotide C·A loop that bound to the ligand with the highest affinity is identical to the consensus site in RNAs that are edited by adenosine deaminases acting on RNA type 2 (ADAR2). These studies provide a detailed investigation of factors affecting small molecule recognition of internal loops containing C·A mismatches, which are present in a variety of RNAs that cause disease.