key: cord-0888945-ulra9loe
authors: Zheng, Ya Ying; Wu, Ying; Begley, Thomas J.; Sheng, Jia
title: Sulfur modification in natural RNA and therapeutic oligonucleotides
date: 2021-04-27
journal: RSC chemical biology
DOI: 10.1039/d1cb00038a
sha: c07fb55dc63ae1e4b7506c9ed1c5dd1235ae8e5c
doc_id: 888945
cord_uid: ulra9loe

Sulfur modifications have been discovered on both DNA and RNA. Sulfur substitution of oxygen atoms at nucleobase or backbone locations in the nucleic acid framework led to a wide variety of sulfur-modified nucleosides and nucleotides. While the discovery, regulation and functions of DNA phosphorothioate (PS) modification, where one of the non-bridging oxygen atoms is replaced by sulfur on the DNA backbone, are important topics, this review focuses on the sulfur modification in natural cellular RNAs and therapeutic nucleic acids. The sulfur modifications on RNAs exhibit diversity in terms of modification location and cellular function, but the various sulfur modifications share common biosynthetic strategies across RNA species, cell types and domains of life. The first section reviews the post-transcriptional sulfur modifications on nucleobases with an emphasis on thiouridine on tRNA and phosphorothioate modification on RNA backbones, as well as the functions of the sulfur modifications on different species of cellular RNAs. The second section reviews the biosynthesis of different types of sulfur modifications and summarizes the general strategy for the biosynthesis of sulfur-containing RNA residues. One of the main goals of investigating sulfur modifications is to aid the genomic drug development pipeline and enhance our understandings of the rapidly growing nucleic acid-based gene therapies. The last section of the review focuses on the current drug development strategies employing sulfur substitution of oxygen atoms in therapeutic RNAs.

Post-transcriptional modifications in cellular RNAs have gained significant attentions over the past few decades. To date, scientists have discovered over 170 chemical modifications in rRNA, tRNA and mRNA in all domains of life. 1 The existence of these naturally occurring modifications added a new dimension beyond the basic building blocks of nucleic acids adenosine (A), uridine (U), cytidine (C) and guanosine (G). The presence of modified nucleosides on the base or ribose or backbone in specific sequence contexts confers unique features that allow for additional layers of gene expression and regulation. [1] [2] [3] Although modifications can provide RNAs with expanded structural diversity, a majority of modifications can affect RNA functionality when working with their cellular targets. Such extraordinary chemical alterations have been observed playing roles in many biological processes such as transcription, translation, transportation, stability and protection of nucleic acids. These chemical modifications are catalytic products of specific writer enzymes that feature different reaction processes such as methylation, demethylation, amination, deamination, isomerization and thiolation, etc. 4 The natural presence of these modifications affords more subtle functions of both DNA and RNA as evidenced by tRNA species with highly evolutionary conserved modifications at specific positions found in all the three domains of life. 5 The development and advancement of new mass spectrometry-based methods serve as the dominant strategies to analyze and quantify the RNA epitranscriptome. 6, 7 Over 170 modifications across RNA species including ribosomal RNA (rRNA), messenger RNA (mRNA), transfer RNA (tRNA), microRNA (miRNA) and other non-coding RNA have been identified. However, the specific roles and the impacts of some modifications in RNA structures and functions have only been appreciated and thoroughly elucidated very recently. 8 Overall speaking, the mapping and characterization for the majority of these modifications remain a big challenge considering the difficulty of studying the effects of a single modification at molecular and cellular levels.

