Transfer RNA structure and coding specificity. II. A D-arm tertiary interaction that restricts coding range

We investigated the structural basis of the kinetic effect on coding specificity by the D-arm mutant (G24 to A) of Escherichia coli tRNATrp. A set of tRNA genes with structural alterations in the D-arm was constructed by site-directed mutagenesis in vitro, and we determined the in vivo translational activities of these tRNAs. Our results suggest that a hydrogen-bond donor in the major groove of the D-helix at position 24 is required for the expansion of tRNA wobble coding specificity. From inspection of tRNA crystal structure, we identified a potential new tertiary pairing of base 24 with the base at position 9 (this base links the acceptor and D-stems). We constructed tRNAs with mutations at position 9 and showed that the phenotypes of position 11-24 D-arm mutants are indeed dependent on the identity of base 9. Our analysis of the effects of these mutations on the interactions of tRNA with the ribosome and with aminoacyl-tRNA synthetase suggests that the conformation or conformational dynamics of the middle of the tRNA molecule alters the kinetics of the interaction with the ribosomal coding site. The 9-23 and putative 9-24 tertiaries, and perhaps other normal tertiary interactions in this region, modulate these kinetics to increase or decrease coding specificity.     

The structure of yeast tRNA(Asp). A model for tRNA interacting with messenger RNA

The anticodon of yeast tRNA(Asp), GUC, presents the peculiarity to be self-complementary, with a slight mismatch at the uridine position. In the orthorhombic crystal lattice, tRNA(Asp) molecules are associated by anticodon-anticodon interactions through a two-fold symmetry axis. The anticodon triplets of symmetrically related molecules are base paired and stacked in a normal helical conformation. A stacking interaction between the anticodon loops of two two-fold related tRNA molecules also exists in the orthorhombic form of yeast tRNA(Phe). In that case however the GAA anticodon cannot be base paired. Two characteristic differences can be correlated with the anticodon-anticodon association: the distribution of temperature factors as determined from the X-ray crystallographic refinements and the interaction between T and D loops. In tRNA(Asp) T and D loops present higher temperature factors than the anticodon loop, in marked contrast to the situation in tRNA(Phe). This variation is a consequence of the anticodon-anticodon base pairing which rigidifies the anticodon loop and stem. A transfer of flexibility to the corner of the tRNA molecule disrupts the G19-C56 tertiary interactions. Chemical mapping of the N3 position of cytosine 56 and analysis of self-splitting patterns of tRNA(Asp) substantiate such a correlation.     

tRNA discrimination by T. thermophilus phenylalanyl-tRNA synthetase at the binding step

The extent of tRNA recognition at the level of binding by Thermus thermophilus phenylalanyl-tRNA synthetase (PheRS), one of the most complex class II synthetases, has been studied by independent measurements of the enzyme association with wild-type and mutant tRNA(Phe)s as well as with non-cognate tRNAs. The data obtained, combined with kinetic data on aminoacylation, clearly show that PheRS exhibits more tRNA selectivity at the level of binding than at the level of catalysis. The anticodon nucleotides involved in base-specific interactions with the enzyme prevail both in the initial binding recognition and in favouring aminoacylation catalysis. Tertiary nucleotides of base pair G19-C56 and base triple U45-G10-C25 contribute primarily to stabilization of the correctly folded tRNA(Phe) structure, which is important for binding. Other nucleotides of the central core (U20, U16 and of the A26-G44 tertiary base pair) are involved in conformational adjustment of the tRNA upon its interaction with the enzyme. The specificity of nucleotide A73, mutation of which slightly reduces the catalytic rate of aminoacylation, is not displayed at the binding step. A few backbone-mediated contacts of PheRS with the acceptor and anticodon stems revealed in the crystal structure do not contribute to tRNA(Phe) discrimination, their role being limited to stabilization of the complex. The highest affinity of T. thermophilus PheRS for cognate tRNA, observed for synthetase-tRNA complexes, results in 100-3000-fold binding discrimination against non-cognate tRNAs.     

