Pseudoknot Structure Research Articles

Gene editing with CRISPR-Cas12a guides possessing ribose-modified pseudoknot handles

CRISPR-Cas12a is a leading technology for development of model organisms, therapeutics, and diagnostics. These applications could benefit from chemical modifications that stabilize or tune enzyme properties. Here we chemically modify ribonucleotides of the AsCas12a CRISPR RNA 5′ handle, a pseudoknot structure that mediates binding to Cas12a. Gene editing in human cells required retention of several native RNA residues corresponding to predicted 2′-hydroxyl contacts. Replacing these RNA residues with a variety of ribose-modified nucleotides revealed 2′-hydroxyl sensitivity. Modified 5′ pseudoknots with as little as six out of nineteen RNA residues, with phosphorothioate linkages at remaining RNA positions, yielded heavily modified pseudoknots with robust cell-based editing. High trans activity was usually preserved with cis activity. We show that the 5′ pseudoknot can tolerate near complete modification when design is guided by structural and chemical compatibility. Rules for modification of the 5′ pseudoknot should accelerate therapeutic development and be valuable for CRISPR-Cas12a diagnostics.

Impacts of single nucleotide deletions from the 3′ end of Senecavirus A 5′ untranslated region on activity of viral IRES and on rescue of recombinant virus

The 5′ untranslated region (UTR) of Senecavirus A (SVA) harbors an internal ribosome entry site (IRES), in which a pseudoknot structure is upstream of start codon AUG. Wild-type SVAs have a highly conserved 13-nt-sequence between the pseudoknot stem II (PKS-II)-forming motif and the AUG. In this study, a single nucleotide was deleted one by one from the 13-nt-sequence within a wild-type SVA minigenome. The result showed that neither mono- nor multi-nucleotide deletions abolished the IRES activity. Furthermore, a single nucleotide was deleted one by one from the 13-nt-sequence within a full-length SVA cDNA clone. The result indicated that nucleotide-deleting SVAs could be rescued from 1- to 5-nt-deleting cDNA clones, whereas only the 1- and 2-nt-deleting viruses were genetically stable during nine serial passages in vitro. Additionally, only the 1-nt-deleting SVA showed similar growth kinetics to that of the wild-type virus, suggesting that the pseudoknot-AUG distance was crucial for SVA replication.

A Novel Frameshifting Inhibitor Having Antiviral Activity against Zoonotic Coronaviruses.

Recent outbreaks of zoonotic coronaviruses, such as Middle East respiratory syndrome coronavirus (MERS-CoV) and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), have caused tremendous casualties and great economic shock. Although some repurposed drugs have shown potential therapeutic efficacy in clinical trials, specific therapeutic agents targeting coronaviruses have not yet been developed. During coronavirus replication, a replicase gene cluster, including RNA-dependent RNA polymerase (RdRp), is alternatively translated via a process called -1 programmed ribosomal frameshift (−1 PRF) by an RNA pseudoknot structure encoded in viral RNAs. The coronavirus frameshifting has been identified previously as a target for antiviral therapy. In this study, the frameshifting efficiencies of MERS-CoV, SARS-CoV and SARS-CoV-2 were determined using an in vitro −1 PRF assay system. Our group has searched approximately 9689 small molecules to identify potential −1 PRF inhibitors. Herein, we found that a novel compound, 2-(5-acetylthiophen-2yl)furo[2,3-b]quinoline (KCB261770), inhibits the frameshifting of MERS-CoV and effectively suppresses viral propagation in MERS-CoV-infected cells. The inhibitory effects of 87 derivatives of furo[2,3-b]quinolines were also examined showing less prominent inhibitory effect when compared to compound KCB261770. We demonstrated that KCB261770 inhibits the frameshifting without suppressing cap-dependent translation. Furthermore, this compound was able to inhibit the frameshifting, to some extent, of SARS-CoV and SARS-CoV-2. Therefore, the novel compound 2-(5-acetylthiophen-2yl)furo[2,3-b]quinoline may serve as a promising drug candidate to interfere with pan-coronavirus frameshifting.

