Search

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.

Initial successes in developing small molecule ligands for non-coding RNAs have underscored their potential as therapeutic targets. More recently, these successes have been aided by advances in biophysical and structural techniques for identification and characterization of more complex RNA structures; these higher-level folds present protein-like binding pockets that offer opportunities to design small molecules that could achieve a degree of selectivity often hard to obtain at the primary and secondary structure level. More specifically, identification of RNA tertiary and quaternary structures and targeting these structures with small molecule ligands have allowed researchers to probe several human diseases and have resulted in promising clinical candidates. In this review we highlight a selection of diverse and exciting successes and the experimental approaches that led to their discovery. These studies include examples of recent developments in RNA-centric assays and ligands that provide insight into the features responsible for the affinity and biological outcome of RNA-targeted chemical probes. This report highlights the potential and emerging opportunities to selectively target RNA tertiary and quaternary structures as a route to better understand and, ultimately, treat many diseases.

With the rapid rate of COVID-19 infections and deaths, treatments and cures besides hand washing, social distancing, masks, isolation, and quarantines are urgently needed. The treatments and vaccines rely on the basic biophysics of the complex viral apparatus. Although proteins are serving as main drug and vaccine targets, therapeutic approaches targeting the 30,000 nucleotide RNA viral genome form important complementary approaches. Indeed, the high conservation of the viral genome, its close evolutionary relationship to other viruses, and the rise of gene editing and RNA-based vaccines all argue for a focus on the RNA agent itself. One of the key steps in the viral replication cycle inside host cells is the ribosomal frameshifting required for translation of overlapping open reading frames. The RNA frameshifting element (FSE), one of three highly conserved regions of coronaviruses, is believed to include a pseudoknot considered essential for this ribosomal switching. In this work, we apply our graph-theory-based framework for representing RNA secondary structures, "RAG (or RNA-As-Graphs)," to alter key structural features of the FSE of the SARS-CoV-2 virus. Specifically, using RAG machinery of genetic algorithms for inverse folding adapted for RNA structures with pseudoknots, we computationally predict minimal mutations that destroy a structurally important stem and/or the pseudoknot of the FSE, potentially dismantling the virus against translation of the polyproteins. Our microsecond molecular dynamics simulations of mutant structures indicate relatively stable secondary structures. These findings not only advance our computational design of RNAs containing pseudoknots, they pinpoint key residues of the SARS-CoV-2 virus as targets for antiviral drugs and gene editing approaches.

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.

Toxoplasma gondii is an obligate intracellular parasite that causes one of the most common parasitic infections in humans and other warm-blooded animals. Currently, there are no effective treatments for inhibiting the formation of chronic tissue cysts in infected hosts. Thus, the development of a vaccine to protect against toxoplasmosis is an attractive option for avoiding infection. The aim of this study was to design an epitope-based vaccine for T. gondii. In the present study, an in silico approach was used to predict and analyze B-cell and T-cell epitopes and the transmembrane domain of proteins SAG1, MIC3, and ROP8. We also predicted the antigenicity, allergenicity, secondary and tertiary structures, and physicochemical characteristics of a chimeric protein. Next, codon optimization and mRNA structure prediction were conducted using bioinformatics tools, and the designed construct was chemically synthesized and cloned into the pET28a vector. SAG1 (amino acid positions 85–235), MIC3 (30–180), and ROP8 (85–185) were found to have several strong immunodominant epitopes that were joined with a rigid linker A(EAAAK)2A. Although the resultant protein called MRS ( M IC3, R OP8, and S AG1) did not turn out to be an allergen, its antigenicity was estimated to be 0.7983. Additionally, MRS was selected as the best vaccine candidate on the basis of its secondary and tertiary structures. The number of amino acids, molecular weight, and numbers of negatively and positively charged residues of MRS were 427 and 45,661.31 Da, 45, and 50, respectively. ΔG of the best-predicted structure was −413.0 kcal/mol, and the first nucleotides at the 5′ end did not form a stable hairpin or pseudoknot. Finally, successful expression and verification of the expressed MRS protein showed that in silico analysis was almost accurate. This vaccine candidate selected by in silico tools should be validated in experimental studies.

Approximately 17 years after the severe acute respiratory syndrome coronavirus (SARS-CoV) epidemic, the world is currently facing the COVID-19 pandemic caused by SARS corona virus 2 (SARS-CoV-2). According to the most optimistic projections, it will take more than a year to develop a vaccine, so the best short-term strategy may lie in identifying virus-specific targets for small molecule-based interventions. All coronaviruses utilize a molecular mechanism called programmed -1 ribosomal frameshift (-1 PRF) to control the relative expression of their proteins. Previous analyses of SARS-CoV have revealed that it employs a structurally unique three-stemmed mRNA pseudoknot that stimulates high -1 PRF rates and that it also harbors a -1 PRF attenuation element. Altering -1 PRF activity impairs virus replication, suggesting that this activity may be therapeutically targeted. Here, we comparatively analyzed the SARS-CoV and SARS-CoV-2 frameshift signals. Structural and functional analyses revealed that both elements promote similar -1 PRF rates and that silent coding mutations in the slippery sites and in all three stems of the pseudoknot strongly ablate -1 PRF activity. We noted that the upstream attenuator hairpin activity is also functionally retained in both viruses, despite differences in the primary sequence in this region. Small-angle X-ray scattering analyses indicated that the pseudoknots in SARS-CoV and SARS-CoV-2 have the same conformation. Finally, a small molecule previously shown to bind the SARS-CoV pseudoknot and inhibit -1 PRF was similarly effective against -1 PRF in SARS-CoV-2, suggesting that such frameshift inhibitors may be promising lead compounds to combat the current COVID-19 pandemic.

