Unpaired Nucleotides Research Articles

Chikungunya Virus RNA Secondary Structures Impact Defective Viral Genome Production.

Chikungunya virus (CHIKV) is a mosquito-borne RNA virus that poses an emerging threat to humans. In a manner similar to other RNA viruses, CHIKV encodes an error-prone RNA polymerase which, in addition to producing full-length genomes, gives rise to truncated, non-functional genomes, which have been coined defective viral genomes (DVGs). DVGs have been intensively studied in the context of therapy, as they can inhibit viral replication and dissemination in their hosts. In this work, we interrogate the influence of viral RNA secondary structures on the production of CHIKV DVGs. We experimentally map RNA secondary structures of the CHIKV genome using selective 2'-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP), which couples chemical labelling with next-generation sequencing. We correlate the inferred secondary structure with preferred deletion sites of CHIKV DVGs. We document an increased probability of DVG generation with truncations at unpaired nucleotides within the secondary structure. We then generated a CHIKV mutant bearing synonymous changes at the nucleotide level to disrupt the existing RNA secondary structure (CHIKV-D2S). We show that CHIKV-D2S presents altered DVG generation compared to wild-type virus, correlating with the change in RNA secondary structure obtained by SHAPE-MaP. Our work thus demonstrates that RNA secondary structure impacts CHIKV DVG production during replication.

A new reagent for in vivo structure probing of RNA G and U residues that improves RNA structure prediction alone and combined with DMS.

A key to understanding the roles of RNA in regulating gene expression is knowing their structures in vivo. One way to obtain this information is through probing the structures of RNA with chemicals. To probe RNA structure directly in cells, membrane-permeable reagents that modify the Watson-Crick (WC) face of unpaired nucleotides can be used. Although dimethyl sulfate (DMS) has led to substantial insight into RNA structure, it has limited nucleotide specificity in vivo, with WC face reactivity only at adenine (A) and cytosine (C) at neutral pH. The reagent 1-ethyl-3-(3-dimethylaminopropyl) carbodiimide (EDC) was recently shown to modify the WC face of guanine (G) and uracil (U). Although useful at lower concentrations in experiments that measure chemical modifications by reverse transcription (RT) stops, at higher concentrations necessary for detection by mutational profiling (MaP), EDC treatment leads to degradation of RNA. Here, we demonstrate EDC-stimulated degradation of RNA in Gram-negative and Gram-positive bacteria. In an attempt to overcome these limitations, we developed a new carbodiimide reagent, 1-ethyl-3-(3-dimethylaminopropyl) carbodiimide methiodide (ETC), which we show specifically modifies unpaired Gs and Us in vivo without substantial degradation of RNA. We establish ETC as a probe for MaP and optimize the RT conditions and computational analysis in Escherichia coli Importantly, we demonstrate the utility of ETC as a probe for improving RNA structure prediction both alone and with DMS.

Plant-specific environmental and developmental signals regulate the mismatch repair protein MSH6 in Arabidopsis thaliana

The DNA mismatch repair (MMR) is a postreplicative system that guarantees genomic stability by correcting mispaired and unpaired nucleotides. In eukaryotic nuclei, MMR is initiated by the binding of heterodimeric MutS homologue (MSH) complexes to the DNA error or lesion. Among these proteins, MSH2-MSH6 is the most abundant heterodimer. Even though the MMR mechanism and proteins are highly conserved throughout evolution, physiological differences between species can lead to different regulatory features. Here, we investigated how light, sugar, and/or hormones modulate Arabidopsis thaliana MSH6 expression pattern. We first characterized the promoter region of MSH6. Phylogenetic shadowing revealed three highly conserved regions. These regions were analyzed by the generation of deletion constructs of the MSH6 full-length promoter fused to the β-glucuronidase (GUS) gene. Combined, our in silico and genetic analyses revealed that a 121-bp promoter fragment was necessary for MSH6 expression and contained potential cis-acting elements involved in light- and hormone-responsive gene expression. Accordingly, light exposure or sugar treatment of four-day old A. thaliana seedlings triggered an upregulation of MSH6 in shoot and root apical meristems. Appropriately, MSH6 was also induced by the stem cell inducer WUSCHEL. Further, the stimulatory effect of light was dependent on the presence of phyA. In addition, treatment of seedlings with auxin or cytokinin also caused an upregulation of MSH6 under darkness. Consistent with auxin signals, MSH6 expression was suppressed in the GATA23 RNAi line compared with the wild type. Our results provide evidence that endogenous factors and environmental signals controlling plant growth and development regulate the MSH6 protein in A. thaliana.