The SARS-CoV-2 Programmed -1 Ribosomal Frameshifting Element Crystal Structure Solved to 2.09 Å Using Chaperone-Assisted RNA Crystallography.

The programmed −1 ribosomal frameshifting element (PFSE) of SARS-CoV-2 is a well conserved structured RNA found in all coronaviruses’ genomes. By adopting a pseudoknot structure in the presence of the ribosome, the PFSE promotes a ribosomal frameshifting event near the stop codon of the first open reading frame Orf1a during translation of the polyprotein pp1a. Frameshifting results in continuation of pp1a via a new open reading frame, Orf1b, that produces the longer pp1ab polyprotein. Polyproteins pp1a and pp1ab produce nonstructural proteins NSPs 1–10 and NSPs 1–16, respectively, which contribute vital functions during the viral life cycle and must be present in the proper stoichiometry. Both drugs and sequence alterations that affect the stability of the −1 programmed ribosomal frameshifting element disrupt the stoichiometry of the NSPs produced, which compromise viral replication. For this reason, the −1 programmed frameshifting element is considered a promising drug target. Using chaperone assisted RNA crystallography, we successfully crystallized and solved the three-dimensional structure of the PFSE. We observe a three-stem H-type pseudoknot structure with the three stems stacked in a vertical orientation stabilized by two triple base pairs at the stem 1/stem 2 and stem 1/stem 3 junctions. This structure provides a new conformation of PFSE distinct from the bent conformations inferred from midresolution cryo-EM models and provides a high-resolution framework for mechanistic investigations and structure-based drug design.

Comparative study on tertiary contacts and folding of RNase P RNAs from a psychrophilic, a mesophilic/radiation-resistant, and a thermophilic bacterium

In most bacterial type A RNase P RNAs (P RNAs), two major loop-helix tertiary contacts (L8–P4 and L18–P8) help to orient the two independently folding S- and C-domains for concerted recognition of precursor tRNA substrates. Here, we analyze the effects of mutations in these tertiary contacts in P RNAs from three different species: (i) the psychrophilic bacterium Pseudoalteromonas translucida (Ptr), (ii) the mesophilic radiation-resistant bacterium Deinococcus radiodurans (Dra), and (iii) the thermophilic bacterium Thermus thermophilus (Tth). We show by UV melting experiments that simultaneous disruption of these two interdomain contacts has a stabilizing effect on all three P RNAs. This can be inferred from reduced RNA unfolding at lower temperatures and a more concerted unfolding at higher temperatures. Thus, when the two domains tightly interact via the tertiary contacts, one domain facilitates structural transitions in the other. P RNA mutants with disrupted interdomain contacts showed severe kinetic defects that were most pronounced upon simultaneous disruption of the L8–P4 and L18–P8 contacts. At 37°C, the mildest effects were observed for the thermostable Tth RNA. A third interdomain contact, L9–P1, makes only a minor contribution to P RNA tertiary folding. Furthermore, D. radiodurans RNase P RNA forms an additional pseudoknot structure between the P9 and P12 of its S-domain. This interaction was found to be particularly crucial for RNase P holoenzyme activity at near-physiological Mg2+ concentrations (2 mM). We further analyzed an exceptionally stable folding trap of the G,C-rich Tth P RNA.

18S rRNA variability maps reveal three highly divergent, conserved motifs within Rotifera