Senecavirus A (SVA), formerly known as Seneca Valley virus, is classified into the genus Senecavirus in the family Picornaviridae. Mature virion harbors an approximately 7 300-nt-long, positive-sense, and single-stranded RNA genome, which contains 5' and 3' untranslated regions (UTRs). Internal ribosome entry site (IRES) is identified in the SVA 5' UTR, and includes a RNA pseudoknot upstream of the start codon. This pseudoknot contains two stem structures, pseudoknot stem I and II (PKS-I and -II). The PKS-I is composed of two base-paired motifs (PKS-Ia and -Ib), between which there is an unpaired spacing (UpS). We reported previously that motif mutation in the PKS-II did not abolish the IRES activity, but interfered with SVA recovery from cDNA clone. In this study, we constructed five SVA minigenomes with point mutations in the PKS-I motif. Dual-luciferase reporter assay showed that motif mutations in PKS-I did not significantly interfere with the IRES activity to initiate protein expression. Correspondingly, we constructed five SVA cDNA clones with point mutations in the PKS-I motif. These genetically modified cDNA clones were separately transfected into BSR-T7/5 cells in attempting to rescue competent SVAs. However, only two viruses, namely PKS-Ia- and UpS-mutated recombinants, could be recovered from their individual cDNA clones. It can be concluded that the PKS-Ib is indispensable for viral growth.

We investigate and compare the refolding and cotranscriptional folding process of RNA pseudoknot including PK5 pseudoknot, human telomerase RNA (hTR) pseudoknot and its mutants by using helix-based transition rate model and master equation. It is found that in the refolding process, the final states for PKWT and two mutants PK179AG/110CU and PK△U177 are all pseudoknots. However, the final state of mutant PKDC is a stable long helix. Furthermore, the cotranscriptional folding intermediates are quite different from refolding ones. The intermediate found in refolding process for PKWT is also shown during cotranscription. For PK179AG/110CU, the two intermediates in refolding are not found in cotranscriptional folding. For PK△U177, only one of the intermediates found in refolding is not observed in cotranscription. For PKDC, the unfolding sequence in the refolding process directly transits to native helix. Two intermediates are found during cotranscription. The populations of intermediates and native structure also depend on the cotranscriptional folding rate. All these findings can offer useful insights into cotranscription folding of hTR pseudoknot with different rates.

Prostate specific antigen (PSA) has become a potential biomarker for detecting prostate cancer (PCa) in the early stage. Herein, we report a target-induced resolution for the detection of PSA sensitively and specifically by amperometric electrochemical measurements. To meet a satisfactory performance, three conformations of pre-design DNA aptamers including two stem-loop structures and a double strand structure have been investigated and compared. All of them are immobilized on gold electrode as capture probes with redox-active molecular. The mechanism of signal transduction depends on molecular recognition events involving aptamer conformational changes, thus influencing the charge transfer. A short, single-stranded DNA (ssDNA) pseudoknot forming two stem-loop structural aptamers with labeled MB at the 3′ -terminus was found to posse the highest signal variation than other structure when induced by PSA due to the strong conformational change. With the optimized capture strand, the aptasensor showed the peak current increase of MB by the binding relationship between PSA and the sensor over a wide concentration range of 4 magnitude orders. The proposed aptasensor exhibited a wide detection range from 10 pg/mL to 500 ng/mL with a low detection limit of 1.24 pg/mL (S/N = 3). Moreover, the electrochemical aptasensor demonstrated good reproducibility, sensitivity, selectivity, and reliability for the detection of PSA. We also found the aptasensor had a good response in the human serum samples, making this device easy to operate for the detection of the PSA physiological concentration.

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.

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.

Because of the potential link between -1 programmed ribosomal frameshifting and response of a pseudoknot (PK) RNA to force, a number of single molecule pulling experiments have been performed on PKs to decipher the mechanism of programmed ribosomal frameshifting. Motivated in part by these experiments, we performed simulations using a coarse-grained model of RNA to describe the response of a PK over a range of mechanical forces ($f$s) and monovalent salt concentrations ($C$s). The coarse-grained simulations quantitatively reproduce the multistep thermal melting observed in experiments, thus validating our model. The free energy changes obtained in simulations are in excellent agreement with experiments. By varying $f$ and $C$, we calculated the phase diagram that shows a sequence of structural transitions, populating distinct intermediate states. As $f$ and $C$ are changed, the stem-loop tertiary interactions rupture first, followed by unfolding of the $3^{\prime}$-end hairpin ($\textrm{I}\rightleftharpoons\textrm{F}$). Finally, the $5^{\prime}$-end hairpin unravels, producing an extended state ($\textrm{E}\rightleftharpoons\textrm{I}$). A theoretical analysis of the phase boundaries shows that the critical force for rupture scales as $\left(\log C_{\textrm{m}}\right)^{\alpha}$ with $\alpha=1\,(0.5)$ for $\textrm{E}\rightleftharpoons\textrm{I}$ ($\textrm{I}\rightleftharpoons\textrm{F}$) transition. This relation is used to obtain the preferential ion-RNA interaction coefficient, which can be quantitatively measured in single-molecule experiments, as done previously for DNA hairpins. A by-product of our work is the suggestion that the frameshift efficiency is likely determined by the stability of the $5^{\prime}$-end hairpin that the ribosome first encounters during translation., Comment: Final draft accepted in Journal of Molecular Biology, 16 pages including Supporting Information

