Mapping the structural landscape of the yeast Ty3 retrotransposon RNA genome.

Long terminal repeat (LTR)-retrotransposons are significant contributors to the evolution and diversity of eukaryotic genomes. Their RNA genomes (gRNA) serve as a template for protein synthesis and reverse transcription to a DNA copy, which can integrate into the host genome. Here, we used the SHAPE-MaP strategy to explore Ty3 retrotransposon gRNA structure in yeast and under cell-free conditions. Our study reveals the structural dynamics of Ty3 gRNA and the well-folded core, formed independently of the cellular environment. Based on the detailed map of Ty3 gRNA structure, we characterized the structural context of cis-acting sequences involved in reverse transcription and frameshifting. We also identified a novel functional sequence as a potential initiator for Ty3 gRNA dimerization. Our data indicate that the dimer is maintained by direct interaction between short palindromic sequences at the 5' ends of the two Ty3 gRNAs, resembling the model characteristic for other retroelements like HIV-1 and Ty1. This work points out a range of cell-dependent and -independent Ty3 gRNA structural changes that provide a solid background for studies on RNA structure-function relationships important for retroelement biology.

The silencing of ets-4 mRNA relies on the functional cooperation between REGE-1/Regnase-1 and RLE-1/Roquin-1.

Regnase-1 is an evolutionarily conserved endoribonuclease. It degrades diverse mRNAs important for many biological processes including immune homeostasis, development and cancer. There are two competing models of Regnase-1-mediated mRNA silencing. One model postulates that Regnase-1 works together with another RNA-binding protein, Roquin-1, which recruits Regnase-1 to specific mRNAs. The other model proposes that the two proteins function separately. Studying REGE-1, the Caenorhabditis elegans ortholog of Regnase-1, we have uncovered its functional relationship with RLE-1, the nematode counterpart of Roquin-1. While both proteins are essential for mRNA silencing, REGE-1 and RLE-1 appear to associate with target mRNA independently of each other. Thus, although the functional interdependence between REGE-1/Regnase-1 and RLE-1/Roquin-1 is conserved, the underlying mechanisms may display species-specific variation, providing a rare perspective on the evolution of this important post-transcriptional regulatory mechanism.

In vivo structure of the Ty1 retrotransposon RNA genome.

Long terminal repeat (LTR)-retrotransposons constitute a significant part of eukaryotic genomes and influence their function and evolution. Like other RNA viruses, LTR-retrotransposons efficiently utilize their RNA genome to interact with host cell machinery during replication. Here, we provide the first genome-wide RNA secondary structure model for a LTR-retrotransposon in living cells. Using SHAPE probing, we explore the secondary structure of the yeast Ty1 retrotransposon RNA genome in its native in vivo state and under defined in vitro conditions. Comparative analyses reveal the strong impact of the cellular environment on folding of Ty1 RNA. In vivo, Ty1 genome RNA is significantly less structured and more dynamic but retains specific well-structured regions harboring functional cis-acting sequences. Ribosomes participate in the unfolding and remodeling of Ty1 RNA, and inhibition of translation initiation stabilizes Ty1 RNA structure. Together, our findings support the dual role of Ty1 genomic RNA as a template for protein synthesis and reverse transcription. This study also contributes to understanding how a complex multifunctional RNA genome folds in vivo, and strengthens the need for studying RNA structure in its natural cellular context.

Computational Pipeline for Reference-Free Comparative Analysis of RNA 3D Structures Applied to SARS-CoV-2 UTR Models.

RNA is a unique biomolecule that is involved in a variety of fundamental biological functions, all of which depend solely on its structure and dynamics. Since the experimental determination of crystal RNA structures is laborious, computational 3D structure prediction methods are experiencing an ongoing and thriving development. Such methods can lead to many models; thus, it is necessary to build comparisons and extract common structural motifs for further medical or biological studies. Here, we introduce a computational pipeline dedicated to reference-free high-throughput comparative analysis of 3D RNA structures. We show its application in the RNA-Puzzles challenge, in which five participating groups attempted to predict the three-dimensional structures of 5'- and 3'-untranslated regions (UTRs) of the SARS-CoV-2 genome. We report the results of this puzzle and discuss the structural motifs obtained from the analysis. All simulated models and tools incorporated into the pipeline are open to scientific and academic use.

