RNA pseudoknots: structure, detection, and prediction

A Hitchhiker's guide to RNA-RNA structure and interaction prediction tools

RNA biology has risen to prominence after a remarkable discovery of diverse functions of noncoding RNA (ncRNA). Most untranslated transcripts often exert their regulatory functions into RNA-RNA complexes via base pairing with complementary sequences in other RNAs. An interplay between RNAs is essential, as it possesses various functional roles in human cells, including genetic translation, RNA splicing, editing, ribosomal RNA maturation, RNA degradation and the regulation of metabolic pathways/riboswitches. Moreover, the pervasive transcription of the human genome allows for the discovery of novel genomic functions via RNA interactome investigation. The advancement of experimental procedures has resulted in an explosion of documented data, necessitating the development of efficient and precise computational tools and algorithms. This review provides an extensive update on RNA-RNA interaction (RRI) analysis via thermodynamic- and comparative-based RNA secondary structure prediction (RSP) and RNA-RNA interaction prediction (RIP) tools and their general functions. We also highlighted the current knowledge of RRIs and the limitations of RNA interactome mapping via experimental data. Then, the gap between RSP and RIP, the importance of RNA homologues, the relationship between pseudoknots, and RNA folding thermodynamics are discussed. It is hoped that these emerging prediction tools will deepen the understanding of RNA-associated interactions in human diseases and hasten treatment processes.

Genomics and pathogenesis of the avian coronavirus infectious bronchitis virus

Infectious bronchitis virus (IBV) is a member of the family Coronaviridae, together with viruses such as SARS-CoV, MERS-CoV and SARS-CoV-2 (the causative agent of the COVID-19 global pandemic). In this family of viruses, interspecies transmission has been reported, so understanding their pathobiology could lead to a better understanding of the emergence of new serotypes. IBV possesses a single-stranded, non-segmented RNA genome about 27.6 kb in length that encodes several non-structural and structural proteins. Most functions of these proteins have been confirmed in IBV, but some other proposed functions have been based on research conducted on other members of the family Coronaviridae. IBV has variable tissue tropism depending on the strain, and can affect the respiratory, reproductive, or urinary tracts; however, IBV can also replicate in other organs. Additionally, the pathogenicity of IBV is also variable, with some strains causing only mild clinical signs, while infection with others results in high mortality rates in chickens. This paper extensively and comprehensibly reviews general aspects of coronaviruses and, more specifically, IBV, with emphasis on protein functions and pathogenesis. The pathogenicity of the Australian strains of IBV is also reviewed, describing the variability between the different groups of strains, from the classical to the novel and recombinant strains. Reverse genetic systems, cloning and cell culture growth techniques applicable to IBV are also reviewed.

LTPConstraint: a transfer learning based end-to-end method for RNA secondary structure prediction

BACKGROUND

RNA secondary structure is very important for deciphering cell's activity and disease occurrence. The first method which was used by the academics to predict this structure is biological experiment, But this method is too expensive, causing the promotion to be affected. Then, computing methods emerged, which has good efficiency and low cost. However, the accuracy of computing methods are not satisfactory. Many machine learning methods have also been applied to this area, but the accuracy has not improved significantly. Deep learning has matured and achieves great success in many areas such as computer vision and natural language processing. It uses neural network which is a kind of structure that has good functionality and versatility, but its effect is highly correlated with the quantity and quality of the data. At present, there is no model with high accuracy, low data dependence and high convenience in predicting RNA secondary structure.

RESULTS

This paper designs a neural network called LTPConstraint to predict RNA secondary structure. The network is based on many network structure such as Bidirectional LSTM, Transformer and generator. It also uses transfer learning to train modelso that the data dependence can be reduced.

CONCLUSIONS

LTPConstraint has achieved high accuracy in RNA secondary structure prediction. Compared with the previous methods, the accuracy improves obviously both in predicting the structure with pseudoknot and the structure without pseudoknot. At the same time, LTPConstraint is easy to operate and can achieve result very quickly.

iTRAQ-based analysis of leaf proteome identifies important proteins in secondary metabolite biosynthesis and defence pathways crucial to cross-protection against TMV

Cross-protection is a phenomenon in which infection with a mild virus strain protects host plants against subsequent infection with a closely related severe virus strain. This study showed that a mild strain mutant virus, Tobacco mosaic virus (TMV)-43A could cross protect Nicotiana benthamiana plants against wild-type TMV. Furthermore, we investigated the host responses at the proteome level to identify important host proteins involved in cross-protection. We used the isobaric tags for relative and absolute quantification (iTRAQ) technique to analyze the proteome profiles of TMV, TMV-43A and cross-protected plants at different time-points. Our results showed that TMV-43A can cross-protect N. benthamiana plants from TMV. In cross-protected plants, photosynthetic activities were augmented, as supported by the increased accumulation of 1-deoxy-D-xylulose-5-phosphate reductoisomerase (DXR) and geranylgeranyl diphosphate synthase (GGPS) enzymes, which are crucial for chlorophyll biosynthesis. The increased abundance of ROS scavenging enzymes like thioredoxins and L-ascorbate peroxidase would prevent oxidative damage in cross-protected plants. Interestingly, the abundance of defence-related proteins (14-3-3 and NbSGT1) decreased, along with a reduction in virus accumulation during cross-protection. In conclusion, we have identified several important host proteins that are crucial in cross-protection to counter TMV infection in N. benthamiana plants. BIOLOGICAL SIGNIFICANCE: TMV is the most studied model for host-virus interaction in plants. It can infect wide varieties of plant species, causing significant economic losses. Cross protection is one of the methods to combat virus infection. A few cross-protection mechanisms have been proposed, including replicase/coat protein-mediated resistance, RNA silencing, and exclusion/spatial separation between virus strains. However, knowledge on host responses at the proteome level during cross protection is limited. To address this knowledge gap, we have leveraged on a global proteomics analysis approach to study cross protection. We discovered that TMV-43A (protector) protects N. benthamiana plants from TMV (challenger) infection through multiple host pathways: secondary metabolite biosynthesis, photosynthesis, defence, carbon metabolism, protein translation and processing and amino acid biosynthesis. In the secondary metabolite biosynthesis pathway, enzymes 1-deoxy-D-xylulose-5-phosphate reductoisomerase (DXR) and geranylgeranyl diphosphate synthase (GGPS) play crucial roles in chlorophyll biosynthesis during cross protection. In addition, accumulation of ROS scavenging enzymes was also found in cross-protected plants, providing rescues from excessive oxidative damage. Reduced abundance of plant defence proteins is correlated to reduced virus accumulation in host plants. These findings have increased our knowledge in host responses during cross-protection.