Metabolite-dependent conformational switching in RNA riboswitches is now widely accepted as a critical regulatory mechanism for gene expression in bacterial systems. More recently, similar gene regulation mechanisms have been found to be important for viral systems as well. One of the most abundant and best-studied systems is the tRNA-like structure (TLS) domain, which has been found to occur in many plant viruses spread across numerous genera. In this work, folding dynamics for the TLS domain of Brome Mosaic Virus have been investigated using single-molecule fluorescence resonance energy transfer techniques. In particular, burst fluorescence methods are exploited to observe metal-ion ([Mn+])-induced folding in freely diffusing RNA constructs resembling the minimal TLS element of brome mosaic virus RNA3. The results of these experiments reveal a complex equilibrium of at least three distinct populations. A stepwise, or consecutive, thermodynamic model for TLS folding is developed, which is in good agreement with the [Mn+]-dependent evolution of conformational populations and existing structural information in the literature. Specifically, this folding pathway explains the metal-ion dependent formation of a functional TLS domain from unfolded RNAs via two consecutive steps: 1) hybridization of a long-range stem interaction, followed by 2) formation of a 3′-terminal pseudoknot. These two conformational transitions are well described by stepwise dissociation constants for [Mg2+] (K1 = 328 ± 30 μM and K2 = 1092 ± 183 μM) and [Na+] (K1 = 74 ± 6 mM and K2 = 243 ± 52 mM)-induced folding. The proposed thermodynamic model is further supported by inhibition studies of the long-range stem interaction using a complementary DNA oligomer, which effectively shifts the dynamic equilibrium toward the unfolded conformation. Implications of this multistep conformational folding mechanism are discussed with regard to regulation of virus replication.

RNA folding dynamics plays important roles in various functions of RNAs. To date, coarse-grained modeling has been successfully employed to simulate RNA folding dynamics on the energy landscape composed of secondary structures. In such a modeling, the energy barrier height between metastable structures is a key parameter that crucially affects the simulation results. Although a number of approaches ranging from the exact method to heuristic ones are available to predict the barrier heights, developing an efficient heuristic for this purpose is still an algorithmic challenge. We developed a novel RNA folding pathway prediction method, ACOfoldpath, based on Ant Colony Optimization (ACO). ACO is a widely used powerful combinatorial optimization algorithm inspired from the food-seeking behavior of ants. In ACOfoldpath, to accelerate the folding pathway prediction, we reduce the search space by utilizing originally devised structure generation rules. To evaluate the performance of the proposed method, we benchmarked ACOfoldpath on the known nineteen conformational RNA switches. As a result, ACOfoldpath successfully predicted folding pathways better than or comparable to the previous heuristics. The results of RNA folding dynamics simulations and pseudoknotted pathway predictions are also presented.

Summary ZTP, the pyrophosphorylated analog of ZMP (5-amino-4-imidazole carboxamide ribose-5′-monophosphate), was identified as an alarmone that senses 10-formyl-tetrahydroflate deficiency in bacteria. Recently, a pfl riboswitch was identified that selectively binds ZMP and regulates genes associated with purine biosynthesis and one-carbon metabolism. We report on the structure of the ZMP-bound Thermosinus carboxydivorans pfl riboswitch sensing domain, thereby defining the pseudoknot-based tertiary RNA fold, the binding-pocket architecture, and principles underlying ligand recognition specificity. Molecular recognition involves shape complementarity, with the ZMP 5-amino and carboxamide groups paired with the Watson-Crick edge of an invariant uracil, and the imidazole ring sandwiched between guanines, while the sugar hydroxyls form intermolecular hydrogen bond contacts. The burial of the ZMP base and ribose moieties, together with unanticipated coordination of the carboxamide by Mg 2+ , contrasts with exposure of the 5′-phosphate to solvent. Our studies highlight the principles underlying RNA-based recognition of ZMP, a master regulator of one-carbon metabolism.

SummaryThe ZTP riboswitch is a widespread family of regulatory RNAs that upregulate de novo purine synthesis in response to increased intracellular levels of ZTP or ZMP. As an important intermediate in purine biosynthesis, ZMP also serves as a proxy for the concentration of N10-formyl-tetrahydrofolate, a key component of one-carbon metabolism. Here, we report the structure of the ZTP riboswitch bound to ZMP at a resolution of 1.80 Å. The RNA contains two subdomains brought together through a long-range pseudoknot further stabilized through helix-helix packing. ZMP is bound at the subdomain interface of the RNA through a set of interactions with the base, ribose sugar, and phosphate moieties of the ligand. Unique to nucleobase recognition by RNAs, the Z base is inner-sphere coordinated to a magnesium cation bound by two backbone phosphates. This interaction, along with steric hindrance by the backbone, imparts specificity over chemically similar compounds such as ATP/AMP.