Background18S rRNA is a major component of the small subunit of the eukaryotic ribosome and an important phylogenetic marker for many groups, often to the point of being the only marker available for some. A core structure across eukaryotes exists for this molecule that can help to inform about its evolution in different groups. Using an alignment of 18S rDNA for Rotifera as traditionally recognized (=Bdelloidea, Monogononta, and Seisonacea, but not Acanthocephala), I fitted sequences for three exemplar species (Adineta vaga, Brachionus plicatilis, and Seison nebaliae, respectively) to the core structure and used these maps to reveal patterns of evolution for the remainder of this diverse group of microscopic animals.ResultsThe obtained variability maps of the 18S rRNA molecule revealed a pattern of high diversity among the three major rotifer clades coupled with strong conservation within each of bdelloids and monogononts. A majority of individual sites (ca. 60%) were constant even across rotifers as a whole with variable sites showing only intermediate rates of evolution. Although the three structural maps each showed good agreement with the inferred core structure for eukaryotic 18S rRNA and so were highly similar to one another at the secondary and tertiary levels, the overall pattern is of three highly distinct, but conserved motifs within the group at the primary sequence level. A novel finding was that of a variably expressed deletion at the 3' end of the V3 hypervariable region among some bdelloid species that occasionally extended into and included the pseudoknot structure following this region as well as the central “square” of the 18S rRNA molecule. Compared to other groups, levels of variation and rates of evolution for 18S rRNA in Rotifera roughly matched those for Gastropoda and Acanthocephala, despite increasing evidence for the latter being a clade within Rotifera.ConclusionsThe lack of comparative data for comparable groups makes interpretation of the results (i.e., very low variation within each of the three major rotifer clades, but high variation between them) and their potential novelty difficult. However, these findings in combination with the high morphological diversity within rotifers potentially help to explain why no clear consensus has been reached to date with regard to the phylogenetic relationships among the major groups.

Microsecond Folding of preQ1 Riboswitch and Its Biological Significance Revealed by Two-Dimensional Fluorescence Lifetime Correlation Spectroscopy.

Riboswitches are regulatory elements of bacterial mRNA which function with conformational switching upon binding of specific cellular metabolites. In particular, transcriptional riboswitches regulate gene expression kinetically through the conformational change of the aptamer domain. In this study, we investigate the conformational dynamics and ligand binding mechanisms of the aptamer domain of a transcriptional prequeuosine (preQ1) riboswitch from Bacillus subtilis using two-dimensional fluorescence lifetime correlation spectroscopy (2D FLCS) with microsecond time resolution. The obtained time-resolved single-molecule data indicate that the aptamer domain undergoes folding/unfolding including three forms, which are attributed to hairpin (O), pseudoknot-like (pF), and H-type pseudoknot (fF) structures. It is found that a cofactor, Mg2+, binds only to the fF form with the conformational selection mechanism. In contrast, it is indicated that the ligand, preQ1, binds to the O form with the induced-fit mechanism and significantly accelerates the microsecond O → pF folding process. It is also shown that the binding with preQ1 substantially stabilizes the fF form that is generated from the pF form with a long time constant (>10 ms). Combining these results with the results of a former smFRET study on the slower time scale, we obtain an overall picture of the folding/unfolding dynamics of the aptamer domain as well as its energy landscape. On the basis of the picture obtained, we discuss the significance of the microsecond folding/unfolding of the aptamer domain for biological function of the riboswitch and propose the molecular mechanism of the gene expression controlled by the structural dynamics of the aptamer domain.

Salt-Dependent RNA Pseudoknot Stability: Effect of Spatial Confinement.

Macromolecules, such as RNAs, reside in crowded cell environments, which could strongly affect the folded structures and stability of RNAs. The emergence of RNA-driven phase separation in biology further stresses the potential functional roles of molecular crowding. In this work, we employed the coarse-grained model that was previously developed by us to predict 3D structures and stability of the mouse mammary tumor virus (MMTV) pseudoknot under different spatial confinements over a wide range of salt concentrations. The results show that spatial confinements can not only enhance the compactness and stability of MMTV pseudoknot structures but also weaken the dependence of the RNA structure compactness and stability on salt concentration. Based on our microscopic analyses, we found that the effect of spatial confinement on the salt-dependent RNA pseudoknot stability mainly comes through the spatial suppression of extended conformations, which are prevalent in the partially/fully unfolded states, especially at low ion concentrations. Furthermore, our comprehensive analyses revealed that the thermally unfolding pathway of the pseudoknot can be significantly modulated by spatial confinements, since the intermediate states with more extended conformations would loss favor when spatial confinements are introduced.

