Unraveling the Structure-Spectrum Relationship of Yeast Phenylalanine Transfer RNA: Insights from Theoretical Modeling of Infrared Spectroscopy.

Yeast phenylalanine tRNA (tRNAphe) is a paradigmatic model in structural biology. In this work, we combine molecular dynamics simulations and spectroscopy modeling to establish a direct link between its structure, conformational dynamics, and infrared (IR) spectra. Employing recently developed vibrational frequency maps and coupling models, we apply a mixed quantum/classical treatment of the line shape theory to simulate the IR spectra of tRNAphe in the 1600-1800 cm-1 region across its folded and unfolded conformations and under varying concentrations of Mg2+ ions. The predicted IR spectra of folded and unfolded tRNAphe are in good agreement with experimental measurements, validating our theoretical framework. We then elucidate how the characteristic L-shaped tertiary structure of the tRNA and its modulation in response to diverse chemical environments give rise to distinct IR absorption peaks and line shapes. These calculations effectively bridge IR spectroscopy experiments and atomistic molecular simulations, unraveling the molecular origins of the observed IR spectra of tRNAphe. This work presents a robust theoretical protocol for modeling the IR spectroscopy of nucleic acids, which will facilitate its application as a sensitive probe for detecting the fluctuating secondary and tertiary structures of these essential biological macromolecules.

High resolution landscape of ribosomal RNA processing and surveillance.

Ribosomal RNAs are processed in a complex pathway. We profiled rRNA processing intermediates in yeast at single-molecule and single-nucleotide levels with circularization, targeted amplification and deep sequencing (CircTA-seq), gaining significant mechanistic insights into rRNA processing and surveillance. The long form of the 5' end of 5.8S rRNA is converted to the short form and represents an intermediate of a unified processing pathway. The initial 3' end processing of 5.8S rRNA involves trimming by Rex1 and Rex2 and Trf4-mediated polyadenylation. The 3' end of 25S rRNA is formed by sequential digestion by four Rex proteins. Intermediates with an extended A1 site are generated during 5' degradation of aberrant 18S rRNA precursors. We determined precise polyadenylation profiles for pre-rRNAs and show that the degradation efficiency of polyadenylated 20S pre-rRNA critically depends on poly(A) lengths and degradation intermediates released from the exosome are often extensively re-polyadenylated.

Kinetic and structural insights into the requirement of fungal tRNA ligase for a 2'-phosphate end.

Fungal RNA ligase (LIG) is an essential tRNA splicing enzyme that joins 3'-OH,2'-PO4 and 5'-PO4 RNA ends to form a 2'-PO4,3'-5' phosphodiester splice junction. Sealing entails three divalent cation-dependent adenylate transfer steps. First, LIG reacts with ATP to form a covalent ligase-(lysyl-Nζ)-AMP intermediate and displace pyrophosphate. Second, LIG transfers AMP to the 5'-PO4 RNA terminus to form an RNA-adenylate intermediate (A5'pp5'RNA). Third, LIG directs the attack of an RNA 3'-OH on AppRNA to form the splice junction and displace AMP. A defining feature of fungal LIG vis-à-vis canonical polynucleotide ligases is the requirement for a 2'-PO4 to synthesize a 3'-5' phosphodiester bond. Fungal LIG consists of an N-terminal adenylyltransferase domain and a unique C-terminal domain. The C-domain of Chaetomium thermophilum LIG (CthLIG) engages a sulfate anion thought to be a mimetic of the terminal 2'-PO4 Here, we interrogated the contributions of the C-domain and the conserved sulfate ligands (His227, Arg334, Arg337) to ligation of a pRNA2'p substrate. We find that the C-domain is essential for end-joining but dispensable for ligase adenylylation. Mutations H227A, R334A, and R337A slowed the rate of step 2 RNA adenylation by 420-fold, 120-fold, and 60-fold, respectively, vis-à-vis wild-type CthLIG. An R334A-R337A double-mutation slowed step 2 by 580-fold. These results fortify the case for the strictly conserved His-Arg-Arg triad as the enforcer of the 2'-PO4 end-specificity of fungal tRNA ligases and as a target for small molecule interdiction of fungal tRNA splicing.

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.

Structural basis of exoribonuclease-mediated mRNA transcription termination.