High-resolution RNA tertiary structures in Zika virus stem-loop A for the development of inhibitory small molecules.

Flaviviruses such as Zika (ZIKV) and dengue virus (DENV) are positive-sense RNA viruses belonging to Flaviviridae The flavivirus genome contains a 5' end stem-loop promoter sequence known as stem-loop A (SLA) that is recognized by the flavivirus polymerase NS5 during viral RNA synthesis and 5' guanosine cap methylation. The crystal structures of ZIKV and DENV SLAs show a well-defined fold, consisting of a bottom stem, side loop, and top stem-loop, providing unique interaction sites for small molecule inhibitors to disrupt the promoter function. To facilitate the identification of small molecule binding sites in flavivirus SLA, we determined high-resolution structures of the bottom and top stems of ZIKV SLA, which contain a single U- or G-bulge, respectively. Both bulge nucleotides exhibit multiple orientations, from folded back on the adjacent nucleotide to flipped out of the helix, and are stabilized by stacking or base triple interactions. These structures suggest that even a single unpaired nucleotide can provide flexibility to RNA structures, and its conformation is mainly determined by the stabilizing chemical environment. To facilitate discovery of small molecule inhibitors that interfere with the functions of ZIKV SLA, we screened and identified compounds that bind to the bottom and top stems of ZIKV SLA.

Expression of immune checkpoints PD-L1, CTLA4, LAG3 in the microenvironment of colon adenocarcinoma depending on MMR status

Study of the features of expression of immune checkpoint proteins PD-L1, CTLA4 and LAG3 in the microenvironment of colon adenocarcinoma depending on MMR status. The study group consisted of 32 patients with a morphologically confirmed diagnosis of colon cancer; all of them underwent surgical treatment in the form of hemicolonectomy or resection. The work assessed samples of tumor tissue obtained as a result of surgery, the study was carried out in 3 stages: morphological examination of histological slides of colon tumors at the light-optical level, immunohistochemistry examination of tumor samples to determine the dMMR/pMMR status of carcinoma using a panel of antibodies to proteins of the unpaired nucleotide repair system MLH1, MSH2, MSH6 and PMS2, multiplex analysis of PD-L1, CTLA4, LAG3, CD3+, CD8+, CD163+ markers using the Vectra 3.0.3 tissue scanning system (Perkin Elmer, USA). Significant differences in the expression of PD-L1, CTLA4, LAG3 in the area of the invasive tumor margin were revealed between the dMMR and pMMR groups of colon adenocarcinomas in patients comparable in clinical and morphological characteristics and treatment. In the group of tumors with dMMR status, an increase in the expression of all studied markers was noted. The number of CD3+ TILs was also significantly higher in the invasive margin of tumors with dMMR status. Similarly, in this group of colon carcinomas, a large number of CD163+ macrophages were noted both in the center and in the invasive margin zone. No statistically significant differences were found in the expression of immune checkpoints and the composition of TILs in the central zone of tumors with different MMR status. A study using multiplex immunohistochemical analysis showed that MMR-deficient colon adenocarcinomas are characterized by more pronounced immune infiltration and increased expression of immune checkpoints in microenvironmental cells, mainly in the area of invasive tumor growth. The data obtained may be important for understanding the mechanisms of immune-mediated control of tumor growth and the choice of immunotherapy tactics depending on MMR status.

Post-transcriptional capping generates coenzyme A-linked RNA

ABSTRACT NAD can be inserted co-transcriptionally via non-canonical initiation to form NAD-RNA. However, that mechanism is unlikely for CoA-linked RNAs due to low intracellular concentration of the required initiator nucleotide, 3’-dephospho-CoA (dpCoA). We report here that phosphopantetheine adenylyltransferase (PPAT), an enzyme of CoA biosynthetic pathway, accepts RNA transcripts as its acceptor substrate and transfers 4’-phosphopantetheine to yield CoA-RNA post-transcriptionally. Synthetic natural (RNAI) and small artificial RNAs were used to identify the features of RNA that are needed for it to serve as PPAT substrate. RNAs with 4–10 unpaired nucleotides at the 5’ terminus served as PPAT substrates, but RNAs having <4 unpaired nucleotides did not undergo capping. No capping was observed when the +1A was changed to G or when 5’ triphosphate was removed by RNA pyrophosphohydrolase (RppH), suggesting the enzyme recognizes pppA-RNA as an ATP analog. PPAT binding affinities were equivalent for transcripts with +1A, +1 G, or 5’OH (+1A), indicating that productive enzymatic recognition is driven more by local positioning effects than by overall binding affinity. Capping rates were independent of the number of unpaired nucleotides in the range of 4–10 nucleotides. Capping was strongly inhibited by ATP, reducing CoA-RNA production ~70% when equimolar ATP and substrate RNA were present. Dual bacterial expression of candidate RNAs with different 5’ structures followed by CoA-RNA CaptureSeq revealed 12-fold enrichment of the better PPAT substrate, consistent with in vivo CoA-capping of RNA transcripts by PPAT. These results suggest post-transcriptional RNA capping as a possible mechanism for the biogenesis of CoA-RNAs in bacteria.