Stem II-disrupting pseudoknot does not abolish ability of Senecavirus A IRES to initiate protein expression, but inhibits recovery of virus from cDNA clone

Senecavirus A (SVA) is classified into the genus Senecavirus in the family Picornaviridae. Its genome is a positive-sense, single-stranded and nonsegmented RNA, approximately 7300 nucleotides in length. A picornaviral genome is essentially an mRNA, which, albeit unmodified with 5′ cap structure, can still initiate protein expression by the internal ribosome entry site (IRES). The SVA genome contains a hepatitis C virus-like IRES, in which a pseudoknot structure plays an important role in initiating protein expression. In this study, we constructed a set of SVA (CH-LX-01-2016 strain) minigenomes with all combinations of point mutations in its pseudoknot stem II (PKS-II). The results showed that any combination of point mutations could not significantly interfere with the IRES to initiate protein expression. Further, we constructed a full-length SVA cDNA clone, in which the PKS-II-forming cDNA motif was subjected to site-directed mutagenesis for totally disrupting the PKS-II formation in IRES. Such a modified SVA cDNA clone was transfected into BSR-T7/5 cells, consequently demonstrating that the PKS-II-disrupting IRES interfered neither with protein expression nor with antigenome replication, whereas a competent SVA could not be rescued from the cDNA clone. It was speculated that the mutated motif possibly disrupted a packaging signal for virion assembly, therefore causing the failure of SVA rescue.

Structural Dynamics of SARS-CoV-2 Frameshift Signal Studied by Single-Molecule Force Spectroscopy Reveal Topologically Distinct Conformers

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causing the COVID-19 pandemic uses −1 programmed ribosomal frameshifting (−1 PRF) to control expression of key viral proteins. −1 PRF involves a shift in the reading frame of the ribosome at a specific location in the RNA message, stimulated by a pseudoknot structure in the mRNA located 5-7 nt downstream of the ‘slippery’ sequence where the reading-frame shift occurs. To understand how this pseudoknot responds to the mechanical tension applied by ribosomes during translation, which is thought to play a key role in stimulating frameshifting, we probed its structural dynamics under force using optical tweezers.

DNA-Tetrahedral-Nanostructure-Based Entropy-Driven Amplifier for High-Performance Photoelectrochemical Biosensing.

In virtue of the inherent molecular recognition and programmability, DNA has recently become the most promising for high-performance biosensors. The rationally engineered nucleic acid architecture will be very advantageous to hybridization efficiency, specificity, and sensitivity. Herein, a robust and split-mode photoelectrochemical (PEC) biosensor for miRNA-196a was developed based on an entropy-driven tetrahedral DNA (EDTD) amplifier coupled with superparamagnetic nanostructures. The DNA tetrahedron structure features in rigidity and structural stability that contribute to obtain precise identification units and specific orientations, improving the hybridization efficiency, sensitivity, and selectivity of the as-designed PEC biosensor. Further, superparamagnetic Fe3O4@SiO2@CdS particles integrated with DNA nanostructures are beneficial for the construction of a split-mode, highly selective, and reliable PEC biosensor. Particularly, the enzyme- and hairpin-free EDTD amplifier eliminates unnecessary interference from the complex secondary structure of pseudoknots or kissing loops in typical hairpin DNAs, significantly lowers the background noise, and improves the detection sensitivity. This PEC biosensor is capable of monitoring miRNA-196a in practical settings with additional advantages of efficient electrode fabrication, stability, and reproducibility. This strategy can be extended to various miRNA assays in complex biological systems with excellent performance.

A novel gray wolf optimizer with RNA crossover operation for tackling the non-parametric modeling problem of FCC process

