N 2-methylguanosine is iso-energetic with guanosine in RNA duplexes and GNRA tetraloops

Challenges with Simulating Modified RNA: Insights into Role and Reciprocity of Experimental and Computational Approaches

RNA is critical to a broad spectrum of biological and viral processes. This functional diversity is a result of their dynamic nature; the variety of three-dimensional structures that they can fold into; and a host of post-transcriptional chemical modifications. While there are many experimental techniques to study the structural dynamics of biomolecules, molecular dynamics simulations (MDS) play a significant role in complementing experimental data and providing mechanistic insights. The accuracy of the results obtained from MDS is determined by the underlying physical models i.e., the force-fields, that steer the simulations. Though RNA force-fields have received a lot of attention in the last decade, they still lag compared to their protein counterparts. The chemical diversity imparted by the RNA modifications adds another layer of complexity to an already challenging problem. Insight into the effect of RNA modifications upon RNA folding and dynamics is lacking due to the insufficiency or absence of relevant experimental data. This review provides an overview of the state of MDS of modified RNA, focusing on the challenges in parameterization of RNA modifications as well as insights into relevant reference experiments necessary for their calibration.

Calculations of pKa Values for a Series of Naturally Occurring Modified Nucleobases

Modified nucleobases are found in functionally important regions of RNA and are often responsible for essential structural roles. Many of these nucleobase modifications are dynamically regulated in nature, with each modification having a different biological role in RNA. Despite the high abundance of modifications, many of their characteristics are still poorly understood. One important property of a nucleobase is its pK_a value, which has been widely studied for unmodified nucleobases, but not for the modified versions. In this study, the pK_a values of modified nucleobases were determined by performing ab initio quantum mechanical calculations using a B3LYP density functional with the 6-31+G(d,p) basis set and a combination of implicit-explicit solvation systems. This method, which was previously employed to determine the pK_a values of unmodified nucleobases, is applicable to a variety of modified nucleobases. Comparisons of the pK_a values of modified nucleobases give insight into their structural and energetic impacts within nucleic acids.

Spliceosomal snRNA Epitranscriptomics

Small nuclear RNAs (snRNAs) are critical components of the spliceosome that catalyze the splicing of pre-mRNA. snRNAs are each complexed with many proteins to form RNA-protein complexes, termed as small nuclear ribonucleoproteins (snRNPs), in the cell nucleus. snRNPs participate in pre-mRNA splicing by recognizing the critical sequence elements present in the introns, thereby forming active spliceosomes. The recognition is achieved primarily by base-pairing interactions (or nucleotide-nucleotide contact) between snRNAs and pre-mRNA. Notably, snRNAs are extensively modified with different RNA modifications, which confer unique properties to the RNAs. Here, we review the current knowledge of the mechanisms and functions of snRNA modifications and their biological relevance in the splicing process.

Naturally occurring modified ribonucleosides

The chemical identity of RNA molecules beyond the four standard ribonucleosides has fascinated scientists since pseudouridine was characterized as the "fifth" ribonucleotide in 1951. Since then, the ever-increasing number and complexity of modified ribonucleosides have been found in viruses and throughout all three domains of life. Such modifications can be as simple as methylations, hydroxylations, or thiolations, complex as ring closures, glycosylations, acylations, or aminoacylations, or unusual as the incorporation of selenium. While initially found in transfer and ribosomal RNAs, modifications also exist in messenger RNAs and noncoding RNAs. Modifications have profound cellular outcomes at various levels, such as altering RNA structure or being essential for cell survival or organism viability. The aberrant presence or absence of RNA modifications can lead to human disease, ranging from cancer to various metabolic and developmental illnesses such as Hoyeraal-Hreidarsson syndrome, Bowen-Conradi syndrome, or Williams-Beuren syndrome. In this review article, we summarize the characterization of all 143 currently known modified ribonucleosides by describing their taxonomic distributions, the enzymes that generate the modifications, and any implications in cellular processes, RNA structure, and disease. We also highlight areas of active research, such as specific RNAs that contain a particular type of modification as well as methodologies used to identify novel RNA modifications. This article is categorized under: RNA Processing > RNA Editing and Modification.

RNA modifications and the link to human disease

Ribonucleic acid (RNA) is involved in translation and transcription, which are the mechanisms in which cells express genes (Alberts et al., 2002). The three classes of RNA discussed are transfer RNA (tRNA), messenger RNA (mRNA), and ribosomal RNA (rRNA). mRNA is the transcript encoded from DNA, rRNA is associated with ribosomes, and tRNA is associated with amino acids and is used to read mRNA transcripts to make proteins (Lodish, Berk, Zipursky, et al., 2000). Interestingly, the function of tRNA, rRNA, and mRNA can be significantly altered by chemical modifications at the co-transcriptional and post-transcriptional levels, and there are over 171 of these modifications identified thus far (Boccaletto et al., 2018; Modomics-Modified bases, 2017). Several of these modifications are linked to diseases such as cancer, diabetes, and neurological disorders. In this review, we will introduce a few RNA modifications with biological functions and how dysregulation of these RNA modifications is linked to human disease.

The chemical diversity of RNA modifications

Nucleic acid modifications in DNA and RNA ubiquitously exist among all the three kingdoms of life. This trait significantly broadens the genome diversity and works as an important means of gene transcription regulation. Although mammalian systems have limited types of DNA modifications, over 150 different RNA modification types have been identified, with a wide variety of chemical diversities. Most modifications occur on transfer RNA and ribosomal RNA, however many of the modifications also occur on other types of RNA species including mammalian mRNA and small nuclear RNA, where they are essential for many biological roles, including developmental processes and stem cell differentiation. These post-transcriptional modifications are enzymatically installed and removed in a site-specific manner by writer and eraser proteins respectively, while reader proteins can interpret modifications and transduce the signal for downstream functions. Dysregulation of mRNA modifications manifests as disease states, including multiple types of human cancer. In this review, we will introduce the chemical features and biological functions of these modifications in the coding and non-coding RNA species.