Altering the Cleaving Effector in Chimeric Molecules that Target RNA Enhances Cellular Selectivity.

Small molecules that target RNA and effect their cleavage are useful chemical probes and potential lead medicines. In this study, we investigate factors affecting degradation of two cancer-associated RNA targets, the mRNA that encodes the transcription factor JUN (c-Jun) and the hairpin precursor to microRNA-372 (pre-miR-372). The two RNA targets harbor the same small-molecule binding site juxtaposed with different neighboring structures. Specifically, pre-miR-372 has AU pairs and contiguous purines on one strand near the small-molecule binding site, making it an ideal substrate for oxidative cleavage via the direct degrader bleomycin A5. In contrast, while JUN mRNA has a similar number of AU pairs near the small-molecule binding site, it lacks contiguous purine nucleotides. Instead, it contains unpaired pyrimidine nucleotides, which are preferred substrates for RNase L, a ribonuclease that can be recruited to RNA with heterobifunctional ribonuclease targeting chimeras (RiboTACs). We hypothesized that structural features surrounding the binding site could be leveraged to program selective small-molecule degradation by alteration of the cleaving module. Indeed, the bleomycin degrader cleaves pre-miR-372 in gastric cancer cells but not JUN mRNA. Conversely, the RiboTAC cleaves JUN mRNA but not pre-miR-372. Thus, the selection of the appropriate cleaving effector moiety for an RNA-binding small molecule can be leveraged to cleave an RNA selectively in a predictable manner, which could have broad implications.

Quantifying RNA structures and interactions with a unified reduced chain representation model

RNA is a pivotal molecule that plays critical roles in various cellular processes. Quantifying RNA structures and interactions is essential to understanding RNA function and developing RNA-based therapeutics. Using a unified five-bead model and a non-redundant database, this paper investigates the structural features and interactions of five commonly occurring RNA motifs, i.e., double-stranded helices, hairpin loops, internal/bulge loops, multi-branched junctions, and single-stranded terminal tails. Analyzing detailed distributions of RNA local structural features and base-base interactions reveals a preference for helical structures in both local backbone structures and base orientations. The interactions between adjacent bases exhibit motif-specific and sequence-dependent characteristics, reflecting the distinct topological constraints imposed by different loop-helix connection modes and the varying pairing and stacking interactions among different sequences. These findings shed light on the stability of RNA helices, emphasizing their significance in providing dominant base pairing and stacking interactions for RNA structures and stability. The four non-helix motifs encompass unpaired nucleotide loops and exhibit diverse base-base interactions, contributing to the structural diversity observed in RNA. Overall, the complexity of RNA structure arises from the intricate interplay of base-base interactions.

DNA Volume, Topology, and Flexibility Dictate Nanopore Current Signals.

Nanopores have developed into powerful single-molecule sensors capable of identifying and characterizing small polymers, such as DNA, by electrophoretically driving them through a nanoscale pore and monitoring temporary blockades in the ionic pore current. However, the relationship between nanopore signals and the physical properties of DNA remains only partly understood. Herein, we introduce a programmable DNA carrier platform to capture carefully designed DNA nanostructures. Controlled translocation experiments through our glass nanopores allowed us to disentangle this relationship. We vary DNA topology by changing the length, strand duplications, sequence, unpaired nucleotides, and rigidity of the analyte DNA and find that the ionic current drop is mainly determined by the volume and flexibility of the DNA nanostructure in the nanopore. Finally, we use our understanding of the role of DNA topology to discriminate circular single-stranded DNA molecules from linear ones with the same number of nucleotides using the nanopore signal.