Sulfur based modifications are among the rapidly growing modifications being reported, with many important roles. Sulfur accounts for one of the most abundant and essential elements in our body and is found in the amino acids methionine and cysteine. 9 Sulfur also plays important roles in the structure and regulation of proteins and coordinates many biomolecules in the cells via the formation of covalent disulfide bonds between cysteine residues in the proteins. Disulfide bond formation sometimes governs the enzymatic activity in many proteins, which points to the regulatory function of sulfur in cysteine. 10 In addition to cysteine and methionine, sulfur is also found in glutathione, the vitamins biotin and thiamine, the cofactor S-adenosyl-Met, coenzyme A, the molybdenum cofactor (MoCo) and many other secondary metabolic products. 10, 11 Glutathione is a tripeptide containing cysteine and is known to detoxify metabolites and other stressors, making it an essential aspect of cellular homeostasis. 12 These crucial sulfur-containing biomolecules and their biochemical functions render sulfur a good natural choice for nucleic acid modification. Several posttranscriptionally installed sulfur modifications in transfer ribonucleic acid (tRNA) have been discovered to maintain structural stability, promote effective codon recognition and enhance translation fidelity. 13 tRNAs are by far the most modified nucleic acid species. More specifically, over 60 chemical modifications are found on uridine (U) in tRNA, of which 16 are thiolated at the C2 position of U. 1, 14 Nucleoside thiolation includes 2-thiouridine (s 2 U), 4-thiouridine (s 4 U), 2-thiocytidine (s 2 C), 2-thioribothymidine (s 2 T), 2-methythioadenosine (ms 2 A), and other derivatives such as 2-selenouridine (se 2 U) and 2-geranyl thiouridine (ges 2 U). 5, [15] [16] [17] The s 2 and s 4 thiolated nucleobases are derived by replacing the keto-oxygen with sulfur. 18 Although the structural role of thiomodifications should promote functions in RNA species across all three domains of life, here we highlight the sulfur decorations on tRNA from well-studied organisms, 19, 20 and also integrate the newly discovered backbone phosphorothioate modification (Table 1) 23 Modification to the wobble 34th position in tRNA, particularly s 2 U, is critical for achieving higher accuracy and efficiency in protein translation. Moreover, modifications located around the anticodon in positions 32 and 37 also have roles in improving reading frame maintenance. 5, 24, 25 The thiolated nucleosides in the wobble position have also shown to be important as identity elements for certain aminoacyl-tRNA synthetase, in addition to ensuring correct codon binding in the ribosome. 23, 25, 26 Sulfur containing amino acids have also been shown to regulate translational capacity and metabolic homeostatic through modulation of tRNA thiolation. In yeast, the availability of sulfur containing amino acids methionine and cysteine can 