The life of U6 small nuclear RNA, from cradle to grave

Removal of introns from precursor messenger RNA (pre-mRNA) and some noncoding transcripts is an essential step in eukaryotic gene expression. In the nucleus, this process of RNA splicing is carried out by the spliceosome, a multi-megaDalton macromolecular machine whose core components are conserved from yeast to humans. In addition to many proteins, the spliceosome contains five uridine-rich small nuclear RNAs (snRNAs) that undergo an elaborate series of conformational changes to correctly recognize the splice sites and catalyze intron removal. Decades of biochemical and genetic data, along with recent cryo-EM structures, unequivocally demonstrate that U6 snRNA forms much of the catalytic core of the spliceosome and is highly dynamic, interacting with three snRNAs, the pre-mRNA substrate, and >25 protein partners throughout the splicing cycle. This review summarizes the current state of knowledge on how U6 snRNA is synthesized, modified, incorporated into snRNPs and spliceosomes, recycled, and degraded.

TRMT1-Catalyzed tRNA Modifications Are Required for Redox Homeostasis To Ensure Proper Cellular Proliferation and Oxidative Stress Survival

Mutations in the tRNA methyltransferase 1 (TRMT1) gene have been identified as the cause of certain forms of autosomal-recessive intellectual disability (ID). However, the molecular pathology underlying ID-associated TRMT1 mutations is unknown, since the biological role of the encoded TRMT1 protein remains to be determined. Here, we have elucidated the molecular targets and function of TRMT1 to uncover the cellular effects of ID-causing TRMT1 mutations. Using human cells that have been rendered deficient in TRMT1, we show that TRMT1 is responsible for catalyzing the dimethylguanosine (m2,2G) base modification in both nucleus- and mitochondrion-encoded tRNAs. TRMT1-deficient cells exhibit decreased proliferation rates, alterations in global protein synthesis, and perturbations in redox homeostasis, including increased endogenous ROS levels and hypersensitivity to oxidizing agents. Notably, ID-causing TRMT1 variants are unable to catalyze the formation of m2,2G due to defects in RNA binding and cannot rescue oxidative stress sensitivity. Our results uncover a biological role for TRMT1-catalyzed tRNA modification in redox metabolism and show that individuals with TRMT1-associated ID are likely to have major perturbations in cellular homeostasis due to the lack of m2,2G modifications.

Comparative Structural Dynamics of tRNA(Phe) with Respect to Hinge Region Methylated Guanosine: A Computational Approach

Transfer RNAs (tRNAs) contain various uniquely modified nucleosides thought to be useful for maintaining the structural stability of tRNAs. However, their significance for upholding the tRNA structure has not been investigated in detail at the atomic level. In this study, molecular dynamic simulations have been performed to assess the effects of methylated nucleic acid bases, N (2)-methylguanosine (m(2)G) and N (2)-N (2)-dimethylguanosine (m 2 (2) G) at position 26, i.e., the hinge region of E. coli tRNA(Phe) on its structure and dynamics. The results revealed that tRNA(Phe) having unmodified guanosine in the hinge region (G26) shows structural rearrangement in the core of the molecule, resulting in lack of base stacking interactions, U-turn feature of the anticodon loop, and TΨC loop. We show that in the presence of the unmodified guanosine, the overall fold of tRNA(Phe) is essentially not the same as that of m(2)G26 and m 2 (2) G26 containing tRNA(Phe). This structural rearrangement arises due to intrinsic factors associated with the weak hydrogen-bonding patterns observed in the base triples of the tRNA(Phe) molecule. The m(2)G26 and m 2 (2) G26 containing tRNA(Phe) retain proper three-dimensional fold through tertiary interactions. Single-point energy and molecular electrostatics potential calculation studies confirmed the structural significance of tRNAs containing m(2)G26 and m 2 (2) G26 compared to tRNA with normal G26, showing that the mono-methylated (m(2)G26) and dimethylated (m 2 (2) G26) modifications are required to provide structural stability not only in the hinge region but also in the other parts of tRNA(Phe). Thus, the present study allows us to better understand the effects of modified nucleosides and ionic environment on tRNA folding.

Intermolecular 'cross-torque': the N4-cytosine propargyl residue is rotated to the 'CH'-edge as a result of Watson-Crick interaction

Propargyl groups are attractive functional groups for labeling purposes, as they allow CuAAC-mediated bioconjugation. Their size minimally exceeds that of a methyl group, the latter being frequent in natural nucleotide modifications. To understand under which circumstances propargyl-containing oligodeoxynucleotides preserve base pairing, we focused on the exocyclic amine of cytidine. Residues attached to the exocyclic N4 may orient away from or toward the Watson-Crick face, ensuing dramatic alteration of base pairing properties. ROESY-NMR experiments suggest a uniform orientation toward the Watson-Crick face of N(4)-propargyl residues in derivatives of both deoxycytidine and 5-methyl-deoxycytidine. In oligodeoxynucleotides, however, UV-melting indicated that N(4)-propargyl-deoxycytidine undergoes standard base pairing. This implies a rotation of the propargyl moiety toward the 'CH'-edge as a result of base pairing on the Watson-Crick face. In oligonucleotides containing the corresponding 5-methyl-deoxycytidine derivative, dramatically reduced melting temperatures indicate impaired Watson-Crick base pairing. This was attributed to a steric clash of the propargyl moiety with the 5-methyl group, which prevents back rotation to the 'CH'-edge, consequently preventing Watson-Crick geometry. Our results emphasize the tendency of an opposing nucleic acid strand to mechanically rotate single N(4)-substituents to make way for Watson-Crick base pairing, providing no steric hindrance is present on the 'CH'-edge.