Efficient termination is required for robust gene transcription. Eukaryotic organisms use a conserved exoribonuclease-mediated mechanism to terminate the mRNA transcription by RNA polymerase II (Pol II)1-5. Here we report two cryogenic electron microscopy structures of Saccharomyces cerevisiae Pol II pre-termination transcription complexes bound to the 5'-to-3' exoribonuclease Rat1 and its partner Rai1. Our structures show that Rat1 displaces the elongation factor Spt5 to dock at the Pol II stalk domain. Rat1 shields the RNA exit channel of Pol II, guides the nascent RNA towards its active centre and stacks three nucleotides at the 5' terminus of the nascent RNA. The structures further show that Rat1 rotates towards Pol II as it shortens RNA. Our results provide the structural mechanism for the Rat1-mediated termination of mRNA transcription by Pol II in yeast and the exoribonuclease-mediated termination of mRNA transcription in other eukaryotes.

Highly purified DNA-containing cell envelopes from fungi for direct use in PCR.

Efficient DNA sample preparation from fungi with the rigid cell walls is still critical for successful polymerase chain reaction (PCR), one of the basic platforms in molecular diagnostics of fungi, especially in medical mycology. Common methods that involve different chaotropes to yield DNA samples have found a limited application for fungi. Here we describe a novel procedure for efficient production of permeable fungal cell envelopes with DNA inside as suitable templates for PCR. This procedure is facile, relies on boiling of fungal cells in aqueous solutions of selected chaotropic agents and additives and enables to remove RNA and proteins from PCR template samples. The use of chaotropic solutions containing 7 M urea, 1% sodium dodecyl sulfate (SDS), up to100 mM ammonia and/or 25 mM sodium citrate was the best option to yield highly purified DNA-containing cell envelopes from all fungal strains under study, including clinical Candida and Cryptococcusisolates. After treatment with the selected chaotropic mixtures, the fungal cell walls had undergone loosening and were no longer a barrier to release DNA in PCR as evident from electron microscopy examinations and successful target gene amplifications. Overall, the developed simple, fast, and low-cost approach to produce PCR-suitable templates in the form of DNA encased by permeable cell walls can find application in molecular diagnostics.

Comprehensive re-analysis of hairpin small RNAs in fungi reveals loci with conserved links.

RNA interference is an ancient mechanism with many regulatory roles in eukaryotic genomes, with small RNAs acting as their functional element. While there is a wide array of classes of small-RNA-producing loci, those resulting from stem-loop structures (hairpins) have received profuse attention. Such is the case of microRNAs (miRNAs), which have distinct roles in plants and animals. Fungi also produce small RNAs, and several publications have identified miRNAs and miRNA-like (mi/milRNA) hairpin RNAs in diverse fungal species using deep sequencing technologies. Despite this relevant source of information, relatively little is known about mi/milRNA features in fungi, mostly due to a lack of established criteria for their annotation. To systematically assess mi/milRNA characteristics and annotation confidence, we searched for publications describing mi/milRNA loci and re-assessed the annotations for 41 fungal species. We extracted and normalized the annotation data for 1727 reported mi/milRNA loci and determined their abundance profiles, concluding that less than half of the reported loci passed basic standards used for hairpin RNA discovery. We found that fungal mi/milRNA are generally more similar in size to animal miRNAs and were frequently associated with protein-coding genes. The compiled genomic analyses identified 25 mi/milRNA loci conserved in multiple species. Our pipeline allowed us to build a general hierarchy of locus quality, identifying more than 150 loci with high-quality annotations. We provide a centralized annotation of identified mi/milRNA hairpin RNAs in fungi which will serve as a resource for future research and advance in understanding the characteristics and functions of mi/milRNAs in fungal organisms.

Chemical Modifications of Ribosomal RNA.