have direct effects on the thiolation status of wobble uridine U34. When those amino acids are limiting the sulfur source, s 2 U levels are downregulated, which serves as a cue to increase biosynthesis enzymes in methionine and cysteine production. Thus, tRNA thiolation is essential for achieving metabolic homeostasis and cell growth. 27 Mutations in the s 2 U biosynthetic gene encoding the 2-thiouridylase TrmU are linked to acute infantile liver failure, 28 and respiratory defects alongside with nonsyndromic deafness. 29 Depletion of s 2 U caused by point mutations in mitochondrial tRNA genes can lead to a variety of human mitochondrial pathologies, including mitochondrial myopathy, encephalopathy, lactic acidosis, and stroke-like episodes (MELAS) and myoclonus epilepsy associated with ragged-red fibers (MERRF). [30] [31] [32] [33] Modifications occurring within the acceptor stem, the TcC-loop and the D-loop serve to govern structural stability, enhance recognition by aminoacyl-tRNA synthetases and fine-tune the overall performance. 13, 34 The s 4 U modification is a thionucleoside usually found at the position 8 in tRNA within the intersection of the acceptor stem and D-loop. This modification is responsible for sensing near-UV light in bacteria contributing to its unique l max at 330 nm. When the cell is challenged with near-UV light, s 4 U crosslinks with cytidine 13, resulting in the destruction of the tRNA tertiary structure to cease translation. 35 In some extreme thermophilic bacteria such as Thermus thermophilus and Pyrococcus furiosus, 5-methyl-2thioridine (m 5 s 2 U) is found at position 54 on the TcC loop. The presence of the 2-thiolation content could strengthen the tRNA structure at elevated temperatures. [36] [37] [38] It is worth mentioning that both thiolated nucleobases s 2 U and s 4 U have impacts on RNA structure, which includes base pairing stability, specificity and an alteration of thermodynamic parameters. As determined by UV-thermal melting and imino-proton NMR experiments, the order of duplex stability for a pentamer RNA sequence containing either s 2 U, s 4 U or U and 2 0 -O-methyl-ribonucleotide to a complementary strand supported that s 2 U is the most stable while s 4 U is the least stable when pairing with A. 39 When comparing the base pairing affinity, s 2 U can stabilize base pairing with A more than G, with vice versa for s 4 U as it can stabilize base pairing with G more than A. 40 Further studies also demonstrated that s 2 U containing oligonucleotides can stabilize both A and U base pairing in RNA duplexes. The thermodynamic stabilizing effect is possibly due to the tendency of s 2 U to pre-orient the RNA strand to single-stranded and facilitate the hybridization with its complementary target. In addition, the enhancement of the base stacking and highly polarizable nature of the sulfur atom could also contribute the stability enhancement. 14 S 2 U has also shown to favor the C3 0 -endo conformation, which substantially stabilizes the nucleotides that likewise impart duplex thermostability. 39, 41 Geranyl-2-thiouridine (ges 2 U) and seleno-2-thiorudine (se 2 U) derivatives on tRNA Wobble tRNA modifications can significantly alter anticodoncodon interactions, which points to their regulatory potential to speed up or slow down translation and affect fidelity. 34, 41 The 2-thiolation of U34 can increase the codon recognition efficiency of tRNA Glu UUC , 42 and plays a critical role in the ribosome binding of tRNA Lys UUU . 43, 44 In addition, the s 2 U can be further modified, forming a s 2 U based modification network that will be described in the next section. Particularly, the sulfur atom in s 2 U can be further replaced with selenium to form 2-selenouridine (se 2 U) by the tRNA modifying enzyme MnmH (also called SelU) using selenophosphate as the cofactor. MnmH has also been shown to use geranyl pyrophosphate as the cofactor to form geranylated 2-thiouridine (ges 2 U) in tRNAs specific for lysine (tRNA Lys UUU ), glutamine (tRNA Gln UUG ), and glutamic acid (tRNA Glu UUC ). Although the mechanisms of this dual function of MnmH remain elusive, MnmH has been known to contain a rhodanese domain, which is the key active working domain for both geranylation and selenation processes, 45, 46 and is ubiquitous protein that catalyzes the transfer of a sulfur atom from thiosulfate (rhodanese) or 3-mercaptopyruvate to e.g. cyanide in vitro. 47, 48 In addition, MnmH also contains a P-loop domain with a Walk-A motif, which is present in many ATPor GTP-binding proteins and also involved in substrate binding. 49, 50 It has also been suggested that ges 2 U is the intermediate product in the se 2 U formation process. 46, 51 The ges 2 U modification is lipid like with a 10-carbon hydrophobic chain on the wobble U. The ges 2 U modification can be further modified to 5-methylaminomethyl-2-gerenythioluridine (mnm 5 ges 2 U) or 5carboxymethylaminomethyl-2-geranylthiouridine (cmnm 5 ges 2 U). 17, 52 This special lipid modification has been found in several bacteria including E. coli, E. aerogenes, P. aeruginosa, and S. typhimurium at a frequency of up to 6.7% in all the U34 residues of the three tRNAs (B400 geranylated nucleotides per cell). 17 Notably, the conjugation of the geranyl group has not been found in any eukaryotic cells, and the lipid modification is distinct in the bacterial kingdom. It has been known that both 2-thiouridine (s 2 U) and 2-selenouridine (se 2 U) enhance the discrimination of U:A pair over U:G pair ( Fig. 1A-D) . 14, 53 In contrast, the geranylated tRNA Glu UUC was observed to increase the codon recognition efficiency to GAG over GAA, 17 indicating that the geranylated uridine has a stronger base pairing preference with G over A at the end of codons. Indeed, by analyzing potential base pairing patterns ( Fig. 1E-H) , one direct consequence of the geranyl group is to switch the N3 of uridine from a hydrogen bond donor to an acceptor, which induces different base pairing patterns. As a result, only G that has two connected hydrogen bonding donors can form relatively stable two hydrogen bonds with the geranylated uridine. Based on this hydrogen bonding pattern and the fact that each of the two codons for glutamic acid, lysine, and glutamine ends in either A or G, it is speculated that geranylated uridine might be involved in translational regulation, by promoting base pairing with the G-ending codon for each amino acid, while restricting pairing with the A-ending codons.

Post-transcriptional modifications of cellular RNAs are indispensable for cellular development and exert influence in

RSC Chemical Biology biological processes in both normal and disease cells. 54, 55 The modifications on RNAs known earlier are either on nitrogenous base or ribose. The discovery of phosphorothioate (PS) modification on DNA has been documented in bacteria as well as in archaea. 56, 57 In a phosphorothioate internucleotide linkage, one of the non-bridging oxygen atoms is replaced with the sulfur, this renders the stereogenic phosphorus center. In both bacteria and archaea, the DNA phosphorothioate modifications have been elucidated to be in Rp configuration, with occurrence frequency of 4-31 PS modifications per 10 4 deoxy-nucleotides. 58 More recently, the discovery of phosphorothioate modification in RNA was also reported in both prokaryotes and eukaryotes. 59 The method of detection involved synthesizing RNA oligonucleotides containing PS modifications and subjecting it to nuclease enzymatic digestion to determine the optimal kinetics in which the normal phosphodiester bonds are hydrolyzed while the phosphorothioate bonds are intact. Different samples with prokaryotes and eukaryotes total RNAs were extracted and digested. The resulting mono-nucleosides and phosphorothioate dinucleotides together with 16 synthetic RNA PS dinucleotide standards were subjected to ultra-performance liquid chromatography coupled with mass spectrometry (UHPLC-MS) for analysis of RNA PS linkages. Quantification analysis estimated the number of phosphorothioate dinucleotide (PS) modification per 10 000 nucleotides in total RNAs of HeLa, mouse liver and DM S2 as 4.7, 7.6, and 6.7 respectively. Additionally, CpsA, GpsC and GpsG are also found in the total RNA of L. lactis. The GpsG occurs at the highest frequency among all the samples tested and exists in stereo-specific Rp configuration. Further experiment showed that this GpsG modification locates on the rRNA in E. coli, L. lactis, and HeLa cells, and it is not detected in the rRNA-depleted total RNAs from these cells. Despite this big leap in the road of RNA modification, however, the cellular functions of RNA PS modification are still under investigation.