Ultrastable pRNA hexameric ring gearing hexameric phi29 DNA-packaging motor by revolving without rotating and coiling

Biomotors have previously been classified into two categories: linear and rotational motors. It has long been popularly believed that viral DNA packaging motors are rotation motors. We have recently found that the DNA-packaging motor of bacteriophage phi29 uses a third mechanism: revolution without rotation. phi29 motor consists of three-coaxial rings of hexameric RNA, a hexameric ATPase, and a dodecameric channel. The motor uses six ATP to revolve one helical turn of dsDNA around the hexameric ring of ATPase gp16. Each dodecameric segment tilts at a 30°-angle and runs anti-parallel to the dsDNA helix to facilitate translation in one direction. The negatively charged phosphate backbone interacts with four positively charged lysine rings, resulting in four steps of transition. This review will discuss how the novel pRNA meets motor requirements for translocation concerning structure, stoichiometry, and thermostability; how pRNA studies have led to the generation of the concept of RNA nanotechnology; and how pRNA is fabricated into nanoparticles to deliver siRNA, miRNA, and ribozymes to cancer and virus-infected cells.

Uniqueness, advantages, challenges, solutions, and perspectives in therapeutics applying RNA nanotechnology

The field of RNA nanotechnology is rapidly emerging. RNA can be manipulated with the simplicity characteristic of DNA to produce nanoparticles with a diversity of quaternary structures by self-assembly. Additionally RNA is tremendously versatile in its function and some RNA molecules display catalytic activities much like proteins. Thus, RNA has the advantage of both worlds. However, the instability of RNA has made many scientists flinch away from RNA nanotechnology. Other concerns that have deterred the progress of RNA therapeutics include the induction of interferons, stimulation of cytokines, and activation of other immune systems, as well as short pharmacokinetic profiles in vivo. This review will provide some solutions and perspectives on the chemical and thermodynamic stability, in vivo half-life and biodistribution, yield and production cost, in vivo toxicity and side effect, specific delivery and targeting, as well as endosomal trapping and escape.

Recent advances in the rational design of silica-based nanoparticles for gene therapy

Gene therapy has attracted much attention in modern society and provides a promising approach for treating genetic disorders, diseases and cancers. Safe and effective vectors are vital tools to deliver genetic molecules to cells. This review summarizes recent advances in the rational design of silica-based nanoparticles and their applications in gene therapy. An overview of different types of genetic agents available for gene therapy is provided. The engineering of various silica nanoparticles is described, which can be used as versatile complexation tools for genetic agents and advanced gene therapy. Several challenges are raised and future research directions in the area of gene therapy using silica-based nanoparticles are proposed.

Uniqueness, advantages, challenges, solutions, and perspectives in therapeutics applying RNA nanotechnology

The field of RNA nanotechnology is rapidly emerging. RNA can be manipulated with the simplicity characteristic of DNA to produce nanoparticles with a diversity of quaternary structures by self-assembly. Additionally RNA is tremendously versatile in its function and some RNA molecules display catalytic activities much like proteins. Thus, RNA has the advantage of both worlds. However, the instability of RNA has made many scientists flinch away from RNA nanotechnology. Other concerns that have deterred the progress of RNA therapeutics include the induction of interferons, stimulation of cytokines, and activation of other immune systems, as well as short pharmacokinetic profiles in vivo. This review will provide some solutions and perspectives on the chemical and thermodynamic stability, in vivo half-life and biodistribution, yield and production cost, in vivo toxicity and side effect, specific delivery and targeting, as well as endosomal trapping and escape.

Ultrastable synergistic tetravalent RNA nanoparticles for targeting to cancers

One of the advantages of nanotechnology is the feasibility to construct therapeutic particles carrying multiple therapeutics with defined structure and stoichiometry. The field of RNA nanotechnology is emerging. However, controlled assembly of stable RNA nanoparticles with multiple functionalities which retain their original role is challenging due to refolding after fusion. Herein, we report the construction of thermodynamically stable X-shaped RNA nanoparticles to carry four therapeutic RNA motifs by self-assembly of reengineered small RNA fragments. We proved that each arm of the four helices in the X-motif can harbor one siRNA, ribozyme, or aptamer without affecting the folding of the central pRNA-X core, and each daughter RNA molecule within the nanoparticle folds into their respective authentic structures and retains their biological and structural function independently. Gene silencing effects were progressively enhanced as the number of the siRNA in each pRNA-X nanoparticles gradually increased from one to two, three, and four. More importantly, systemic injection of ligand-containing nanoparticles into the tail-vein of mice revealed that the RNA nanoparticles remained intact and strongly bound to cancers without entering the liver, lung or any other organs or tissues, while remaining in cancer tissue for more than 8 h.