Structure of transfer RNA molecules containing the long variable loop

A structure is proposed for the type II tRNA molecules containing the long variable loop and the tertiary base interactions here are compared with type I tRNAs having the short variable loop. The type II tRNAs are similar to the type I tRNAs in their tertiary base pairing interactions but differ from them generally by not having the tertiary base triples. The long variable loop, which is comprised of a helical stem and a loop at the end of it, emerges from the deep groove side of the dihydrouridine helix, and is tilted roughly 30° to the plane formed by the amino acid-pseudouridine and anticodon-dihydrouridine helices found in yeast tRNA^Phe. The fact that many of the type I tRNAs also lack the full compliment of base triples suggests that the tertiary base pairs may alone suffice to sustain the tRNA fold required for its biological function. The base triples and the variable loop appear to have little functional significance. The base type at position 9 is correlated with the number of base triples and G-C base pairs in the dihydrouridine stem.     

Quantum chemical studies of structures and binding in noncanonical RNA base pairs: the trans Watson-Crick:Watson-Crick family

The trans Watson-Crick/Watson-Crick family of base pairs represent a geometric class that play important structural and possible functional roles in the ribosome, tRNA, and other functional RNA molecules. They nucleate base triplets and quartets, participate as loop closing terminal base pairs in hair pin motifs and are also responsible for several tertiary interactions that enable sequentially distant regions to interact with each other in RNA molecules. Eleven representative examples spanning nine systems belonging to this geometric family of RNA base pairs, having widely different occurrence statistics in the PDB database, were studied at the HF/6-31G (d, p) level using Morokuma decomposition, Atoms in Molecules as well as Natural Bond Orbital methods in the optimized gas phase geometries and in their crystal structure geometries, respectively. The BSSE and deformation energy corrected interaction energy values for the optimized geometries are compared with the corresponding values in the crystal geometries of the base pairs. For non protonated base pairs in their optimized geometry, these values ranged from -8.19 kcal/mol to -21.84 kcal/mol and compared favorably with those of canonical base pairs. The interaction energies of these base pairs, in their respective crystal geometries, were, however, lesser to varying extents and in one case, that of A:A W:W trans, it was actually found to be positive. The variation in RMSD between the two geometries was also large and ranged from 0.32-2.19 A. Our analysis shows that the hydrogen bonding characteristics and interaction energies obtained, correlated with the nature and type of hydrogen bonds between base pairs; but the occurrence frequencies, interaction energies, and geometric variabilities were conspicuous by the absence of any apparent correlation. Instead, the nature of local interaction energy hyperspace of different base pairs as inferred from the degree of their respective geometric variability could be correlated with the identities of free and bound hydrogen bond donor/acceptor groups present in interacting bases in conjunction with their tertiary and neighboring group interaction potentials in the global context. It also suggests that the concept of isostericity alone may not always determine covariation potentials for base pairs, particularly for those which may be important for RNA dynamics. These considerations are more important than the absolute values of the interaction energies in their respective optimized geometries in rationalizing their occurrences in functional RNAs. They highlight the importance of revising some of the existing DNA based structure analysis approaches and may have significant implications for RNA structure and dynamics, especially in the context of structure prediction algorithms.     

Functional importance of RNA interactions in selection of translation initiation codons

RNA base pairing between the initiation codon and anticodon loop of initiator tRNA is essential but not sufficient for the selection of the 'correct' mRNA translational start site by ribosomes. In prokaryotes, additional RNA interactions between small ribosomal subunit RNA and mRNA sequences just upstream of the start codon can efficiently direct the ribosome to the initiation site. Although there is presently no proof for a similar important ribosomal RNA interaction in eukaryotes, the 5' non-coding regions of their mRNAs and 'consensus sequences' surrounding initiation codons have been shown to be strong determinants for initiation-site selection, but the exact mechanisms are not yet understood. Intramolecular base pairing in mRNA and participation of translation initiation factors can strongly influence the formation of mRNA–small ribosomal subunit–initiator tRNA complexes and modulate translational activities in both prokaryotes and eukaryotes. Only recently has it been appreciated that alternative mechanisms may also contribute to the selection of initiation codons in all organisms. Although direct proof is currently lacking, there is accumulating evidence that additional cis -acting mRNA elements and trans -acting proteins may form specific 'bridging' interactions with ribosomes during translation initiation.     