Microarrays for identifying binding sites and probing structure of RNAs

Oligonucleotide microarrays are widely used in various biological studies. In this review, application of oligonucleotide microarrays for identifying binding sites and probing structure of RNAs is described. Deep sequencing allows fast determination of DNA and RNA sequence. High-throughput methods for determination of secondary structures of RNAs have also been developed. Those methods, however, do not reveal binding sites for oligonucleotides. In contrast, microarrays directly determine binding sites while also providing structural insights. Microarray mapping can be used over a wide range of experimental conditions, including temperature, pH, various cations at different concentrations and the presence of other molecules. Moreover, it is possible to make universal microarrays suitable for investigations of many different RNAs, and readout of results is rapid. Thus, microarrays are used to provide insight into oligonucleotide sequences potentially able to interfere with biological function. Better understanding of structure-function relationships of RNA can be facilitated by using microarrays to find RNA regions capable to bind oligonucleotides. That information is extremely important to design optimal sequences for antisense oligonucleotides and siRNA because both bind to single-stranded regions of target RNAs.

Counterion and polythymidine loop-length-dependent folding and thermodynamic stability of DNA hairpins reveal the unusual counterion-dependent stability of tetraloop hairpins

Stem-loop DNA hairpins containing a 5-base-pair (bp) stem and single-stranded polythymidine loop were investigated using thermodynamic melting analysis and stopped-flow kinetics. These studies revealed the thermodynamic stability and folding kinetics as a function of loop length and counterion concentration. Our results show the unusually high thermodynamic stability for tetraloop or 4 poly(dT) loop hairpin as compared with longer loop length hairpins. Furthermore, this exceptional stability is highly counterion-dependent. For example, in the higher counterion concentration regime of 50 mM NaCl and above, the tetraloop hairpin displays enhanced stability as compared with longer loop length hairpins. However, at lower counterion concentration of 25 mM NaCl and below, the thermal stability of tetraloop hairpin is consistent with the longer loop hairpins. The enhanced stability of tetraloop hairpins at higher counterion concentration can be explained on the basis of the combined entropic effect of loop closure as well as base stacking in the loop regions. The stability of longer loop length hairpins at all counterion concentrations as well as tetraloop hairpin at lower counterion concentration can be explained on the basis of entropic effect of loop closure alone. The thermodynamic parameters at lower and higher counterion concentrations were determined to quantify the enhanced stability of base-stacking effects occurring at higher counterion concentrations. For example, for 100 mM NaCl, excess Gibbs energy and enthalpy due to base stacking within the tetraloops were measured to be -1.2 ± 0.14 and -3.28 ± 0.32 kcal/mol, respectively, whereas, no excess of Gibbs energy and enthalpy was observed for 0, 5, 10, and 25 mM NaCl. These findings suggest significant base-stacking interactions occurring in the loop region of the tetraloop hairpins at higher counterion concentration and less significant base-stacking interactions in the lower counterion concentration regime. We suggest that at higher counterion concentrations, hydrophobic collapse of the nucleotides in the loop may be enhanced due to the increased polarity of the solvent, thereby enhancing base-stacking interactions that contribute to unusually high stability.

MD simulation studies to investigate iso-energetic conformational behaviour of modified nucleosides m(2)G and m(2) 2G present in tRNA

Modified nucleic acid bases are most commonly found in tRNA. These may contain modifications from simple methylation to addition of bulky groups. Methylation of the four canonical nucleotide bases at a wide variety of positions is particularly prominent among the known modification. Methylation of N2 group of guanine is a relatively common modification in tRNA and rRNA. N2-methylguanosine (m²G) is the second most often encountered nucleoside in E. coli tRNAs. N2, N2- dimethylguanosine (m²₂G) is found in the majority of eukaryotic tRNAs and involved in forming base pair interactions with adjacent bases. Hence, in order to understand the structural significance of these methylated nucleic acid bases we have carried out molecular dynamics simulation to see the salvation effect. The results obtained shows iso-energetic conformational behaviors for m²G and m²₂G. The simulation trajectory of m²G shows regular periodical fluctuations suggesting that m²G is equally stable as either s-cis or s-trans rotamers. The two rotamers of m²G may interact canonically or non-canonically with opposite base as s-trans m²G26:C/A/U44 and s-cis m²G26:A/U44. The free rotations around the C-N bond could be the possible reason for these iso-energetic conformations. Dimethylation of G has almost no influence on base pairing with either A or U. Thus, these results reveal that modified nucleosides m²G and m²₂G may play an important role to prevent tRNA from adopting the unusual mitochondrial like conformation.

Conformational preferences of modified nucleoside N(2)-methylguanosine (m(2)G) and its derivative N(2), N(2)-dimethylguanosine (m(2)(2)G) occur at 26th position (hinge region) in tRNA