In this section, we summarize the biosynthesis of sulfur modifications located on nucleobases emphasising on 2-thiouridine, 4-thioruidine and geranyl derivatives on tRNA. Sulfur-modified nucleosides in eukaryotes and bacteria all utilize free L-cysteine in the cell as primary sulfur source. In general, the biosynthetic pathways for sulfur modifications in tRNA begin with the cysteine desulfurase protein, which activates L-cysteine residue to generate persulfide-SSH molecule. Enzyme linked persulfide is then delivered to downstream sulfur acceptor proteins and eventually incorporated into tRNA. 60, 61 Biosynthesis of 2-Thiouridine (s 2 U)

The 2-thiouridine modification at the wobble position 34 of the anticodon stem loop in the tRNAs for glutamate acid, glutamine, and lysine is critical for enhancing translation due to its functional role in reducing frameshift and promoting codon recognition during protein synthesis. 5 The biosynthesis of these thionucleosides has been well elucidated in eukaryote such as Saccharomyces cerevisiae and bacteria as in E. coli. The synthetic pathway can be categorized into two groups depending on the need for iron-sulfur [Fe-S] clusters. The overall pathway is summarized in Fig. 2 .

In Gram-negative bacteria such as E. coli, the biosynthesis of 2-thiouridine requires two conserved enzymes, IscS and MnmA along with other enzymes: TusA, TusBCD heterohexameric complex and TusE enzymes, which serve as sulfur-relay system for the flow of persulfide sulfur in the pathway. IscS, the cysteine desulfurase, liberates sulfur from L-cysteine to form IscS-SSH enzyme bound persulfide intermediate. The binding of IscS with TusA stimulates the transfer of sulfur persulfide to TusA. The small sulfur-carrier protein TusA then transfers the sulfur to TusD aided by the interaction of TusE with the In Gram-positive bacterium such as Bacillus subtilis, the biosynthesis of s 2 U is distinct from that of E. coli. Bacillus subtilis does not contain Tus proteins and cysteine desulfurase IscS. However, genomic analysis showed its genome interestingly encodes an active IscS-like cysteine desulfurase YrvO. It is worth noting that YrvO coding sequence is 31 bases upstream of MnmA gene and in the same genomic region as IscS gene. In vitro and in vivo studies showed that both of the YrvO and MnmA proteins are the sole requirement for s 2 U synthesis in Bacillus subtilis. In this two-component pathway, YrvO transfers persulfide sulfur to the cysteine residue within MnmA at the expense of ATP. Although the exact mechanism awaits further experimental investigation, the function of YrvO and MnmA in s 2 U tRNA formation and the lack of Tus sulfur relay proteincsystem suggested the presence of a devoted cysteine desulfurase that can perform direct sulfur transfer from cysteine to MnmA. 64 Bacteria contains at least two biosynthesis mechanisms of s 2 U tRNA: the 5-methoxycarbonylmethyl-2-thioruidne (mnm 5 s 2 U34) in cytoplasm and the 5-carboxymethylaminomethyl-2-thiouridine (cmnm 5 s 2 U34) in mitochondria. The 5-mnm and 5-cmnm groups are introduced through the MnmCDEG pathways 65 and are also important for efficient codon recognition. 42, 66 The overall 2-thiouridine centered modification pathways are illustrated in Fig. 3 with several key enzymes still missed in this network. The differences in the biosynthesis pathways between bacteria and yeast is due to the involvement of [Fe-S] clusters and the usage of protein-thiocarboxylate as sulfur donor. 67, 68 The mechanistic pathway for the 2-thiolation of cy-tRNAs in S. cerevisiae begins with cysteine desulfurase Nfs1 catalyzing the first step, which is the transfer of persulfide to the rhodanese domain (Rhd) on Uba4 via Rhd of Tum1. Rhodanese serves as a sulfur carrier biomolecule and catalyzes the transfer reaction to downstream sulfur acceptor enzymes. Ubiquitin-related modifier 1 or Urm1 acts as a sulfur carrier while acyl-adenylate forms Urm1-COAMP intermediate and then thiocarboxylate Urm1-COSH at the C-terminus by receiving persulfide sulfur from the Rhd of Uba4. The sulfur from activated thiocarboxylate is then transferred to tRNA via Ncs6/Ncs2 heterodimer complex which catalyzes the final step of 2-thiolation U34 cy-tRNAs. 5 