A two-level particle swarm optimization (TL-PSO) algorithm is proposed for training stochastic context-sensitive hidden Markov model (cs-HMM), that addresses a thrust area of bioinformatics i.e. structural alignment of pseudoknotted non-coding RNAs (ncRNAs). Due to the well-conserved sequences and corresponding secondary structures of ncRNAs, the structural information becomes imperative for performing their alignments. Proposed approach is unique in the sense: it is the first idea so far which works on optimization of the model length; also it is the first swarm intelligence technique that is proposed for training csHMM. The two-level strategy with training and cross training sets helps in increasing the diversity of the particles so as to avoid trapping in local optima, yields more accurate estimation parameters, preserves the structure of the model and provides the best compression from the model. TL-PSO yields a trained stochastic model with position-dependent probabilities that achieves high prediction ratios than the compared non-stochastic scoring matrix based csHMM approaches. TL-PSO is also tested solely for sequence alignment of proteins, by training the conventional HMMs. TLPSO-HMM produces an effective framework for sequence alignment in terms of alignment quality and prediction accuracy than the competitive state-of-the-art and family of PSO based algorithms. Conjointly, TLPSO-csHMM finds higher prediction measures than competitive RNA structural alignment techniques for pseudoknotted and non-pseudoknotted RNA structures of diverse complexities.

RNA interactions are fundamental in many cellular processes, where two or more RNA molecules can be involved. Multiple RNA interactions are also believed to be much more complex than pairwise interactions. Recently, multiple RNA interaction prediction has been formulated as a maximization problem. Here we extensively examine this optimization problem under several biologically meaningful interaction models. We present a polynomial time algorithm for the problem when the order of interacting RNAs is known and pseudoknot interactions are allowed; for the general problem without an assumed RNA order, we prove the NP-hardness for both variants (allowing and disallowing pseudoknot interactions), and present a constant ratio approximation algorithm for each of them.

Summary CRISPR-Cas12a (Cpf1), a type V CRISPR-associated nuclease, provides bacterial immunity against bacteriophages and plasmids but also serves as a tool for genome editing. Foreign nucleic acids are integrated into the CRISPR locus, prompting transcription of CRISPR RNAs (crRNAs) that guide Cas12a cleavage of foreign complementary DNA. However, mobile genetic elements counteract Cas12a with inhibitors, notably type V-A anti-CRISPRs (AcrVAs). We present cryoelectron microscopy structures of Cas12a-crRNA bound to AcrVA1 and AcrVA4 at 3.5 and 3.3 A resolutions, respectively. AcrVA1 is sandwiched between the recognition (REC) and nuclease (NUC) lobes of Cas12a and inserts into the binding pocket for the protospacer-adjacent motif (PAM), a short DNA sequence guiding Cas12a targeting. AcrVA1 cleaves crRNA in a Cas12a-dependent manner, inactivating Cas12a-crRNA complexes. The AcrVA4 dimer is anchored around the crRNA pseudoknot of Cas12a-crRNA, preventing required conformational changes for crRNA-DNA heteroduplex formation. These results uncover molecular mechanisms for CRISPR-Cas12a inhibition, providing insights into bacteria-phage dynamics.

Pseudoknot formation in the core region of the telomerase RNA has been demonstrated to be important for telomerase activity in vertebrates, ciliates, and yeast. Characterization of the Saccharomyces cerevisiae telomerase RNA (TLC1) pseudoknot identified tertiary structural interactions that are also important for telomerase activity, as previously observed for the Kluyveromyces lactis and human telomerase RNA pseudoknots. In addition, the contributions of backbone ribose 2'-OH groups in the pseudoknot to telomerase catalysis were investigated previously, using 2'-OH (ribose) to 2'-H (deoxyribose) or 2'-O-methyl substitutions in the stem 2 helix, and it was proposed that one or more 2'-OH groups from the stem 2 sequences at or near the triple helix participate in telomerase catalysis. Based on these studies and investigations of the structural and thermodynamic properties of the TLC1 RNA pseudoknot region, we have examined the structural and thermodynamic perturbations of the 2'-O-methyl and 2'-H substituted pseudoknots, using UV-monitored thermal denaturation, native gel electrophoresis, and circular dichroism spectroscopy. Our results demonstrate the presence of A-form helical geometry perturbations in the backbone sugar substituted pseudoknots, show a correlation between thermodynamic stability and telomerase activity, and are consistent with the identification of the U809 ribose 2'-OH as a potential contributor to telomerase activity.

Mechanical resistance of pseudoknot (PK) in mRNA deemed relevant in stimulating programmed ribosomal frameshifting (PRF) has lead to a number of single molecule pulling experiments in order to decipher the PRF mechanism. Motivated by these studies, we performed simulations to describe the response of a Beet Western Yellow Virus (BWYV) PK over a range of mechanical forces (fs) and monovalent salt concentrations (Cs) using a coarse-grained model. The simulations based on the model quantitatively reproduces the three-stage thermal melting of BWYV PK and the stability of the folded state (F) with respect to the unfolded state (U) inferred in ensemble experiments setting the stage for making testable predictions. The predicted phase diagram in the [C, f] plane shows that a sequence of structural transitions populating intermediates occur as f and C are changed. The least stable tertiary interactions rupture first followed by unfolding of the 3′-end hairpin. Finally, the most stable 5′-end hairpin unravels producing a stretched state. A theoretical analysis of the phase boundaries, separating the F and an intermediate state (I), and I and U states shows that the critical force for rupture scales as (log Cm)α with α=1 (=0.5) for U−I (I−F) transition. This relation is used to obtain the ion-preferential coefficients, which suggests that ion-RNA interactions can be quantitatively measured using laser optical tweezers experiments. We predict that −1 PRF efficiency should be determined by the stability of the 5′ hairpin that the ribosome first encounters. A good correlation between the increase in −1 PRF efficiency and increase in the stability of 5′ hairpin supports this prediction.

