Structural effects of modified ribonucleotides and magnesium in transfer RNAs

Modified nucleotides are ubiquitous and important to tRNA structure and function. To understand their effect on tRNA conformation, we performed a series of molecular dynamics simulations on yeast tRNA^Phe and tRNA_init, Escherichia coli tRNA_init and HIV tRNA^Lys. Simulations were performed with the wild type modified nucleotides, using the recently developed CHARMM compatible force field parameter set for modified nucleotides (J. Comput. Chem. 2016, 37, 896), or with the corresponding unmodified nucleotides, and in the presence or absence of Mg²⁺. Results showed a stabilizing effect associated with the presence of the modifications and Mg²⁺ for some important positions, such as modified guanosine in position 37 and dihydrouridines in 16/17 including both structural properties and base interactions. Some other modifications were also found to make subtle contributions to the structural properties of local domains. While we were not able to investigate the effect of adenosine 37 in tRNA_init and limitations were observed in the conformation of E. coli tRNA_init, the presence of the modified nucleotides and of Mg²⁺ better maintained the structural features and base interactions of the tRNA systems than in their absence indicating the utility of incorporating the modified nucleotides in simulations of tRNA and other RNAs.     

Solution conformations of unmodified and A(37)N(6)-dimethylallyl modified anticodon stem-loops of Escherichia coli tRNA(Phe)

The modification of RNA nucleotide bases, a fundamental process in all cells, alters the chemical and physical properties of RNA molecules and broadly impacts the physiological properties of cells. tRNA molecules are by far the most diverse-modified RNA species within cells, containing as a group >80% of the known 96 chemically unique nucleic acid modifications. The greatest varieties of modifications are located on residue 37 and play a role in ensuring fidelity and efficiency of protein synthesis. The enzyme dimethylallyl (Delta(2)-isopentenyl) diphosphate:tRNA transferase catalyzes the addition of a dimethylallyl group to the exocyclic amine nitrogen (N6) of A(37) in several tRNA species. Using a 17 residue oligoribonucleotide corresponding to the anticodon arm of Escherichia coli tRNA(Phe), we have investigated the structural and dynamic changes introduced by the dimethylallyl group. The unmodified RNA molecule adopts stem-loop conformation composed of seven base-pairs and a compact three nucleotide loop. This conformation is distinctly different from the U-turn motif that characterizes the anticodon arm in the X-ray crystal structure of the fully modified yeast tRNA(Phe). The adoption of the tri-nucleotide loop by the purine-rich unmodified tRNA(Phe) anticodon arm suggests that other anticodon sequences, especially those containing pyrimidine bases, also may favor a tri-loop conformation. Introduction of the dimethylallyl modification increases the mobility of nucleotides of the loop region but does not dramatically alter the RNA conformation. The dimethylallyl modification may enhance ribosome binding through multiple mechanisms including destabilization of the closed anticodon loop and stabilization of the codon-anticodon helix.     

Role of modified nucleosides of yeast tRNAPhe in ribosomal binding

Naturally occurring nucleoside modifications are an intrinsic feature of transfer RNA (tRNA), and have been implicated in the efficiency, as well as accuracy-of codon recognition. The structural and functional contributions of the modified nucleosides in the yeast tRNA(Phe) anticodon domain were examined. Modified nucleosides were site-selectively incorporated, individually and in combinations, into the heptadecamer anticodon stem and loop domain, (ASL(Phe)). The stem modification, 5-methylcytidine, improved RNA thermal stability, but had a deleterious effect on ribosomal binding. In contrast, the loop modification, 1-methylguanosine, enhanced ribosome binding, but dramatically decreased thermal stability. With multiple modifications present, the global ASL stability was mostly the result of the individual contributions to the stem plus that to the loop. The effect of modification on ribosomal binding was not predictable from thermodynamic contributions or location in the stem or loop. With 4/5 modifications in the ASL, ribosomal binding was comparable to that of the unmodified ASL. Therefore, modifications of the yeast tRNA(Phe) anticodon domain may have more to do with accuracy of codon reading than with affinity of this tRNA for the ribosomal P-site. In addition, we have used the approach of site-selective incorporation of specific nucleoside modifications to identify 2'O-methylation of guanosine at wobble position 34 (Gm34) as being responsible for the characteristically enhanced chemical reactivity of C1400 in Escherichia coli 16S rRNA upon ribosomal footprinting of yeast tRNA(Phe). Thus, effective ribosome binding of tRNA(Phe) is a combination of anticodon stem stability and the correct architecture and dynamics of the anticodon loop. Correct tRNA binding to the ribosomal P-site probably includes interaction of Gm34 with 16S rRNA C1400.     

Modified nucleotides in tRNA(Lys) and tRNA(Val) are important for translocation