Conformational preferences of the modified nucleosides N(2)-methylguanosine (m(2)G) and N(2), N(2)-dimethylguanosine (m(2)(2)G) have been studied theoretically by using quantum chemical perturbative configuration interaction with localized orbitals (PCILO) method. Automated complete geometry optimization using semiempirical quantum chemical RM1, along with ab initio molecular orbital Hartree-Fock (HF-SCF), and density functional theory (DFT) calculations has also been made to compare the salient features. Single-point energy calculation studies have been made on various models of m(2)G26:C/A/U44 and m(2)(2)G26:C/A/U44. The glycosyl torsion angle prefers "syn" (χ = 286°) conformation for m(2)G and m(2)(2)G molecules. These conformations are stabilized by N(3)-HC2' and N(3)-HC3' by replacing weak interaction between O5'-HC(8). The N(2)-methyl substituent of (m(2)G26) prefers "proximal" or s-trans conformation. It may also prefer "distal" or s-cis conformation that allows base pairing with A/U44 instead of C at the hinge region. Thus, N(2)-methyl group of m(2)G may have energetically two stable s-trans m(2)G:C/A/U or s-cis m(2)G:A/U rotamers. This could be because of free rotations around C-N bond. Similarly, N(2), N(2)-dimethyl substituent of (m(2)(2)G) prefers "distal" conformation that may allow base pairing with A/U instead of C at 44th position. Such orientations of m(2)G and m(2)(2)G could play an important role in base-stacking interactions at the hinge region of tRNA during protein biosynthesis process.

Identification of a tertiary interaction important for cooperative ligand binding by the glycine riboswitch

The glycine riboswitch has a tandem dual aptamer configuration, where each aptamer is a separate ligand-binding domain, but the aptamers function together to bind glycine cooperatively. We sought to understand the molecular basis of glycine riboswitch cooperativity by comparing sites of tertiary contacts in a series of cooperative and noncooperative glycine riboswitch mutants using hydroxyl radical footprinting, in-line probing, and native gel-shift studies. The results illustrate the importance of a direct or indirect interaction between the P3b hairpin of aptamer 2 and the P1 helix of aptamer 1 in cooperative glycine binding. Furthermore, our data support a model in which glycine binding is sequential; where the binding of glycine to the second aptamer allows tertiary interactions to be made that facilitate binding of a second glycine molecule to the first aptamer. These results provide insight into cooperative ligand binding in RNA macromolecules.

Modification of 16S ribosomal RNA by the KsgA methyltransferase restructures the 30S subunit to optimize ribosome function

All organisms incorporate post-transcriptional modifications into ribosomal RNA, influencing ribosome assembly and function in ways that are poorly understood. The most highly conserved modification is the dimethylation of two adenosines near the 3' end of the small subunit rRNA. Lack of these methylations due to deficiency in the KsgA methyltransferase stimulates translational errors during both the initiation and elongation phases of protein synthesis and confers resistance to the antibiotic kasugamycin. Here, we present the X-ray crystal structure of the Thermus thermophilus 30S ribosomal subunit lacking these dimethylations. Our data indicate that the KsgA-directed methylations facilitate structural rearrangements in order to establish a functionally optimum subunit conformation during the final stages of ribosome assembly.

Bioinformatics analysis suggests base modifications of tRNAs and miRNAs in Arabidopsis thaliana

Background

Modifications of RNA bases have been found in some mRNAs and non-coding RNAs including rRNAs, tRNAs, and snRNAs, where modified bases are important for RNA function. Little is known about RNA base modifications in Arabidopsis thaliana.

Results

In the current work, we carried out a bioinformatics analysis of RNA base modifications in tRNAs and miRNAs using large numbers of cDNA sequences of small RNAs (sRNAs) generated with the 454 technology and the massively parallel signature sequencing (MPSS) method. We looked for sRNAs that map to the genome sequence with one-base mismatch (OMM), which indicate candidate modified nucleotides. We obtained 1,187 sites with possible RNA base modifications supported by both 454 and MPSS sequences. Seven hundred and three of these sites were within tRNA loci. Nucleotide substitutions were frequently located in the T arm (substitutions from A to U or G), upstream of the D arm (from G to C, U, or A), and downstream of the D arm (from G to U). The positions of major substitution sites corresponded with the following known RNA base modifications in tRNAs: N1-methyladenosine (m¹A), N2-methylguanosine (m²G), and N2-N2-methylguanosine (m²₂G).

Conclusion

These results indicate that our bioinformatics method successfully detected modified nucleotides in tRNAs. Using this method, we also found 147 substitution sites in miRNA loci. As with tRNAs, substitutions from A to U or G and from G to C, U, or A were common, suggesting that base modifications might be similar in tRNAs and miRNAs. We suggest that miRNAs contain modified bases and such modifications might be important for miRNA maturation and/or function.

Structural rearrangements in the active site of the Thermus thermophilus 16S rRNA methyltransferase KsgA in a binary complex with 5'-methylthioadenosine

Posttranscriptional modification of ribosomal RNA (rRNA) occurs in all kingdoms of life. The S-adenosyl-L-methionine-dependent methyltransferase KsgA introduces the most highly conserved rRNA modification, the dimethylation of A1518 and A1519 of 16S rRNA. Loss of this dimethylation confers resistance to the antibiotic kasugamycin. Here, we report biochemical studies and high-resolution crystal structures of KsgA from Thermus thermophilus. Methylation of 30S ribosomal subunits by T. thermophilus KsgA is more efficient at low concentrations of magnesium ions, suggesting that partially unfolded RNA is the preferred substrate. The overall structure is similar to that of other methyltransferases but contains an additional alpha-helix in a novel N-terminal extension. Comparison of the apoenzyme with complex structures with 5'-methylthioadenosine or adenosine bound in the cofactor-binding site reveals novel features when compared with related enzymes. Several mobile loop regions that restrict access to the cofactor-binding site are observed. In addition, the orientation of residues in the substrate-binding site indicates that conformational changes are required for binding two adjacent residues of the substrate rRNA.

Nucleotide analog interference mapping