In thermophilic bacterium such as T. thermophilus, the biosynthesis of m 5 s 2 U at position 54 in tRNA involves the TtuA/TtuB pathway. Cysteine desulfurase IscS or SufS begins the persulfide sulfur relay to the subsequent enzyme TtuB, a ubiquitin like sulfur carrier. TtuC activates TtuB at the C-terminal of Gly to form acyl-adenylated intermediate while being thiocarboxylated by accepting the sulfur from IscS/SufS. The final transfer of sulfur to tRNA in thiouridine formation is mediated by TtuA. 5

4-Thiouridine at position 8 of tRNA serves as a photo-sensor for near-ultraviolet radiation found in bacteria. In Gram-negative bacterium such as E. coli, the synthesis of 4-thio modification required both IscS and ThiI enzymes. Earlier studies have reported that both enzymes also play roles in the synthesis of thiamine. [73] [74] [75] IscS is a pyridoxal 5 0 -phosphate (PLP) dependent cysteine desulfurase that converts L-cysteine to L-alanine by catalyzing the cleavage of the carbon-sulfur bond in cysteine while liberating the sulfur atom to generate IscS linked persulfide intermediate IscS-SSH at its active site Cys-328. 76, 77 The subsequent step is the transfer of this activated persulfide sulfur to the acceptor molecule ThiI on Cys-456. ThiI contains the rhodanese-like sulfurtransferase domain which serves as the sulfur acceptor site. This persulfide on ThiI adenylates on U8 tRNA and donates sulfur for the synthesis of s 4 U. 78 Characterization of ThiI reveals that this enzyme is common in the biosynthetic pathway for both thiamin and 4-thiouridine in tRNA. The overall synthesis pathway of s 4 U is shown in Fig. 4 .

The E. coli ThiI contains four domains: the C-terminal rhodanese like domain (RLD), the PP-loop pyrophosphatase domain (PP-loop domain), the N-terminal ferredoxin-like domain (NFLD), and the thiouridine synthases, RNA methyltransferase and pseudouridine synthases domain (THUMP domain). 79 The NFLD and THUMP domains are involved in tRNA binding at the acceptor stem region. The proposed mechanism for the transfer of persulfide on the ThiI to form 4-thiouridine in tRNA involves the adenylation of 4-carbonyl group at the U8 tRNA in the PP-loop domain of ThiI utilizing the Mg-ATP. This activates U8 for the incorporation of sulfur atom. Sulfur is delivered from IscS to the Cys-456 in ThiI, forming ThiI persulfide. The transfer of terminal sulfur atom of persulfide on ThiI to U8 tRNA occurs in the rhodanese like domain (Rhd). Two mechanisms have been proposed for sulfur transfer mediated by ThiI. 80 The nucleophilic attacking of persulfide sulfur on Cys-456 of ThiI active site to the C4 of uridine-8 in tRNA could initiate this process in concomitant with the release of AMP. In addition, this mechanism involves the enzymatic disulfide bond formation between Cys-456 and Cys-344 located in the PP-loop domain and generates s 4 U tRNA as the product. The assembly of disulfide bond serves to regenerate the enzyme for subsequent cycle which requires the reductive cleavage of the disulfide bond by a reductase through either a thiol group or another cysteine residue of ThiI. The regenerated cysteine active site is ready for another sulfur Fig. 3 Chemical illustration of the biosynthetic pathway of s 2 U derivatives on both position 5 and 2. Compound 1 is uridine (U) with no substitution at the N5 position. When N5 is substituted with different moieties of X to produce compound 1a: 5-carboxymethylaminomethyluridine (cmnm5U), 1b: 5aminomethyluridine (nm5U) both through enzyme MnmE, G, and 1c: 5-methylaminomethyluridine (mnm5U) through enzyme MnmC,D. These products can be further modified at the N2 position by the enzyme MnmA to generate compound 2 derivatives containing 2-thiouridine, with the products depending upon the X moieties. For example, U will be converted to s2U, cmnm5U to cmnm5s 2 U, nm5U to nm5s 2 U and mnm5U to mnm5s 2 U. Moreover, enzyme MnmH will further modified compound 2 derivatives eventually to compound 4 via common geranyl intermediate compound 3. The final products of compound 4 derivatives are Se 2 U via ges 2 U, cmnm5se 2 U via cmnm5ges 2 U, nm5se 2 U via nm5ges 2 U and mnm5se 2 U via mnm5ges 2 U intermediates. It is worth noting that enzyme MnmE, G can also use s2U and se 2 U as substrates and MnmC,D can also use nm5s 2 U and nm5se 2 U as substrates.