Ribosomes translate genetic information encoded by messenger RNAs (mRNAs) into proteins. Accurate decoding by the ribosome depends on the proper interaction between the mRNA codon and the anticodon of transfer RNA (tRNA). tRNAs from all kingdoms of life are enzymatically modified at distinct sites, particularly in and near the anticodon. Yet, the role of these naturally occurring tRNA modifications in translation is not fully understood. Here we show that modified nucleosides at the first, or wobble, position of the anticodon and 3'-adjacent to the anticodon are important for translocation of tRNA from the ribosome's aminoacyl site (A site) to the peptidyl site (P site). Thus, naturally occurring modifications in tRNA contribute functional groups and conformational dynamics that are critical for accurate decoding of mRNA and for translocation to the P site during protein synthesis.     

28 The study of tRNA modifications by molecular dynamics

Modified nucleosides are prevalent in tRNA. Experimental studies reveal that modifications play an important role in tuning tRNA activity. In this study, molecular dynamics (MD) simulations were used to investigate how modifications alter tRNA structure and dynamics. The X-ray crystal structures of tRNA-Asp, tRNA-Phe, and tRNA-iMet, both with and without modifications, were used as initial structures for 333-ns time-scale MD trajectories with AMBER. For each tRNA molecule, three independent trajectory calculations were performed. Force field parameters were built using the RESP procedure of Cieplak et al. for 17 nonstandard tRNA residues. The global root-mean-square deviations (RMSDs) of atomic positions show that modifications only introduce significant rigidity to tRNA-Phe’s global structure. Interestingly, regional RMSDs of anticodon stem-loop suggest that modified tRNA has more rigid structure compared to the unmodified tRNA in this domain. The anticodon RMSDs of the modified tRNAs, however, are higher than those of corresponding unmodified tRNAs. These findings suggest that rigidity of the anticodon arm is essential for tRNA translocation in the ribosome complex, and, on the other hand, flexibility of anticodon might be critical for anticodon–codon recognition. We also measure the angle between the 3D L-shaped arms of tRNA; backbone atoms of acceptor stem and TψC stem loop are selected to indicate one vector, and backbone atoms of anticodon stem and D stem loop are selected to indicate the other vector. By measuring the angle between two vectors, we find that the initiator tRNA has a narrower range of hinge motion compared to tRNA-Asp and tRNA-Phe, which are elongator tRNA. This suggests that elongator tRNAs, which might require significant flexibility in this hinge to transition from the A–to-P site in the ribosome, have evolved to specifically accommodate this need.     

The conformations of locked nucleic acids (LNA)

We have used 2D NMR spectroscopy to study the sugar conformations of oligonucleotides containing a conformationally restricted nucleotide (LNA) with a 2'-O, 4'-C-methylene bridge. We have investigated a modified 9-mer single stranded oligonucleotide as well as three 9- and 10-mer modified oligonucleotides hybridized to unmodified DNA. The single-stranded LNA contained three modifications whereas the duplexes contained one, three and four modifications, respectively. The LNA:DNA duplexes have normal Watson-Crick base-pairing with all the nucleotides in anti-conformation. By use of selective DQF-COSY spectra we determined the ratio between the N-type (C3'-endo) and S-type (C2'-endo) sugar conformations of the nucleotides. In contrast to the corresponding single-stranded DNA (ssDNA), we found that the sugar conformations of the single-stranded LNA oligonucleotide (ssLNA) cannot be described by a major S-type conformer of all the nucleotides. The nucleotides flanking an LNA nucleotide have sugar conformations with a significant population of the N-type conformer. Similarly, the sugar conformations of the nucleotides in the LNA:DNA duplexes flanking a modification were also shown to have significant contributions from the N-type conformation. In all cases, the sugar conformations of the nucleotides in the complementary DNA strand in the duplex remain in the S-type conformation. We found that the locked conformation of the LNA nucleotides both in ssLNA and in the duplexes organize the phosphate backbone in such a way as to introduce higher population of the N-type conformation. These conformational changes are associated with an improved stacking of the nucleobases. Based on the results reported herein, we propose that the exceptional stability of the LNA modified duplexes is caused by a quenching of concerted local backbone motions (preorganization) by the LNA nucleotides in ssLNA so as to decrease the entropy loss on duplex formation combined with a more efficient stacking of the nucleobases.     

Phylogenetic comparative chemical footprint analysis of the interaction between ribonuclease P RNA and tRNA.