Nucleotide analog interference mapping (NAIM) is a powerful chemogenetic technique that rapidly identifies chemical groups essential for RNA function. Using a series of phosphorothioate-tagged nucleotide analogs, each carrying different modifications of nucleobase or backbone functionalities, it is possible to simultaneously, yet individually, assess the contribution of particular functional groups to an RNA's activity at every position within the molecule. In contrast to traditional mutagenesis, which modifies RNA on the nucleobase level, the smallest mutable unit in a NAIM analysis is a single atom, providing a detailed description of interactions at critical nucleotides. Because the method introduces modified nucleotides by in vitro transcription, NAIM offers a straightforward and efficient approach to study any RNA that has a selectable function, and it can be applied to RNAs of nearly any length.

Effects of N2,N2-dimethylguanosine on RNA structure and stability: crystal structure of an RNA duplex with tandem m2 2G:A pairs

Methylation of the exocyclic amino group of guanine is a relatively common modification in rRNA and tRNA. Single methylation (N(2)-methylguanosine, m(2)G) is the second most frequently encountered nucleoside analog in Escherichia coli rRNAs. The most prominent case of dual methylation (N(2),N(2)-dimethylguanosine, m(2) (2)G) is found in the majority of eukaryotic tRNAs at base pair m(2) (2)G26:A44. The latter modification eliminates the ability of the N(2) function to donate in hydrogen bonds and alters its pairing behavior, notably vis-à-vis C. Perhaps a less obvious consequence of the N(2),N(2)-dimethyl modification is its role in controlling the pairing modes between G and A. We have determined the crystal structure of a 13-mer RNA duplex with central tandem m(2) (2)G:A pairs. In the structure both pairs adopt an imino-hydrogen bonded, pseudo-Watson-Crick conformation. Thus, the sheared conformation frequently seen in tandem G:A pairs is avoided due to a potential steric clash between an N(2)-methyl group and the major groove edge of A. Additionally, for a series of G:A containing self-complementary RNAs we investigated how methylation affects competitive hairpin versus duplex formation based on UV melting profile analysis.

Probing RNA structure and function by nucleotide analog interference mapping

Nucleotide analog interference mapping (NAIM) can be used to simultaneously, yet individually, identify structurally or catalytically important functional groups within an RNA molecule. Phosphorothioate-tagged nucleotides and nucleotide analogs are randomly incorporated into an RNA of interest by in vitro transcription. The phosphorothioate tag marks the site of substitution and identifies sites at which the modification affects the structure or function of the RNA molecule. This technique has been expanded to include identification of hydrogen bonding pairs (NAIS), ionizable functional groups, metal ion ligands, and the energetics of protein binding (QNAIM). The analogs, techniques, and data analysis used in NAIM are described here.

Combined Approaches to Site-Specific Modification of RNA

Both natural and unnatural modifications in RNA are of interest to biologists and chemists. More than 100 different analogues of the four standard RNA nucleosides have been identified in nature. Unnatural modifications are useful for structure and mechanistic studies of RNA. This Review highlights chemical, enzymatic, and combined (semisynthesis) approaches to generate site specifically modified RNAs. The availability of these methods for site-specific modifications of RNAs of all sizes is important in order to study the relationships between RNA chemical composition, structure, and function.

Chemical basis of glycine riboswitch cooperativity

The glycine binding riboswitch forms a unique tandem aptamer structure that binds glycine cooperatively. We employed nucleotide analog interference mapping (NAIM) and mutagenesis to explore the chemical basis of glycine riboswitch cooperativity. Based on the interference pattern, at least two sites appear to facilitate cooperative tertiary interactions, namely, the minor groove of the P1 helix from aptamer 1 and the major groove of the P3a helix from both aptamers. Mutation of these residues altered both the cooperativity and binding affinity of the riboswitch. The data support a model in which the P1 helix of the first aptamer participates in a tertiary interaction important for cooperativity, while nucleotides in the P1 helix of the second aptamer interface with the expression platform. These data have direct analogy to well-characterized mutations in hemoglobin, which provides a framework for considering cooperativity in this RNA-based system.

Identification and role of functionally important motifs in the 970 loop of Escherichia coli 16S ribosomal RNA

The 970 loop (helix 31) of Escherichia coli 16S ribosomal RNA contains two modified nucleotides, m(2)G966 and m(5)C967. Positions A964, A969, and C970 are conserved among the Bacteria, Archaea, and Eukarya. The nucleotides present at positions 965, 966, 967, 968, and 971, however, are only conserved and unique within each domain. All organisms contain a modified nucleoside at position 966, but the type of the modification is domain specific. Biochemical and structure studies have placed this loop near the P site and have shown it to be involved in the decoding process and in binding the antibiotic tetracycline. To identify the functional components of this ribosomal RNA hairpin, the eight nucleotides of the 970 loop of helix 31 were subjected to saturation mutagenesis and 107 unique functional mutants were isolated and analyzed. Nonrandom nucleotide distributions were observed at each mutated position among the functional isolates. Nucleotide identity at positions 966 and 969 significantly affects ribosome function. Ribosomes with single mutations of m(2)G966 or m(5)C967 produce more protein in vivo than do wild-type ribosomes. Overexpression of initiation factor 3 specifically restored wild-type levels of protein synthesis to the 966 and 967 mutants, suggesting that modification of these residues is important for initiation factor 3 binding and for the proper initiation of protein synthesis.

Expanding the nucleotide repertoire of the ribosome with post-transcriptional modifications

In all kingdoms of life, RNAs undergo specific post-transcriptional modifications. More than 100 different analogues of the four standard RNA nucleosides have been identified. Modifications in ribosomal RNAs (rRNAs) are highly prevalent and cluster in regions of the ribosome that have functional importance, have a high level of nucleotide conservation, and typically lack proteins. Modifications also play roles in determining antibiotic resistance or sensitivity. A wide spectrum of chemical diversity from the modifications provides the ribosome with a broader range of possible interactions between rRNA regions, transfer RNA, messenger RNA, proteins, or ligands by influencing local rRNA folds and fine-tuning the translation process. The collective importance of the modified nucleosides in ribosome function has been demonstrated for a number of organisms, and further studies may reveal how the individual players regulate these functions through synergistic or cooperative effects.