transfer from IscS and the cycle starts again. 79 Alternatively, Cys-344 can also attack the bridging sulfur of Cys-456 followed by the release of sulfide that can feed directly into tRNA for the synthesis of s 4 U. The key feature of both mechanisms is the formation of disulfide bond between Cys-456 and Cys-344 as mark for the completion of the catalytic cycle. 80, 81 The s 4 U synthesis pathway in Gram-positive bacteria and archaea diverges from those of Gram-negative bacteria as in E. coli, although the exact mechanism is still unknown. Experimental data have revealed that majority of Grampositive bacteria such as Bacillus subtilis do not contain IscS protein but actually encode another cysteine desulfurase enzyme NifZ. The sequence of NifZ exists four bases directly upstream of ThiI. The short sequence of ThiI produces protein that lacks the rhodanese domain but retains NFLD, THUMP and PP-loop domains. PP-loop domain reserves the Cys-344, the equivalent cysteine residue position in E. coli ThiI that is important for adenylation and thiolation of U8 tRNA. 82 Despite lacking the rhodanese domain, B. subtilis ThiI and NifZ are sufficient to generate thiolation in tRNA as ThiI can directly receive persulfide sulfur from NifZ. In E. coli, Cys-344 involves in the final sulfur transferase step accompanied by the formation of disulfide bond with Cys-456 in the rhodanese domain. The missing of this important Cys-456 counterpart suggests the possibility of other unidentified ThiI cysteine for completing the thiolation. 82 Similar to Gram-positive bacteria, the mechanism of s 4 U thiolation needs further elucidation since majority of archaea also lack cysteine desulfurase IscS and the rhodanese (Rhd) domain in ThiI protein. Study in Methanococcus maripaludis reveals that ThiI contains three conserved Cys residues (two come from a CXXC motif) in the putative catalytic region. A single mutation study showed that all three Cys residues are critical for persulfide transferase and replacing any Cys with Ala abolishes the function of ThiI as persulfide carrier. ThiI also contains [3Fe-4S] clusters which are essential for thiolation functionality in tRNA. 83 Biosynthesis of RNA phosphorothioate (PS) backbone modification RNA phosphorothioate (PS) modification was recently reported in both prokaryotes and eukaryotes. 59 However, its corresponding regulatory enzymes and biological functions remain elusive. It has been demonstrated that the PS modification in RNA shares some common features with DNA counterpart such 

RSC Chemical Biology as the stereospecificity and relative abundance of the linked dinucleotides. For example, the GpsG, the most frequent RNA PS dinucleotide found so far, is only detected in the Rp stereo form. The abundance of RNA PS modification is also comparable to that of DNA PS modification found in bacteria. Therefore it is possible that this modification is regulated by the same set of enzymes or enzymes of similar functions. The DNA phosphorothioate backbone modification with one of the non-bridging oxygens replaced by the sulfur atom is the product of the DndA-E gene clusters together with DndFGH genes which constitute a defense mechanism that resemble restriction-modification system and protect bacteria from foreign invaders. 84 The DndA-E gene clusters that regulate DNA PS modifications in both bacteria and archaeal are conserved in all domains of life. DndA actually shares significant nucleotide sequence identity and protein structural similarity with the pyridoxal 5-phosphate (PLP) dependent cysteine desulfurases IscS and NifS in eukaryotes. 59, 84 In addition, DndA was also reported as the sulfur source provider in the biosynthesis of DNA PS modification in bacteria. To test the possibility of DndA being involved in the biosynthesis of RNA PS modification as a writer, this gene was knocked out in E. coli BUN21/pML300 strain and the total RNAs were then extracted from both wild type and knock out strains. Quantification of GpsG with LC-MS/MS shows the modification level is reduced by half in DndA-knocked out strain comparing to the wild type BUN21. This result supports that DndA plays important roles in RNA PS biosynthesis pathway. Moreover, bacteria strains lacking the DndA gene were previously reported to still have detected DNA PS modification, which implies that there might be other devoted cysteine desulfurase such as IscS that serves as alternative sulfur source in this backbone thiolation process. Similar in DNA PS modification, there might be other enzymes in the Dnd gene clusters worthy of further investigation, probably guided by bioinformatic studies. 59 