Most cellular RNAs engage in intrastrand base-pairing that gives rise to complex three-dimensional folds. This self-pairing presents an impediment toward binding of the RNA by nucleic acid-based ligands. An important step in the discovery of RNA-targeting ligands is therefore to identify those regions in a folded RNA that are accessible toward the nucleic acid-based ligand. Because the folding of RNA targets can involve interactions between nonadjacent regions and employ both Watson-Crick and non-Watson-Crick base-pairing, screening of candidate binder ensembles is typically necessary. Microarray-based screening approaches have shown great promise in this regard and have suggested that achieving complete sequence coverage would be a valuable attribute of a next generation system. Here, we report a custom microarray displaying a library of RNA-interacting polynucleotides comprising all possible 2'-OMe RNA sequences from 4- to 8-nucleotides in length. We demonstrate the utility of this array in identifying RNA-interacting polynucleotides that bind tightly and specifically to the highly conserved, functionally essential template/pseudoknot domain of human telomerase RNA and that inhibit telomerase function in vitro.

Kissing hairpin interactions form when the loop residues of two hairpins have Watson-Crick complementarity. In a unimolecular context, kissing interactions are important for tertiary folding and pseudoknot formation, whereas in a bimolecular context, they provide a basis for molecular recognition. In some cases, kissing complexes can be a prelude to strand displacement reactions where the two hairpins resolve to form a stable extended intermolecular duplex. The kinetics and thermodynamics of kissing-complex formation and their subsequent strand-displacement reactions are poorly understood. Here, biophysical techniques including isothermal titration calorimetry, surface plasmon resonance, and single-molecule fluorescence have been employed to probe the factors that govern the stability of kissing complexes and their subsequent structural rearrangements. We show that the general understanding of RNA duplex formation can be extended to kissing complexes but that kissing complexes display an unusual level of stability relative to simple duplexes of the same sequence. These interactions form and break many times at room temperature before becoming committed to a slow, irreversible forward transition to the strand-displaced form. Furthermore, using smFRET we show that the primary difference between stable and labile kissing complexes is based almost completely on their off rates. Both stable and labile complexes form at the same rate within error, but less stable species dissociate rapidly, allowing us to understand how these complexes can help generate specificity along a folding pathway or during a gene regulation event.

RNAs can fold into various conformations in their biological mechanisms. Elucidation of RNA folding processes is important to the understanding and regulation of RNA biological functions. Here we developed a nanopore electric snapshot technology for probing the folding pathway of RNA tertiary structures at room temperature. Our target is the gene 32 mRNA pseudoknot of bacteriophage T2. Driven by a transmembrane voltage applied, the target RNA can be captured by the nanopore from one side, and immobilized in the pore via the attached streptavidin.

In a recent issue of Science, Akiyama et al. (2016) prove the existence of a pseudoknot that stabilizes a nuclease-resistant RNA structure in the 3’ untranslated region of Zika virus. This reinforced structure blocks the 5’→3’ exonuclease Xrn1 for the production of pathogenic subgenomic flaviviral RNAs.

We present a new algorithm to predict RNA secondary structures with pseudoknots. The method is based on a classification of RNA structures according to their topological genus. The algorithm utilizes a simplified parametrization of the free energies for pair stacking, loop penalties, etc. and in addition a free energy penalty proportional to the topological genus of the pairing graph. Our method can take into account all pseudoknot topologies and achieves high success rates compared to state-of-the-art methods. This shows that the genus is a promising concept to classify pseudoknots.

Turnip yellow mosaic virus (TYMV) contains a tRNA-like structure (TLS) in its 3′ untranslated region (3′ UTR). This highly structured element induces valylation of the viral RNA by host cell enzymes and is important for virus proliferation. Directly upstream of the TYMV TLS is an upstream pseudoknot domain (UPD) that has been considered to be structurally distinct from the TLS. However, using a combination of functional, biochemical, and biophysical assays, we show that the entire 3′ UTR of the viral genome is a single structured element in the absence of cellular protein. This packing architecture stabilizes the RNA structure and creates a better substrate for aminoacylation, and thus the UPD and TLS are functionally and structurally coupled. It has been proposed that the TYMV TLS acts as a molecular switch between translation and replication. Our results suggest that this putative switch could be based on structural changes within the global architecture of the UTR induced by interactions with the ribosome. The TYMV TLS·UPD might demonstrate how RNA structural plasticity can play a role in regulation of biological processes.

In this paper, we enumerate k-noncrossing RNA pseudoknot structures with a given minimum arc- and stack-length. That is, we study the numbers of RNA pseudoknot structures with arc-length >=3, stack-length >=@s and in which there are at most k-1 mutually crossing bonds, denoted by T"k","@s^[^3^](n). We prove that the numbers ofk-noncrossing RNA structures with arc-length >=3 and stack-length >=2 satisfy T"k","2^[^3^](n)^~C"kn^-^(^k^-^1^)^^^2^-^k^-^1^2(@c"k","2^[^3^])^-^n. In the case k=3,4,5, we derive T"3","2^[^3^](n)^~C"3n^-^52.5721^n, T"4","2^[^3^](n)~C"4n^-^2^1^23.0306^n, and T"5","2^[^3^](n)~C"5n^-^1^83.4092^n, respectively, where C"3,C"4,C"5 are constants. Our results are of importance for prediction algorithms for RNA pseudoknot structures.