Identity elements required for enzymatic formation of N2,N2-dimethylguanosine from N2-monomethylated derivative and its possible role in avoiding alternative conformations in archaeal tRNA

Here, we have investigated the specificity of purified recombinant tRNA:m(2)(2)G10 methyltransferase of Pyrococcus abyssi ((Pab)Trm-m(2)(2)G10 enzyme). This archaeal enzyme catalyses mono- and dimethylation of the N(2)-exocyclic amino group of guanine at position 10 of several tRNA species. Our results indicate that only few identity elements are required for the efficient formation of m(2)(2)G10. They are composed of a G10.U25 wobble base-pair in the dihydrouridine arm (D-arm) and a four nucleotide variable loop (V-loop) within a canonical three-dimensional (3D) structure. The types of base-pairs in the D-arm or amino acid acceptor stem are also important for the enzymatic reaction, but appear to affect only the rate of tRNA methylation. However, in tRNA species harbouring a G10-C25 Watson-Crick base-pair and/or five nucleotide V-loop, only m(2)G10 is produced. To impair the monomethylation reaction, drastic amputation in the T-arm is required. Our observations contrast with those reported earlier for the identity elements required for a remotely related Pyrococcus furiosus Trm-m(2)(2)G26 enzyme (alias (Pfu)Trm1) that also catalyses the two step formation of m(2)(2)G but at position 26 in several tRNA species. In this case, a G10-C25 base-pair together with the five nucleotide V-loop were shown to be required for efficient formation of m(2)(2)G26. Thus, in the Pyrococcus genus, the major identity elements that preclude formation of m(2)(2)G at positions 10 or 26 in tRNA are mutually exclusive. Therefore, the Trm-m(2)(2)G10 and Trm-m(2)(2)G26 enzymes have evolved independently towards different specificities. In addition, identity elements for m(2)/m(2)(2)G10 formation in archaeal tRNA are different from the ones required for m(2)G10 formation in eukaryal tRNA. We propose that archaeal tRNA:m(2)(2)G10 methyltransferases, unlike the orthologous eukaryal tRNA:m(2)G10 methyltransferases, evolved towards m(2)(2)G10 specificity due to the possible requirement of preventing formation of alternative structures in G/C rich archaeal tRNA species.

N2-Methylation of Guanosine at Position 10 in tRNA Is Catalyzed by a THUMP Domain-containing, S-Adenosylmethionine-dependent Methyltransferase, Conserved in Archaea and Eukaryota

In sequenced genomes, genes belonging to the cluster of orthologous group COG1041 are exclusively, and almost ubiquitously, found in Eukaryota and Archaea but never in Bacteria. The corresponding gene products exhibit a characteristic Rossmann fold, S-adenosylmethionine-dependent methyltransferase domain in the C terminus and a predicted RNA-binding THUMP (thiouridine synthases, RNA methyltransferases, and pseudouridine synthases) domain in the N terminus. Recombinant PAB1283 protein from the archaeon Pyrococcus abyssi GE5, a member of COG1041, was purified and shown to behave as a monomeric 39-kDa entity. This protein (EC 2.1.1.32), now renamed (Pab)Trm-G10, which is extremely thermostable, forms a 1:1 complex with tRNA and catalyzes the adenosylmethionine-dependent methylation of the exocyclic amino group (N(2)) of guanosine located at position 10. Depending on the experimental conditions used, as well as the tRNA substrate tested, the enzymatic reaction leads to the formation of either N(2)-monomethyl (m(2)G) or N(2)-dimethylguanosine (m(2)(2)G). Interestingly, (Pab)Trm-G10 exhibits different domain organization and different catalytic site architecture from another, earlier characterized, tRNA-dimethyltransferase from Pyrococcus furiosus ((Pfu)Trm-G26, also known as (Pfu)Trm1, a member of COG1867) that catalyzes an identical two-step dimethylation of guanosine but at position 26 in tRNAs and is also conserved among all sequenced Eukaryota and Archaea. The co-occurrence of these two guanosine dimethyltransferases in both Archaea and Eukaryota but not in Bacteria is a hallmark of distinct tRNAs maturation strategies between these domains of life.

Effects of base substitutions in an RNA hairpin from molecular dynamics and free energy simulations

Contributions of individual interactions in the GGCGCAAGCC hairpin containing a GCAA tetraloop were studied by computer simulations using base substitutions. The G in the first tetraloop position was replaced by inosine (I) or adenosine (A), and the G in the C-G basepair closing the tetraloop was replaced by I. These substitutions eliminate particular hydrogen bonds proposed in the nuclear magnetic resonance model of the GCAA tetraloop. Molecular dynamics simulations of the GCAA tetraloop in aqueous solvent displayed a well-defined hydrogen pattern between the first and last loop nucleotides (G and A) stabilized by a bridging water molecule. Substitution of G-->I in the basepair closing the tetraloop did not significantly influence the loop structure and dynamics. The ICAA loop maintained the overall structure, but displayed variation in the hydrogen-bond network within the tetraloop itself. Molecular dynamics simulations of the ACAA loop led to conformational heterogeneity of the resulting structures. Changes of hairpin formation free energy associated with substitutions of individual bases were calculated by the free energy perturbation method. The calculated decrease of the hairpin stability upon G-->I substitution in the C-G basepair closing the tetraloop was in good agreement with experimental thermodynamic data. Our theoretical estimates for G-->I and G-->A mutations located in the tetraloop suggest larger loop destabilization than corresponding experimental results. The extent of conformational sampling of the structures resulting from base substitutions and its impact on the calculated free energy was discussed.