Advances in biotechnology and genomics have led to rapid advances in nucleic acid based therapeutics, all of which have great potential in clinical applications for the treatment of cancer, infectious diseases, cardiovascular and neurological disorders. Currently, nucleic acids based therapeutic strategies such as antisense oligonucleotides (ASOs), RNA interference (RNAi) or small interfering RNA (siRNA) and aptamers all have been demonstrated to have promising results for biomedical applications since they are able to alter gene expression, thereby preventing or alleviating disease development. 85 These synthetic nucleic acids are designed to modulate gene expression with high specificity via target sequence recognition and hybridization. Antisense oligonucleotides (ASOs) can hybridize target mRNA by Watson-Crick base pairing and silence the gene expression by employing one of the three antisense mechanisms (Fig. 5) : steric inhibition of mRNA translation; mRNA degradation caused by RNase H cleaving the RNA strand of the DNA-RNA duplex; 86, 87 and induction of splicing switch by exon skipping which restores the production of a protein. 88 siRNA mediates mRNA degradation in the 

Review process of RNA interference (RNAi). Ago2, a RNase H like endonuclease can cleave targeted mRNA that is complementary to the guide RNA in siRNA bound to Ago2. 89, 90 Aptamers are single stranded nucleic acid that resemble monoclonal antibodies by self-assembling into 3D structure with the capability of binding to a variety of targets ranging from small molecules to proteins to whole cells with high affinity and specificity. 91 It is also noteworthy the potential therapeutic applications of CRISPR-Cas9 technologies in genome editing. 92 The clustered regularly interspaced short palindromic repeat (CRISPR) with CRISPR associated protein together constitute adaptive immunity possessed by prokaryotes, as a defense mechanism against foreign DNA or RNA-DNA interference. The single guide RNA (sgRNA) is versatile, its sequence can be programmed to guide Cas protein to cause precise gene editing at any target DNA loci. 93 To improve in vivo activity and specificity, the engineered sgRNA with chemical modification can be obtained by synthesis. 94 One of the main restrictions or challenges in the development of nucleic acid based drugs is the nuclease degradation before reaching the target sites. To solve this problem, different chemical modifications have been used to increase the biostability of these oligonucleotides such as base methylation, 2 0 -substitutions like 2 0 -fluoro (2 0 -F), 2 0 -methoxy (2 0 -OMe), 2 0 -O-methoxyethyl (2 0 -MOE), locked nucleic acid (LNA), phosphorodiamidate morpholino (PMO), N-acetylgalactosamine (GalNAc) and phosphorothioate (PS), etc. The incorporation of these modifications have proven to drastically stabilize the synthetic agents. 95 In the case of CRISPR-Cas9 system, sgRNA with PS along or in conjunction with 2 0 -modifications have reported to reduce off target and significantly enhanced the editing ability by improving target binding of Cas9 protein. 94, 96 Table 2 shows the FDA approved oligo therapies. [97] [98] [99] [100] [101] [102] [103] [104] [105] [106] These modified RNA-based drugs have shown resistance toward nuclease degradation, improving target selectivity, increasing binding affinity, as well as enhancing metabolic uptake. In particular, nusinersen, mipomersen, inotersen and volanesorsen are all modified with 2 0 -methoxy and phosphorothioate backbone linkages and fomivirsen, a DNA based ASO, have all PS modifications. 106 PS modifications, where a sulfur atom replaces one of the non-bridging oxygens on the oligonucleotide strand, are one of the most important and widely used strategies to improve the drug-like properties of nucleic acid based therapeutics to enhance both cellular uptake and biological stability towards nuclease degradation. 107 This characteristic is essential for the intended therapeutic effect as the molecules will remain in circulation until reaching to the target location for gene silencing. 108 Out of a series of chemical modifications either on the backbone, nucleobase, or ribose sugar, the internucleotide PS linkage exhibits most resistance to RNase H cleavage. 109 The chemical differences between the phosphorothioate and phosphate group arise from a larger van der Waals radius of 1.85 Å of sulfur atom versus 1.44 Å of oxygen, which resulted in the lengthened P-S bond. Moreover, phosphorothioate occurs as the mixture of stereo-isomers (Rp and Sp) since its phosphate center is connected with four different ligands. This chirality renders the isomers with different binding affinity and interaction modes with enzymes. In general, the Rp stereo-isomers have a stronger binding to the target mRNA and form a more stable complex with higher melting temperatures. 110 However, Rp configuration has shown to be more susceptible to RNase H degradation compared to the Sp linkage that has shown to possess a higher resistance to nuclease cleavage activity. 111 For example, Sp isomer confers greater stability for 3 0 exonuclease by RNase A and snake venom phosphodiesterase than Rp isomer. 112, 113 Therefore, stereo-control synthesis of oligonucleotides must be attainable to examine the fundamental properties of individual isomers. Instead of using the regular nucleoside phosphoramidites for solid phase synthesis, using the bicyclic oxazaphospholidine derivatives as monomer building block has demonstrated the success of obtaining a high quality of stereo-pure Rp or Sp diastereomers. The 3 0 -O-bicyclic oxazaphospholidine derivatives were generated by the reaction between the 3 0 -OH of the corresponding protected nucleosides and 2-chloro-1,3,2-oxazaphospholidine derivatives, and the resultant trans-oxazaphospholidine isomers were configurationally very stable without any loss of diastereopurity even under acidic condition. 114 

