Yeast tRNA ligase (Trl1) is an essential trifunctional enzyme that catalyzes exon-exon ligation during tRNA biogenesis and the non-conventional splicing of HAC1 mRNA during the unfolded protein response (UPR). The UPR regulates the protein folding capacity of the endoplasmic reticulum (ER). ER stress activates Ire1, an ER-resident kinase/RNase, which excises an intron from HAC1 mRNA followed by exon-exon ligation by Trl1. The spliced product encodes for a potent transcription factor that drives the UPR. Here we report the crystal structure of Trl1 RNA ligase domain from Chaetomium thermophilum at 1.9 Å resolution. Structure-based mutational analyses uncovered kinetic competition between RNA ligation and degradation during HAC1 mRNA splicing. Incompletely processed HAC1 mRNA is degraded by Xrn1 and the Ski/exosome complex. We establish cleaved HAC1 mRNA as endogenous substrate for ribosome-associated quality control. We conclude that mRNA decay and surveillance mechanisms collaborate in achieving fidelity of non-conventional mRNA splicing during the UPR.https://doi.org/10.7554/eLife.44199.001
RNA ligases are found in all domains of life. They catalyze the ligation of RNA molecules via phosphodiester bonds during different RNA processing events, such as repair, editing and splicing (Popow et al., 2012). The fungal tRNA ligase Trl1 (previously named Rlg1) is encoded by an essential gene and is involved in tRNA splicing and the unfolded protein response (UPR). Trl1 is a tripartite enzyme (Figure 1A), consisting of an N-terminal adenylyltransferase domain (ligase; LIG), which belongs to the nucleotidyltransferase superfamily alongside DNA ligases and RNA capping enzymes, a C-terminal cyclic phosphodiesterase domain (CPD) and a central polynucleotide kinase (KIN) domain (Phizicky et al., 1986; Xu et al., 1990). Trl1 substrates initially contain a 2’,3’ cyclic phosphate and a 5’-OH group at the RNA termini (Greer et al., 1983). The ligation reaction progresses via three enzymatic steps (Figure 1B). First, the CPD activity opens the 2’,3’ cyclic phosphate by hydrolysis to form a 3’-OH/2’ phosphate terminus, and, second, the KIN activity phosphorylates the 5’-OH in an NTP-dependent reaction – preferring GTP over ATP. These first two steps ‘heal’ (i.e., modify) the RNA termini in preparation for the ligation reaction. Third, the LIG activity ‘seals’ the healed ends through ATP-dependent phosphodiester bond formation. This final reaction occurs in three nucleotidyl transfer steps: (1) Trl1-LIG reacts with ATP to form a covalent LIG-(lysyl-N)-AMP intermediate; (2) the bound AMP is transferred to the 5’ phosphate end to form a 5’-to-5’ RNA-adenylate; (3) Trl1-LIG catalyzes the attack by the 3’-OH on the RNA-adenylate to form a phosphodiester bond, releasing AMP (Greer et al., 1983; Phizicky et al., 1986). The remaining 2’ phosphate at the splice junction is removed by an additional enzyme, Tpt1, which is a 2’ phosphotransferase that, like Trl1, is also essential for cell viability (Banerjee et al., 2019; Culver et al., 1997; Culver et al., 1993).
Trl1 ligates tRNA halves after intron excision by the tRNA splicing endonuclease (SEN) complex (Greer et al., 1983; Peebles et al., 1983). In addition, it is a key component of the UPR, a major intracellular stress signaling pathway (Sidrauski et al., 1996). All eukaryotic cells monitor and adjust the protein folding capacity of their endoplasmic reticulum (ER) through the UPR signaling network. The evolutionarily most conserved – and in fungi sole – branch of the UPR signals via a unique, non-conventional mRNA splicing reaction: Protein folding stress in the ER activates the cytosolic endonuclease domain of the transmembrane stress sensor Ire1 (Cox et al., 1993; Mori et al., 1993). The activated RNase cleaves HAC1 mRNA in fungi or XBP1 mRNA in metazoans at the splice sites (Calfon et al., 2002; Cox and Walter, 1996; Yoshida et al., 2001). Next, a conformational change within the RNA actively ejects the intron and coordinates the two exons (Gonzalez et al., 1999; Peschek et al., 2015), which are ligated by Trl1 in fungi (Sidrauski et al., 1996) or the RTCB ligase complex in metazoans (Jurkin et al., 2014; Kosmaczewski et al., 2014; Lu et al., 2014). The spliced HAC1 and XBP1 mRNAs are then translated to produce the active transcription factors Hac1 and XBP1 that drive expression of UPR target genes in yeast and metazoan cells, respectively (Cox and Walter, 1996; Mori et al., 1996; Nikawa et al., 1996).
Here, we report the crystal structure of the Trl1 ligase domain from Chaetomium thermophilum. Guided by our structural data, we provide evidence that the non-conventional splicing of HAC1 mRNA competes with RNA decay in the cell. We further establish cleaved HAC1 mRNA as an endogenous substrate for ribosome-associated mRNA quality control.
Earlier studies identified several residues within yeast tRNA ligase Trl1 that are critical for its essential function in tRNA splicing. Intriguingly, a yeast genetic screen in Saccharomyces cerevisiae (sc) revealed a single His(148)Tyr point mutation lying within Trl1-LIG that in vivo abolished the UPR signaling yet did not impair tRNA splicing. This mutant allele (henceforth referred to as trl1-H148Y, originally named rlg1-100) paved the way to discovery of Trl1 as the RNA ligase during Ire1-mediated splicing of HAC1 mRNA (Sidrauski et al., 1996). To understand the structural basis for the two distinct roles of Trl1, we sought to crystallize the full-length enzyme, as well as functional modules of it. While our efforts failed using S. cerevisiae-derived protein, we were successful using protein derived from Chaetomium thermophilum (ct), a thermophilic fungus that previously proved invaluable for structural studies (Amlacher et al., 2011; Bock et al., 2014). ctTrl1 has the same tripartite domain structure as scTrl1 (Figure 1A) with 39% overall sequence identity (Figure 1—figure supplement 1). Moreover, expression of ctTrl1 rescued S. cerevisiae trl1∆ cells from the lethal effects of the deletion, both under normal growth conditions as well as during ER stress induced by tunicamycin (Tm), an inhibitor of N-linked glycosylation (Figure 2A). Hence, ctTrl1 catalyzes in S. cerevisiae cells the ligation of both tRNA halves, thus sustaining normal growth, and HAC1 mRNA exons during the UPR. In addition, recombinantly expressed ctTrl1 efficiently ligated the exons of a HAC1-derived RNA substrate, HAC1U-508 (‘U’ for unspliced, that is intron containing), after its endonucleolytic cleavage by Ire1 completing the splicing reaction in vitro (Figure 2B).
Aiming toward structural analysis, we expressed ctTrl1 constructs of varying length with incorporated selenomethionine in Escherichia coli and purified the proteins. While the full-length protein resisted crystallization, we obtained diffracting crystals belonging to space group P212121 of ctTrl1-LIG (residues 1 to 414) bound to the α,β-non-hydrolyzable ATP analog α,β-methyleneadenosine 5′-triphosphate (AMPcPP). The structure was solved by SeMet-SAD phasing to 1.9 Å resolution (see Table 1 for data collection and refinement statistics). The overall structure delineates two subdomains: an N-terminal (LIG-N; aa 13–326) and a C-terminal (LIG-C; aa 327–407) subdomain (Figure 3A, colored blue and yellow, respectively). LIG-N is composed of three central antiparallel β-sheets surrounded by seven α-helices. The overall architecture of this domain resembles the bacteriophage T4 RNA ligase (T4Rnl1) structure consistent with it carrying the adenylyltransferase activity (El Omari et al., 2006). The most N-terminal residues (~35 aa), which wrap around the domain, are absent in T4Rnl1. LIG-C has an all-helical fold that is unrelated, based on structural comparison by Dali (Holm and Laakso, 2016), to all other structures in the PDB from the same superfamily, which besides ATP-dependent RNA ligases also includes DNA ligases and mRNA capping enzymes.
LIG’s active site is buried in the center of a positively charged cleft that spans LIG-N and is extended by LIG-C (Figure 3B). The adenosine nucleotide binding pocket is contained in LIG-N. In our structure, Nζ of Lys148 in the highly conserved nucleotide binding motif I [Kx(D/H/N)G] (Figure 1—figure supplement 1), which is the site of covalent AMP attachment, is more than 4 Å away from the α-phosphonate (Figure 3C), suggestive of conformational rearrangements necessary for the formation of the Trl1-(lysyl-N)-AMP intermediate. Further analysis of the nucleotide binding site identified the residues interacting with AMPcPP. Thr146, Leu147, Glu149, Arg305 and Lys325 form contacts with the adenine base via hydrogen bonding, and His241 via π-stacking. The adenosine nucleoside base of AMPcPP is in syn conformation. Arg99 functions as hydrogen bond donor with the 3’ oxygen of the ribose ring and together with Lys169 forms salt-bridges to the γ-phosphate (Figure 3C).
Two conserved residues, Arg334 and Arg337, within an α-helix (α11) of the LIG-C subdomain coordinate a single well-ordered sulfate ion via salt bridges near the active site pocket. In addition, His227 is within distance to form a third salt bridge, and Asn150 forms a hydrogen bond with the sulfate ion (Figure 3D). Notably, Arg337 and His227 have been previously identified as essential residues by alanine scanning mutagenesis (Wang and Shuman, 2005). We surmise that the sulfate ion represents a surrogate for one of the phosphate groups of Trl1’s RNA substrate.
Based on the strong homology of S. cerevisiae and C. thermophilum Trl1 (Figure 1—figure supplement 1), we found all residues previously identified to be essential for the RNA ligase activity of scTrl1 conserved between both species. When we mapped these residues on the structure of ctTrl1-LIG (Figure 4A), we found that they all cluster around the active site. In addition, these residues are all surface exposed, which points to their importance to binding and correct positioning of RNA substrates. Interestingly, His148 was also identified as an essential residue as its change to Ala, Asn, or Gln scored as lethal in vivo (Wang and Shuman, 2005). By contrast, cells bearing the sc-trl1-H148Y allele were viable and shown previously to block UPR-induced nonconventional mRNA splicing (Sidrauski et al., 1996). This conserved histidine (His182 in ctTrl1) is surface-exposed and located on the periphery of the active-site groove in the vicinity of other essential residues (Figure 4A, colored red). Since the positioning of His182 suggests its participation in substrate binding, we performed in vitro ligation assays with recombinant, purified wild-type (WT) and scTrl1-LIG-H148Y. As substrate we used Ire1-cleaved HAC1U-508 RNA whose ends were enzymatically ‘healed’ (modified) to yield an appropriate substrate for Trl1-LIG (see Materials and methods for details). Equivalent concentrations of the WT enzyme yielded more spliced (i.e. ligated) product of the HAC1U-508 RNA substrate compared to scTrl1-LIG-H148Y (Figure 4B). These results indicate that scTrl1-H148Y is catalytically impaired towards a HAC1-derived RNA substrate but not inactive. To test if reduced substrate affinities cause the inhibiting effect of the H148Y mutation on RNA ligation, we determined the dissociation constants (KD) of scTrl1-LIG-WT and scTrl1-LIG-H148Y towards HAC1-derived 5’ and 3’ exon RNA oligonucleotides using fluorescence titration experiments. Despite the reduced ligation efficiency of the H148Y mutant, we only observed a small impact on binding affinity towards the RNA substrates (Figure 4—figure supplement 1). To obtain a quantitative comparison of the differences in ligation kinetics, we tested both enzymes in an RNA oligonucleotide ligation assay. To this end, we monitored ligation of a fluorescein (FAM)-labeled HAC1-derived 5’ exon RNA oligonucleotide (last 10 nucleotides before the splice site; enzymatically modified 2’ phosphate/3’ OH end) to an unlabeled 3’ exon RNA oligonucleotide (first 20 nucleotides after the splice site; harboring a 5’ phosphate). The isolated ligase domain (scTrl1-LIG-WT or -H148Y) was pre-incubated with the 3’ exon fragment to improve the reaction efficiency. The pre-formed LIG/3’ exon complex was used in excess over the FAM-labeled 5’ exon fragment and the time-course of ligation was analyzed by denaturing urea-PAGE (Figure 4C). The results were fitted to a first-order model (Figure 4D) and revealed a 2.2-fold longer half-life for the H148Y mutant (t1/2 = 9.2 min) compared to WT ligase (t1/2 = 4.1 min). These kinetic data confirmed the impeding effect of the His-to-Tyr mutation on HAC1 exon-exon ligation. Interestingly, the mutant ligase did not only display slower ligation kinetics but also plateaued at a lower amount of ligation product (8%) compared to wild-type enzyme (50%).
To explain how the reduced ligase activity of Trl1-H148Y would impair HAC1 mRNA splicing so severely in cells as to entirely block UPR signaling, we surmised that, in the presence of RNA decay machinery in vivo, the balance between exon ligation and RNA decay might be tipped towards decay. To test this notion directly, we expressed Trl1-H148Y from the strong constitutive alcohol dehydrogenase 1 (ADH1) promoter. Both scTrl1-H148Y and ctTrl1-H182Y allowed growth on Tm-containing medium, as did expression of the corresponding LIG only construct (Figure 5A, compare sector 1 to 2 and 4, and sector 1 to 3 and 5, respectively). The suppression of the mutant phenotype emphasizes that the His148-to-Tyr mutation does not fully inhibit Trl1-LIG functionality and that it is indeed the reduced enzymatic activity of the LIG domain only that becomes limiting, causing the UPR-deficient phenotype of trl1-H148Y cells.
Second, we tested whether unligated, Ire1-cleaved fragments of the HAC1 mRNA are substrates for RNA decay pathways using yeast growth assays (Figure 5B). While the capped 5’ exon represents a possible substrate for 3’-to-5’ exonucleases, the 3’ exon, which still has a poly(A) tail, hence, would likely be degraded in the 5’-to-3’ direction (Figure 5C). The two major conserved exonucleolytic RNA decay enzymes are Xrn1 (5’-to-3’) and the multi-protein RNA exosome complex (3’-to-5’). Since deletion of both pathways is synthetically lethal in yeast, we tested their contribution individually by deletion of XRN1 and SKI2 (encoding for an RNA helicase in the exosome-assisting Ski complex). Indeed, individual deletion of either XRN1 or SKI2 rescued growth of trl1-H148Y cells upon ER stress (Figure 5B). While growth of ski2Δ cells was indistinguishable from that of WT cells, xrn1Δ cells exhibited a reduced colony size, which might be explained by the annotated slow growth phenotype of XRN1 deletion being exacerbated in the presence of ER stress. To ensure that the observed rescue indeed resulted from restored UPR signaling, we generated HAC1 deletions in the respective strains and tested their growth. Neither the xrn1Δ hac1Δ nor the ski2Δ hac1Δ strain were viable in the presence of ER stress (Figure 5—figure supplement 1). The rescue by deletion of the major cytosolic RNA degradation pathway in either the 5’-to-3’ or the 3’-to-5’ direction substantiates that non-conventional splicing and RNA decay are competing outcomes after Ire1 cleavage of HAC1 mRNA.
As mentioned above, the cleaved HAC1 exons should only be susceptible from one end to the attack by exonucleases (Figure 5C); yet, deletion of RNA decay in only one direction was sufficient to allow a HAC1-dependent rescue of the trl-H148Y phenotype. To understand the mechanism of rescue, we used reverse transcription (RT-)PCR to monitor HAC1 mRNA splicing in these cells. To this end, we used oligonucleotide primers in both exons flanking the intron to obtain differently sized amplification products for unspliced and spliced HAC1 mRNA based on the presence or absence of the 252 nt intron. As expected, HAC1 mRNA was spliced upon Tm-induced ER stress in the WT strain and not affected by deletion of XRN1 or SKI2 (Figure 5D, lanes 1–6). In agreement with previous results, upon ER stress we detected an overall reduction of unspliced HAC1 mRNA and no HAC1 mRNA splicing in trl1-H148Y cells (Figure 5D, lanes 7 and 8). Deletion of XRN1 in trl1-H148Y cells restored splicing, albeit not to the same degree as in the WT strain (Figure 5D, lanes 11 and 12). We surmise that stabilization of the severed HAC1 3’ exon in the absence of the major 5’-to-3’ exonuclease Xrn1 allows enzymatically impaired Trl1-H148Y to catch up and ligate the HAC1 exons to restore the UPR.
To our surprise, despite the complete phenotypic rescue shown above (Figure 5B), we did not detect spliced HAC1 mRNA in ski2Δ trl1-H148Y cells upon ER stress (Figure 5D, lanes 9 and 10), but rather a reduction of unspliced HAC1 mRNA, reminiscent of the parental trl1-H148Y mutant strain. We confirmed this result using RT-PCR to amplify the 5’ exon only (Figure 5D, row 2, compare lanes 8 and 10). Both experiments did not account for RNA lacking a poly(A) tail since we used oligo(dT) primers to generate the cDNA templates. Interestingly, when we performed the same RT-PCR on the 5’ exon with cDNA templates from the entire cellular RNA (using random hexamers for priming), we now observed an increase in the amount of the HAC1 5’ exon upon SKI2 deletion, when compared to the parental trl1-H148Y strain (Figure 5D, row 3, compare lanes 8 and 10). Taken together, the RT-PCR analyses suggest that deletion of either XRN1 or SKI2 rescued the UPR-deficiency of the trl1-H148Y allele by different mechanisms: The genetic ablation of the major 5’-to-3’ exonuclease Xrn1 restored HAC1 mRNA splicing in the UPR-deficient yeast cells. By contrast, as we did not detect spliced HAC1 mRNA in ski2Δ trl1-H148Y cells, abrogation of cytosolic RNA exosome by disruption of the Ski complex must have suppressed ER stress sensitivity solely via stabilization of the severed 5’ exon.
Previous studies showed that a truncated form of Hac1 produced from an mRNA bearing a stop codon immediately after the 5’ exon and hence lacking the entire C-terminal portion encoded by the 3’ exon (Hac1trunc) is a fundamentally still functional transcription factor that can induce the UPR (Cox and Walter, 1996; Di Santo et al., 2016). Since the stabilized, cleaved 5’ exon in the ski2Δ trl1-H148Y strain reported here does not contain a stop codon, it should be subject to co-translational mRNA surveillance by no-go decay (NGD). This ribosome recycling pathway hinges on the splitting of ribosomal subunits by the non-canonical release factors Dom34 (Pelota in mammals) and Hbs1. We generated DOM34 deletion strains to test for the possible involvement of NGD in translation of the severed HAC1 5’ exon. Growth analysis of these strains showed that deletion of DOM34 re-sensitized ski2Δ trl1-H148Y cells to ER stress (Figure 6A). Thus, rescue of trl1-H148Y by depletion of Ski2 (i.e. cytosolic RNA exosome activity) is dependent on successful ribosome rescue at the 5’ non-conventional splice site of the cleaved HAC1 mRNA.
To further corroborate the findings from the growth and HAC1 mRNA RT-PCR assays, we analyzed Hac1 levels by immunoblotting of extracts from cells expressing epitope tagged 3x-HA-Hac1. We confirmed the production of Hac1 upon ER stress in WT cells and the lack of its expression in trl1-H148Y cells (Figure 6B, lanes 1–2 and 9–10). Two distinct mechanisms prevent the accumulation of translation product from unspliced HAC1 mRNA. A long-range base pairing interaction between the HAC1U intron and its 5’ untranslated region (UTR) represses Hac1 protein production (Rüegsegger et al., 2001) by inhibiting translation initiation (Sathe et al., 2015). In addition, any translation products from the unspliced HAC1 mRNA are efficiently degraded due to an in-frame degron encoded in the intron (Di Santo et al., 2016). As expected from the results shown above, upon ER stress induction we detected full-length Hac1 in xrn1Δ trl1-H148Y cells (Figure 6B, lane 12). The expression levels were lower than in WT cells, in accordance with the reduced amount of spliced HAC1 mRNA produced in this strain (Figure 5D, compare lane 2 with lane 12). By contrast, while we detected no Hac1 in ski2Δ trl1-H148Y cells exposed to ER stress, we saw instead a faster migrating band (Figure 6B, lane 14), consistent with the expression of a C-terminally truncated form of Hac1 in these cells. In the same strain, we also detected a series of higher molecular mass smeared out bands, presumably a banding pattern resulting from ubiquitination. Interestingly, deletion of SKI2 in the wild-type strain also led to an accumulation of truncated Hac1 and higher molecular mass laddering, albeit to a lesser extent than in the respective trl1-H148Y strain. This result indicated that RNA decay and ribosome-associated mRNA surveillance, while their effects are exacerbated in the UPR mutant strain, play an important role in maintaining the fidelity of HAC1 mRNA splicing in wild-type cells.
Since the discovery of Trl1 as the tRNA ligase in S. cerevisiae by Abelson and co-workers, its atomic structure has remained elusive. Trl1 was later identified in a yeast genetic screen as an indispensable component of the UPR. The deficiency of the underlying H148Y mutation within Trl1-LIG to splice HAC1 mRNA in response to ER stress posed an intriguing mystery. The lack of high-resolution structural data rendered the quest for mechanistic understanding of this phenotype challenging and thus the problem remained unsolved for more than two decades. In the present study, we present the crystal structure of Trl1-LIG and provide new insights into the intricate kinetic competition between HAC1 mRNA splicing and mRNA decay.
The crystal structure of Trl1-LIG from C. thermophilum revealed a two-domain architecture comprised of a canonical adenylyltransferase domain and a C-terminal domain with a unique all-helical fold within the same enzyme superfamily. Previous structural and enzymatic data on T4Rnl1 showed that its C-terminal domain is not required to catalyze RNA ligation but confers specificity for tRNA repair (El Omari et al., 2006). Based on these observations and the structural similarities between the adenylyltransferase domains of ctTrl1-LIG and T4Rnl1, we suspect a similar function, that is conferring substrate specificity, of the C-terminal domain of ctTrl1-LIG. It will be interesting to see if the requirement for a 2’ phosphate on the 5’ exon is mediated by this unique domain.
A major motivation for our efforts towards an atomic structure of Trl1-LIG was to understand the separation of function resulting from the UPR mutant allele trl1-H148Y. The reduced ligation kinetics of Trl1-H148Y tip the balance almost completely in favor of degrading cleaved HAC1 mRNA. Both major exonucleolytic decay pathways of the cytosol, Xrn1 and the RNA exosome with the Ski complex, are part of this kinetic competition of RNA processing enzymes, albeit with different roles.
Our data suggest that Xrn1 degrades the 3’ exon fragment of cleaved HAC1 and plays an important role in allowing efficient translation from the spliced mRNA. The HAC1 intron in S. cerevisiae blocks translation initiation by forming complementary base pairing with the 5’ UTR of the mRNA (Rüegsegger et al., 2001; Sathe et al., 2015). Degradation of the intron after splicing was reported as required for lifting the translational block and the production of Hac1 (Mori et al., 2010). Deletion of XRN1, even in the presence of WT Trl1, resulted in lower levels of Hac1 upon ER stress. Interestingly, Xrn1 also degrades tRNA introns during tRNA splicing, a process with equivalent biochemistry as HAC1 mRNA splicing and shared functionality of Trl1 (Wu and Hopper, 2014). It should be noted that in all cases phosphorylation of the 5’ end, as catalyzed by the central polynucleotide kinase part of Trl1, is required to render the RNA molecules substrates for Xrn1. In consequence, this kinase-mediated decay poses a permanent competition for Trl1-mediated ligation of any RNA substrate in the cell. Degradation of tRNA introns in yeast has been previously reported to be catalyzed by Trl1-KIN and Xrn1in an analogous two-step process (Wu and Hopper, 2014).
We showed that abrogation of the RNA exosome by deletion of SKI2 rescued ER stress sensitivity of trl1-H148Y. However, splicing was not restored in ski2Δ; instead, stabilization of the isolated 5’ exon fragment of HAC1 allowed UPR signaling. Previous studies have shown by adding a stop codon after the 5’ exon of HAC1 that the 18 amino acids encoded in the 3’ exon are dispensable for the transactivation function of Hac1p. In trl1-H148Y ski2Δ, the truncated form of Hac1p would need to be translated and released from an mRNA that lacks a stop codon. Indeed, we could demonstrate, that Hac1trunc is produced from the isolated HAC1 5’ exon in sufficient amounts to initiate the transcriptional response. The nascent chain release in this context is dependent on ribosome splitting by the NGD pathway. This finding is in accordance with a previous study that identified HAC1 mRNA as the strongest Dom34 target by ribosome profiling (Guydosh and Green, 2014). A recent study demonstrated that RIDD in Schizosaccharomyces pombe also depends on NGD-mediated ribosome rescue to allow degradation of cleaved mRNAs by the Ski complex/exosome (Guydosh et al., 2017; Kimmig et al., 2012). Together, these findings establish a crucial role of NGD in both functional outputs of the Ire1 branch of the UPR. Interestingly, production of Hac1trunc coincides with a high-molecular weight laddering that we interpret as ubiquitination of the truncated form. We surmise that the HAC1 mRNA is an endogenous substrate of the recently characterized ribosome quality control (RQC) complex (Brandman et al., 2012; Shen et al., 2015).
The UPR is of crucial importance for cellular homeostasis. In yeast, Hac1 initiates a major transcriptional response that re-shapes the entire ER and upregulates hundreds of target genes comprising ~5% of the genome (Travers et al., 2000). Thus, tight regulation of the response is of great importance for cell physiology. Our data show, that deletion of SKI2 leads to the accumulation of Hac1trunc – even in the context of WT Trl1. It is very likely that single cleavage by Ire1 occurs frequently and that RNA decay pathways assist in degrading these HAC1 mRNA fragments to prevent premature UPR signaling. We conclude that tight regulation of the UPR in yeast is achieved by the interplay of several control mechanisms to which this study adds RNA decay and surveillance by mRNA quality control factors.
Trl1 homologs are present in all human fungal pathogens. They are essential enzymes for them due to their role in tRNA splicing. They also represent promising targets for antifungal drug discovery because their structure and mechanism are distinct from the RTCB-type tRNA ligases in metazoans (as well as archaea and many bacteria). The absence of human homologs of Trl1-LIG could therefore be explored as a promising new drug target in antifungal therapy, whose development has largely stalled in recent years. Together with the recently published structure of Candida albicans Trl1 kinase domain Trl1-KIN (Remus et al., 2017), our study presents new opportunities for inhibitor screens guided by atomic-resolution Trl1 structures.
Full-length ctTrl1 (NCBI Entrez Gene ID: 18257519, CTHT_0034810 tRNA ligase-like protein) was cloned from C. thermophilum cDNA into pET15b (EMD Millipore; pPW3206) as described previously (Peschek et al., 2015). The expression construct for ctTrl1-LIG (residues 1–414; pPW3207) was generated by PCR-based introduction of a stop codon (oligo: 5’-GCTAGCGTGTAGCGTGACATCATTC-3’) into the pET15b vector using site-directed mutagenesis. The expression construct for ctTrl1-CPD (residues 576–846; pPW3208) was generated by cloning the target sequence into pET15b using PCR amplification, followed by restriction digestion and DNA ligation. Similarly, the expression constructs for WT (pPW3209) and H148Y (pPW3210) scTrl1-LIG (residues 1–388) were cloned into pET47b using the restriction enzymes SmaI and BamHI. The expression construct for E. coli RtcA (pPW3420) was generated using the In-Fusion HD Cloning Plus kit (Takara Bio). For the In-Fusion reaction, pET28a (EMD Millipore) was linearized by NdeI (New England BioLabs) and the RtcA insert was amplified from E. coli DNA by PCR. All described constructs were confirmed by DNA sequencing. See Supplementary file 1 for a complete list of plasmids used in this study.
All recombinant proteins in this study were expressed in E. coli BL21-CodonPlus (DE3)-RIL (Agilent Technologies). The cells were grown in Luria broth medium including the appropriate antibiotic at 37°C, expression was induced by 1 mM isopropyl-1-thio-β-D-galactopyranoside (IPTG), cells were harvested by centrifugation and lysed using an EmulsiFlex-C3 (Avestin) high-pressure homogenizer (exceptions to the above as noted).
Recombinant, His6-tagged ctTrl1 was purified by Ni2+ affinity chromatography using a HisTrap FF column (GE Healthcare Life Sciences), followed by size-exclusion chromatography using a HiLoad 16/60 Superdex200 pg column (GE Healthcare Life Sciences) in 20 mM Tris/HCl pH 7.1, 300 mM NaCl, and 1 mM MgCl2 (Peschek et al., 2015).
Selenomethionine (SeMet)-substituted His6-tagged ctTrl1-LIG (residues 1–414) was expressed using feedback inhibition of methionine biosynthesis by adding selected amino acids to the E. coli culture prior to induction. In detail, cells were grown in standard M9 minimal medium supplemented with thiamine (0.5% w/v) and trace elements. When the OD600 reached 0.6, the temperature was lowered to 28°C and the feedback-inhibition amino acids mix (1 g/l of lysine, threonine and phenylalanine, 0.5 g/l leucine, isoleucine and valine, and 0.5 g/l L(+)-SeMet) was added and, after 30 min, the cells were induced with IPTG. After 16 hr expression, cells were harvested and disrupted in lysis buffer (40 mM NaH2PO4/Na2HPO4 pH 7.4, 1.2 M NaCl, 25 mM imidazole, 5 mM dithiothreitol [DTT]) plus protease inhibitor (cOmplete EDTA-free protease inhibitor cocktail, Roche). The cleared supernatant was applied to a Ni2+ affinity chromatography column (HisTrap FF 5 ml, GE Healthcare Life Sciences) and ctTrl1-LIG eluted with a 25–500 mM imidazole gradient (same as lysis buffer). LIG-containing fractions were dialysed against AEX buffer (20 mM HEPES/NaOH pH 8, 10 mM NaCl, 5 mM DTT) and further purified by anion exchange chromatography (HiTrap Q HP 5 ml, GE Healthcare Life Sciences) using a 10–1000 mM NaCl gradient in AEX buffer. Finally, the protein was applied to size-exclusion chromatography using a HiLoad 16/60 Superdex200 pg column (GE Healthcare Life Sciences) in SEC buffer (10 mM HEPES/NaOH pH 7.5, 100 mM NaCl, 5 mM DTT), concentrated and flash frozen in N2 (l).
His6-tagged ctTrl1-CPD (residues 576–846) was expressed at 30°C for 5 hr. Cells were harvested and disrupted in lysis buffer (50 mM Tris/HCl pH 7.5, 1 M NaCl, 20 mM imidazole, 2 mM DTT) plus protease inhibitor (cOmplete EDTA-free protease inhibitor cocktail, Roche). The cleared supernatant was applied to a Ni2+ affinity chromatography column (HisTrap FF 5 ml, GE Healthcare Life Sciences) and ctTrl1-CPD eluted with a 20–500 mM imidazole gradient (same as lysis buffer). CPD-containing fractions were further purified by size-exclusion chromatography using a HiLoad 16/60 Superdex200 pg column (GE Healthcare Life Sciences) in SEC buffer (10 mM Tris/HCl pH 7.5, 100 mM NaCl, 20 mM imidazole, 1 mM tris(2-carboxyethyl)phosphine hydrochloride [TCEP]), concentrated and flash frozen in N2 (l).
His6-tagged scTrl1-LIG-WT and -H148Y (residues 1–388) were expressed at 20°C for 16 hr. Cells were harvested and disrupted in lysis buffer (40 mM HEPES/NaOH pH 7.5, 1 M NaCl, 25 mM imidazole, 2 mM DTT, 5% glycerol) plus protease inhibitor (cOmplete EDTA-free protease inhibitor cocktail, Roche). The cleared supernatant was applied to a Ni2+ affinity chromatography column (HisTrap FF 5 ml, GE Healthcare Life Sciences) and scTrl1-LIG eluted with a 25–500 mM imidazole gradient (same as lysis buffer). LIG-containing fractions were dialysed against 3C cleavage buffer (25 mM HEPES/NaOH pH 8, 75 mM NaCl, 2 mM DTT, 5% glycerol) followed by proteolytic removal of the His6-tag by incubation with recombinant His10-tagged human rhinovirus 3C (HRV3C) protease. To separate the cleaved target protein from the uncleaved fraction and HRV3C protease, another round of Ni2+ affinity chromatography was performed. The flow-through was dialysed against AEX buffer (25 mM HEPES/NaOH pH 8, 25 mM NaCl, 2 mM DTT, 5% glycerol) and further purified by anion exchange chromatography (HiTrap Q HP 5 ml, GE Healthcare Life Sciences) using a 25–1000 mM NaCl gradient in AEX buffer. Finally, the protein was applied to size-exclusion chromatography using a HiLoad 16/60 Superdex200 pg column (GE Healthcare Life Sciences) in SEC buffer (25 mM HEPES/NaOH pH 7.5, 150 mM NaCl, 1 mM TCEP, 5% glycerol), concentrated and flash frozen in N2 (l).
The cytosolic kinase/RNase portion of yeast Ire1 including 32 residues of the linker region (Ire1-KR32) was recombinantly expressed in E. coli as glutathione S-transferase (GST) fusion protein and purified as described previously (Korennykh et al., 2009).
His6-tagged E. coli RtcA was expressed at 37°C for 4 hr. Cells were harvested and disrupted in lysis buffer (50 mM Tris/HCl pH 7.5, 300 mM NaCl, 25 mM imidazole, 2 mM DTT, 10% glycerol) plus protease inhibitor (cOmplete EDTA-free protease inhibitor cocktail, Roche). The cleared supernatant was applied to a Ni2+ affinity chromatography column (HisTrap FF 5 ml, GE Healthcare Life Sciences) and RtcA eluted with a 20–500 mM imidazole gradient (same as lysis buffer). Protein-containing fractions were diluted to 50 mM NaCl and further purified by anion exchange chromatography (HiTrap Q HP 5 ml, GE Healthcare Life Sciences) using a 50–1000 mM NaCl gradient. Finally, the protein was applied to size-exclusion chromatography using a HiLoad 16/60 Superdex75 pg column (GE Healthcare Life Sciences) in SEC buffer (20 mM Tris/HCl pH 7.5, 100 mM NaCl, 1 mM DTT, 5% glycerol), concentrated and flash frozen in N2 (l).
For endonucleolytic cleavage, in vitro transcribed, polyacrylamide gel electrophoresis (PAGE)-purified, refolded HAC1U-508 RNA (Gonzalez et al., 1999) was incubated at 20 ng/μl with Ire1KR32 at 1 μM in cleavage buffer (20 mM HEPES/NaOH pH 7.5, 70 mM NaCl, 2 mM Mg(OAc)2, 1 mM TCEP, and 5% glycerol). For splicing, ligation of the Ire1-cleaved HAC1U-508 RNA was initiated with 500 nM ctTrl1 (1 mM ATP, 1 mM GTP) or the indicated amounts of scTrl1-LIG. To allow ligation by Trl1-LIG, the ends of the cleaved RNA were first modified by 1 μM ctTrl1-CPD and 667 units/ml T4 polynucleotide kinase (3’ phosphatase minus, New England BioLabs) plus 4 mM ATP. Stop solution (10 M urea, 0.1% SDS, 1 mM EDTA) was added at five-fold excess to stop the reactions. Samples were then denatured by heating at 80°C for 3 min and analyzed by 6% TBE–urea gels (Thermo Fisher Scientific) stained with SYBR Gold nucleic acid stain (Thermo Fisher Scientific).
Oligonucleotide-guided RNase H digestion was used to identify RNA products of HAC1 in vitro splicing. The assay was performed similarly as described previously (Mori et al., 2010). After completion of the ligation reaction, 250 pmol of antisense DNA oligonucleotides were added to a total volume of 7.5 μL. The oligonucleotides were designed to hybridize with the HAC1 intron (oligo a: TAATCACGGCGGACAGTA; oligo b: TTGAAGGTACTTTAACCG). Samples were heated to 75°C, incubated at 43°C for 10 min to allow hybridization, then slowly cooled to 37°C and incubated with 5 U (addition of 1 μL) of RNase H (5,000 U/ml; New England BioLabs) for 30 min. After proteinase K treatment, the samples were analyzed by PAGE as described in the previous section.
All HAC1-derived RNA oligonucleotides were purchased from Integrated DNA Technologies. The 5’ exon and 3’ exon RNA substrates were synthesized with 3’ and 5’ phosphate groups at their respective splice site ends. RNA ligation kinetics were measured using RNA oligonucleotides based on the exon sequences closest to the splice sites. We used a 10 nt long 5’ exon fragment with a terminal fluorescein (FAM) fluorophore (5’-FAM-CGUAAUCCAG-3’-PO4) and a 20 nt long 3’ exon fragment (5’-PO4-AAGCGCAGUCAGGUUUGAAU-3’). The 5’exon fragment was modified by RtcA and ctTrl1-CPD in the presence of ATP to yield 2’ phosphate/3’ hydroxyl ends. Both RNA substrates were PAGE-purified prior to the assay. 250 nM scTrl1-LIG (WT or H148Y) were incubated with 500 nM of 3’ exon oligonucleotide and 0.5 mM ATP in assay buffer (20 mM Tris/HCl pH 7.5, 150 mM NaCl, 10 mM MgCl2, 1 mM TCEP, 5% glycerol) for 10 min at 30°C. The ligation kinetics were initiated by addition of 50 nM FAM-labeled 5’ exon fragment. Aliquots were taken at the indicated timepoints and mixed with a 3-fold excess of Stop solution. The ligation reaction was analyzed by 15% TBE–urea gels (Thermo Fisher Scientific) and imaged for FAM fluorescence using a Typhoon 9400 Variable Mode Imager (GE Lifesciences). The relative intensity of the bands was determined by densitometry using ImageJ (Schindelin et al., 2012; Schneider et al., 2012). Kinetic data were fitted to a first-order model using Origin (Version 2019, OriginLab Corporation).
Affinity measurements were performed using RNA oligonucleotides based on the HAC1 exons’ most proximal 10 nucleotides to the splice site. The 5’ exon and 3’ exon RNA substrates were modified with Cy5 fluorophores at their distal termini (5’ exon RNA: 5’-Cy5-CGUAAUCCAG-3’-PO4; 3’ exon RNA: 5’-PO4-AAGCGCAGUC-Cy5). The 5’exon fragment was modified by RtcA and ctTrl1-CPD in the presence of ATP to yield 2’ phosphate/3’ hydroxyl ends. Both RNA substrates were PAGE-purified prior to the binding assay. Binding was measured by changes in fluorescence using a Monolith NT.115 (NanoTemper). The RNA substrate concentration was held constant at 10 nM and the concentration of scTrl1-LIG. RNA binding data were fitted to a Hill model using Origin (Version 2019, OriginLab Corporation).
Crystals of SeMet-ctTrl1-LIG were grown at 20°C by hanging drop vapor diffusion. Two volumes of SeMet-ctTrl1-LIG (12.5 mg/ml), 2 mM AMPcPP (Sigma, M6517) and 2 mM MgCl2 were mixed with one volume of reservoir solution containing 1.7 M ammonium sulfate, 200 mM NaCl, 100 mM HEPES (pH 7.3) and 5 mM DTT. Obtained crystals were harvested directly from the crystallization drops, cryoprotected by soaking in reservoir solution including 25% glycerol and frozen in liquid nitrogen. X-ray diffraction data were collected at the Advanced Photon Source beamline 23ID-D and processed using MOSFLM (Battye et al., 2011) and AIMLESS (Evans and Murshudov, 2013) as implemented in the CCP4 software suite (Winn et al., 2011).
Phasing and initial model building were carried out using CRANK2 for automated structure solution (Skubák and Pannu, 2013). The phases were obtained using single-wavelength anomalous dispersion (SAD) data from a single crystal with SHELX (Sheldrick, 2010). The initial model was built using BUCCANEER (Cowtan, 2006) and used to solve the structure from a second high-resolution dataset. Interactive model building was performed using COOT (Emsley et al., 2010) and the structure was refined with REFMAC5 (Murshudov et al., 2011) including TLS refinement. The model optimization and quality assessment was carried out with PDB-REDO (Joosten et al., 2014) and MOLPROBITY (Davis et al., 2004). All structure figures were generated using PyMOL (The PyMOL Molecular Graphics System, Version 2.0 Schrödinger, LLC) and electrostatic surface potential calculations were performed using PDB2PQR (Dolinsky et al., 2004) and APBS (Baker et al., 2001). Ligand interactions were analyzed using the protein-ligand interaction profiler (PLIP) web service (Salentin et al., 2015). A composite omit map was generated in Phenix (Adams et al., 2010) to exclude model bias and to verify ligand density.
All yeast strains in this study were based on S. cerevisiae W303 that was made TRP1 and ADE2 by repairing the endogenous auxotrophy. All strains are listed in Supplementary file 2. Deletion of non-essential genes was performed in haploid cells, deletion of TRL1 was performed in diploid cells. Each deletion mutant was generated using a PCR-based homologous recombination method as described previously (Janke et al., 2004). Yeast transformations were performed by the standard LiOAc/single stranded carrier DNA/PEG method, cells were recovered on YPD medium and grown on the appropriate selection medium. Functional complementation of trl1Δ and trl1-H148Y strains was performed by transforming the strain of interest with different Trl1 constructs cloned into the p416ADH (ATCC 87376; pPW3211) yeast expression vector (pPW3212-3215; see Supplementary file 1 for a complete list of plasmids).
For growth tests under ER stress, 5-fold serial dilutions were spotted onto either standard YPD plates (-Tm) or YPD plates containing the indicated concentration of tunicamycin (Sigma, T7765) dissolved in DMSO (+Tm). All growth tests were performed as three biological replicates; representative images are shown in the respective figures.
For immunoblot analysis, endogenous HAC1 was N-terminally tagged with a 3x hemagglutinin (HA) tag (n-YPYDVPDYAGYPYDVPDYAGSYPYDVPDYA-c) in the strains of interest using the Delitto perfetto method (Storici et al., 2001).
Total RNA was prepared from yeast using phenol-chloroform extraction followed by ethanol precipitation. The extracted RNA was resolubilized in water, the quality was assessed by agarose gel electrophoresis and the concentration determined by absorbance spectroscopy. For RT-PCR, cDNA was generated using SuperScript II reverse transcriptase (Thermo Fisher Scientific) with oligo-(dT)12-18 or random hexamers primers (both Thermo Fisher Scientific) from 1 μg of total RNA. The cDNAs of interest were amplified from the reverse transcription reactions (1:100 dilution) with Phusion polymerase (Thermo Scientific) using the appropriate primer sets and optimal cycle numbers (see Supplementary file 3 for primer sequences and cycle numbers). RT-PCR reactions were analyzed by agarose gel electrophoresis.
Extraction of yeast proteins was performed as described elsewhere (Kushnirov, 2000). Proteins were separated by SDS-PAGE (any kD Criterion TGX stain-free gels, Bio-Rad) and transferred onto Protran 0.2 μm pore nitrocellulose membrane (Perkin Elmar). Blots were probed with the following antibodies diluted in TBS-T buffer containing 5% nonfat dry milk: mouse anti-HA (1:2,000; Sigma H3663), rabbit anti GAPDH (1:1,000; abcam ab9485), HRP conjugated goat anti-mouse IgG (1:10,000; Promega W4021), HRP conjugated goat anti-rabbit IgG (1:10,000; Promega W4011). Blots were developed using SuperSignal West Dura Extended Duration Substrate (Pierce, Life Technologies) and chemiluminescence was detected on a ChemiDoc XRS + imager (Bio-Rad).
Structural coordinates for the ctTrl1-LIG structure have been deposited in the Protein Data Bank under accession code 6N67.
PHENIX: a comprehensive Python-based system for macromolecular structure solutionActa Crystallographica Section D Biological Crystallography 66:213–221.https://doi.org/10.1107/S0907444909052925
iMOSFLM: a new graphical interface for diffraction-image processing with MOSFLMActa Crystallographica. Section D, Biological Crystallography 67:271–281.https://doi.org/10.1107/S0907444910048675
An integrated approach for genome annotation of the eukaryotic thermophile chaetomium thermophilumNucleic Acids Research 42:13525–13533.https://doi.org/10.1093/nar/gku1147
A 2'-phosphotransferase implicated in tRNA splicing is essential in Saccharomyces cerevisiaeJournal of Biological Chemistry 272:13203–13210.https://doi.org/10.1074/jbc.272.20.13203
MOLPROBITY: structure validation and all-atom contact analysis for nucleic acids and their complexesNucleic Acids Research 32:W615–W619.https://doi.org/10.1093/nar/gkh398
PDB2PQR: an automated pipeline for the setup of Poisson-Boltzmann electrostatics calculationsNucleic Acids Research 32:W665–W667.https://doi.org/10.1093/nar/gkh381
Molecular architecture and ligand recognition determinants for T4 RNA ligaseJournal of Biological Chemistry 281:1573–1579.https://doi.org/10.1074/jbc.M509658200
REFMAC5 for the refinement of macromolecular crystal structuresActa Crystallographica Section D Biological Crystallography 67:355–367.https://doi.org/10.1107/S0907444911001314
Saccharomyces cerevisiae IRE2/HAC1 is involved in IRE1-mediated KAR2 expressionNucleic Acids Research 24:4222–4226.https://doi.org/10.1093/nar/24.21.4222
Structural basis for the GTP specificity of the RNA kinase domain of fungal tRNA ligaseNucleic Acids Research 45:12945–12953.https://doi.org/10.1093/nar/gkx1159
PLIP: fully automated protein-ligand interaction profilerNucleic Acids Research 43:W443–W447.https://doi.org/10.1093/nar/gkv315
Evidence that Base-pairing interaction between intron and mRNA leader sequences inhibits initiation of HAC1 mRNA translation in yeastJournal of Biological Chemistry 290:21821–21832.https://doi.org/10.1074/jbc.M115.649335
Fiji: an open-source platform for biological-image analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019
Experimental phasing with SHELXC/D/E: combining chain tracing with density modificationActa Crystallographica. Section D, Biological Crystallography 66:479–485.https://doi.org/10.1107/S0907444909038360
Automatic protein structure solution from weak X-ray dataNature Communications 4:2777.https://doi.org/10.1038/ncomms3777
In vivo site-directed mutagenesis using oligonucleotidesNature Biotechnology 19:773–776.https://doi.org/10.1038/90837
Overview of the CCP4 suite and current developmentsActa Crystallographica. Section D, Biological Crystallography 67:235–242.https://doi.org/10.1107/S0907444910045749
Nahum SonenbergReviewing Editor; McGill University, Canada
Cynthia WolbergerSenior Editor; Johns Hopkins University School of Medicine, United States
Javier MartinezReviewer; IMBA-Autrian Academy of Sciences, Austria
In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.
Thank you for submitting your article "tRNA ligase structure reveals kinetic competition between non-conventional mRNA splicing and mRNA decay" for consideration by eLife. Your article has been reviewed by two peer reviewers, and the evaluation has been overseen by Nahum Sonenberg, the Reviewing Editor and Cynthia Wolberger as the Senior Editor. The following individual involved in review of your submission has agreed to reveal their identity: Javier Martinez (Reviewer #2).
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
The authors studied the mechanism by which HAC1 pre-mRNA is spliced in the cytoplasm of yeast cells, which is triggered by UPR. To this end, the authors solved the crystal structure of the Saccharomyces cerevisiae Trl1 homolog from Chaetomium thermophilum. They exploited in particular the His-to-Tyr mutation at Trl1 position 148 (within the ligase domain), which abolishes UPR signaling but not splicing, gaining important insight into the HAC1 pre-mRNA splicing mechanism by demonstrating kinetic competition between non-conventional mRNA splicing and RNA degradation, which effects the outcome of Ire1-initiated mRNA processing. Thus, the mRNA decay and surveillance mechanisms interact to ensure fidelity of non-conventional mRNA splicing during the UPR.
Both reviewers agree that the paper presents an interesting study, which would warrant publication in eLife if amended in response to the reviewers' critiques. The Reviewing Editor concurs with the reviewers' opinions. As you can see in the comments below both reviewers believe that the key results are shown in Figure 3, but they argue that some issues related to these results should be experimentally addressed.
1) What is the basis of the reduction in splicing efficiency in the trl1-H148Y enzyme? Does the enzyme have a lower affinity for specific RNA fragments? An answer to this question might provide a link to the understanding of the enhanced sensitivity to Xrn1 or the exosome. This experiment could be readily done.
2) What is the band marked with an asterisk in Figure 1D, lane 3?
3) The results in Figure 4 should be quantified.
4) Discussion: Where is it shown that Xrn1 degrades the 3' exon?https://doi.org/10.7554/eLife.44199.022
Both reviewers agree that the paper presents an interesting study, which would warrant publication in eLife if amended in response to the reviewers' critiques. The Reviewing Editor concurs with the reviewer' opinions. As you can see in the comments below both reviewers believe that the key results are shown in Figure 3, but they argue that some issues related to these results should be experimentally addressed.
1) What is the basis of the reduction in splicing efficiency in the trl1-H148Y enzyme? Does the enzyme have a lower affinity for specific RNA fragments? An answer to this question might provide a link to the understanding of the enhanced sensitivity to Xrn1 or the exosome. This experiment could be readily done.
We agree with the reviewers’ assessment that a detailed mechanistic understanding of the reduced splicing (ligation) efficiency of the H148Y mutant will add to our understanding of the Trl1 ligation mechanism. We tested several different in vitro approaches to determine binding affinities for both exons. We measured affinities by monitoring fluorescence changes of labeled RNA oligonucleotides upon binding to scTrl1-LIG. Unexpectedly, these experiments yielded very similar affinities for WT and mutant Trl1-LIG towards 5’ and 3’ exon fragment (differing by only ~1.4 and ~1.1-fold, respectively), respectively). We included these fluorescence binding measurements as new Figure 4—figure supplement 1.
We conclude from this work that, if binding affinities are not significantly impeded, the reduced activity of the mutant enzyme must result from a changed catalytic activity. To confirm the reduced ligation kinetics (efficiency) of the Trl1-H148Y enzyme (Figure 4B), we performed time-resolved in vitro ligation experiments with labeled RNA oligonucleotides. These results confirmed the effect of the His-to-Tyr mutation, causing a ~2-fold reduced apparent rate constants. Interestingly, the mutant ligase is not only a 2x slower enzyme but also plateaued at a >5-fold lower amount of ligation product compared to the wild-type enzyme. These results are depicted as new Figures 4C (PAGE gel) and 4D (quantification and fitting). Taken together, these results strongly support our kinetic competition model.
2) What is the band marked with an asterisk in Figure 1D, lane 3?
The * demarks a circularized tandem HAC1 intron. We now added experimental evidence that this is indeed the case. To this end, we used antisense DNA oligos complementary to either the presumptive ligation (circularization) site or the middle of the intron. After annealing to the RNA products of the in vitro splicing reaction and treatment with RNase H, which specifically hydrolyzes the phosphodiester bonds of RNA in RNA/DNA hybrids, the band in question disappeared and a band at the exact position of the severed intron (252 nt) re-appeared. We conclude that * represents two cleaved introns that were circularized by two intermolecular ligation events at their 5’ and 3’ ends. The results are depicted in Figure 2—figure supplement 1. We added the RNase H assay to the Materials and methods.
3) The results in Figure 4 should be quantified.
We quantified the HAC1 splicing results and added the results to Figure 5 (previously Figure 4).
4) Discussion: Where is it shown that Xrn1 degrades the 3' exon?
We softened the statement by changing “Our data show…” to “Our data suggest…”. Based on our RT-PCR results (Figure 5D), we could show that levels of the 5’ exon are not affected by the presence or absence of Xrn1. Although we did not directly “show” the degradation of the 3’ exon by Xrn1, it is the only HAC1 mRNA component left in the non-conventional splicing reaction. Our RT-PCR data demonstrate that deletion of XRN1 restores splicing, which requires a fully intact 3’ exon as substrate.https://doi.org/10.7554/eLife.44199.023
- Jirka Peschek
- Peter Walter
- Peter Walter
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
We thank Christof Osman, Voytek Okreglak, Weihan Li, and members of the Walter Lab for their valuable help and insightful discussion; Damian Ekiert, Oren Rosenberg and Lan Wang for helpful advice regarding X-ray crystallography. We would like to thank the 2016 CCP4 School at Argonne including all its instructors for help and guidance to JP during data collection, initial data analysis and structure determination. We thank the entire staff at GM/CA@APS, which has been funded in whole or in part with Federal funds from the National Cancer Institute (ACB-12002) and the National Institute of General Medical Sciences (AGM-12006). This research used resources of the Advanced Photon Source, a US Department of Energy (DOE) Office of Science User Facility operated for the DOE Office of Science by Argonne National Laboratory under Contract No. DE-AC02-06CH11357. JP acknowledges a long-term fellowship from the Human Frontier Science Program. This work was supported by NIH grant GM032384 to PW. PW is an Investigator of the Howard Hughes Medical Institute.
- Cynthia Wolberger, Johns Hopkins University School of Medicine, United States
- Nahum Sonenberg, McGill University, Canada
- Javier Martinez, IMBA-Autrian Academy of Sciences, Austria
- Received: December 6, 2018
- Accepted: June 11, 2019
- Version of Record published: June 25, 2019 (version 1)
© 2019, Peschek and Walter
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.