Cellular RNAs in all three kingdoms of life are modified with diverse chemical modifications. These chemical modifications expand the topological repertoire of RNAs, and fine-tune their functions. Ribosomal RNA in yeast contains more than 100 chemically modified residues in the functionally crucial and evolutionary conserved regions. The chemical modifications in the rRNA are of three types-methylation of the ribose sugars at the C2-positionAbstract (Nm), isomerization of uridines to pseudouridines (Ψ), and base modifications such as (methylation (mN), acetylation (acN), and aminocarboxypropylation (acpN)). The modifications profile of the yeast rRNA has been recently completed, providing an excellent platform to analyze the function of these modifications in RNA metabolism and in cellular physiology. Remarkably, majority of the rRNA modifications and the enzymatic machineries discovered in yeast are highly conserved in eukaryotes including humans. Mutations in factors involved in rRNA modification are linked to several rare severe human diseases (e.g., X-linked Dyskeratosis congenita, the Bowen-Conradi syndrome and the William-Beuren disease). In this chapter, we summarize all rRNA modifications and the corresponding enzymatic machineries of the budding yeast.

Transcriptional neighborhoods regulate transcript isoform lengths and expression levels.

Sequence features of genes and their flanking regulatory regions are determinants of RNA transcript isoform expression and have been used as context-independent plug-and-play modules in synthetic biology. However, genetic context-including the adjacent transcriptional environment-also influences transcript isoform expression levels and boundaries. We used synthetic yeast strains with stochastically repositioned genes to systematically disentangle the effects of sequence and context. Profiling 120 million full-length transcript molecules across 612 genomic perturbations, we observed sequence-independent alterations to gene expression levels and transcript isoform boundaries that were influenced by neighboring transcription. We identified features of transcriptional context that could predict these alterations and used these features to engineer a synthetic circuit where transcript length was controlled by neighboring transcription. This demonstrates how positional context can be leveraged in synthetic genome engineering.

Pseudouridine synthase 7 is an opportunistic enzyme that binds and modifies substrates with diverse sequences and structures.

Pseudouridine (Ψ) is a ubiquitous RNA modification incorporated by pseudouridine synthase (Pus) enzymes into hundreds of noncoding and protein-coding RNA substrates. Here, we determined the contributions of substrate structure and protein sequence to binding and catalysis by pseudouridine synthase 7 (Pus7), one of the principal messenger RNA (mRNA) modifying enzymes. Pus7 is distinct among the eukaryotic Pus proteins because it modifies a wider variety of substrates and shares limited homology with other Pus family members. We solved the crystal structure of Saccharomyces cerevisiae Pus7, detailing the architecture of the eukaryotic-specific insertions thought to be responsible for the expanded substrate scope of Pus7. Additionally, we identified an insertion domain in the protein that fine-tunes Pus7 activity both in vitro and in cells. These data demonstrate that Pus7 preferentially binds substrates possessing the previously identified UGUAR (R = purine) consensus sequence and that RNA secondary structure is not a strong requirement for Pus7-binding. In contrast, the rate constants and extent of Ψ incorporation are more influenced by RNA structure, with Pus7 modifying UGUAR sequences in less-structured contexts more efficiently both in vitro and in cells. Although less-structured substrates were preferred, Pus7 fully modified every transfer RNA, mRNA, and nonnatural RNA containing the consensus recognition sequence that we tested. Our findings suggest that Pus7 is a promiscuous enzyme and lead us to propose that factors beyond inherent enzyme properties (e.g., enzyme localization, RNA structure, and competition with other RNA-binding proteins) largely dictate Pus7 substrate selection.

Single-Base Resolution Mapping Reveals Distinct 5-Formylcytidine in Saccharomyces cerevisiae mRNAs.

5-Formylcytidine (f5C) is one type of post-transcriptional RNA modification, which is known at the wobble position of tRNA in mitochondria and essential for mitochondrial protein synthesis. Here, we show a method to detect f5C modifications in RNA and a transcriptome-wide f5C mapping technique, named f5C-seq. It is developed based on the treatment of pyridine borane, which can reduce f5C to 5,6-dihydrouracil, thus inducing C-to-T transition in f5C sites during PCR to achieve single-base resolution detection. More than 1000 f5C sites were identified after mapping in Saccharomyces cerevisiae by f5C-seq. Moreover, codon composition demonstrated a preference for f5C within wobble sites in mRNA, suggesting the potential role in regulation of translation. These findings expand the scope of the understanding of cytosine modifications in mRNA.

Biochemical approach for isolation of polyadenylated RNAs with bound proteins from yeast.