RNA Binding Properties of the Ty1 LTR-Retrotransposon Gag Protein.

A universal feature of retroelement propagation is the formation of distinct nucleoprotein complexes mediated by the Gag capsid protein. The Ty1 retrotransposon Gag protein from Saccharomyces cerevisiae lacks sequence homology with retroviral Gag, but is functionally related. In addition to capsid assembly functions, Ty1 Gag promotes Ty1 RNA dimerization and cyclization and initiation of reverse transcription. Direct interactions between Gag and retrotransposon genomic RNA (gRNA) are needed for Ty1 replication, and mutations in the RNA-binding domain disrupt nucleation of retrosomes and assembly of functional virus-like particles (VLPs). Unlike retroviral Gag, the specificity of Ty1 Gag-RNA interactions remain poorly understood. Here we use microscale thermophoresis (MST) and electrophoretic mobility shift assays (EMSA) to analyze interactions of immature and mature Ty1 Gag with RNAs. The salt-dependent experiments showed that Ty1 Gag binds with high and similar affinity to different RNAs. However, we observed a preferential interaction between Ty1 Gag and Ty1 RNA containing a packaging signal (Psi) in RNA competition analyses. We also uncover a relationship between Ty1 RNA structure and Gag binding involving the pseudoknot present on Ty1 gRNA. In all likelihood, the differences in Gag binding affinity detected in vitro only partially explain selective Ty1 RNA packaging into VLPs in vivo.

On the Way to Understanding the Interplay between the RNA Structure and Functions in Cells: A Genome-Wide Perspective.

RNAs adopt specific structures in order to perform their biological activities. The structure of RNA is an important layer of gene expression regulation, and can impact a plethora of cellular processes, starting with transcription, RNA processing, and translation, and ending with RNA turnover. The development of high-throughput technologies has enabled a deeper insight into the sophisticated interplay between the structure of the cellular transcriptome and the living cells environment. In this review, we present the current view on the RNA structure in vivo resulting from the most recent transcriptome-wide studies in different organisms, including mammalians, yeast, plants, and bacteria. We focus on the relationship between the mRNA structure and translation, mRNA stability and degradation, protein binding, and RNA posttranscriptional modifications.

RNA-Puzzles Round III: 3D RNA structure prediction of five riboswitches and one ribozyme.