Ribonuclease P RNA is the catalytic moiety of the ribonucleoprotein enzyme that endonucleolytically cleaves precursor sequences from the 5' ends of pre-tRNAs. The bacterial RNase P RNA-tRNA complex was examined with a footprinting approach, utilizing chemical modification to determine RNase P RNA nucleotides that potentially contact tRNA. RNase P RNA was modified with dimethylsulfate or kethoxal in the presence or absence of tRNA, and sites of modification were detected by primer extension. Comparison of the results reveals RNase P bases that are protected from modification upon binding tRNA. Analyses were carried out with RNase P RNAs from three different bacteria: Escherichia coli, Chromatium vinosum and Bacillus subtilis. Discrete bases of these RNAs that lie within conserved, homologous portions of the secondary structures are similarly protected. One protection among all three RNAs was attributed to the precursor segment of pre-tRNA. Experiments using pre-tRNAs containing precursor segments of variable length demonstrate that a precursor segment of only 2-4 nucleotides is sufficient to confer this protection. Deletion of the 3'-terminal CCA sequence of tRNA correlates with loss of protection of a particular loop in the RNase P RNA secondary structure. Analysis of mutant tRNAs containing sequential 3'-terminal deletions suggests a relative orientation of the bound tRNA CCA to that loop.     

Rp-phosphorothioate modifications in RNase P RNA that interfere with tRNA binding.

We have used Rp-phosphorothioate modifications and a binding interference assay to analyse the role of phosphate oxygens in tRNA recognition by Escherichia coli ribonuclease P (RNase P) RNA. Total (100%) Rp-phosphorothioate modification at A, C or G positions of RNase P RNA strongly impaired tRNA binding and pre-tRNA processing, while effects were less pronounced at U positions. Partially modified E. coli RNase P RNAs were separated into tRNA binding and non-binding fractions by gel retardation. Rp-phosphorothioate modifications that interfered with tRNA binding were found 5' of nucleotides A67, G68, U69, C70, C71, G72, A130, A132, A248, A249, G300, A317, A330, A352, C353 and C354. Manganese rescue at positions U69, C70, A130 and A132 identified, for the first time, sites of direct metal ion coordination in RNase P RNA. Most sites of interference are at strongly conserved nucleotides and nine reside within a long-range base-pairing interaction present in all known RNase P RNAs. In contrast to RNase P RNA, 100% Rp-phosphorothioate substitutions in tRNA showed only moderate effects on binding to RNase P RNAs from E. coli, Bacillus subtilis and Chromatium vinosum, suggesting that pro-Rp phosphate oxygens of mature tRNA contribute relatively little to the formation of the tRNA-RNase P RNA complex.     

A systematic, ligation-based approach to study RNA modifications

Over 100 different chemical types of modifications have been identified in thousands of sites in tRNAs, rRNAs, mRNAs, small nuclear RNAs, and other RNAs. Some modifications are highly conserved, while others are more specialized. They include methylation of bases and the ribose backbone, rotation, and reduction of uridine, base deamination, elaborate addition of ring structures, carbohydrate moieties, and more. We have developed a systematic approach to detect and quantify the extent of known RNA modifications. The method is based on the enzymatic ligation of oligonucleotides using the modified or unmodified RNA as the template. The efficiency of ligation is very sensitive to the presence and the type of modifications. First, two oligo pairs for each type of modification are identified. One pair greatly prefers ligation using the unmodified RNA template over the modified RNA template or vice versa. The other pair has equal reactivity with unmodified and modified RNA. Second, separate ligations with each of the two oligo pairs and the total RNA mixture are performed to detect the presence or absence of modifications. Multiple modification sites can be examined in the same ligation reaction. The feasibility of this method is demonstrated for three 2'O-methyl modification sites in yeast rRNA.     

Conformation effects of base modification on the anticodon stem-loop of Bacillus subtilis tRNA(Tyr)

tRNA molecules contain 93 chemically unique nucleotide base modifications that expand the chemical and biophysical diversity of RNA and contribute to the overall fitness of the cell. Nucleotide modifications of tRNA confer fidelity and efficiency to translation and are important in tRNA-dependent RNA-mediated regulatory processes. The three-dimensional structure of the anticodon is crucial to tRNA-mRNA specificity, and the diverse modifications of nucleotide bases in the anticodon region modulate this specificity. We have determined the solution structures and thermodynamic properties of Bacillus subtilis tRNA^Tyr anticodon arms containing the natural base modifications N⁶-dimethylallyl adenine (i⁶A₃₇) and pseudouridine (ψ₃₉). UV melting and differential scanning calorimetry indicate that the modifications stabilize the stem and may enhance base stacking in the loop. The i⁶A₃₇ modification disrupts the hydrogen bond network of the unmodified anticodon loop including a C₃₂-A₃₈⁺ base pair and an A₃₇-U₃₃ base-base interaction. Although the i⁶A₃₇ modification increases the dynamic nature of the loop nucleotides, metal ion coordination reestablishes conformational homogeneity. Interestingly, the i⁶A₃₇ modification and Mg²⁺ are sufficient to promote the U-turn fold of the anticodon loop of Escherichia coli tRNA^Phe, but these elements do not result in this signature feature of the anticodon loop in tRNA^Tyr.