This work proposes a novel gray wolf optimizer (GWO) with RNA crossover-operation and adaptive control parameter scheme (named as RNA-GWO). The RNA crossover-operation can better enhance the population diversity of RNA-GWO, it is designed according to the special pseudoknot structure of RNA molecule. The adaptive control parameter scheme is proposed to replace the linear one to balance the exploration and exploitation capacities during the optimization process of RNA-GWO. The effectiveness of RNA-GWO is verified by a suite of IEEE CEC 2017 benchmark functions. The experimental results prove that RNA-GWO can obtain higher accuracy solutions than the other six meta-heuristic algorithms. The RNA-GWO is also employed to solve the parameter problem of the wavelet neural network (WNN) non-parametric modeling method and applied to model the FCC process. The simulation results demonstrate that the RNA-GWO can provide favorable parameters for WNN, the model outputs of RNA-GWO optimized WNN can better agree with the experimental data of FCC process.

Full-length genome characterization of a new fusagravirus hosted by the spring orange peel fungus Caloscypha fulgens.

Mycoviruses widely exist in diverse lineages of fungi, yet there are only a few studies on mycovirus infection in uncultivated fungi. We here report the presence of a dsRNA mycovirus in saprotrophic spring orange peel fungus Caloscypha fulgens. A novel dsRNA virus, named "Caloscypha fulgens fusagravirus 1" (CfFV1), was isolated from a single ascocarp of C. fulgens, and its molecular features were revealed. The full-length cDNA of CfFV1 comprises 9,548 nucleotides with a calculated GC content of 47.9% and two discontinuous open reading frames (ORF 1 and 2). A-1 ribosomal frameshift region with two distinctive elements, including a canonical slippery heptanucleotide (AAAAAAC) and a pseudoknot structure, predicted as a Recoding Stimulatory Element, was detected in the junction region of ORF1 and ORF2. The deduced amino acid sequence of ORF1 and ORF2 showed the highest similarity to the putative structural protein and RNA-dependent RNA polymerase (RdRp) of Rosellinia necatrix fusagravirus 4 (RnFGV4). Genome organization, sequence similarity, and phylogenetic analysis indicate that this virus belongs to a new member of the proposed family Fusagraviridae. This is the first report of the presence of a mycovirus in the spring orange peel fungus C. fulgens. Keywords: mycovirus; dsRNA; proposed Fusagraviridae; uncultivated fungi; Caloscypha fulgens.

ATTfold: RNA Secondary Structure Prediction With Pseudoknots Based on Attention Mechanism.

Accurate RNA secondary structure information is the cornerstone of gene function research and RNA tertiary structure prediction. However, most traditional RNA secondary structure prediction algorithms are based on the dynamic programming (DP) algorithm, according to the minimum free energy theory, with both hard and soft constraints. The accuracy is particularly dependent on the accuracy of soft constraints (from experimental data like chemical and enzyme detection). With the elongation of the RNA sequence, the time complexity of DP-based algorithms will increase geometrically, as a result, they are not good at coping with relatively long sequences. Furthermore, due to the complexity of the pseudoknots structure, the secondary structure prediction method, based on traditional algorithms, has great defects which cannot predict the secondary structure with pseudoknots well. Therefore, few algorithms have been available for pseudoknots prediction in the past. The ATTfold algorithm proposed in this article is a deep learning algorithm based on an attention mechanism. It analyzes the global information of the RNA sequence via the characteristics of the attention mechanism, focuses on the correlation between paired bases, and solves the problem of long sequence prediction. Moreover, this algorithm also extracts the effective multi-dimensional features from a great number of RNA sequences and structure information, by combining the exclusive hard constraints of RNA secondary structure. Hence, it accurately determines the pairing position of each base, and obtains the real and effective RNA secondary structure, including pseudoknots. Finally, after training the ATTfold algorithm model through tens of thousands of RNA sequences and their real secondary structures, this algorithm was compared with four classic RNA secondary structure prediction algorithms. The results show that our algorithm significantly outperforms others and more accurately showed the secondary structure of RNA. As the data in RNA sequence databases increase, our deep learning-based algorithm will have superior performance. In the future, this kind of algorithm will be more indispensable.

Disruption of Zika Virus xrRNA1-Dependent sfRNA1 Production Results in Tissue-Specific Attenuated Viral Replication.