RNA is known to exhibit a jerky dynamics, as intramolecular thermal motion, on0.1 micros time scales, is punctuated by infrequent structural rearrangements on much longer time scales, i.e. from10 micros up to a few minutes or even hours. These rare stochastic events correspond to the formation or dissociation of entire stems through cooperative base pairing/unpairing transitions. Such a clear separation of time scales in RNA dynamics has made it possible to implement coarse grained RNA simulations, which predict RNA folding and unfolding pathways including kinetically trapped structures on biologically relevant time scales of seconds to minutes. RNA folding simulations also enable to predict the formation of pseudoknots, that is, helices interior to loops, which mechanically restrain the relative orientations of other non-nested helices. But beyond static structural constraints, pseudoknots can also strongly affect the folding and unfolding dynamics of RNA, as the order by which successive helices are formed and dissociated can lead to topologically blocked transition intermediates. The resulting knotty dynamics can enhance the stability of RNA switches, improve the efficacy of co-transcriptional folding pathways and lead to unusual self-assembly properties of RNA.

Summary The modified nucleotide queuosine (Q) is almost universally found in the anticodon wobble position of specific tRNAs. In many bacteria, biosynthesis of Q is modulated by a class of regulatory mRNA elements called riboswitches. The preQ 1 riboswitch, found in the 5′UTR of bacterial genes involved in synthesis of the Q precursors preQ 0 and preQ 1 , contains the smallest known aptamer domain. We report the solution structure of the preQ 1 riboswitch aptamer domain from Bacillus subtilis bound to preQ 1 , which is a unique compact pseudoknot with three loops and two stems that encapsulates preQ 1 at the junction between the two stems. The pseudoknot only forms in the presence of preQ 1 , and the 3′ A-rich tail of the aptamer domain is an integral part of the pseudoknot. In the absence of preQ 1 , the A-rich tail forms part of the antiterminator. These structural studies provide insight into riboswitch transcriptional control of preQ 1 biosynthesis.

RNA molecules are crucial in different levels of cellular function, ranging from translation and regulating genes to coding for proteins. Additionally, nucleic acids (RNA and DNA molecules) are designed for novel applications in biotechnology. Understanding the structure of a molecule is important in inferring its function, and computational methods for structure prediction have captured the interest of many researchers.Some functions of RNA molecules in cells, such as gene regulation, result from the binding of one RNA molecule to another, so-called target RNA molecule. This has led to recent interest in prediction of the secondary structure formed from interacting molecules. In this paper, we provide a brief overview of methods, applications, and challenges in computational prediction of nucleic acid secondary structure, both for single strands and for interacting strands.

In this paper we enumerate k-noncrossing RNA pseudoknot structures with given minimum stack-length, @s. One main result of the paper is the asymptotic formula for their number: T"k","@s(n)~n^-^(^(^k^-^1^)^^^2^+^(^k^-^1^)^/^2^)@c"k","@s^-^n, where @c"k","@s is explicitly known. Our results show that the number of k-noncrossing structures without isolated base pairs is significantly smaller than the number of all k-noncrossing structures. In particular we prove that, for large n, the number of 3- and 4-noncrossing RNA structures with stack-length >=2 is given by 311.24704!n(n-1)...(n-4)2.5881^n and 1.217@?10^7n^-^2^1^23.0382^n, respectively. Our results are of importance for prediction algorithms and provide evidence for the existence of neutral networks of RNA pseudoknot structures.

Telomerase is a ribonucleoprotein complex that replicates the 3’ ends of linear chromosomes by successive additions of telomere repeat DNA. The telomerase holoenzyme contains two essential components for catalysis, a telomerase reverse transcriptase (TERT) and telomerase RNA (TER). The TER includes a template for telomere repeat synthesis as well as other domains required for function. We report the solution structure of the wild type minimal conserved human TER pseudoknot (PKWT) refined with an extensive set of RDCs, and a detailed analysis of the effect of the bulge U177 on pseudoknot structure, dynamics analyzed by RDC and 13C relaxation measurements, and base pair stability. The overall structure of PKWT is highly similar to the previously reported ΔU177 pseudoknot (PKDU) that has a deletion of a conserved bulge U important for catalytic activity. For direct comparison to PKWT, the structure of PKDU was re-refined with a comparable set of RDCs. Both pseudoknots contain a catalytically essential triple helix at the junction of the two stems, including two stem 1-loop 2 minor groove triples, a junction loop 1-loop 2 Hoogsteen base pair, and stem 2-loop 1 major groove U•A-U Watson-Crick-Hoogsteen triples located directly above the bulge U177. However there are significant differences in stabilities of base pairs near the bulge and the dynamics of some nucleotides. The stability of the base pairs in stem 2 surrounding the bulge U177 is greatly decreased, with the result that the Watson-Crick pairs in the triple helix begin to unfold before the Hoogsteen pairs, which may affect telomerase assembly and activity. The bulge U is positioned in the minor groove on the opposite face from the triple helical interactions, and sterically blocks the A176 2’OH which has recently been proposed to play a role in catalysis. The bulge U may serve as a hinge to provide backbone flexibility in this region.

The comparative sequence analysis is the most reliable method for RNA secondary structure prediction, and many algorithms based on it have been developed in last several decades. This paper considers RNA structure prediction as a 2-classes classification problem: given a sequence alignment, to decide whether or not two columns of alignment form a base pair. We employed Support Vector Machine (SVM) to predict potential paired sites, and selected co-variation information, thermodynamic information and the fraction of complementary bases as feature vectors. Considering the effect of sequence similarity upon co-variation score, we introduced a similarity weight factor, which could adjust the contribution of co-variation and thermodynamic information toward prediction according to sequence similarity. The test on 49 Rfam-seed alignments showed the effectiveness of our method, and the accuracy was better than many similar algorithms. Furthermore, this method could predict simple pseudoknot.