It is also worth mentioning the thio-containing methyl-7guanosine triphosphate cap at the 5 0 -end of mRNA (5 0 m7G) connected via a 5 0 ,5 0 -triphosphate linkage with the first 

nucleotide, which is an evolutionarily conserved across eukaryotic species. This modification confers unique functions such as protection against degradation, pre-mRNA processing, nuclear exportation and modulating protein synthesis. 115 In cancer immunotherapy, phosphorothioate modified cap analogs have been shown to profoundly increase the biostability and translational efficiency of RNA vaccines in immature dendritic cells and induce superior immune responses in vivo. 116 In another study where the mRNA bears the phosphorothioate modification at either the a, b, or g position of the 5 0 ,5 0 -triphosphate chain, it has been demonstrated that the thio-modified cap generally binds tighter to eIF4E, the eukaryotic translation initiation factor 4E, and are more resistant to the decapping pyrophosphatase DcpS. Particularly, the b-substituted analog m 7 GPPsPG has the strongest binding affinity toward eIF4E, and g-substituted analog shows strong resistant to hydrolysis by DcpS. 117 The 5 0 -m 7 G cap analog containing 6-thioguanosine can be applied to study nucleic acid and cap dependent protein interaction through photoinducible crosslinking experiments. 2 0 -O-methylation within the m 7 G moiety ensures correct cap incorporation during mRNA synthesis, which is known as anti-reverse cap analogs (ARCAs). 118 The 6-thioguanosine serves as a photo-activable nucleoside which absorbs higher wavelength light compared with the natural nucleic acids and proteins. This allows for identification of the binding site between nucleic acid and protein when the system is selectively activated with certain wavelength. 119 Cap analogs modified with 1,2-dithiodiodiphosphate containing a sulfur substitute at the two neighboring phosphate moieties have shown to protect mRNA from decapping and enhance its overall translational efficiency. 120 The thio-modified cap has also been tried in formulating COVID-19 mRNA vaccines to increase their biostability and overall efficiency, which represents a very useful strategy for future vaccine development. 121 

Overall, the biosynthetic pathways of nucleosides containing sulfur are quite complicated since these processes usually involve different sulfur-transporting enzymes rather than a direct single reaction. In addition, these pathways could vary significantly in different species and organisms. More detailed functions and working mechanisms in these biological reactions are still largely unknown. In general, the unique functions of sulfur modifications are based on the biochemical properties of sulfur-containing motifs in the nucleic acid structures. Sulfur atoms are significantly less electronegative (EN = 2.58) than oxygen (EN = 3.44), leading to changes in noncovalent interactions such as hydrogen bonding properties when oxygen atoms are replaced by sulfurs. These reduced or enhanced non-covalent interactions rendered by sulfur substitution alter a wide range of cellular RNAs' functions from base pairing specificity to three dimensional structures and RNA stability toward nuclease degradation. The sulfur residues can also respond to cellular stress and environmental stimuli to induce a cascade of biological processes. The phosphorothioate bond containing a chiral center with two stereo configurations does not alter the RNA structure but substantially reduce their susceptibility towards nuclease degradation. Advancement in chemical synthesis and the huge progress in understanding genetic mutations and molecular drivers of diseases have led to rapid advancement of modified oligonucleotide-based therapeutics and provided solutions to solve the RNA degradation issues. As we have seen in COVID-19 vaccines, mRNA holds great promise in future drug development and the thiomodifications including the thio-caps on mRNA will play important roles in achieving more effective therapeutics for targeted diseases.

There are no conflicts to declare.

Fine-tuning of RNA functions by modification and editing

The arabidopsis book

Drug Des

We are grateful to NSF (CHE-1845486 and MCB-1715234) and the NERF grant from the University at Albany, State University of New York, for the financial support.