The difference in the type of codon-anticodon base pairing at the ribosomal P-site is one of the determinants of the translational rate

By utilizing an enzymatically reconstructed tRNA variant containing an altered anticodon sequence, we have examined the different biochemical behavior of translation between the Watson-Crick type and the wobble type base pair interactions at the first anticodon position. We have found that the Watson-Crick type base pair has an advantage in translation in contrast to the wobble type base pair by comparing the efficiency of transpeptidation of native tRNA(Phe) (anticodon; GmAA) with its variant tRNA (anticodon; AAA) in the poly(U)-programmed ribosome system. Thomas et al. [Proc. Natl. Acad. Sci. U.S. (1988) 85, 4242-4246] showed that the wobble codon at the ribosomal A-site accepted its cognate tRNA less efficiently than the Watson-Crick base pairing codon. We report here that the wobble interaction at the ribosomal P-site also affected the rate of translation. This variable translational rate may be a mechanism of gene regulation through preferential codon usage.     

Thermodynamics of RNA-RNA duplexes with 2- or 4-thiouridines: implications for antisense design and targeting a group I intron

Antisense compounds are designed to optimize selective hybridization of an exogenous oligonucleotide to a cellular target. Typically, Watson-Crick base pairing between the antisense compound and target provides the key recognition element. Uridine (U), however, not only stably base pairs with adenosine (A) but also with guanosine (G), thus reducing specificity. Studies of duplex formation by oligonucleotides with either an internal or a terminal 2- or 4-thiouridine (s(2)U or s(4)U) show that s(2)U can increase the stability of base pairing with A more than with G, while s(4)U can increase the stability of base pairing with G more than with A. The latter may be useful when binding can be enhanced by tertiary interactions with a s(4)U-G pair. To test the effects of s(2)U and s(4)U substitutions on tertiary interactions, binding to a group I intron ribozyme from mouse-derived Pneumocystis carinii was measured for the hexamers, r(AUGACU), r(AUGACs(2)U), and r(AUGACs(4)U), which mimic the 3' end of the 5' exon. The results suggest that at least one of the carbonyl groups of the 3' terminal U of r(AUGACU) is involved in tertiary interactions with the catalytic core of the ribozyme and/or thio groups change the orientation of a terminal U-G base pair. Thus thio substitutions may affect tertiary interactions. Studies of trans-splicing of 5' exon mimics to a truncated rRNA precursor, however, indicate that thio substitutions have negligible effects on overall reactivity. Therefore, modified bases can enhance the specificity of base pairing while retaining other activities and, thus, increase the specificity of antisense compounds targeting cellular RNA.     

Correlation between three-dimensional structure and chemical reactivity of transfer RNA

The bases of yeast tRNA^Phe which react with carbodiimide and methoxyamine have been determined and this information has been combined with chemical modification studies of other workers to produce a composite picture of base accessibility in this tRNA. The results are compared with the three-dimensional structure which we have recently determined. The bases which react chemically lie in exposed positions in the three-dimensional model and those which do not are either in the double helical stem regions or else are involved in maintaining the tertiary structure through pairing or stacking interactions.     

The structure of yeast tRNAAsp. A model for tRNA interacting with messenger RNA

Abstract