Triplex-Forming Peptide Nucleic Acid Controls Dynamic Conformations of RNA Bulges.

RNA folding is driven by the formation of double-helical segments interspaced by loops of unpaired nucleotides. Among the latter, bulges formed by one or several unpaired nucleotides are one of the most common structural motifs that play an important role in stabilizing RNA-RNA, RNA-protein, and RNA-small molecule interactions. Single-nucleotide bulges can fold in alternative structures where the unpaired nucleobase is either looped-out (flexible) in a solvent or stacked-in (intercalated) between the base pairs. In the present study, we discovered that triplex-forming peptide nucleic acids (PNAs) had unusually high affinity for single-purine-nucleotide bulges in double-helical RNA. Depending on the PNA's sequence, the triplex formation shifted the equilibrium between looped-out and stacked-in conformations. The ability to control the dynamic equilibria of RNA's structure will be an important tool for studying structure-function relationships in RNA biology and may have potential in novel therapeutic approaches targeting disease-related RNAs.

Deficiency of the Arabidopsis mismatch repair MSH6 attenuates Pseudomonas syringae invasion

The MutS homolog 6 (MSH6) is a nuclear DNA mismatch repair (MMR) gene that encodes the MSH6 protein. MSH6 interacts with MSH2 to form the MutSα heterodimer. MutSα corrects DNA mismatches and unpaired nucleotides arising during DNA replication, deamination of 5-methylcytosine, and recombination between non-identical DNA sequences. In addition to correcting DNA biosynthetic errors, MutSα also recognizes chemically damaged DNA bases. Here, we show that inactivation of MSH6 affects the basal susceptibility of Arabidopsis thaliana to Pseudomonas syringae pv tomato DC3000. The msh6 T-DNA insertional mutant exhibited a reduced susceptibility to the bacterial invasion. This heightened basal resistance of msh6 mutants appears to be dependent on an increased stomatal closure, an accumulation of H2O2 and double-strand breaks (DSBs) and a constitutive expression of pathogenesis-related (NPR1 and PR1) and DNA damage response (RAD51D and SOG1) genes. Complementation of this mutant with the MSH6 wild type allele under the control of its own promoter resulted in reversal of the basal bacterial resistance phenotype and the stomatal closure back to wild type levels. Taken together, these results demonstrate that inactivation of MSH6 increases Arabidopsis basal susceptibility to the bacterial pathogen and suggests a link between DNA repair and stress signaling in plants.

Sulfonylation of RNA 2'-OH groups.

The nucleophilic reactivity of RNA 2'-OH groups in water has proven broadly useful in probing, labeling, and conjugating RNA. To date, reactions selective to ribose 2'-OH have been limited to bond formation with short-lived carbonyl electrophiles. Here we report that many activated small-molecule sulfonyl species can exhibit extended lifetimes in water and retain 2'-OH reactivity. The data establish favorable aqueous solubility for selected reagents and successful RNA-selective reactions at stoichiometric and superstoichiometric yields, particularly for aryl sulfonyltriazole species. We report that the latter are considerably more stable than most prior carbon electrophiles in aqueous environments and tolerate silica chromatography. Furthermore, an azide-substituted sulfonyltriazole reagent is developed to introduce labels into RNA via click chemistry. In addition to high-yield reactions, we find that RNA sulfonylation can also be performed under conditions that give trace yields necessary for structure mapping. Like acylation, the reaction occurs with selectivity for unpaired nucleotides over those in the duplex structure, and a sulfonate adduct causes reverse transcriptase stops, suggesting potential use in RNA structure analysis. Probing of rRNA is demonstrated in human cells, indicating possible cell permeability. The sulfonyl reagent class enables a new level of control, selectivity, versatility, and ease of preparation for RNA applications.

Microsatellite instability and DNA mismatch repair deficiency detection in tumors of various sites

Microsatellite instability, which is caused by a deficiency in the DNA unpaired nucleotide repair system, is an important pathogenetic event for some tumors. In addition, the detection of this molecular feature becomes an independent prognostic factor in the course of the disease and a predictor for the appointment of therapy with immune checkpoint inhibitors. Immunohistochemistry is a reliable and available method for detecting a deficiency in the DNA mismatch repair system, and it has recommended as a screening for hereditary syndromes associated with microsatellite instability. This article discusses the advantages and disadvantages of this research method from the point of view of the practitioner.

Diverse Reagent Scaffolds Provide Differential Selectivity of 2'-OH Acylation in RNA.