The functional anatomy of an intrinsic transcription terminator

To induce dissociation of the transcription elongation complex, a typical intrinsic terminator forms a G.C-rich hairpin structure upstream from a U-rich run of approximately eight nucleotides that define the transcript 3' end. Here, we have adapted the nucleotide analog interference mapping (NAIM) approach to identify the critical RNA atoms and functional groups of an intrinsic terminator during transcription with T7 RNA polymerase. The results show that discrete components within the lower half of the hairpin stem form transient termination-specific contacts with the RNA polymerase. Moreover, disruption of interactions with backbone components of the transcript region hybridized to the DNA template favors termination. Importantly, comparative NAIM of termination events occurring at consecutive positions revealed overlapping but distinct sets of functionally important residues. Altogether, the data identify a collection of RNA terminator components, interactions and spacing constraints that govern efficient transcript release. The results also suggest specific architectural rearrangements of the transcription complex that may participate in allosteric control of intrinsic transcription termination.

Synthesis of a 3-methyluridine phosphoramidite to investigate the role of methylation in a ribosomal RNA hairpin

The synthesis of a 5'-O-BzH-2'-O-ACE-protected-3-methyluridine phosphoramidite is reported [BzH, benzhydryloxy-bis(trimethylsilyloxy)silyl; ACE, bis(2-acetoxyethoxy)methyl]. The phosphoramidite was employed in solid-phase RNA synthesis to generate a series of RNA hairpins containing single or multiple modifications, including the common nucleoside pseudouridine. Three 19-nucleotide hairpin RNAs that represent the 1920-loop region (G(1906)-C(1924)) of Escherichia coli 23S ribosomal RNA were generated. Modifications were present at positions 1911, 1915, and 1917. The stabilities and structures of the three RNAs were examined by using thermal melting, circular dichroism, and NMR spectroscopy

Methylation of the nucleobases in RNA oligonucleotides mediates duplex-hairpin conversion

We have systematically investigated the duplex to hairpin conversion of oligoribonucleotides under the aspect of nucleobase methylation. The first part of our study refers to the self-complementary sequence rCGCGAAUUCGCGA, which forms a stable Watson-Crick base paired duplex under various buffer conditions. It is shown that this sequence is forced to adopt a hairpin conformation if one of the central 6 nt is replaced by the corresponding methylated nucleotide, such as 1-methylguanosine N(2),N(2)-dimethylguanosine, N(6),N(6)-dimethyladenosine (m(6)(2)A) or 3-methyluridine. On the other hand, the duplex structure is retained and even stabilized by replacement of a central nucleotide with N(2)-methylguanosine (m(2)G) or N(4)-methylcytidine. A borderline case is represented by N(6)-methyladenosine (m(6)A). Although generally a duplex-preserving modification, our data indicate that m(6)A in specific strand positions and at low strand concentrations is able to effectuate duplex-hairpin conversion. Our studies also include the ssu ribosomal helix 45 sequence motif, rGACCm(2)GGm(6)(2)Am(6)(2)AGGUC. In analogy, it is demonstrated that the tandem m(6)(2)A nucleobases of this oligoribonucleotide prevent duplex formation with complementary strands. Therefore, it can be concluded that nucleobase methylations at the Watson-Crick base pairing site provide the potential not only to modulate but to substantially affect RNA structure by formation of different secondary structure motifs.

Influence of N6-isopentenyladenosine (i(6)A) on thermal stability of RNA duplexes

The thermodynamic stability of self-complementary oligoribonucleotides containing N6-isopentenyladenosine (i(6)A) or N6-isopentanyladenosine (p(6)A) was determined. The base pairs i(6)A.U and p(6)A.U were placed in either an internal (separated and tandem) and a terminal position within the duplex, or unpaired i(6)A and p(6)A as a 3'-dangling ends. The thermal unfolding of the oligomers was determined by means of UV melting profiles and the thermodynamic parameters: enthalpy (DeltaH degrees ), entropy (DeltaS degrees) and free energy (DeltaG degrees (37)) as well as the melting temperature (T(m)) were calculated. Both modified nucleosides destabilized the duplexes, however, the effect depended on the position of the modified adenosine within the duplex. The similarity of the behavior of oligomers containing i(6)A and p(6)A suggests a negligible effect of the double bond on the thermal stability. The largest destabilization was observed when derivatives of adenosine were placed in an internal position. The effect of 3'-dangling ends suggests that the presence of the N6-isopentenyl- or N6-isopentanyl substitutent affects hydrogen bonding rather than stacking within duplex.

Translesional synthesis on a DNA template containing N2-methyl-2'-deoxyguanosine catalyzed by the Klenow fragment of Escherichia coli DNA polymerase I

Formaldehyde is produced in most living systems and is present in the environment. Evidence that formaldehyde causes cancer in experimental animals infers that it may be a carcinogenic hazard to humans. Formaldehyde reacts with the exocyclic amino group of deoxyguanosine, resulting in the formation of N2-methyl-2'-deoxyguanosine (N2-Me-dG) via reduction of the Schiff base. The same reaction is likely to occur in living cells, because cells contain endogenous reductants such as ascorbic acid and gluthathione. To explore the miscoding properties of formaldehyde-derived DNA adducts a site-specifically modified oligodeoxynucleotide containing a N2-Me-dG was prepared and used as the template in primer extension reactions catalyzed by the Klenow fragment of Escherichia coli DNA polymerase I. The primer extension reaction was slightly stalled one base before the N2-Me-dG lesion, but DNA synthesis past this lesion was readily completed. The fully extended products were analyzed to quantify the miscoding specificities of N2-Me-dG. Preferential incorporation of dCMP, the correct base, opposite the lesion was observed, along with small amounts of misincorporation of dTMP (9.4%). No deletions were detected. Steady-state kinetic studies indicated that the frequency of nucleotide insertion for dTMP was only 1.2 times lower than for dCMP and the frequency of chain extension from the 3'-terminus of a dT:N2-Me-dG pair was only 2.1 times lower than from a dC:N2-Me-dG pair. We conclude that N2-Me-dG is a miscoding lesion capable of generating G-->A transition mutations.

