1
|
Genome-Wide RNA Secondary Structure Prediction. Methods Mol Biol 2023; 2586:35-48. [PMID: 36705897 DOI: 10.1007/978-1-0716-2768-6_3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
The information of RNA secondary structure has been widely applied to the inference of RNA function. However, a classical prediction method is not feasible to long RNAs such as mRNA due to the problems of computational time and numerical errors. To overcome those problems, sliding window methods have been applied while their results are not directly comparable to global RNA structure prediction. In this chapter, we introduce ParasoR, a method designed for parallel computation of genome-wide RNA secondary structures. To enable genome-wide prediction, ParasoR distributes dynamic programming (DP) matrices required for structure prediction to multiple computational nodes. Using the database of not the original DP variable but the ratio of variables, ParasoR can locally compute the structure scores such as stem probability or accessibility on demand. A comprehensive analysis of local secondary structures by ParasoR is expected to be a promising way to detect the statistical constraints on long RNAs.
Collapse
|
2
|
Zhang N, Shi S, Wang X, Ni W, Yuan X, Duan J, Jia TZ, Yoo B, Ziegler A, Russo JJ, Li W, Zhang S. Direct Sequencing of tRNA by 2D-HELS-AA MS Seq Reveals Its Different Isoforms and Dynamic Base Modifications. ACS Chem Biol 2020; 15:1464-1472. [PMID: 32364699 PMCID: PMC7902080 DOI: 10.1021/acschembio.0c00119] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
Post-transcriptional modifications are intrinsic to RNA structure and function. However, methods to sequence RNA typically require a cDNA intermediate and are either not able to sequence these modifications or are tailored to sequence one specific nucleotide modification only. Interestingly, some of these modifications occur with <100% frequency at their particular sites, and site-specific quantification of their stoichiometries is another challenge. Here, we report a direct method for sequencing tRNAPhe without cDNA by integrating a two-dimensional hydrophobic RNA end-labeling strategy with an anchor-based algorithm in mass spectrometry-based sequencing (2D-HELS-AA MS Seq). The entire tRNAPhe was sequenced and the identity, location, and stoichiometry of all eleven different RNA modifications was determined, five of which were not 100% modified, including a 2'-O-methylated G (Gm) in the wobble anticodon position as well as an N2, N2-dimethylguanosine (m22G), a 7-methylguanosine (m7G), a 1-methyladenosine (m1A), and a wybutosine (Y), suggesting numerous post-transcriptional regulations in tRNA. Two truncated isoforms at the 3'-CCA tail of the tRNAPhe (75 nt with a 3'-CC tail (80% abundance) and 74 nt with a 3'-C tail (3% abundance)) were identified in addition to the full-length 3'-CCA-tailed tRNAPhe (76 nt, 17% abundance). We discovered a new isoform with A-G transitions/editing at the 44 and 45 positions in the tRNAPhe variable loop, and discuss possible mechanisms related to the emergence and functions of the isoforms with these base transitions or editing. Our method revealed new isoforms, base modifications, and RNA editing as well as their stoichiometries in the tRNA that cannot be determined by current cDNA-based methods, opening new opportunities in the field of epitranscriptomics.
Collapse
Affiliation(s)
- Ning Zhang
- Department of Biological and Chemical Sciences, New York Institute of Technology, New York, New York 10023, United States
- Department of Chemical Engineering, Columbia University, New York, New York 10027, United States
| | - Shundi Shi
- Department of Chemical Engineering, Columbia University, New York, New York 10027, United States
| | - Xuanting Wang
- Department of Chemical Engineering, Columbia University, New York, New York 10027, United States
| | - Wenhao Ni
- Department of Biological and Chemical Sciences, New York Institute of Technology, New York, New York 10023, United States
| | - Xiaohong Yuan
- Department of Biological and Chemical Sciences, New York Institute of Technology, New York, New York 10023, United States
| | - Jiachen Duan
- Department of Biological and Chemical Sciences, New York Institute of Technology, New York, New York 10023, United States
| | - Tony Z Jia
- Earth-Life Science Institute, Tokyo Institute of Technology, Meguro-ku, Tokyo 152-8550, Japan
- Blue Marble Space Institute of Science, Seattle, Washington 98154, United States
| | - Barney Yoo
- Department of Chemistry, Hunter College, City University of New York, New York, New York 10065, United States
| | - Ashley Ziegler
- Department of Biological and Chemical Sciences, New York Institute of Technology, New York, New York 10023, United States
| | - James J Russo
- Department of Chemical Engineering, Columbia University, New York, New York 10027, United States
| | - Wenjia Li
- Department of Computer Science, New York Institute of Technology, New York, New York 10023, United States
| | - Shenglong Zhang
- Department of Biological and Chemical Sciences, New York Institute of Technology, New York, New York 10023, United States
| |
Collapse
|
3
|
Wang Z, Gaba A, Sachs MS. A highly conserved mechanism of regulated ribosome stalling mediated by fungal arginine attenuator peptides that appears independent of the charging status of arginyl-tRNAs. J Biol Chem 1999; 274:37565-74. [PMID: 10608810 DOI: 10.1074/jbc.274.53.37565] [Citation(s) in RCA: 55] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
The Arg attenuator peptide (AAP) is an evolutionarily conserved peptide involved in Arg-specific negative translational control. It is encoded as an upstream open reading frame (uORF) in fungal mRNAs specifying the small subunit of Arg-specific carbamoyl phosphate synthetase. We examined the functions of the Saccharomyces cerevisiae CPA1 and Neurospora crassa arg-2 AAPs using translation extracts from S. cerevisiae, N. crassa, and wheat germ. Synthetic RNA containing AAP and firefly luciferase (LUC) sequences were used to program translation; analyses of LUC activity indicated that the AAPs conferred Arg-specific negative regulation in each system. The AAPs functioned either as uORFs or fused in-frame at the N terminus of LUC. Mutant AAPs lacking function in vivo did not function in vitro. Therefore, trans-acting factors conferring AAP-mediated regulation are in both fungal and plant systems. Analyses of ribosome stalling in the fungal extracts by primer extension inhibition (toeprint) assays showed that these AAPs acted similarly to stall ribosomes in the region immediately distal to the AAP coding region in response to Arg. The regulatory effect increased as the Arg concentration increased; all of the arginyl-tRNAs examined appeared maximally charged at low Arg concentrations. Therefore, AAP-mediated Arg-specific regulation appeared independent of the charging status of arginyl-tRNA.
Collapse
Affiliation(s)
- Z Wang
- Department of Biochemistry, Oregon Graduate Institute of Science and Technology, Beaverton, Oregon 97006-8921, USA
| | | | | |
Collapse
|
6
|
Selker E, Yanofsky C. A phenylalanine tRNA gene from Neurospora crassa: conservation of secondary structure involving an intervening sequence. Nucleic Acids Res 1980; 8:1033-42. [PMID: 6449692 PMCID: PMC323971 DOI: 10.1093/nar/8.5.1033] [Citation(s) in RCA: 19] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
We have isolated and sequenced a tRNAPhe gene from Neurospora crassa. Hybridization analyses suggest that trnaPhe is the only tRNA encoded on the cloned 5 kb DNA fragment. The tRNAPhe gene contains an intervening sequence 16 nucleotides in length located one nucleotide 3' to the anticodon position. The tRNAPhe coding region of Neurospora and yeast are 91% conserved, whereas their intervening sequences are only 50% identical. The pattern of sequence conservation is consistent with a proposed secondary structure for the tRNA precursor in which the anticodon is base paired with the middle of the intervening sequence and the splice points are located in adjacent single-stranded loops. The DNA sequence following the tRNAPhe coding region is similar to sequences following other genes transcribed by RNA polymerase III in that it is AT-rich and includes a tract of A residues in the coding strand. In contrast, the sequence preceding the Neurospora tRNAPhe coding region does not resemble sequences preceding other sequenced tRNA genes.
Collapse
|