The anticodon of yeast tRNA^Asp, GUC, presents the peculiarity to be self-complementary, with a slight mismatch at the uridine position. In the orthorhombic crystal lattice, tRNA^Asp molecules are associated by anticodon-anticodon interactions through a two-fold symmetry axis. The anticodon triplets of symmetrically related molecules are base paired and stacked in a normal helical conformation. A stacking interaction between the anticodon loops of two two-fold related tRNA molecules also exists in the orthorhombic form of yeast tRNA^Phe. In that case however the GAA anticodon cannot be base paired. Two characteristic differences can be correlated with the anticodon-anticodon association: the distribution of temperature factors as determined from the X-ray crystallographic refinements and the interaction between T and D loops. In tRNA^Asp T and D loops present higher temperature factors than the anticodon loop, in marked contrast to the situation in tRNA^Phe. This variation is a consequence of the anticodon-anticodon base pairing which rigidities the anticodon loop and stem. A transfer of flexibility to the corner of the tRNA molecule disrupts the G19-C56 tertiary interactions. Chemical mapping of the N3 position of cytosine 56 and analysis of self-splitting patterns of tRNA^Asp substantiate such a correlation.     

The SL1-SL2 (stem-loop) domain is the primary determinant for stability of the gamma retroviral genomic RNA dimer

Retroviral genomes are assembled from two sense-strand RNAs by noncovalent interactions at their 5' ends, forming a dimer. The RNA dimerization domain is a potential target for antiretroviral therapy and represents a compelling RNA folding problem. The fundamental dimerization unit for the Moloney murine sarcoma gamma retrovirus spans a 170-nucleotide minimal dimerization active sequence. In the dimer, two self-complementary sequences, PAL1 and PAL2, form intermolecular duplexes, and an SL1-SL2 (stem-loop) domain forms loop-loop base pairs, mediated by GACG tetraloops, and extensive tertiary interactions. To develop a framework for assembly of the retroviral RNA dimer, we quantified the stability of and established nucleotide resolution secondary structure models for sequence variants in which each motif was compromised. Base pairing and tertiary interactions between SL1-SL2 domains contribute a large free energy increment of -10 kcal/mol. In contrast, even though the PAL1 and PAL2 intermolecular duplexes span 10 and 16 bp in the dimer, respectively, they contribute only -2.5 kcal/mol to stability, roughly equal to a single new base pair. First, these results emphasize that the energetic costs for disrupting interactions in the monomer state nearly balance the PAL1 and PAL2 base pairing interactions that form in the dimer. Second, intermolecular duplex formation plays a biological role distinct from simply stabilizing the structure of the retroviral genomic RNA dimer.     

Nonsense suppressor and antisuppressor mutations at the 1409-1491 base pair in the decoding region of Escherichia coli 16S rRNA.

Using a genetic selection for suppressors of a UGA nonsense mutation in trpA, we have isolated a G to A transition mutation at position 1491 in the decoding region of 16S rRNA. This suppressor displayed no codon specificity, suppressing UGA, UAG and UAA nonsense mutations and +1 and -1 frameshift mutations in lacZ. Subsequent examination of a series of mutations at G1491 and its base-pairing partner C1409 revealed various effects on nonsense suppression and frameshifting. Mutations that prevented Watson-Crick base pairing between these residues were observed to increase misreading and frameshifting. However, double mutations that retained pairing potential produced an antisuppressor or hyperaccurate phenotype. Previous studies of antibiotic resistance mutations and antibiotic and tRNA footprints have placed G1491 and C1409 near the site of codon-anticodon pairing. The results of this study demonstrate that the nature of the interaction of these two residues influences the fidelity of tRNA selection.     

Rapid ribosomal translocation depends on the conserved 18-55 base pair in P-site transfer RNA

The L shape of tRNA is stabilized by the 'tertiary core' region, which contains base-pairing interactions between the D and T loops. Distortions of the L shape accompany tRNA movement across the ribosomal surface. Here, using single-turnover rapid kinetics assays, we determine the effects of mutations within the tertiary core of P site-bound tRNA(fMet) on three measures of the rate of translocation, the part of the elongation cycle involving the most extensive tRNA movement. Mutations in the strictly conserved G18.U55 base pair result in as much as an 80-fold decrease in the rate of translocation, demonstrating the importance of the 18-55 interaction for rapid translocation. This implicates the core region as a locus for functionally important dynamic interactions with the ribosome and leads to the proposal that translocation of ribosome-bound tRNAs may be sequential rather than concerted.     