RNA-Puzzles is a collective experiment in blind 3D RNA structure prediction. We report here a third round of RNA-Puzzles. Five puzzles, 4, 8, 12, 13, 14, all structures of riboswitch aptamers and puzzle 7, a ribozyme structure, are included in this round of the experiment. The riboswitch structures include biological binding sites for small molecules (S-adenosyl methionine, cyclic diadenosine monophosphate, 5-amino 4-imidazole carboxamide riboside 5'-triphosphate, glutamine) and proteins (YbxF), and one set describes large conformational changes between ligand-free and ligand-bound states. The Varkud satellite ribozyme is the most recently solved structure of a known large ribozyme. All puzzles have established biological functions and require structural understanding to appreciate their molecular mechanisms. Through the use of fast-track experimental data, including multidimensional chemical mapping, and accurate prediction of RNA secondary structure, a large portion of the contacts in 3D have been predicted correctly leading to similar topologies for the top ranking predictions. Template-based and homology-derived predictions could predict structures to particularly high accuracies. However, achieving biological insights from de novo prediction of RNA 3D structures still depends on the size and complexity of the RNA. Blind computational predictions of RNA structures already appear to provide useful structural information in many cases. Similar to the previous RNA-Puzzles Round II experiment, the prediction of non-Watson-Crick interactions and the observed high atomic clash scores reveal a notable need for an algorithm of improvement. All prediction models and assessment results are available at http://ahsoka.u-strasbg.fr/rnapuzzles/.

Retroviral-like determinants and functions required for dimerization of Ty1 retrotransposon RNA.

During replication of long terminal repeat (LTR)-retrotransposons, their proteins and genome (g) RNA assemble into virus-like particles (VLPs) that are not infectious but functionally related to retroviral virions. Both virions and VLPs contain gRNA in a dimeric form, but contrary to retroviruses, little is known about how gRNA dimerization and packaging occurs in LTR-retrotransposons. The LTR-retrotransposon Ty1 from Saccharomyces cerevisiae is an informative model for studying LTR-retrotransposon and retrovirus replication. Using structural, mutational and functional analyses, we explored dimerization of Ty1 genomic RNA. We provide direct evidence that interactions of self-complementary PAL1 and PAL2 palindromic sequences localized within the 5'UTR are essential for Ty1 gRNA dimer formation. Mutations disrupting PAL1-PAL2 complementarity restricted RNA dimerization in vitro and Ty1 mobility in vivo. Although dimer formation and mobility of these mutants was inhibited, our work suggests that Ty1 RNA can dimerize via alternative contact points. In contrast to previous studies, we cannot confirm a role for PAL3, tRNAiMet as well as recently proposed initial kissing-loop interactions in dimer formation. Our data also supports the critical role of Ty1 Gag in RNA dimerization. Mature Ty1 Gag binds in the proximity of sequences involved in RNA dimerization and tRNAiMet annealing, but the 5' pseudoknot in Ty1 RNA may constitute a preferred Gag-binding site. Taken together, these results expand our understanding of genome dimerization and packaging strategies utilized by LTR-retroelements.

For good and for bad: the role of endogenous retroelements in humans / Na dobre i na zle: rola endogennych retroelementów u czlowieka.

Transposable elements (TEs) are the sequences that are able to "jump" across the genome. They are found in virtually all organisms including human. Although in human, the majority of TEs lost their ability to autonomous transposition, they make up almost half of our genome, and played important roles in genome evolution. Fast progress in deep sequencing and functional analysis has revealed the importance of domes­ticated copies of transposable elements, including their regulatory sequences, transcripts and proteins in normal cells functioning. However, a growing numer of evidence suggest the involvment of TEs in development and progression of autoimmune and neurodegenerative disaeses as well as in many types of cancer. In this review we summarize the current state of knowledge about the LTR retroelements: endogenous retroviruses (ERVs) and Ty3/Gypsy retrotransposons, and their role in human organism.

RNAComposer and RNA 3D structure prediction for nanotechnology.

RNAs adopt specific, stable tertiary architectures to perform their activities. Knowledge of RNA tertiary structure is fundamental to understand RNA functions beginning with transcription and ending with turnover. Contrary to advanced RNA secondary structure prediction algorithms, which allow good accuracy when experimental data are integrated into the prediction, tertiary structure prediction of large RNAs still remains a significant challenge. However, the field of RNA tertiary structure prediction is rapidly developing and new computational methods based on different strategies are emerging. RNAComposer is a user-friendly and freely available server for 3D structure prediction of RNA up to 500 nucleotide residues. RNAComposer employs fully automated fragment assembly based on RNA secondary structure specified by the user. Importantly, this method allows incorporation of distance restraints derived from the experimental data to strengthen the 3D predictions. The potential and limitations of RNAComposer are discussed and an application to RNA design for nanotechnology is presented.

Ty1 retrovirus-like element Gag contains overlapping restriction factor and nucleic acid chaperone functions.

Ty1 Gag comprises the capsid of virus-like particles and provides nucleic acid chaperone (NAC) functions during retrotransposition in budding yeast. A subgenomic Ty1 mRNA encodes a truncated Gag protein (p22) that is cleaved by Ty1 protease to form p18. p22/p18 strongly inhibits transposition and can be considered an element-encoded restriction factor. Here, we show that only p22 and its short derivatives restrict Ty1 mobility whereas other regions of GAG inhibit mobility weakly if at all. Mutational analyses suggest that p22/p18 is synthesized from either of two closely spaced AUG codons. Interestingly, AUG1p18 and AUG2p18 proteins display different properties, even though both contain a region crucial for RNA binding and NAC activity. AUG1p18 shows highly reduced NAC activity but specific binding to Ty1 RNA, whereas AUG2p18 shows the converse behavior. p22/p18 affects RNA encapsidation and a mutant derivative defective for RNA binding inhibits the RNA chaperone activity of the C-terminal region (CTR) of Gag-p45. Moreover, affinity pulldowns show that p18 and the CTR interact. These results support the idea that one aspect of Ty1 restriction involves inhibition of Gag-p45 NAC functions by p22/p18-Gag interactions.

The matrix domain contributes to the nucleic acid chaperone activity of HIV-2 Gag.

BACKGROUND: The Gag polyprotein is a multifunctional regulator of retroviral replication and major structural component of immature virions. The nucleic acid chaperone (NAC) activity is considered necessary to retroviral Gag functions, but so far, NAC activity has only been confirmed for HIV-1 and RSV Gag polyproteins. The nucleocapsid (NC) domain of Gag is proposed to be crucial for interactions with nucleic acids and NAC activity. The major function of matrix (MA) domain is targeting and binding of Gag to the plasma membrane but MA can also interact with RNA and influence NAC activity of Gag. Here, we characterize RNA binding properties and NAC activity of HIV-2 MA and Gag, lacking p6 domain (GagΔp6) and discuss potential contribution of NC and MA domains to HIV-2 GagΔp6 functions and interactions with RNA. RESULTS: We found that HIV-2 GagΔp6 is a robust nucleic acid chaperone. HIV-2 MA protein promotes nucleic acids aggregation and tRNA(Lys3) annealing in vitro. The NAC activity of HIV-2 NC is affected by salt which is in contrast to HIV-2 GagΔp6 and MA. At a physiological NaCl concentration the tRNA(Lys3) annealing activity of HIV-2 GagΔp6 or MA is higher than HIV-2 NC. The HIV-2 NC and GagΔp6 show strong binding to the packaging signal (Ψ) of HIV-2 RNA and preference for the purine-rich sequences, while MA protein binds mainly to G residues without favouring Ψ RNA. Moreover, HIV-2 GagΔp6 and NC promote HIV-2 RNA dimerization while our data do not support MA domain participation in this process in vitro. CONCLUSIONS: We present that contrary to HIV-1 MA, HIV-2 MA displays NAC activity and we propose that MA domain may enhance the activity of HIV-2 GagΔp6. The role of the MA domain in the NAC activity of Gag may differ significantly between HIV-1 and HIV-2. The HIV-2 NC and MA interactions with RNA are not equivalent. Even though both NC and MA can facilitate tRNA(Lys3) annealing, MA does not participate in RNA dimerization in vitro. Our data on HIV-2 indicate that the role of the MA domain in the NAC activity of Gag differs not only between, but also within, retroviral genera.

A self-encoded capsid derivative restricts Ty1 retrotransposition in Saccharomyces.

Retrotransposons and retroviral insertions have molded the genomes of many eukaryotes. Since retroelements transpose via an RNA intermediate, the additive nature of the replication cycle can result in massive increases in copy number if left unchecked. Host organisms have countered with several defense systems, including domestication of retroelement genes that now act as restriction factors to minimize propagation. We discovered a novel truncated form of the Saccharomyces Ty1 retrotransposon capsid protein, dubbed p22 that inhibits virus-like particle (VLP) assembly and function. The p22 restriction factor expands the repertoire of defense proteins targeting the capsid and highlights a novel host-parasite strategy. Instead of inhibiting all transposition by domesticating the restriction gene as a distinct locus, Ty1 and budding yeast may have coevolved a relationship that allows high levels of transposition when Ty1 copy numbers are low and progressively less transposition as copy numbers rise. Here, we offer a perspective on p22 restriction, including its mode of expression, effect on VLP functions, interactions with its target, properties as a nucleic acid chaperone, similarities to other restriction factors, and future directions.

Similarities and differences in the nucleic acid chaperone activity of HIV-2 and HIV-1 nucleocapsid proteins in vitro.

BACKGROUND: The nucleocapsid domain of Gag and mature nucleocapsid protein (NC) act as nucleic acid chaperones and facilitate folding of nucleic acids at critical steps of retroviral replication cycle. The basic N-terminus of HIV-1 NC protein was shown most important for the chaperone activity. The HIV-2 NC (NCp8) and HIV-1 NC (NCp7) proteins possess two highly conserved zinc fingers, flanked by basic residues. However, the NCp8 N-terminal domain is significantly shorter and contains less positively charged residues. This study characterizes previously unknown, nucleic acid chaperone activity of the HIV-2 NC protein. RESULTS: We have comparatively investigated the in vitro nucleic acid chaperone properties of the HIV-2 and HIV-1 NC proteins. Using substrates derived from the HIV-1 and HIV-2 genomes, we determined the ability of both proteins to chaperone nucleic acid aggregation, annealing and strand exchange in duplex structures. Both NC proteins displayed comparable, high annealing activity of HIV-1 TAR DNA and its complementary nucleic acid. Interesting differences between the two NC proteins were discovered when longer HIV substrates, particularly those derived from the HIV-2 genome, were used in chaperone assays. In contrast to NCp7, NCp8 weakly facilitates annealing of HIV-2 TAR RNA to its complementary TAR (-) DNA. NCp8 is also unable to efficiently stimulate tRNALys3 annealing to its respective HIV-2 PBS motif. Using truncated NCp8 peptide, we demonstrated that despite the fact that the N-terminus of NCp8 differs from that of NCp7, this domain is essential for NCp8 activity. CONCLUSION: Our data demonstrate that the HIV-2 NC protein displays reduced nucleic acid chaperone activity compared to that of HIV-1 NC. We found that NCp8 activity is limited by substrate length and stability to a greater degree than that of NCp7. This is especially interesting in light of the fact that the HIV-2 5'UTR is more structured than that of HIV-1. The reduced chaperone activity observed with NCp8 may influence the efficiency of reverse transcription and other key steps of the HIV-2 replication cycle.

The in vitro loose dimer structure and rearrangements of the HIV-2 leader RNA.

RNA dimerization is an essential step in the retroviral life cycle. Dimerization and encapsidation signals, closely linked in HIV-2, are located in the leader RNA region. The SL1 motif and nucleocapsid protein are considered important for both processes. In this study, we show the structure of the HIV-2 leader RNA (+1-560) captured as a loose dimer. Potential structural rearrangements within the leader RNA were studied. In the loose dimer form, the HIV-2 leader RNA strand exists in vitro as a single global fold. Two kissing loop interfaces within the loose dimer were identified: SL1/SL1 and TAR/TAR. Evidence for these findings is provided by RNA probing using SHAPE, chemical reagents, enzymes, non-denaturing PAGE mobility assays, antisense oligonucleotides hybridization and analysis of an RNA mutant. Both TAR and SL1 as isolated domains are bound by recombinant NCp8 protein with high affinity, contrary to the hairpins downstream of SL1. Foot-printing of the SL1/NCp8 complex indicates that the major binding site maps to the SL1 upper stem. Taken together, these data suggest a model in which TAR hairpin III, the segment of SL1 proximal to the loop and the PAL palindromic sequence play specific roles in the initiation of dimerization.

Cell Compartment-Specific Folding of Ty1 Long Terminal Repeat Retrotransposon RNA Genome.

The structural transitions RNAs undergo during trafficking are not well understood. Here, we used the well-developed yeast Ty1 retrotransposon to provide the first structural model of genome (g) RNA in the nucleus from a retrovirus-like transposon. Through a detailed comparison of nuclear Ty1 gRNA structure with those established in the cytoplasm, virus-like particles (VLPs), and those synthesized in vitro, we detected Ty1 gRNA structural alterations that occur during retrotransposition. Full-length Ty1 gRNA serves as the mRNA for Gag and Gag-Pol proteins and as the genome that is reverse transcribed within VLPs. We show that about 60% of base pairs predicted for the nuclear Ty1 gRNA appear in the cytoplasm, and active translation does not account for such structural differences. Most of the shared base pairs are represented by short-range interactions, whereas the long-distance pairings seem unique for each compartment. Highly structured motifs tend to be preserved after nuclear export of Ty1 gRNA. In addition, our study highlights the important role of Ty1 Gag in mediating critical RNA-RNA interactions required for retrotransposition.

RNAthor - fast, accurate normalization, visualization and statistical analysis of RNA probing data resolved by capillary electrophoresis.

RNAs adopt specific structures to perform their functions, which are critical to fundamental cellular processes. For decades, these structures have been determined and modeled with strong support from computational methods. Still, the accuracy of the latter ones depends on the availability of experimental data, for example, chemical probing information that can define pseudo-energy constraints for RNA folding algorithms. At the same time, diverse computational tools have been developed to facilitate analysis and visualization of data from RNA structure probing experiments followed by capillary electrophoresis or next-generation sequencing. RNAthor, a new software tool for the fully automated normalization of SHAPE and DMS probing data resolved by capillary electrophoresis, has recently joined this collection. RNAthor automatically identifies unreliable probing data. It normalizes the reactivity information to a uniform scale and uses it in the RNA secondary structure prediction. Our web server also provides tools for fast and easy RNA probing data visualization and statistical analysis that facilitates the comparison of multiple data sets. RNAthor is freely available at http://rnathor.cs.put.poznan.pl/.

New, extended hairpin form of the TAR-2 RNA domain points to the structural polymorphism at the 5' end of the HIV-2 leader RNA.

The HIV-2 TAR RNA domain (TAR-2) plays a key role in the trans-activation of HIV-2 transcription as it is the target for the Tat-2 protein and several cell factors. Here, we show that the TAR-2 domain exists in vitro in two global, alternative forms: a new, extended hairpin form with two conformers and the already proposed branched hairpins form. This points strongly to the structural polymorphism of the 5' end of the HIV-2 leader RNA. The evidence comes from the non-denaturing PAGE mobility assay, 2D structure prediction, enzymatic and Pb2+- or Mg2+-induced RNA cleavages. Existence of the TAR-2 extended form was further proved by the examination of engineered TAR-2 mutants stabilized either in the branched or extended structure. The TAR-2 extended form predominates with an increasing magnesium concentration. Gel retardation assays reveal that both TAR-2 wt and its mutant, unable to form branched structure, bind Tat-2 protein with comparable, high affinity, while RNA hairpins I and II, derived from TAR-2 branched structure model, show much less protein binding. We propose that an internal loop region of the TAR-2 extended hairpin form is a potential Tat-2 binding site.

Evolutionary plasticity of the NHL domain underlies distinct solutions to RNA recognition.

RNA-binding proteins regulate all aspects of RNA metabolism. Their association with RNA is mediated by RNA-binding domains, of which many remain uncharacterized. A recently reported example is the NHL domain, found in prominent regulators of cellular plasticity like the C. elegans LIN-41. Here we employ an integrative approach to dissect the RNA specificity of LIN-41. Using computational analysis, structural biology, and in vivo studies in worms and human cells, we find that a positively charged pocket, specific to the NHL domain of LIN-41 and its homologs (collectively LIN41), recognizes a stem-loop RNA element, whose shape determines the binding specificity. Surprisingly, the mechanism of RNA recognition by LIN41 is drastically different from that of its more distant relative, the fly Brat. Our phylogenetic analysis suggests that this reflects a rapid evolution of the domain, presenting an interesting example of a conserved protein fold that acquired completely different solutions to RNA recognition.

A normal genetic variation modulates synaptic MMP-9 protein levels and the severity of schizophrenia symptoms.

Matrix metalloproteinase 9 (MMP-9) has recently emerged as a molecule that contributes to pathological synaptic plasticity in schizophrenia, but explanation of the underlying mechanisms has been missing. In the present study, we performed a phenotype-based genetic association study (PGAS) in > 1,000 schizophrenia patients from the Göttingen Research Association for Schizophrenia (GRAS) data collection and found an association between the MMP-9 rs20544 C/T single-nucleotide polymorphism (SNP) located in the 3'untranslated region (UTR) and the severity of a chronic delusional syndrome. In cultured neurons, the rs20544 SNP influenced synaptic MMP-9 activity and the morphology of dendritic spines. We demonstrated that Fragile X mental retardation protein (FMRP) bound the MMP-9 3'UTR We also found dramatic changes in RNA structure folding and alterations in the affinity of FMRP for MMP-9 RNA, depending on the SNP variant. Finally, we observed greater sensitivity to psychosis-related locomotor hyperactivity in Mmp-9 heterozygous mice. We propose a novel mechanism that involves MMP-9-dependent changes in dendritic spine morphology and the pathophysiology of schizophrenia, providing the first mechanistic insights into the way in which the single base change in the MMP-9 gene (rs20544) influences gene function and results in phenotypic changes observed in schizophrenia patients.