In vivo characterization of RNA-protein interactions is the key for understanding RNA regulatory mechanisms. Herein, we describe a protocol for detection of proteins interacting with polyadenylated RNAs in the yeast Saccharomyces cerevisiae. Proteins are crosslinked to nucleic acids in vivo by ultraviolet (UV) irradiation of cells, and poly(A)-containing RNAs with bound proteins are isolated from cell lysates using oligo[dT]25 beads. RBPs can be detected by immunoblot analysis or with mass spectrometry to define the mRNA-binding proteome (mRBPome) and its changes under stress. For complete details on the use and execution of this protocol, please refer to Matia-González et al. (2021, 2015).

Harmful R-loops are prevented via different cell cycle-specific mechanisms.

Identifying how R-loops are generated is crucial to know how transcription compromises genome integrity. We show by genome-wide analysis of conditional yeast mutants that the THO transcription complex, prevents R-loop formation in G1 and S-phase, whereas the Sen1 DNA-RNA helicase prevents them only in S-phase. Interestingly, damage accumulates asymmetrically downstream of the replication fork in sen1 cells but symmetrically in the hpr1 THO mutant. Our results indicate that: R-loops form co-transcriptionally independently of DNA replication; that THO is a general and cell-cycle independent safeguard against R-loops, and that Sen1, in contrast to previously believed, is an S-phase-specific R-loop resolvase. These conclusions have important implications for the mechanism of R-loop formation and the role of other factors reported to affect on R-loop homeostasis.

A budding yeast model for human disease mutations in the EXOSC2 cap subunit of the RNA exosome complex.

RNA exosomopathies, a growing family of diseases, are linked to missense mutations in genes encoding structural subunits of the evolutionarily conserved, 10-subunit exoribonuclease complex, the RNA exosome. This complex consists of a three-subunit cap, a six-subunit, barrel-shaped core, and a catalytic base subunit. While a number of mutations in RNA exosome genes cause pontocerebellar hypoplasia, mutations in the cap subunit gene EXOSC2 cause an apparently distinct clinical presentation that has been defined as a novel syndrome SHRF (short stature, hearing loss, retinitis pigmentosa, and distinctive facies). We generated the first in vivo model of the SHRF pathogenic amino acid substitutions using budding yeast by modeling pathogenic EXOSC2 missense mutations (p.Gly30Val and p.Gly198Asp) in the orthologous S. cerevisiae gene RRP4 The resulting rrp4 mutant cells show defects in cell growth and RNA exosome function. Consistent with altered RNA exosome function, we detect significant transcriptomic changes in both coding and noncoding RNAs in rrp4-G226D cells that model EXOSC2 p.Gly198Asp, suggesting defects in nuclear surveillance. Biochemical and genetic analyses suggest that the Rrp4 G226D variant subunit shows impaired interactions with key RNA exosome cofactors that modulate the function of the complex. These results provide the first in vivo evidence that pathogenic missense mutations present in EXOSC2 impair the function of the RNA exosome. This study also sets the stage to compare exosomopathy models to understand how defects in RNA exosome function underlie distinct pathologies.

Overlapping regions of Caf20 mediate its interactions with the mRNA-5'cap-binding protein eIF4E and with ribosomes.

By interacting with the mRNA 5' cap, the translation initiation factor eIF4E plays a critical role in selecting mRNAs for protein synthesis in eukaryotic cells. Caf20 is a member of the family of proteins found across eukaryotes termed 4E-BPs, which compete with eIF4G for interaction with eIF4E. Caf20 independently interacts with ribosomes. Thus, Caf20 modulates the mRNA selection process via poorly understood mechanisms. Here we performed unbiased mutagenesis across Caf20 to characterise which regions of Caf20 are important for interaction with eIF4E and with ribosomes. Caf20 binding to eIF4E is entirely dependent on a canonical motif shared with other 4E-BPs. However, binding to ribosomes is weakened by mutations throughout the protein, suggesting an extended binding interface that partially overlaps with the eIF4E-interaction region. By using chemical crosslinking, we identify a potential ribosome interaction region on the ribosome surface that spans both small and large subunits and is close to a known interaction site of eIF3. The function of ribosome binding by Caf20 remains unclear.

Discovery of 17 conserved structural RNAs in fungi.