A boost for the emerging field of RNA nanotechnology

This Nano Focus article highlights recent advances in RNA nanotechnology as presented at the First International Conference of RNA Nanotechnology and Therapeutics, which took place in Cleveland, OH, USA (October 23-25, 2010) ( http://www.eng.uc.edu/nanomedicine/RNA2010/ ), chaired by Peixuan Guo and co-chaired by David Rueda and Scott Tenenbaum. The conference was the first of its kind to bring together more than 30 invited speakers in the frontier of RNA nanotechnology from France, Sweden, South Korea, China, and throughout the United States to discuss RNA nanotechnology and its applications. It provided a platform for researchers from academia, government, and the pharmaceutical industry to share existing knowledge, vision, technology, and challenges in the field and promoted collaborations among researchers interested in advancing this emerging scientific discipline. The meeting covered a range of topics, including biophysical and single-molecule approaches for characterization of RNA nanostructures; structure studies on RNA nanoparticles by chemical or biochemical approaches, computation, prediction, and modeling of RNA nanoparticle structures; methods for the assembly of RNA nanoparticles; chemistry for RNA synthesis, conjugation, and labeling; and application of RNA nanoparticles in therapeutics. A special invited talk on the well-established principles of DNA nanotechnology was arranged to provide models for RNA nanotechnology. An Administrator from National Institutes of Health (NIH) National Cancer Institute (NCI) Alliance for Nanotechnology in Cancer discussed the current nanocancer research directions and future funding opportunities at NCI. As indicated by the feedback received from the invited speakers and the meeting participants, this meeting was extremely successful, exciting, and informative, covering many groundbreaking findings, pioneering ideas, and novel discoveries.

Engineering RNA for targeted siRNA delivery and medical application

RNA engineering for nanotechnology and medical applications is an exciting emerging research field. RNA has intrinsically defined features on the nanometre scale and is a particularly interesting candidate for such applications due to its amazing diversity, flexibility and versatility in structure and function. Specifically, the current use of siRNA to silence target genes involved in disease has generated much excitement in the scientific community. The intrinsic ability to sequence-specifically downregulate gene expression in a temporally- and spatially controlled fashion has led to heightened interest and rapid development of siRNA-based therapeutics. Although methods for gene silencing have been achieved with high efficacy and specificity in vitro, the effective delivery of nucleic acids to specific cells in vivo has been a hurdle for RNA therapeutics. This article covers different RNA-based approaches for diagnosis, prevention and treatment of human disease, with a focus on the latest developments of non-viral carriers of siRNA for delivery in vivo. The applications and challenges of siRNA therapy, as well as potential solutions to these problems, the approaches for using phi29 pRNA-based vectors as polyvalent vehicles for specific delivery of siRNA, ribozymes, drugs or other therapeutic agents to specific cells for therapy will also be addressed.

On quantitative effects of RNA shape abstraction

Over the last few decades, much effort has been taken to develop approaches for identifying good predictions of RNA secondary structure. This is due to the fact that most computational prediction methods based on free energy minimization compute a number of suboptimal foldings and we have to identify the native folding among all these possible secondary structures. Using the abstract shapes approach as introduced by Giegerich et al. (Nucleic Acids Res 32(16):4843-4851, 2004), each class of similar secondary structures is represented by one shape and the native structures can be found among the top shape representatives. In this article, we derive some interesting results answering enumeration problems for abstract shapes and secondary structures of RNA. We compute precise asymptotics for the number of different shape representations of size n and for the number of different shapes showing up when abstracting from secondary structures of size n under a combinatorial point of view. A more realistic model taking primary structures into account remains an open challenge. We give some arguments why the present techniques cannot be applied in this case.

Structural and functional analyses of stem-loop 1 of the Sindbis virus genome

Alphavirus genome function is controlled by elements at both the 5' and 3' ends. The 5' 220 nt of the Sindbis virus genome is predicted to consist of four stem-loop structures the first of which has been demonstrated to be required for efficient minus-strand RNA synthesis. To understand the role of the structure of the first stem-loop (SL1) in regulating genome function, we performed enzymatic and chemical probing analyses. There were significant differences between the computer-predicted structures and our experimental data. In the 5' terminus, two loop regions appear to be interacting in a complex and interdependent fashion with non-Watson-Crick interactions involving multiple adenosine residues playing a critical role in determining the overall structure. Some of the mutations that disrupted these interactions had significant affects, both positive and negative, on minus-strand synthesis, and translational efficiency was generally increased. In the context of full-length virus, these structural changes resulted in reduced virus growth kinetics particularly in mosquito cells suggesting host-specific effects of mutations in this region of the viral genome. Possible SL1 structures based on our experimental data are discussed.

Predicting RNA pseudoknot folding thermodynamics

Based on the experimentally determined atomic coordinates for RNA helices and the self-avoiding walks of the P (phosphate) and C₄ (carbon) atoms in the diamond lattice for the polynucleotide loop conformations, we derive a set of conformational entropy parameters for RNA pseudoknots. Based on the entropy parameters, we develop a folding thermodynamics model that enables us to compute the sequence-specific RNA pseudoknot folding free energy landscape and thermodynamics. The model is validated through extensive experimental tests both for the native structures and for the folding thermodynamics. The model predicts strong sequence-dependent helix-loop competitions in the pseudoknot stability and the resultant conformational switches between different hairpin and pseudoknot structures. For instance, for the pseudoknot domain of human telomerase RNA, a native-like and a misfolded hairpin intermediates are found to coexist on the (equilibrium) folding pathways, and the interplay between the stabilities of these intermediates causes the conformational switch that may underlie a human telomerase disease.

Crystal structure of a luteoviral RNA pseudoknot and model for a minimal ribosomal frameshifting motif

To understand the role of structural elements of RNA pseudoknots in controlling the extent of -1-type ribosomal frameshifting, we determined the crystal structure of a high-efficiency frameshifting mutant of the pseudoknot from potato leaf roll virus (PLRV). Correlations of the structure with available in vitro frameshifting data for PLRV pseudoknot mutants implicate sequence and length of a stem-loop linker as modulators of frameshifting efficiency. Although the sequences and overall structures of the RNA pseudoknots from PLRV and beet western yellow virus (BWYV) are similar, nucleotide deletions in the linker and adjacent minor groove loop abolish frameshifting only with the latter. Conversely, mutant PLRV pseudoknots with up to four nucleotides deleted in this region exhibit nearly wild-type frameshifting efficiencies. The crystal structure helps rationalize the different tolerances for deletions in the PLRV and BWYV RNAs, and we have used it to build a three-dimensional model of the PRLV pseudoknot with a four-nucleotide deletion. The resulting structure defines a minimal RNA pseudoknot motif composed of 22 nucleotides capable of stimulating -1-type ribosomal frameshifts.

A heuristic approach for detecting RNA H-type pseudoknots

MOTIVATION

RNA H-type pseudoknots are ubiquitous pseudoknots that are found in almost all classes of RNA and thought to play very important roles in a variety of biological processes. Detection of these RNA H-type pseudoknots can improve our understanding of RNA structures and their associated functions. However, the currently existing programs for detecting such RNA H-type pseudoknots are still time consuming and sometimes even ineffective. Therefore, efficient and effective tools for detecting the RNA H-type pseudoknots are needed.

RESULTS

In this paper, we have adopted a heuristic approach to develop a novel tool, called HPknotter, for efficiently and accurately detecting H-type pseudoknots in an RNA sequence. In addition, we have demonstrated the applicability and effectiveness of HPknotter by testing on some sequences with known H-type pseudoknots. Our approach can be easily extended and applied to other classes of more general pseudoknots.

AVAILABILITY

The web server of our HPknotter is available for online analysis at http://bioalgorithm.life.nctu.edu.tw/HPKNOTTER/ CONTACT: cllu@mail.nctu.edu.tw, chiu@cc.nctu.edu.tw

In vitro characterization of a base pairing interaction between the primer binding site and the minimal packaging signal of avian leukosis virus genomic RNA

The 5' leader region of avian sarcoma-leukosis viruses (ASLVs) folds into a series of RNA secondary structures which are involved in key steps in the viral replication cycle such as reverse transcription, dimerization and packaging of genomic RNA. The O3 stem and three stem-loops (O3SLa, O3SLb and O3SLc) form the minimal packaging signal that is located downstream of the primer binding site (PBS). The U5-PBS region contributes to packaging via a mechanism that remains unknown. In this in vitro study, we have investigated the possibility of interactions between the R-U5-PBS region and the minimal packaging signal using chemical and enzymatic probing, antisense oligonucleotides and site-directed mutagenesis. We have identified a base pairing interaction between the PBS sequence and the terminal loop of O3SLa. It was found that the PBS/O3SLa interaction was intramolecular since it occurred not only in dimeric RNA but also in monomeric RNA. This interaction probably corresponds to a pseudoknot interaction. The PBS/O3SLa interaction may be formed in vivo since the sequences are highly conserved in ASLV strains. The PBS/O3SLa interaction may regulate the processes of primer tRNA annealing, packaging and initiation of Gag translation through its involvement in leader tertiary structure. Interestingly, we found that in other retroviruses the PBS sequence can also base pair with a terminal loop of the stem-loops involved in RNA packaging.

A point mutation in bioactive RNA results in the failure of mutant heart correction in Mexican axolotls

Ambystoma mexicanum is an intriguing animal model for studying heart development because it carries a mutation in gene c. Hearts of homozygous recessive (c/c) mutant embryos do not contain organized myofibrils and fail to beat. The defect can be corrected by organ-culturing the mutant heart in the presence of RNA from anterior endoderm or endoderm/mesoderm-conditioned medium. By screening a cDNA library made of total conditioned medium RNA from normal axolotl embryonic endoderm, we isolated a single clone (MIR), the synthetic RNA from which corrects the mutant heart defect by promoting myofibrillogenesis and thus was named MIR (myofibrillogenesis inducing RNA). In the present study, we have examined MIR gene expression in mutant axolotl hearts at early pre-heart-beat developmental stages and found its quantitative expression, as detected by RT-PCR, to be the same as in normal hearts. However, careful analysis of sequence data revealed a G-->U point mutation in the mutant MIR RNA. Further computational analyses, using GENEBEE software to compare normal and mutant MIR RNAs show a significant alteration in RNA secondary structure of the point-mutated MIR RNA. The results from bioassay and confocal microscopy immunofluorescent studies demonstrate that, unlike MIR RNA derived from normal embryos, the mutated MIR RNA does not promote myofibrillogenesis in mutant embryonic hearts and fails to rescue/correct the mutant heart defect.

Eukaryotic elongation factor 1A interacts with the upstream pseudoknot domain in the 3' untranslated region of tobacco mosaic virus RNA

The genomic RNA of tobacco mosaic virus (TMV), like that of other positive-strand RNA viruses, acts as a template for both translation and replication. The highly structured 3' untranslated region (UTR) of TMV RNAs plays an important role in both processes; it is not polyadenylated but ends with a tRNA-like structure (TLS) preceded by a conserved upstream pseudoknot domain (UPD). The TLS of tobamoviral RNAs can be specifically aminoacylated and, in this state, can interact with eukaryotic elongation factor 1A (eEF1A)/GTP with high affinity. Using a UV cross-linking assay, we detected another specific binding site for eEF1A/GTP, within the UPDs of TMV and crucifer-infecting tobamovirus (crTMV), that does not require aminoacylation. A mutational analysis revealed that UPD pseudoknot conformation and some conserved primary sequence elements are required for this interaction. Its possible role in the regulation of tobamovirus gene expression and replication is discussed.

Three-dimensional interaction of Phi29 pRNA dimer probed by chemical modification interference, cryo-AFM, and cross-linking

Six pRNAs (p for packaging) of bacterial virus phi29 form a hexamer complex that is an essential component of the viral DNA translocating motor. Dimers, the building block of pRNA hexamer, assemble in the order of dimer --> tetramer --> hexamer. The two-dimensional structure of the pRNA monomer has been investigated extensively; however, the three-dimensional structure concerning the distance constraints of the three stems and loops are unknown. In this report, we probed the three-dimensional structure of pRNA monomer and dimer by photo affinity cross-linking with azidophenacyl. Bases 75-81 of the left stem were found to be oriented toward the head loop and proximate to bases 26-31 in a parallel orientation. Chemical modification interference indicates the involvement of bases 45-71 and 82-91 in dimer formation. Dimer was formed via hand-in-hand contact, a novel RNA dimerization that in some aspects is similar to the kissing loops of the human immunodeficiency virus. The covalently linked dimers were found to be biologically active. Both the native dimer and the covalently linked dimer were found by cryo-atomic force microscopy to be similar in global conformation and size.

A dimer as a building block in assembling RNA. A hexamer that gears bacterial virus phi29 DNA-translocating machinery

Six RNA (pRNA) molecules form a hexamer, via hand-in-hand interaction, to gear bacterial virus phi29 DNA translocation machinery. Here we report the pathway and the conditions for the hexamer formation. Stable pRNA dimers and trimers were assembled in solution, isolated from native gels, and separated by sedimentation, providing a model system for the study of RNA dimers and trimers in a protein-free environment. Cryo-atomic force microscopy revealed that monomers displayed a check mark outline, dimers exhibited an elongated shape, and trimers formed a triangle. Dimerization of pRNA was promoted by a variety of cations including spermidine, whereas procapsid binding and DNA packaging required specific divalent cations, including Mg(2+), Ca(2+), and Mn(2+). Both the tandem and fused pRNA dimers with complementary loops designed to form even-numbered rings were active in DNA packaging, whereas those without complementary loops were inactive. We conclude that dimers are the building blocks of the hexamer, and the pathway of building a hexamer is: dimer --> tetramer --> hexamer. The Hill coefficient of 2.5 suggests that there are three binding sites with cooperative binding on the surface of the procapsid. The two interacting loops played a key role in recruiting the incoming dimer, whereas the procapsid served as the foundation for hexamer assembly.

Gill-associated virus of Penaeus monodon prawns: an invertebrate virus with ORF1a and ORF1b genes related to arteri- and coronaviruses

A 20089 nucleotide (nt) sequence was determined for the 5' end of the (+)-ssRNA genome of gill-associated virus (GAV), a yellow head-like virus infecting Penaeus monodon prawns. Clones were generated from a approximately 22 kb dsRNA purified from lymphoid organ total RNA of GAV-infected prawns. The region contains a single gene comprising two long overlapping open reading frames, ORF1a and ORF1b, of 4060 and 2646 amino acids, respectively. The ORFs are structurally related to the ORF1a and ORF1ab polyproteins of coronaviruses and arteriviruses. The 99 nt overlap between ORF1a and ORF1b contains a putative AAAUUUU 'slippery' sequence associated with -1 ribosomal frameshifting. A 131 nt stem-loop with the potential to form a complex pseudoknot resides 3 nt downstream of this sequence. Although different to the G/UUUAAAC frameshift sites and 'H-type' pseudoknots of nidoviruses, in vitro transcription/translation analysis demonstrated that the GAV element also facilitates read-through of the ORF1a/1b junction. As in coronaviruses, GAV ORF1a encodes a 3C-like cysteine protease domain located between two hydrophobic regions. However, its sequence suggests some structural relationship to the chymotrypsin-like serine proteases of arteriviruses. ORF1b encodes homologues of the 'SDD' polymerase, which among (+)-RNA viruses is unique to nidoviruses, as well as metal-ion-binding and helicase domains. The presence of a dsRNA replicative intermediate and ORF1a and ORF1ab polyproteins translated by a-1 frameshift suggests that GAV represents the first invertebrate member of the Order NIDOVIRALES:

Structure, stability and function of RNA pseudoknots involved in stimulating ribosomal frameshifting

Programmed -1 ribosomal frameshifting has become the subject of increasing interest over the last several years, due in part to the ubiquitous nature of this translational recoding mechanism in pathogenic animal and plant viruses. All cis-acting frameshift signals encoded in mRNAs are minimally composed of two functional elements: a heptanucleotide "slippery sequence" conforming to the general form X XXY YYZ, followed by an RNA structural element, usually an H-type RNA pseudoknot, positioned an optimal number of nucleotides (5 to 9) downstream. The slippery sequence itself promotes a low level ( approximately 1 %) of frameshifting; however, downstream pseudoknots stimulate this process significantly, in some cases up to 30 to 50 %. Although the precise molecular mechanism of stimulation of frameshifting remains poorly understood, significant advances have been made in our knowledge of the three-dimensional structures, thermodynamics of folding, and functional determinants of stimulatory RNA pseudoknots derived from the study of several well-characterized frameshift signals. These studies are summarized here and provide new insights into the structural requirements and mechanism of programmed -1 ribosomal frameshifting.

Secondary structure of vertebrate telomerase RNA

Telomerase is a ribonucleoprotein enzyme that maintains telomere length by adding telomeric sequence repeats onto chromosome ends. The essential RNA component of telomerase provides the template for telomeric repeat synthesis. To determine the secondary structure of vertebrate telomerase RNA, 32 new telomerase RNA genes were cloned and sequenced from a variety of vertebrate species including 18 mammals, 2 birds, 1 reptile, 7 amphibians, and 4 fishes. Using phylogenetic comparative analysis, we propose a secondary structure that contains four structural domains conserved in all vertebrates. Ten helical regions of the RNA are universally conserved while other regions vary significantly in length and sequence between different classes of vertebrates. The proposed vertebrate telomerase RNA structure displays a strikingly similar topology to the previously determined ciliate telomerase RNA structure, implying an evolutionary conservation of the global architecture of telomerase RNA.

Genomic structure of three phenotypically different isolates of peach latent mosaic viroid: implications of the existence of constraints limiting the heterogeneity of viroid quasispecies

The peach latent mosaic viroid (PLMVd) is used to study the interactions between a viroid containing hammerhead ribozymes and its natural host, peach. To gain insight into the molecular basis of the phenotypic effects observed upon viroid infection, sequence variants from three PLMVd isolates that differ in symptom expression on the peach indicator GF-305 have been characterized. Analysis of the primary structures of a total of 29 different sequence variants derived from a severe and two latent isolates has revealed a large number of polymorphic positions in the viroid molecule. The variability pattern indicates that preservation of the stability of both hammerhead structures and conservation of a branched secondary structure of the viroid molecule may be factors limiting sequence heterogeneity in PLMVd. Moreover, compensatory mutations in two hairpin loops of the proposed secondary structure, suggesting that a pseudoknot-like interaction may exist between them, have also been observed. Phylogenetic analysis has allowed the allocation of PLMVd molecules into three major groups. This clustering does not strictly correlate with the source isolate from which the variants were obtained, providing insights into the complex mixture of molecules which make up each isolate. Bioassays of individual PLMVd sequence variants on GF-305 peach seedlings have shown that the biological properties of the PLMVd isolates may be correlated with both the complexity of their viroid populations and the presence of specific sequence variants.

The structure of HIV-1 reverse transcriptase complexed with an RNA pseudoknot inhibitor

Small RNA pseudoknots, selected to bind human immunodeficiency virus type 1 (HIV-1) reverse transcriptase tightly, are potent inhibitors of reverse transcriptase. The co-crystal structure of reverse transcriptase complexed with a 33 nucleotide RNA pseudoknot has been determined by fitting the ligand into a high quality, 4-fold averaged 4.8 A resolution electron density map. The RNA is kinked between stems S1 and S2, thereby optimizing its contacts with subunits of the heterodimer. Its binding site extends along the cleft that lies between the polymerase and RNase H active sites, partially overlaps with that observed for duplex DNA and presumably overlaps some portion of the tRNA site. Stem S2 and loop L1 stabilize the 'closed' conformation of the polymerase through extensive electrostatic interactions with several basic residues in helix I of the p66 thumb and in the p66 fingers domain. Presumably, this RNA ligand inhibits reverse transcriptase by binding to a site that partly overlaps the primer-template binding site.

Intracisternal A-type particles express their proteinase in a separate reading frame by translational frameshifting, similar to D-type retroviruses

Intracisternal A-type particles (IAP) are defective endogenous retroviruses that accumulate in the endoplasmic reticulum of rodent cells. IAP genomes share extensive sequence homologies with D-type retroviruses, but were presumed to express the viral proteinase (PR) as part of the gag open reading frame (ORF) while D-type retroviruses express PR in a separate ORF. Here we show that expression of the murine IAP element MIA14 yields three major translation products, corresponding to the Gag, Gag-PR, and Gag-PR-Pol polyproteins. Sequence analysis revealed that MIA14 PR is encoded in its own reading frame, separate from gag and pol. Frameshifting occurred with an efficiency of approximately 25% between the gag and pro ORFs and 35% between pro and pol. The region containing the putative gag-pro frameshift signal consists of a heptanucleotide slippery sequence (A6C) and a stem-loop structure probably forming a pseudoknot. Deletion of this structure element almost completely abolished frameshifting. Insertion of an additional base next to the frameshift signal placed gag and pro in the same ORF and resulted in predominant formation of Gag-PR and Gag-PR-Pol polyproteins which were not processed following in vitro translation. Expression of a similar construct in tissue culture cells, on the other hand, led to efficient intracellular processing of the mutant polyproteins.

An NMR and mutational study of the pseudoknot within the gene 32 mRNA of bacteriophage T2: insights into a family of structurally related RNA pseudoknots

NMR methods were used to investigate a series of mutants of the pseudoknot within the gene 32 messenger RNA of bacteriophage T2, for the purpose of investigating the range of sequences, stem and loop lengths that can form a similar pseudoknot structure. This information is of particular relevance since the T2 pseudoknot has been considered a representative of a large family of RNA pseudoknots related by a common structural motif, previously referred to as 'common pseudoknot motif 1' or CPK1. In the work presented here, a mutated sequence with the potential to form a pseudoknot with a 6 bp stem2 was shown to adopt a pseudoknot structure similar to that of the wild-type sequence. This result is significant in that it demonstrates that pseudoknots with 6 bp in stem2 and a single nucleotide in loop1 are indeed feasible. Mutated sequences with the potential to form pseudoknots with either 5 or 8 bp in stem2 yielded NMR spectra that could not confirm the formation of a pseudoknot structure. Replacing the adenosine nucleotide in loop1 of the wild-type pseudoknot with any one of G, C or U did not significantly alter the pseudoknot structure. Taken together, the results of this study provide support for the existence of a family of similarly structured pseudoknots with two coaxially stacked stems, either 6 or 7 bp in stem2, and a single nucleotide in loop1. This family includes many of the pseudoknots predicted to occur downstream of the frameshift or readthrough sites in a significant number of viral RNAs.

Structural and functional studies of retroviral RNA pseudoknots involved in ribosomal frameshifting: nucleotides at the junction of the two stems are important for efficient ribosomal frameshifting

Ribosomal frameshifting, a translational mechanism used during retroviral replication, involves a directed change in reading frame at a specific site at a defined frequency. Such programmed frameshifting at the mouse mammary tumor virus (MMTV) gag-pro shift site requires two mRNA signals: a heptanucleotide shifty sequence and a pseudoknot structure positioned downstream. Using in vitro translation assays and enzymatic and chemical probes for RNA structure, we have defined features of the pseudoknot that promote efficient frameshifting. Heterologous RNA structures, e.g. a hairpin, a tRNA or a synthetic pseudoknot, substituted downstream of the shifty site fail to promote frameshifting, suggesting that specific features of the MMTV pseudoknot are important for function. Site-directed mutations of the MMTV pseudoknot indicate that the pseudoknot junction, including an unpaired adenine nucleotide between the two stems, provides a specific structural determinant for efficient frameshifting. Pseudoknots derived from other retroviruses (i.e. the feline immunodeficiency virus and the simian retrovirus type 1) also promote frameshifting at the MMTV gag-pro shift site, dependent on the same structure at the junction of the two stems.

The human astrovirus RNA-dependent RNA polymerase coding region is expressed by ribosomal frameshifting

The genomic RNA of human astrovirus serotype 1 (HAst-1) contains three open reading frames (ORFs), 1a, 1b, and 2. ORF 1b is located downstream of, and overlaps, 1a, and it has been suggested on the basis of sequence analysis that expression of ORF 1b is mediated through -1 ribosomal frameshifting. To examine this possibility, a cDNA fragment containing the 1a-1b overlap region was cloned within a reporter gene and placed under the control of the bacteriophage SP6 promoter in a recombinant plasmid. Synthetic transcripts derived from this plasmid, when translated in the rabbit reticulocyte lysate cell-free system, specified the synthesis of polypeptides whose size and antibody reactivity were consistent with an efficient -1 ribosomal frameshift event at the overlap region. The HAst-1 frameshift signal has two essential components, a heptanucleotide slippery sequence, A6C, and a stem-loop structure in the RNA. The presence of this structure was confirmed by complementary and compensatory mutation analysis and by direct structure probing with single- and double-stranded RNA-specific reagents. The HAst-1 frameshift signal, like that present at the overlap of the gag and pro genes of the retrovirus human T-cell lymphotrophic virus type II, does not involve the formation of an RNA pseudoknot.

Identification of unusual RNA folding patterns encoded by bacteriophage T4 gene 60

A 50-nucleotide (nt) untranslated region (coding gap sequence) that interrupts the amino acid coding sequence in T4 gene 60, plus an additional 5 nt upstream and another 3 nt downstream from the gap sequence, shows unusual folding patterns according to RNA structure prediction. A predicted highly stable and significant hairpin structure in the 5' half of the gap sequence and a plausible tertiary structural element computed in the 3' part of the gap sequence seem significant by statistical tests on the wild-type (wt) sequence. This feature is absent in insertion, deletion and substitution variants of the gap sequence, in which template activities are markedly lower than that of the wt. The proposed feature is consistent with currently available data showing that the translational bypass of the coding gap is correlated with a stop codon involved in a stem-loop structure folded in the gap sequence. We suggest that the role of this segment in 'ribosomal bypass' of a portion of the mRNA sequence is a property of its special folded structure.

Mutational analysis of the "slippery-sequence" component of a coronavirus ribosomal frameshifting signal

The ribosomal frameshift signal in the genomic RNA of the coronavirus IBV is composed of two elements, a heptanucleotide "slippery-sequence" and a downstream RNA pseudoknot. We have investigated the kinds of slippery sequence that can function at the IBV frameshift site by analysing the frameshifting properties of a series of slippery-sequence mutants. We firstly confirmed that the site of frameshifting in IBV was at the heptanucleotide stretch UUUAAAC, and then used our knowledge of the pseudoknot structure and a suitable reporter gene to prepare an expression construct that allowed both the magnitude and direction of ribosomal frameshifting to be determined for candidate slippery sequences. Our results show that in almost all of the sequences tested, frameshifting is strictly into the -1 reading frame. Monotonous runs of nucleotides, however, gave detectable levels of a -2/+1 frameshift product, and U stretches in particular gave significant levels (2% to 21%). Preliminary evidence suggests that the RNA pseudoknot may play a role in influencing frameshift direction. The spectrum of slip-sequences tested in this analysis included all those known or suspected to be utilized in vivo. Our results indicate that triplets of A, C, G and U are functional when decoded in the ribosomal P-site following slippage (XXXYYYN) although C triplets were the least effective. In the A-site (XXYYYYN), triplets of C and G were non-functional. The identity of the nucleotide at position 7 of the slippery sequence (XXXYYYN) was found to be a critical determinant of frameshift efficiency and we show that a hierarchy of frameshifting exists for A-site codons. These observations lead us to suggest that ribosomal frameshifting at a particular site is determined, at least in part, by the strength of the interaction of normal cellular tRNAs with the A-site codon and does not necessarily involve specialized "shifty" tRNAs.

Pseudo--half-knot formation with RNA

A pseudo--half-knot can be formed by binding an oligonucleotide asymmetrically to an RNA hairpin loop. This binding motif was used to target the human immunodeficiency virus TAR element, an important viral RNA structure that is the receptor for Tat, the major viral transactivator protein. Oligonucleotides complementary to different halves of the TAR structure bound with greater affinity than molecules designed to bind symmetrically around the hairpin. The pseudo--half-knot--forming oligonucleotides altered the TAR structure so that specific recognition and binding of a Tat-derived peptide was disrupted. This general binding motif may be used to disrupt the structure of regulatory RNA hairpins.

Mutational analysis of the pseudoknot in the tRNA-like structure of turnip yellow mosaic virus RNA. Aminoacylation efficiency and RNA pseudoknot stability

Site-directed mutations were introduced in the connecting loops and one of the two stem regions of the RNA pseudoknot in the tRNA-like structure of turnip yellow mosaic virus RNA. The kinetic parameters of valylation for each mutated RNA were determined in a cell-free extract from wheat germ. Structure mapping was performed on most mutants with enzymic probes, like RNase T1, nuclease S1 and cobra venom ribonuclease. An insertion of four A residues in the four-membered connecting loop L1 that crosses the deep groove of the pseudoknot reduces aminoacylation efficiency. Deletions up to three nucleotides do not affect aminoacylation or RNA pseudoknot formation. Deletion of the entire loop abolishes aminoacylation. Although elimination of the pseudoknot is presumed, this could not be demonstrated. Unlike the mutations in loop L1, all mutations in the three-membered connecting loop L2 that crosses the shallow groove of the RNA pseudoknot decrease the aminoacylation efficiency considerably. Nonetheless, the RNA pseudoknot is still present in most mutated RNAs. These results indicate that a number of mutations can be introduced in both loops without abolishing aminoacylation. Results obtained with the introduction of mismatches and A.U base-pairs in stem S1 of the pseudoknot, containing three G.C base-pairs in wild-type RNA, indicate that the pseudoknot is only marginally stable. Our estimation of the gain of free energy due to the pseudoknot formation is at most 2.0 kcal/mol. The pseudoknot structure can, however, be stabilized upon binding the valyl-tRNA synthetase.

Mutational analysis of the RNA pseudoknot component of a coronavirus ribosomal frameshifting signal

The genomic RNA of the coronavirus IBV contains an efficient ribosomal frameshift signal at the junction of the overlapping 1a and 1b open reading frames. The signal is comprised of two elements, a heptanucleotide "slip-site" and a downstream tertiary RNA structure in the form of an RNA pseudoknot. We have investigated the structure of the pseudoknot and its contribution to the frameshift process by analysing the frameshifting properties of a series of pseudoknot mutants. Our results show that the pseudoknot structure closely resembles that which can be predicted from current building rules, although base-pair formation at the region where the two pseudoknot stems are thought to stack co-axially is not a pre-requisite for efficient frameshifting. The stems, however, must be in close proximity to generate a functional structure. In general, the removal of a single base-pair contact in either stem is sufficient to reduce or abolish frameshifting. No primary sequence determinants in the stems or loops appear to be involved in the frameshift process; as long as the overall structure is maintained, frameshifting is highly efficient. Thus, small insertions into the pseudoknot loops and a deletion in loop 2 that reduced its length to the predicted functional minimum did not influence frameshifting. However, a large insertion (467 nucleotides) into loop 2 abolished frameshifting. A simple stem-loop structure with a base-paired stem of the same length and nucleotide composition as the stacked stems of the pseudoknot could not functionally replace the pseudoknot, suggesting that some particular conformational feature of the pseudoknot determines its ability to promote frameshifting.

RNA pseudoknots: translational frameshifting and readthrough on viral RNAs

Ribosomal frameshifting on retroviral RNAs has been proposed to be mediated by slippage of two adjacent tRNAs into the -1 direction at a specific heptanucleotide sequence. Here we report a computer-aided analysis of the structure around the established or putative frameshift sites in a number of retroviral, coronaviral, toroviral, and luteoviral RNAs and two dsRNA yeast viruses. In almost all cases a stable hairpin was predicted four to nine nucleotides downstream of the shifty heptanucleotide. More than half of the resulting hairpin loops give rise to potential pseudoknotting with sequences downstream of this hairpin. Especially in the case of the shifty heptanucleotides U UUA AAC and G GGA AAC, stable downstream pseudoknots are present. Indications were also found for the presence of pseudoknots downstream of amber stop condons at readthrough sites in some retroviral RNAs.

Prediction of RNA secondary structure, including pseudoknotting, by computer simulation

A computer program is presented which determines the secondary structure of linear RNA molecules by simulating a hypothetical process of folding. This process implies the concept of 'nucleation centres', regions in RNA which locally trigger the folding. During the simulation, the RNA is allowed to fold into pseudoknotted structures, unlike all other programs predicting RNA secondary structure. The simulation uses published, experimentally determined free energy values for nearest neighbour base pair stackings and loop regions, except for new extrapolated values for loops larger than seven nucleotides. The free energy value for a loop arising from pseudoknot formation is set to a single, estimated value of 4.2 kcal/mole. Especially in the case of long RNA sequences, our program appears superior to other secondary structure predicting programs described so far, as tests on tRNAs, the LSU intron of Tetrahymena thermophila and a number of plant viral RNAs show. In addition, pseudoknotted structures are often predicted successfully. The program is written in mainframe APL and is adapted to run on IBM compatible PCs, Atari ST and Macintosh personal computers. On an 8 MHz 8088 standard PC without coprocessor, using STSC APL, it folds a sequence of 700 nucleotides in one and a half hour.