RNA 2'-OH acylation is widely used both for mapping structure and for conjugating RNA, generally relying on selective reactions with unpaired nucleotides over paired ones. Common reagents for this acylation have been chiefly restricted to two similar aryl scaffolds, leaving open the question of how more broadly varied reagent structure might affect selectivity. Here, we prepared a set of 10 structurally diverse acylimidazole reagents and employed deep sequencing to profile their reactivity and selectivity in an RNA library of systematically varied structure. We show that structure-directed reactivity profiles vary significantly with the reagent scaffold, and we document new acylating agents that have altered selectivity profiles, including reagents that show elevated selectivity within loops, as well as compounds with reduced off-target reactivity in loop closing base pairs. Interestingly, we also show that the simplest reagent (acetylimidazole) is cell permeable and is small enough to map RNA structure in the presence of protein contacts that block other reagents. Finally, we describe reagents that show elevated selectivity within small loops, with applications in site-selective labeling. The results provide new tools for improved conjugation and mapping of RNA.

A novel approach for a joint analysis of isomiR and mRNA expression data reveals features of isomiR targeting in breast cancer.

A widely used procedure for selecting significant miRNA-mRNA or isomiR-mRNA pairs out of predicted interactions involves calculating the correlation between expression levels of miRNAs/isomiRs and mRNAs in a series of samples. In this manuscript, we aimed to assess the validity of this procedure by comparing isomiR-mRNA correlation profiles in sets of sequence-based predicted target mRNAs and non-target mRNAs (negative controls). Target prediction was carried out using RNA22 and TargetScan algorithms. Spearman's correlation analysis was conducted using miRNA and mRNA sequencing data of The Cancer Genome Atlas Breast Invasive Carcinoma (TCGA-BRCA) project. Luminal A, luminal B, basal-like breast cancer subtypes, and adjacent normal tissue samples were analyzed separately. Using the sets of putative targets and non-targets, we introduced adjusted isomiR targeting activity (ITA)-the number of negatively correlated potential isomiR targets adjusted by the background (estimated using non-target mRNAs). We found that for most isomiRs a significant negative correlation between isomiR-mRNA expression levels appeared more often in a set of predicted targets compared to the non-targets. This trend was detected for both classical seed region binding types (8mer, 7mer-m8, 7mer-A1, 6mer) predicted by TargetScan and the non-classical ones (G:U wobbles and up to one mismatch or unpaired nucleotide within seed sequence) predicted by RNA22. Adjusted ITA distributions were similar for target sites located in 3'-UTRs and coding mRNA sequences, while 5'-UTRs had much lower scores. Finally, we observed strong cancer subtype-specific patterns of isomiR activity, highlighting the differences between breast cancer molecular subtypes and normal tissues. Surprisingly, our target prediction- and correlation-based estimates of isomiR activities were practically non-correlated with the average isomiR expression levels neither in cancerous nor in normal samples.

DNA Penetration into a Lysozyme Layer at the Surface of Aqueous Solutions.

The interactions of DNA with lysozyme in the surface layer were studied by performing infrared reflection–absorption spectroscopy (IRRAS), ellipsometry, surface tensiometry, surface dilational rheology, and atomic force microscopy (AFM). A concentrated DNA solution was injected into an aqueous subphase underneath a spread lysozyme layer. While the optical properties of the surface layer changed fast after DNA injection, the dynamic dilational surface elasticity almost did not change, thereby indicating no continuous network formation of DNA/lysozyme complexes, unlike the case of DNA interactions with a monolayer of a cationic synthetic polyelectrolyte. A relatively fast increase in optical signals after a DNA injection under a lysozyme layer indicates that DNA penetration is controlled by diffusion. At low surface pressures, the AFM images show the formation of long strands in the surface layer. Increased surface compression does not lead to the formation of a network of DNA/lysozyme aggregates as in the case of a mixed layer of DNA and synthetic polyelectrolytes, but to the appearance of some folds and ridges in the layer. The formation of more disordered aggregates is presumably a consequence of weaker interactions of lysozyme with duplex DNA and the stabilization, at the same time, of loops of unpaired nucleotides at high local lysozyme concentrations in the surface layer.

Evaluation of weak interactions of proteins and organic cations with DNA duplex structures