Nearly all members of a widespread family of bacterial transposable elements related to insertion sequence 3 (IS3), therefore called the IS3 family, very likely use programmed -1 ribosomal frameshifting to produce their transposase, a protein required for mobility. Comparative analysis of the potential frameshift signals in this family suggested that most of the insertion sequences from the IS51 group contain in their mRNA an elaborate pseudoknot that could act as a recoding stimulator. It results from a specific intramolecular interaction between an apical loop and an internal loop from two stem-loop structures. Directed mutagenesis, chemical probing, and gel mobility assays of the frameshift region of one element from the IS51 group, IS3411, provided clear evidences of the existence of the predicted structure. Modeling was used to generate a three-dimensional molecular representation of the apical loop-internal loop complex. We could demonstrate that mutations affecting the stability of the structure reduce both frameshifting and transposition, thus establishing the biological importance of this new type of RNA structure for the control of transposition level.

International audience; The translation of genetic information into proteins is a fundamental process of life. Stepwise addition of amino acids to the growing polypeptide chain requires the coordinated movement of mRNA and tRNAs through the ribosome, a process known as translocation. Here, we review current understanding of the kinetics and mechanics of translocation, with particular emphasis on the structure of a functional mammalian ribosome stalled during translocation by an mRNA pseudoknot. In the context of a pseudoknot-stalled complex, the translocase EF-2 is seen to compress a hybrid-state tRNA into a strained conformation. We propose that this strain energy helps overcome the kinetic barrier to translocation and drives tRNA into the P-site, with EF-2 biasing this relaxation in one direction. The tRNA can thus be considered a molecular spring and EF-2 a Brownian ratchet in a "spring-and-ratchet" system within the translocation process.

Graphical abstract An RNA pseudoknot is one of the essential components of the viral −1 ribosomal frameshifting (−1RF) system. The control of pseudoknot stability is important in regulating the efficiency of −1RF. We conducted the virtual screening of a chemical database targeting the X-ray crystal structure of biotin-bound RNA pseudoknot to identify ligands that may regulate an −1RF system containing biotin-aptamer as an RNA pseudoknot component. After docking screening of about 80,000 compounds, 58 high-ranked hits were selected and their activities were examined by in vitro and cell-based −1 frameshifting assays. Six compounds increased the efficiency of −1 frameshifting, and these are novel small molecule compounds that regulate the −1RF., In many viruses, −1 ribosomal frameshifting (−1RF) regulates synthesis of proteins and is crucial for virus production. An RNA pseudoknot is one of the essential components of the viral −1RF system. Thermodynamic or kinetic control of pseudoknot folding may be important in regulating the efficiency of −1RF. Thus, small molecules that interact with viral RNA pseudoknots may disrupt the −1RF system and show antiviral activity. In this study, we conducted virtual screening of a chemical database targeting the X-ray crystal structure of RNA pseudoknot complexed with biotin to identify ligands that may regulate an −1RF system containing biotin-aptamer as an RNA pseudoknot component. After docking screening of about 80,000 compounds, 58 high-ranked hits were selected and their activities were examined by in vitro and cell-based −1 frameshifting assays. Six compounds increased the efficiency of −1 frameshifting, and these are novel small molecule compounds that regulate the −1RF.

Simian Immunodeficiency Virus (SIV), like its human homologues (HIV-1, HIV-2), requires a −1 translational frameshift event to properly synthesize all of the proteins required for viral replication. The frameshift mechanism is dependent upon a seven nucleotide slippery sequence and a downstream RNA structure. In SIV, the downstream RNA structure has been proposed to be either a stem-loop or a pseudoknot. Here, we report the functional, structural and thermodynamic characterization of the SIV frameshift site RNA. Translational frameshift assays indicate that a stem-loop structure is sufficient to promote efficient frameshifting in vitro. NMR and thermodynamic studies of SIV RNA constructs of varying length further support the absence of any pseudoknot interaction and indicate the presence of a stable stem-loop structure. We determined the structure of the SIV frameshift-inducing RNA by NMR. The structure reveals a highly ordered 12 nucleotide loop containing a sheared G-A pair, cross-strand adenine stacking, two G-C base-pairs, and a novel CCC triloop turn. The loop structure and its high thermostability preclude pseudoknot formation. Sequence conservation and modeling studies suggest that HIV-2 RNA forms the same structure. We conclude that, like the main sub-groups of HIV-1, SIV and HIV-2 utilize stable stem-loop structures to function as a thermodynamic barrier to translation, thereby inducing ribosomal pausing and frameshifting.

Nucleotide sequences of a broad range of Peach Latent Mosaic Viroid (PLMVd) variants were determined. The variants were isolated from peach, pear, and almond tree samples collected in Tunisia. Sequence analysis confirmed the high variability of PLMVd, as no less than 119 new variants were identified. Variations included new polymorphic positions, insertions of 11 to 14 nucleotides, and new mutations within the hammerhead self-cleavage motifs. We provide the first covariation-based evidence for certain stems within the proposed secondary structure. Our covariation analysis also strengthens the view that a pseudoknot closes the replication domain. On the basis of phylogenetic tree studies and informative positions, PLMVd variants are proposed to cluster into groups and subgroups likely to have resulted from recombination events. PLMVd thus emerges as a suitable viroid for retracing the evolution of an RNA genome.