Quantum Chemical Studies of Structures and Binding in Noncanonical RNA Base pairs: The Trans Watson-Crick:Watson-Crick Family

Abstract

The trans Watson-Crick/Watson-Crick family of base pairs represent a geometric class that play important structural and possible functional roles in the ribosome, tRNA, and other functional RNA molecules. They nucleate base triplets and quartets, participate as loop closing terminal base pairs in hair pin motifs and are also responsible for several tertiary interactions that enable sequentially distant regions to interact with each other in RNA molecules. Eleven representative examples spanning nine systems belonging to this geometric family of RNA base pairs, having widely different occurrence statistics in the PDB database, were studied at the HF/6–31G (d, p) level using Morokuma decomposition, Atoms in Molecules as well as Natural Bond Orbital methods in the optimized gas phase geometries and in their crystal structure geometries, respectively. The BSSE and deformation energy corrected interaction energy values for the optimized geometries are compared with the corresponding values in the crystal geometries of the base pairs. For non protonated base pairs in their optimized geometry, these values ranged from ?8.19 kcal/mol to ?21.84 kcal/mol and compared favorably with those of canonical base pairs. The interaction energies of these base pairs, in their respective crystal geometries, were, however, lesser to varying extents and in one case, that of A:A W:W trans, it was actually found to be positive. The variation in RMSD between the two geometries was also large and ranged from 0.32–2.19 Å. Our analysis shows that the hydrogen bonding characteristics and interaction energies obtained, correlated with the nature and type of hydrogen bonds between base pairs; but the occurrence frequencies, interaction energies, and geometric variabilities were conspicuous by the absence of any apparent correlation. Instead, the nature of local interaction energy hyperspace of different base pairs as inferred from the degree of their respective geometric variability could be correlated with the identities of free and bound hydrogen bond donor/acceptor groups present in interacting bases in conjunction with their tertiary and neighboring group interaction potentials in the global context. It also suggests that the concept of isostericity alone may not always determine covariation potentials for base pairs, particularly for those which may be important for RNA dynamics. These considerations are more important than the absolute values of the interaction energies in their respective optimized geometries in rationalizing their occurrences in functional RNAs. They highlight the importance of revising some of the existing DNA based structure analysis approaches and may have significant implications for RNA structure and dynamics, especially in the context of structure prediction algorithms.     

Fidelity of secondary and tertiary interactions in tRNA.

Contingency table analysis has previously been used to detect sequence correlations in RNAs caused by either secondary or tertiary interactions. An approach known as matrix reduction is developed here as an alternative to the usual Chi square test. This approach is especially sensitive to covariation between equivalent positions and is effective at detecting tertiary interactions that exhibit such covariation. Matrix reduction was also effective at detecting Watson-Crick base-pairs that exhibit a high degree of pairing fidelity. The method was applied to five closely related structural classes of tRNA and a number of tertiary interactions were detected in each class.     

High-resolution nuclear magnetic resonance determination of transfer RNA tertiary base pairs in solution. 1. Species containing a small variable loop.

Eight class I tRNA species have been purified to homogeneity and their proton nuclear magnetic resonance (NMR) spectra in the low-field region (-11 to -15 ppm) have been studied at 360 MHz. The low-field spectra contain only one low-field resonance from each base pair (the ring NH hydrogen bond) and hence directly monitor the number of long-lived secondary and tertiary base pairs in solution. The tRNA species were chosen on the basis of their sequence homology with yeast phenylalanine tRNA in the regions which form tertiary base pairs in the crystal structure of this tRNA. All of the spectra show 26 or 27 low-field resonances approximately 7 of which are derived from tertiary base pairs. These results are contrary to previous claims that the NMR spectra indicate the presence of resonances from secondary base pairs only, as well as more recent claims of only 1-3 tertiary resonances, but are in good agreement with the number of tertiary base pairs expected in solution based on the crystal structure. The tertiary base pair resonances are stable up to at least 46 degrees C. Removal of magnesium ions causes structural changes in the tRNA but does not result in the loss of any secondary or tertiary base pairs.