The Zika virus (ZIKV), like other flaviviruses, produces several species of sub-genomic RNAs (sfRNAs) during infection, corresponding to noncoding RNA fragments of different lengths that result from the exonuclease degradation of the viral 3′ untranslated region (UTR). Over the course of infection, these sfRNAs accumulate in the cell as a result of an incomplete viral genome degradation of the 3′ UTR by the host 5′ to 3′ exoribonuclease, Xrn1. The halting of Xrn1 in the 3′ UTR is due to two RNA pseudoknot structures in the 3′ UTR, termed exoribonuclease-resistant RNA1 and 2 (xrRNA1&2). Studies with related flaviviruses have shown that sfRNAs are important for pathogenicity and inhibiting both mosquito and mammalian host defense mechanisms. However, these investigations have not included ZIKV and there is very limited data addressing how sfRNAs impact infection in a whole animal model or specific tissues. In this study, we generate a sfRNA1-deficient ZIKV (X1) by targeted mutation in the xrRNA1 3′ UTR structure. We find that the X1 virus lacks the production of the largest ZIKV sfRNA species, sfRNA1. Using the X1 virus to infect adult Ifnar1−/− mice, we find that while the lack of sfRNA1 does not alter ZIKV replication in the spleen, there is a significant reduction of ZIKV genome replication in the brain and placenta compared to wild-type ZIKV infection. Despite the attenuated phenotype of the X1 ZIKV, mice develop a robust neutralizing antibody response. We conclude that the targeted disruption of xrRNA1 results in tissue-specific attenuation while still supporting robust neutralizing antibody responses. Future studies will need to investigate the tissue-specific mechanisms by which ZIKV sfRNAs influence infection and may utilize targeted xrRNA mutations to develop novel attenuated flavivirus vaccine approaches.

A combinatorial method to isolate short ribozymes from complex ribozyme libraries.

In vitro selections are the only known methods to generate catalytic RNAs (ribozymes) that do not exist in nature. Such new ribozymes are used as biochemical tools, or to address questions on early stages of life. In both cases, it is helpful to identify the shortest possible ribozymes since they are easier to deploy as a tool, and because they are more likely to have emerged in a prebiotic environment. One of our previous selection experiments led to a library containing hundreds of different ribozyme clusters that catalyze the triphosphorylation of their 5′-terminus. This selection showed that RNA systems can use the prebiotically plausible molecule cyclic trimetaphosphate as an energy source. From this selected ribozyme library, the shortest ribozyme that was previously identified had a length of 67 nucleotides. Here we describe a combinatorial method to identify short ribozymes from libraries containing many ribozymes. Using this protocol on the library of triphosphorylation ribozymes, we identified a 17-nucleotide sequence motif embedded in a 44-nucleotide pseudoknot structure. The described combinatorial approach can be used to analyze libraries obtained by different in vitro selection experiments.

Crystal structure and ligand-induced folding of the SAM/SAH riboswitch.

While most SAM riboswitches strongly discriminate between SAM and SAH, the SAM/SAH riboswitch responds to both ligands with similar apparent affinities. We have determined crystal structures of the SAM/SAH riboswitch bound to SAH, SAM and other variant ligands at high resolution. The riboswitch forms an H-type pseudoknot structure with coaxial alignment of the stem–loop helix (P1) and the pseudoknot helix (PK). An additional three base pairs form at the non-open end of P1, and the ligand is bound at the interface between the P1 extension and the PK helix. The adenine nucleobase is stacked into the helix and forms a trans Hoogsteen–Watson–Crick base pair with a uridine, thus becoming an integral part of the helical structure. The majority of the specific interactions are formed with the adenosine. The methionine or homocysteine chain lies in the groove making a single hydrogen bond, and there is no discrimination between the sulfonium of SAM or the thioether of SAH. Single-molecule FRET analysis reveals that the riboswitch exists in two distinct conformations, and that addition of SAM or SAH shifts the population into a stable state that likely corresponds to the form observed in the crystal. A model for translational regulation is presented whereby in the absence of ligand the riboswitch is largely unfolded, lacking the PK helix so that translation can be initiated at the ribosome binding site. But the presence of ligand stabilizes the folded conformation that includes the PK helix, so occluding the ribosome binding site and thus preventing the initiation of translation.