Many non-coding RNAs with known functions are structurally conserved: their intramolecular secondary and tertiary interactions are maintained across evolutionary time. Consequently, the presence of conserved structure in multiple sequence alignments can be used to identify candidate functional non-coding RNAs. Here, we present a bioinformatics method that couples iterative homology search with covariation analysis to assess whether a genomic region has evidence of conserved RNA structure. We used this method to examine all unannotated regions of five well-studied fungal genomes (Saccharomyces cerevisiae, Candida albicans, Neurospora crassa, Aspergillus fumigatus, and Schizosaccharomyces pombe). We identified 17 novel structurally conserved non-coding RNA candidates, which include four H/ACA box small nucleolar RNAs, four intergenic RNAs and nine RNA structures located within the introns and untranslated regions (UTRs) of mRNAs. For the two structures in the 3' UTRs of the metabolic genes GLY1 and MET13, we performed experiments that provide evidence against them being eukaryotic riboswitches.

RNA Precipitation.

Precipitation is a critical step to recover RNA of high purity. This chapter describes the principles of alcoholic precipitation as well as a standard, basic protocol with key advices to observe, but numerous variations on the theme are discussed. Indeed, several important parameters, such as the choice of salt, alcohol, or carrier, have to be considered to improve the efficiency of precipitation and the yield of RNA recovery.

Application of NMR Spectroscopy to Determine the 3D Structure of Small Non-Coding RNAs.

Many RNA architectures were discovered to be involved in a wide range of essential biological processes in all organisms from carrying genetic information to gene expression regulation. The remarkable ability of RNAs to adopt various architectures depending on their environment enables the achievement of their myriads of biological functions. Nuclear Magnetic Resonance (NMR) is a powerful technique to investigate both their structure and dynamics. NMR is also a key tool for studying interactions between RNAs and their numerous partners such as small molecules, ions, proteins, or other nucleic acids.In this chapter, to illustrate the use of NMR for 3D structure determination of small noncoding RNA, we describe detailed methods that we used for the yeast C/D box small nucleolar RNA U14 from sample preparation to 3D structure calculation.

Structural and mechanistic insight into stem-loop RNA processing by yeast Pichia stipitis Dicer.

Dicer is a member of the ribonuclease III enzyme family and processes double-stranded RNA into small functional RNAs. The variation in the domain architecture of Dicer among different species whilst preserving its biological dicing function is intriguing. Here, we describe the structure and function of a novel catalytically active RNase III protein, a non-canonical Dicer (PsDCR1), found in budding yeast Pichia stipitis. The structure of the catalytically active region (the catalytic RNase III domain and double-stranded RNA-binding domain 1 [dsRBD1]) of DCR1 showed that RNaseIII domain is structurally similar to yeast RNase III (Rnt1p) but uniquely presents dsRBD1 in a diagonal orientation, forming a catalytic core made of homodimer and large RNA-binding surface. The second dsRNA binding domain at C-terminus, which is absent in Rnt1, enhances the RNA cleavage activity. Although the cleavage pattern of PsDCR1 anchors an apical loop similar to Rnt1, the cleavage activity depended on the sequence motif at the lower stem, not the apical loop, of hairpin RNA. Through RNA sequencing and RNA mutations, we showed that RNA cleavage by PsDCR1 is determined by the stem-loop structure of the RNA substrate, suggesting the possibility that stem-loop RNA-guided gene silencing pathway exists in budding yeast.

Structural basis for conformational equilibrium of the catalytic spliceosome.

The ATPase Prp16 governs equilibrium between the branching (B∗/C) and exon ligation (C∗/P) conformations of the spliceosome. Here, we present the electron cryomicroscopy reconstruction of the Saccharomyces cerevisiae C-complex spliceosome at 2.8 Å resolution and identify a novel C-complex intermediate (Ci) that elucidates the molecular basis for this equilibrium. The exon-ligation factors Prp18 and Slu7 bind to Ci before ATP hydrolysis by Prp16 can destabilize the branching conformation. Biochemical assays suggest that these pre-bound factors prime the C complex for conversion to C∗ by Prp16. A complete model of the Prp19 complex (NTC) reveals how the branching factors Yju2 and Isy1 are recruited by the NTC before branching. Prp16 remodels Yju2 binding after branching, allowing Yju2 to remain tethered to the NTC in the C∗ complex to promote exon ligation. Our results explain how Prp16 action modulates the dynamic binding of step-specific factors to alternatively stabilize the C or C∗ conformation and establish equilibrium of the catalytic spliceosome.