tmRNA rescues stalled ribosomes in eubacteria by forcing the ribosome to abandon its mRNA template and resume translation with tmRNA itself as a template. Pseudoknot 1 (pk1), immediately upstream of this coding region in tmRNA, is a structural element that is considered essential for tmRNA function based on the analysis of pk1 mutants in vitro. pk1 binds near the ribosomal decoding site and may make base-specific contacts with tmRNA ligands. To study pk1 structure and function in vivo, we have developed a genetic selection that ties the life of Escherichia coli cells to tmRNA activity. Mutation of pk1 at 20% per base and selection for tmRNA activity yielded sequences that retain the same pseudoknot fold. In contrast, selection of active mutants from 10(6) completely random sequences identified hairpin structures that functionally replace pk1. Rational design of a hairpin with increased stability using an unrelated sequence yielded a tmRNA mutant with nearly wild-type activity. We conclude that the role of pk1 in tmRNA function is purely structural and that it can be replaced with a variety of hairpin structures. Our results demonstrate that in the study of functional RNAs, the inactivity of a mutant designed to destroy a given structure should not be interpreted as proof that the structure is necessary for RNA function. Such mutations may only destabilize a global fold that could be formed equally well by an entirely different, stable structure.

A bioinformatics approach to finding new cases of -1 frameshifting in the expression of human genes revealed a classical retrovirus-like heptanucleotide shift site followed by a potential structural stimulator in the paraneoplastic antigen Ma3 and Ma5 genes. Analysis of the sequence 3' of the shift site demonstrated that an RNA pseudoknot in Ma3 is important for promoting efficient -1 frame-shifting. Ma3 is a member of a family of six genes in humans whose protein products contain homology to retroviral Gag proteins. The -1 frameshift site and pseudoknot structure are conserved in other mammals, but there are some sequence differences. Although the functions of the Ma genes are unknown, the serious neurological effects of ectopic expression in tumor cells indicate their importance in the brain.

In this paper, we present the results from a comprehensive study of nanosecond-scale implicit and explicit solvent molecular dynamics simulations of the wild-type telomerase RNA hairpin. The effects of various mutations on telomerase RNA dynamics are also investigated. Overall, we found that the human telomerase hairpin is a very flexible molecule. In particular, periodically the molecule exhibits dramatic structural fluctuations represented by the opening and closing of a non-canonical base-pair region. These structural deviations correspond to significant disruptions of the direct hydrogen bonding network in the helix, widening of the major groove of the hairpin structure, and causing several U and C nucleotides to protrude into the major groove from the helix permitting them to hydrogen bond with, for example, the P3 domain of the telomerase RNA. We suggest that these structural fluctuations expose a nucleation point for pseudoknot formation. We also found that mutations in the pentaloop and non-canonical region stabilize the hairpin. Moreover, our results show that the hairpin with dyskeratosis congenita mutations is more stable and less flexible than the wild-type hairpin due to base stacking in the pentaloop. The results from our molecular dynamics simulations are in agreement with experimental observations. In addition, they suggest a possible mechanism for pseudoknot formation based on the dynamics of the hairpin structure and also may explain the mutational aspects of dyskeratosis congenita.

Human telomerase contains a 451 nt RNA (hTR) and several proteins, including a specialized reverse transcriptase (hTERT). The 5' half of hTR comprises the pseudoknot (core) domain, which includes the RNA template for telomere synthesis and a highly conserved pseudoknot that is required for telomerase activity. The solution structure of this essential pseudoknot, presented here, reveals an extended triple helix surrounding the helical junction. The network of tertiary interactions explains the phylogenetic sequence conservation and existing human and mouse TR functional studies as well as mutations linked to disease. Thermodynamic stability, dimerization potential, and telomerase activity of mutant RNAs that alter the tertiary contacts were investigated. Telomerase activity is strongly correlated with tertiary structure stability, whereas there is no correlation with dimerization potential of the pseudoknot. These studies reveal that a conserved pseudoknot tertiary structure is required for telomerase activity.

Here, we demonstrate that a series of naphthyridine derivatives, naphthyridine carbamate tetramer (NCTn), can bind to the RNA CGG/CGG triad comprised of two single-stranded regions of a hairpin loop and a tail. Complete suppression of the binding by a single mutation indicated simultaneous binding of NCTn between hairpin loop and single stranded tail, leading to the formation of NCTn-induced pseudoknot.

Bacteria contain transfer-messenger RNA (tmRNA), a molecule that during trans-translation tags incompletely translated proteins with a small peptide to signal the proteolytic destruction of defective polypeptides. TmRNA is composed of tRNA- and mRNA-like domains connected by several pseudoknots. Using truncated ribosomal protein L27 as a reporter for tagging in vitro and in vivo, we have developed exceptionally sensitive assays to study the role of Escherichia coli tmRNA in trans-translation. Site-directed mutagenesis experiments showed that pseudoknot 2 and the abutting helix 5 were particularly important for the binding of ribosomal protein S1 to tmRNA. Pseudoknot 4 not only facilitated tmRNA maturation but also promoted tagging. In addition, the three pseudoknots (pk2 to pk4) were shown to play a significant role in the proper folding of the tRNA-like domain. Protein SmpB enhanced tmRNA processing, suggesting a new role for SmpB in trans-translation. Taken together, these results provide unanticipated insights into the functions of the pseudoknots and protein SmpB during tmRNA folding, maturation, and protein synthesis.