Hierarchical Conformational Dynamics Confers Thermal Adaptability to preQ1 RNA Riboswitches

Single-stranded noncoding regulatory RNAs, as exemplified by bacterial riboswitches, are highly dynamic. The conformational dynamics allow the riboswitch to reach maximum switching efficiency under appropriate conditions. Here we characterize the conformational dynamics of preQ1 riboswitches from mesophilic and thermophilic bacterial species at various temperatures. With the integrative use of small-angle X-ray scattering, NMR, and molecular dynamics simulations, we model the ensemble-structures of the preQ1 riboswitch aptamers without or with a ligand bound. We show that the preQ1 riboswitch is sufficiently dynamic and fluctuating among multiple folding intermediates only near the physiological temperature of the microorganism. The hierarchical folding dynamics of the RNA involves the docking of 3′-tail to form a second RNA helix and the helical stacking to form an H-type pseudoknot structure. Further, we show that RNA secondary and tertiary dynamics can be modulated by temperature and by the length of an internal loop. The coupled equilibria between RNA folding intermediates are essential for preQ1 binding, and a four-state exchange model can account for the change of ligand-triggered switching efficiency with temperature. Together, we have established a relationship between the hierarchical dynamics and riboswitch function, and illustrated how the RNA adapts to high temperature.

The conformational landscape of transcription intermediates involved in the regulation of the ZMP-sensing riboswitch from Thermosinus carboxydivorans.

Recently, prokaryotic riboswitches have been identified that regulate transcription in response to change of the concentration of secondary messengers. The ZMP (5-Aminoimidazole-4-carboxamide ribonucleotide (AICAR))-sensing riboswitch from Thermosinus carboxydivorans is a transcriptional ON-switch that is involved in purine and carbon-1 metabolic cycles. Its aptamer domain includes the pfl motif, which features a pseudoknot, impeding rho-independent terminator formation upon stabilization by ZMP interaction. We herein investigate the conformational landscape of transcriptional intermediates including the expression platform of this riboswitch and characterize the formation and unfolding of the important pseudoknot structure in the context of increasing length of RNA transcripts. NMR spectroscopic data show that even surprisingly short pre-terminator stems are able to disrupt ligand binding and thus metabolite sensing. We further show that the pseudoknot structure, a prerequisite for ligand binding, is preformed in transcription intermediates up to a certain length. Our results describe the conformational changes of 13 transcription intermediates of increasing length to delineate the change in structure as mRNA is elongated during transcription. We thus determine the length of the key transcription intermediate to which addition of a single nucleotide leads to a drastic drop in ZMP affinity.

Comparison of Pseudoknotted RNA Secondary Structures by Topological Centroid Identification and Tree Edit Distance

Comparison of RNA structures is one of the most crucial analysis for elucidating their individual functions and promoting medical applications. Because it is widely accepted that their functions and structures are strongly correlated, various methods for RNA secondary structure analysis have been proposed owing to the difficulty in predicting RNA three-dimensional structure directly from its sequence. However, there are few methods dealing with RNA secondary structures with a specific and complex partial structure called pseudoknot despite its significance to biological process, which is a big obstacle for analyzing their functions. In this study, we propose a novel tree representation of pseudoknotted RNA secondary structures by topological centroid identification and their comparison methods based on the tree edit distance. In the proposed method, a given graph representing an RNA secondary structure is transformed to a tree rooted at one of the vertices constituting the topological centroid that is identified by removing cycles with peeling processing for the graph. When comparing tree-represented RNA secondary structures collected from a public database using the tree edit distance and functional gene groups defined by Gene Ontology (GO), the proposed method showed better clustering results according to their GOs than canonical RNA sequence-based comparison. In addition, we also report a case that the combination of the tree edit distance and the sequence edit distance shows a better classification of the pseudoknotted RNA secondary structures.