Biologically active oligodeoxyribonucleotides. Part 12: N2-methylation of 2'-deoxyguanosines enhances stability of parallel G-quadruplex and anti-HIV-1 activity

2'-Deoxyguanosine residues of a 3',5'-end-modified hexadeoxyribonucleotide (R-95288) with anti-HIV-1 activity were substituted with N2-methyl-2'-deoxyguanosine (m2dG). These modified oligodeoxyribonucleotides (ODNs) showed a 2-fold higher activity than R-95288. Also, the CD spectra of these ODNs indicated that the m2dG modification stabilized the tertiary structure of the G-quadruplex.

A chemical phylogeny of group I introns based upon interference mapping of a bacterial ribozyme

Despite its small size, the 205 nt group I intron from Azoarcus tRNA(Ile) is an exceptionally stable self-splicing RNA. This IC3 class intron retains the conserved secondary structural elements common to group I ribozymes, but lacks several peripheral helices. These features make it an ideal system to establish the conserved chemical basis of group I intron activity. We collected nucleotide analog interference mapping (NAIM) data of the Azoarcus intron using 14 analogs that modified the phosphate backbone, the ribose sugar, or the purine base functional groups. In conjunction with a complete interference set collected on the Tetrahymena group I intron (IC1 class), these data define a "chemical phylogeny" of functional groups that are important for the activity of both introns and that may be common chemical features of group I intron catalysts. The data identify the functional moieties most likely to play a conserved role as ligands for catalytic metal ions, the substrate helix, and the guanosine cofactor. These include backbone functional groups whose nucleotide identity is not conserved, and hence are difficult to identify by standard phylogenetic sequence comparisons. The data suggest that both introns utilize an equivalent set of long range tertiary interactions for 5'-splice site selection between the P1 substrate helix and its receptor in the J4/5 asymmetric bulge, as well as an equivalent set of 2'-OH groups for P1 helix docking into most of the single stranded segment J8/7. However, the Azoarcus intron appears to make an alternative set of interactions at the base of the P1 helix and at the 5'-end of the J8/7. Extensive differences were observed within the intron peripheral domains, particularly in P2 and P8 where the Azoarcus data strongly support the proposed formation of a tetraloop-tetraloop receptor interaction. This chemical phylogeny for group I intron catalysis helps to refine structural models of the RNA active site and identifies functional groups that should be carefully investigated for their role in transition state stabilization.

Ribozyme chemogenetics

In this review I will outline several chemogenetic approaches used to determine the chemical basis of large ribozyme function and structure. The term chemogenetics was first used to describe site-specific functional group modification experiments in the analysis of DNA-protein interactions. Within the past few years equivalent experiments have been performed on large catalytic RNAs using both single-site substitution and interference mapping techniques with nucleotide analogues. While functional group mutagenesis is an important aspect of a chemogenetic approach, chemical correlates to genetic revertants and suppressors must also be realized for the genetic analogy to be intellectually valid and experimentally useful. Several examples of functional group revertants and suppressors have now been obtained within the Tetrahymena group I ribozyme. These experiments define an ensemble of tertiary hydrogen bonds that have made it possible to construct a detailed model of the ribozyme catalytic core. The model includes a functionally important monovalent metal ion binding site, a wobble-wobble receptor motif for helix-helix packing interactions, and a minor groove triple helix.

A chemogenetic approach to RNA function/structure analysis

Almost two dozen nucleotide analogs have been synthesized with alpha-phosphorothioate-tagged triphosphates and utilized in an interference modification approach termed Nucleotide Analog Interference Mapping. This method has made it possible to determine the chemical basis of RNA function and structure, including the identification of new rules for RNA helix packing, the functional analysis of a binding site for monovalent metal ions within RNA and the characterization of the catalytic mechanism of RNA enzymes.

Nucleotide analog interference mapping

Single-atom substitution experiments provide atomic resolution biochemical information concerning RNA structure and function. Traditionally, these experiments are performed using chimeric RNAs generated by reassembly of full-length RNA from a synthetic substituted oligonucleotide and a truncated RNA transcript. Unfortunately, this technique is limited by the technical difficulty of assembling and measuring the effect of each singly substituted molecule in a given RNA. Here we review an alternate method for rapidly screening the effect of chemical group substitutions on RNA function. Nucleotide analog interference mapping is a chemogenetic approach that utilizes a series 5'-O-(1-thio)-nucleoside analog triphosphates to simultaneously, yet individually, probe the contribution of a functional group at every nucleotide position in an RNA molecule. A population of randomly substituted RNAs is prepared by including phosphorothioate-tagged nucleotide analogs in an in vitro transcription reaction. The active molecules in the RNA population are selected by an activity assay, and the location of the analog substitution detrimental to activity is identified by cleavage at the phosphorothioate tag with iodine and resolution of the cleavage fragments by gel electrophoresis. This method, which is as easy as RNA sequencing, is applicable to any RNA that can be transcribed in vitro and has an assayable function. Here we describe protocols for the synthesis of phosphorothioate-tagged analogs and their incorporation into RNA transcripts. The incorporation properties and unique biochemical signatures of each individual analog are discussed.