Molecular interactions and reactions in living cells occur with high background concentrations of organic compounds including proteins. Uncharged water-soluble polymers are commonly used cosolutes in studies on molecular crowding, and most studies argue about the effects of intracellular crowding based on results obtained using polymer cosolutes. Further investigations using protein crowders and organic cations are important in understanding the effects of cellular environments on nucleic acids with negatively charged surfaces. We assessed the effects of using model globular proteins, serum proteins, histone proteins, structurally flexible polypeptides, di- and polyamines, and uncharged polymers. Thermal stability analysis of DNA oligonucleotide structures revealed that unlike conventional polymer cosolutes, basic globular proteins (lysozyme and cytochrome c) at high concentrations stabilized long internal and bulge loop structures but not fully matched duplexes. The selective stabilization of long loop structures suggests preferential binding to unpaired nucleotides in loops through weak electrostatic interactions. Furthermore, the ability of the proteins to stabilize the loop structures was enhanced under macromolecular crowding conditions. Remarkably, the effects of basic proteins on the stability of fully matched duplexes were dissimilar to those of basic amino-acid-rich polypeptides and polyamines. This study provides new insights into the interaction of nucleic acid structures with organic cations.

Unpaired nucleotides on the stem of microRNA precursor are important for precise cleavage by Dicer-like 1 in Arabidopsis.

Dicer-like 1 (DCL1) is a core component of the plant microRNA (miRNA) biogenesis machinery. MiRNA is transcribed as a precursor RNA, termed primary miRNA (pri-miRNA), which is cleaved by DCL1 in two steps to generate miRNA/miRNA* duplex. Pri-miRNA is a single-stranded RNA that forms a hairpin structure with a number of unpaired bases, hereafter called mismatches, on its stem. In the present study, by using purified recombinant Arabidopsis DCL1, we presented evidence that mismatches on the stem of pri-miRNA are important for precise DCL1cleavage. We showed that a mismatch at the loop-distal side of the end of miRNA/miRNA* duplex is important for efficient cleavage of pri-miRNA in vitro, as previously suggested in planta. On the contrary, mismatches distant from the miRNA/miRNA* duplex region are important for determining the cleavage position by DCL1. The purified DCL1 proteins cleaved mutant pri-miRNA variants without such mismatches at a position at which wild-type pri-miRNA variants are not usually cleaved, resulting in an increased accumulation of small RNA different from miRNA. Therefore, our results suggest that, in addition to the distance from the ssRNA-dsRNA junction, mismatches on the stem of pri-miRNA function as a determinant for precise processing of pri-miRNA by DCL1 in plants.

The modulation of DNA helical stability and geometry by the sequence and asymmetry of internal loops

Localized domains of unpaired nucleotides can arise in DNA and RNA through several pathways, including inherent sequence, processing errors, and protein interactions. For catalytic nucleic acids and DNA nanostructures, these domains have critical roles that include serving as functional sites and positioning helical domains; however, these unpaired domains can impart genomic instability as in the trinucleotide repeat disorders. Here, we present findings from experimental and computational work that address the role of DNA internal loop characteristics (sequence, size, and asymmetry) on the energetics and structural dynamics of loop-containing helical domains.

Determination of RNA Structure with In Vitro SHAPE Experiments.

The structure of an RNA molecule is often critical for its correct functioning, post-transcriptional regulation, and/or translation. Predicting RNA secondary structures with in silico tools is relatively straightforward with the large array of software and webservers available. However, for long RNAs and RNA at high temperatures, in silico predictions are less accurate and require experimental validation. To this end, a variety of structural probing reagents are commonly used, both for in vitro and in vivo mapping of RNA structure. Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) experiments make use of a nonbase-specific modifying reagent, acylating all conformationally flexible (mainly single-stranded or unpaired) nucleotides and have been employed both for in vitro and in vivo modification of RNA. Here, we describe a protocol for an in vitro SHAPE experiment, starting from in vitro transcribed RNA using a T7 polymerase system. This RNA is folded and subsequently modified in vitro with the SHAPE-reagent N-methyl isatoic anhydride (NMIA). Primer extension employing a radioactive 32P-labeled primer that binds the RNA downstream of the structure of interest will generate cDNA until the reverse transcriptase enzyme is halted by the introduced SHAPE modifications. Denaturing acrylamide gel electrophoresis of the pool of 32P-labeled cDNAs and the corresponding sequencing ladders, followed by autoradiography, will expose these stops in reverse transcription (RT) and will therefore enable to identify single-stranded nucleotides in the RNA of interest. These RT stops and NMIA-modification efficiencies can be quantified with ImageJ software and can be used to validate or increase the accuracy of RNA secondary structure predictions.