Next Article in Journal
The Role of Metal Binding in the Amyotrophic Lateral Sclerosis-Related Aggregation of Copper-Zinc Superoxide Dismutase
Next Article in Special Issue
Synthesis of Disaccharide Nucleosides Utilizing the Temporary Protection of the 2′,3′-cis-Diol of Ribonucleosides by a Boronic Ester
Previous Article in Journal
Chemical Composition of Pinus roxburghii Bark Volatile Oil and Validation of Its Anti-Inflammatory Activity Using Molecular Modelling and Bleomycin-Induced Inflammation in Albino Mice
Previous Article in Special Issue
Identification of Optically Active Pyrimidine Derivatives as Selective 5-HT2C Modulators
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stringent Nucleotide Recognition by the Ribosome at the Middle Codon Position

1
Department of Chemistry, University of Pennsylvania, Philadelphia, PA 19104, USA
2
Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, CA 92093, USA
3
Theoretical Biology and Biophysics Group, Los Alamos National Laboratory, Los Alamos, NM 87545, USA
*
Author to whom correspondence should be addressed.
Molecules 2017, 22(9), 1427; https://doi.org/10.3390/molecules22091427
Submission received: 25 July 2017 / Revised: 15 August 2017 / Accepted: 22 August 2017 / Published: 29 August 2017
(This article belongs to the Special Issue Nucleoside and Nucleotide Analogues)

Abstract

:
Accurate translation of the genetic code depends on mRNA:tRNA codon:anticodon base pairing. Here we exploit an emissive, isosteric adenosine surrogate that allows direct measurement of the kinetics of codon:anticodon base formation during protein synthesis. Our results suggest that codon:anticodon base pairing is subject to tighter constraints at the middle position than at the 5′- and 3′-positions, and further suggest a sequential mechanism of formation of the three base pairs in the codon:anticodon helix.

Graphical Abstract

1. Introduction

Great strides have been made in recent years toward understanding the complex process by which the ribosome utilizes tRNAs to mediate translation of the information encoded within mRNA as consecutive three-nucleotide codons into a polypeptide sequence, benefiting from results of structural, biochemical, biophysical and computational studies [1,2,3,4,5,6,7,8,9,10,11]. Detailed kinetic and dynamic studies have led to the formulation of quantitative schemes for the two major steps within the elongation cycle of peptide synthesis: (i) ternary complex (TC, consisting of aminoacyl(aa)-tRNA.EF-Tu.GTP) binding to a posttranslocation (POST) or 70S initiation complex (70SIC), followed by accommodation of aa-tRNA within the A-site of the ribosome and peptide bond formation, leading to pretranslocation (PRE) complex formation [5,9,10,12,13,14,15,16] and (ii) EF-G.GTP-catalyzed conversion of PRE to POST complex, by which tRNAs bound in the A- and P-sites are translocated, along with the mRNA codons to which they are bound, to the P- and E-sites, respectively [16,17,18,19,20,21,22,23]. These studies rely heavily on functional derivatives of tRNAs, ribosomal proteins, or translational factors, site-specifically labeled with either a fluorophore or a fluorescent quencher [4,16,20,21,22,23,24]. However, direct information about mRNA dynamics during elongation is quite limited, leaving many questions about the mechanism of codon:anticodon base pairing unanswered.
Emissive and responsive isosteric mimetics for purine nucleosides, denoted thG and thA (Figure 1A), which can participate in normal Watson-Crick base pairing, have recently been described [25,26,27]. Previously we demonstrated that changes in the fluorescence of mRNAs containing thG at each of the three positions within the first elongation codon [28] could be used to monitor rates of PRE and POST complex formation. Like thG, thA is strongly emissive (Ф = 0.21) (Figure S1) and its fluorescence is responsive to environmental changes [25]. In the present work we report experiments utilizing mRNAs encoding the tetrapeptide fMetLysPheArg, in which each of the adenosines in the first two codons, AUG and AAA encoding tRNAfMet and tRNALys, respectively, is substituted, one at a time, by thA, yielding labeled mRNAs denoted thA1-mRNA, thA4-mRNA, thA5-mRNA, and thA6-mRNA (Figure 1B). Our results allow comparison of the effects on both PRE and POST complex formation of substituting the more bulky, less polar, thiophene C-S-C grouping for the N7-C8-N9 purine grouping in adenosine, providing an example of using “atomic mutagenesis” [29,30] to elucidate structure:function relationships. Our results demonstrate the tight structural constraints to which A:U base pairing at the middle position is subject, and support a mechanism for the formation of the three base pairs in the codon:anticodon helix that proceeds sequentially rather than simultaneously.

2. Results

2.1. thA Fluorescence Intensity Changes Accompanying PRE and POST Complex Formation

Ribosomes programmed with mRNAs thA1, thA4, thA5, and thA6 displayed substantial functionality in forming 70SIC (Figure 2A), PRE complex (Figure 2B), and POST complex (Figure 2C). Stoichiometries of complex formation, measured at a single standard condition, were generally ≥50% relative to ribosomes programmed with control mRNA lacking thA, denoted fMKFR, although POST complex formation with thA5 (30%) was clearly lower.
For 70SICs programmed with thA4, thA5 or thA6, PRE complex formation on addition of cognate EF-Tu.GTP.Lys-tRNALys (Lys-TC) leads, in each case, to significant decreases in fluorescence intensity (Figure 3A–C). These changes arise as a consequence of codon:anticodon interaction, since almost no changes are seen with thA1 (Figure 3D) or on replacing Lys-TC with non-cognate Arg-TC (Figure S2A–C). Translocation of PRE complexes to POST complexes on addition of EF-G.GTP proceeds with a large decrease in the fluorescence intensity of thA1 as tRNAfMet dissociates from the ribosome (Figure 3D), but with moderate to negligible increases in the fluorescence intensities of thA4, thA5 and thA6 (Figure 3A–C). These intensity changes, summarized in Figure 3E, provide spectroscopic signatures for monitoring the kinetics of PRE and POST complex formation during the polypeptide elongation cycle, as described below. Since thA fluorescence responds to alterations in several factors, including solvation, base-pairing, stacking and interactions with proximal chromophores, it is not currently possible to ascribe fluorescence intensity changes to specific molecular perturbations in the thA microenvironment.

2.2. Kinetics of PRE Complex Formation

We and others previously have used changes in the proflavin (prf) fluorescence intensities of aa-tRNAaa(prf) (as part of a TC) and fMet-tRNAfMet(prf) to monitor, respectively, the initial rates of TC binding to the 70SIC [12] and of aa-tRNA accommodation into the A-site prior to dipeptide formation [31]. In this work we compare the rates of prf fluorescence change of Lys-tRNALys(prf) and fMet-tRNAfMet(prf) with those obtained with thA present in thA4-mRNA (Figure 4A), thA5-mRNA (Figure 4B), or thA6-mRNA (Figure 4C) to determine the rates of different substeps within overall PRE complex formation. Reaction mixtures contained either a single fluorescent component, thA alone, or two fluorescent components, thA plus either Lys-tRNALys(prf) or fMet-tRNAfMet(prf). In the latter cases, changes in both thA and prf fluorescence were monitored simultaneously.
Apparent rate constants, collected in Table 1, were evaluated using Equation (1) or Equation (2) (Section 4). In all three cases, the Lys-TC(prf) fluorescence change was the most rapid. The apparent first order constants for 70SICs programmed with thA4-mRNA (Figure 4A), thA6-mRNA (Figure 4C), and fMKFR-mRNA, 44–62 s−1 (Table 1), correspond closely to previously published values [12], whereas the value with thA5-mRNA is ~3-fold slower (18 s−1) (Table 1). For both thA4-mRNA and thA6-mRNA, rates of fMet-tRNAfMet(prf) fluorescence change, which also agree with previously published values [31], are somewhat slower than those for thA fluorescence change. These results lead to the conclusion that, for these labeled mRNAs, the change in thA fluorescence that depends on codon:anticodon interaction follows TC binding and precedes Lys-tRNALys accommodation, a kinetic mechanism that is consistent with the accepted view of the sequence of events involved in PRE complex formation [5,9].
This order of reaction is captured by globally fitting the results for thA4-mRNA and thA6-mRNA (Figure 4A,B) to Scheme 1 for PRE complex formation (more detailed schemes are presented in references [9] and [32]), in which steps 1, 2 and 3 correspond to initial binding of TC to 70SIC (Complex C1), formation of the codon:anticodon helix (Complex C2) and Lys-tRNALys accommodation into the A-site followed by rapid peptide bond formation [5,12,32,33] to complete PRE complex formation, respectively. The fitted microscopic rate constants and relative fluorescence intensities for this Scheme are presented in Table 2.
The situation is more complicated for thA5-mRNA (Figure 4C), for which the thA fluorescence change has clear biphasic character, with a minor, fluorescence decrease that occurs more rapidly than the monophasic change seen with thA4-mRNA or thA6-mRNA and a major fluorescence change that occurs much less rapidly, at a rate that is quite similar to that observed for fMet-tRNAfMet(prf) fluorescence change (Table 1). These differences suggest that, although we are able to globally fit the results obtained with thA5-mRNA to Scheme 1, giving the values for k1, k2, and k3 displayed in Table 2, the identity of C2 differs from that assigned in the case of thA4-mRNA or thA6-mRNA (see Section 3).

2.3. Kinetics of POST Complex Formation

Apparent rate constants for EF-G.GTP-dependent translocation of PRE complexes programmed with thA-containing mRNAs to form POST complexes were determined from the time-dependent changes in fluorescence intensity of either thA (kA) or fMet-aa-tRNAaa(prf) (kprf) [17,18] (Figure 5, Table 3; the lack of a spectral change precluded kA determination for PRE complexes programmed with thA6).
The relative values of kA vs. kprf were different for PRE complexes programmed with thA4, thA1, or thA5. For thA4-programmed ribosomes kA is indistinguishable from kprf, indicating that the movements of tRNA and mRNA are strongly coupled. In contrast, kprf for thA1-programmed PRE complexes is 3-fold larger than kA. This most likely reflects a more dramatic change in thA1 fluorescence on tRNAfMet dissociation from the ribosomal E-site, a step following translocation. Movements of tRNA and mRNA appear quite decoupled for thA5-programmed PRE complexes, with the apparent rate constants for the major and minor changes in prf fluorescence intensity, kprf1 and kprf2, being 6-fold larger and 2-fold smaller, respectively, than kA. The rate constant kprf1 for thA5 is also considerably larger (13 s−1) than the values of kprf for thA1, thA4, and thA6 (2.6–4.4 s−1), an acceleration that might arise from weaker codon:anticodon interaction (see Section 3). The similarity in the values of kA for thA1 and of kprf2 for thA5 (Table 2) suggests that the latter may also correspond to tRNAfMet dissociation.

3. Discussion

The results presented here show that thA, like thG [24], can be incorporated into functional mRNAs, giving rise to fluorescence intensity changes on conversion of 70SICs to PRE complexes and PRE complexes to POST complexes. These changes make possible the monitoring of microenvironment changes of specific codon nucleotides during a cycle of polypeptide elongation, allowing measurement of the kinetics of both PRE complex formation on cognate TC binding to the ribosome and of mRNA translocation catalyzed by EF-G.GTP. Our most interesting results were obtained for PRE complex formation on the first AAA elongation codon, and it is on these results that we will focus in this section. The rate constants presented in Table 2 demonstrate that whereas thA substitution for adenosine at the first and third positions (thA4 and thA6) show little or no interference with normal mRNA function both in TC binding (k1, 120–170 µM−1·s−1) and in the overall rate constant for aa-tRNA accommodation following TC binding [k2k3/(k2 + k3), 6 s−1], major perturbations in mRNA function are observed for thA substitution in the middle position (thA5), including a 3–4-fold decrease in the rate of TC binding (40 µM−1·s−1) and a much larger 15-fold decrease in the overall rate of aa-tRNA accommodation (0.4 s−1).
We attribute these thA5-induced perturbations to the greater sensitivity on the middle codon position to structural changes arising from substitution of the bulkier, less polar, thiophene C-S-C grouping for the N7-C8-N9 purine grouping. Although we posit that TC binding (step 1, Scheme 1) does not involve codon:anticodon base pairing, the more extensive backbone-rRNA H-bonding network at the middle codon position [1,34] could make it more sensitive to structural change than at the 5′ and 3′ codon positions. This interpretation is consistent with results measuring the effects of disruption of the H-bonding network on TC binding, which are more pronounced at the 5-position than at the 4- or 6-positions [35]. The larger effect on the overall rate constant for aa-tRNA accommodation is also reflected in the large differences in the time courses of thA fluorescence change that directly monitors codon:anticodon interaction, which is biphasic for thA5 with rate constants of 40 s−1 (minor phase) and 0.4 s−1 (major phase) but monophasic for thA4 and thA6, with rate constants of ~10–11 s−1. These differences raise two major questions. First, why does the major fluorescence change with thA5 proceed so much more slowly than the monophasic changes observed with thA4 and thA6? Second, why is the thA5 fluorescence change biphasic?
In addressing the first question we note that codon:anticodon base pairing at the middle position is subject tighter structural constraints within the ribosomal decoding center than at the 5′ and 3′ codon positions [34,36,37], in accord with the observation that misreading by the ribosome occurs least frequently at the middle codon nucleotide [38,39]. Specifically, A5:U is stacked between two adjacent base pairs, which form a tight cavern. In contrast, neither 5′-A4:U nor 3′-A6:U, which each have only one stacking interaction, are faced with such a cavern. The “solvent volume” of a base-pair, defined as the volume within a cutoff radius (5 Å) from the base-pair that can be occupied by a solvent molecule [40], provides a quantitative measure of the steric constraint imposed by this cavern and is considerably reduced for A5:U as compared with either A4:U or A6:U (Table 4). Higher constraint at the middle base-pair position provides at least a partial rationale, based on steric crowding, for the dramatically reduced rate of base-pairing. Other factors may also be in play. Thus, PRE complex formation using mRNAs containing the guanosine surrogate thG (Figure 1A), showed a much less marked reduction in the rate of thG fluorescence change when thG is placed at the middle base pair as compared with the 5′ and 3′ base-pairs [28], even though solvent volumes calculated for G:C base pairs also show a much reduced value for the middle base pair (Table 4). A clear rationale for the apparent context dependence of the rate reduction phenomenon is currently lacking, but might be related to the much more stable stacking and hydrogen bonding interactions generally seen for G-C as opposed to A-U base pairs [41,42].
With regard to the biphasic nature of thA5 fluorescence change, we speculate that the answer lies in the detailed kinetic mechanism by which the three base pairs in the codon:anticodon helix are formed, whether simultaneously or sequentially. Although such helix formation in solution proceeds with a single apparent rate constant, it has been assumed that multiple base pairings occur in a series of rapid, unresolved steps [43], in accord with simulation studies of RNA helix formation [44,45]. The same scenario would explain why virtually all of the fluorescence change on conversion of C1 to PRE complex with either thA4-mRNA or thA6-mRNA occurs in Step 2 (Table 2). However, the putative slow formation of the middle base pair with thA5-mRNA affords partial kinetic resolution of the base pairing steps. Accordingly, we interpret the rapid, minor phase for thA5 as corresponding to the formation a first base pair to either A4 or A6, and the very slow, major phase to formation of the second base pair to thA5, which is followed by much the more rapid formation of the third base pair. Interestingly, the 3.7 ± 1.5-fold higher k2 value for thA5-mRNA vs. thA4- or thA6-mRNA (Table 2), which we attribute to formation of the first base pair, is the approximate value one would expect if individual base pairs to thA4- or thA6-mRNA are formed sequentially with each having a first-order rate constant equal to the k2 value found for thA5-mRNA. Very slow formation of the base pair with thA5 would make that step rate-determining for all subsequent steps leading to PRE complex formation, accounting for the similarity in the apparent rates of the major decreases in the fluorescence intensities of both thA and tRNAfMet(prf) (Table 1).
In summary, our results indicate that emissive isosteric purine nucleoside mimetics can be used not only to measure rates of specific steps in ribosome-catalyzed polypeptide elongation, but also to probe position-dependent constraints that modulate the rate of codon:anticodon helix formation, particularly at A-U base pairs. In this connection, it would be interesting to extend these studies to the recently synthesized isothiazolo-based nucleosides, in which the CH group of thA and thG is replaced by N, in a position corresponding to the purines’ N7 [46]. Extension of this approach to examination of the helix formation within other RNA structures, such as ribozymes, IRESs and riboswitches, should be relatively straightforward.

4. Materials and Methods

GTP and puromycin were obtained from Sigma-Aldrich (St. Louis, MO, USA). Escherichia coli tRNAfMet and tRNALys were obtained from Chemical Block (Moscow, Russia). Previously reported procedures were used to prepare tight-coupled 70S ribosomes from E. coli MRE600 cells [31], cloned E. coli His-tagged proteins IF1, IF2, IF3 and EF-G [31], fMet-tRNAfMet and Lys-tRNALys [18]. Proflavin (prf) labeled tRNAs were prepared similarly to described procedures [18] except that tRNAs, reduced with NaBH4, were aminoacylated (and, in the case of fMet-tRNAfMet, formylated) and purified on FPLC (MonoQ) before being labeled with prf [24]. Proflavin labeling stoichiometries were 1.0 ± 0.1 prf/tRNA and 1.1 ± 0.1 prf/tRNA for fMet-tRNAfMet(prf), labeled at position 20, and Lys-tRNALys(prf), labeled at position 16, 17, or 20, respectively.

4.1. Synthesis of Modified mRNAs

mRNAs were prepared using standard solid-phase phosphoramidite synthesis on an Expedite 8909 synthesizer on a 500 Å CPG (1 μmol scale) column. Commercially available reagents and phosphoramidites (Glen Research, Sterling, VA, USA), and a thA-phosphoramidite, prepared as described [25], were employed. Cleavage from the solid support and deprotection were accomplished with 50:50 mixture of MeNH2 in water (40 wt %) and in ethanol (33 wt %) at 35 °C for 6 h. The 2′-TOM group was removed by TEA·3HF at 65 °C for 3 h and the residue was desalted by precipitation (Glen Report 19–20, Glen Research). All oligonucleotides were purified by preparative polyacrylamide gel electrophoresis (PAGE) using the crush and soak method; the desired band was cut out, pulverized, extracted with 50 mM TEAA (pH 7.0) for a minimum of 12 h (while shaking) and decanted. The buffer containing the purified oligonucleotide was lyophilized and the residue was taken up in 0.2 M TEAB (pH 7.0) buffer and desalted on a Sep-pak C-18 (Waters, Milford, MA, USA). The oligonucleotides were eluted with 40% acetonitrile in water. All oligonucleotides were >98% pure as determined by analytical high resolution PAGE (Figure S3). The structures were confirmed by MALDI-TOF mass spectrometry (thA1, calculated M 9567.6, found 9588.0 [M + Na − H]+; thA4, calculated M 9567.6, found 9586.1 [M + Na − H]+; thA5, calculated M 9567.6, found 9591.9 [M + Na − H]+; thA6, calculated M 9567.2, found 9589.6 [M + Na − H]+). Purified oligonucleotides were quantified by UV absorbance at 260 nm at 70 °C.

4.2. Complex Preparation and Quantification

70SIC was prepared, purified and quantified following a previously published procedure [28]. Modified mRNA was used as half of 70S concentration for experiments presented in Figure 2 or 2-fold excess over 70S for all other experiments. All experiments reported in this paper used purified 70SIC except for those shown in Figure 2B, for which the ultracentrifugation step described in [28] was omitted. TC was formed as described [28]. PRE complex was formed, purified and quantified by the amount of co-sedimenting [3H]-Lys [28]. All experiments reported in this paper used purified PRE complexes. POST complexes were prepared and quantified as described [28].

4.3. Fluorescence Experiments

Fluorescence experiments were carried out in Buffer A. Steady-state fluorescence was measured on a Fluorolog-3 spectrofluorometer (Horiba Jobin Yvon, Edison, NJ, USA). Emission and excitation spectra for thA-containing mRNAs were determined using excitation at 360 nm and monitoring emission at 410 nm, respectively. PRE complex traces shown in Figure 3 are the average of two independently determined traces: (i) PRE complex prepared in situ by addition of TC to 70SIC; (ii) purified PRE complex used in the preparation of POST complex. To account for small concentration differences arising from the two preparation methods, the traces for purified PRE complexes were normalized to those for in situ generated PRE complexes, using an 8 nm window covering the fluorescence maximum to generate a normalization factor. The resulting two very similar traces were then averaged. POST complexes were prepared in situ from purified PRE complexes by incubation with EF-G.GTP (2.5 µM for 20 s), and traces were normalized using the same normalization factors as above. For all traces background fluorescence of the corresponding complex made with unlabeled mRNA was subtracted from the observed fluorescence of each labeled ribosome complex prior to normalization to the peak intensity value of 70SIC for each mRNA. For Figure 3E, multiple determinations of single experiments indicate that fluorescence intensity differences of less than 10% are unreliable because of small errors associated with both pipetting and background fluorescence subtraction.
Kinetic experiments were carried out at 25 °C on a KinTek SF-300X stopped-flow spectrofluorometer (KinTek, Snow Shoe, PA, USA). The excitation wavelengths and band-pass filters used were: 360 nm and 410 ± 10-nm for thA-containing mRNAs; 462 nm and 515 ± 15-nm for proflavin. In each independent kinetic experiment, traces of fluorescence intensity changes were obtained as an average of at least 4 shots. In stopped-flow experiments containing both thA and a prf-labeled tRNA, thA and prf signals were collected on successive pulses from the same syringes. Each experiment was performed at least twice.

4.4. Solvent Volume Calculations

Structural models of codon-anticodon-ribosome interactions were based on structures of Yusupova and co-workers (protein databank accession code 3TVF [34]). Bases in the codon and anticodon were changed via molecular replacement using the swapna routine of chimera, which replaces base moieties. Structures were subsequently checked for clashes [47]. For each individual codon-anticodon base pair considered, an environment consisting of residues of the tRNA, mRNA and ribosome within 5.5 Å was used. To extract the solvent volume of these environments, a two-probe volume voxelator method was used, with a minor probe radius of 1 Å, major probe radius of 10 Å, and high grid resolution [40]. Chimera was used to obtain the subset of these solvent volumes within distance rbp of the codon-anticodon base pair of interest, with rbp = 5 Å. The volumes (in Å3) and areas (in Å2) of these solvent volumes were calculated with Chimera.

4.5. Apparent Rate Constants

Apparent rate constants (kapps) were obtained by either single (Equation (1)) or double (Equation (2)) exponential fitting of each independent experiment using Origin (OriginLab), respectively.
F = F0 + F1ek1appt
F = F0 + F1ekapp1t + F2ekapp2t

Supplementary Materials

Supplementary materials are available online.

Acknowledgments

This work was supported by NIH grants GM080376, B.S.C.; GM069773, Y.T.; GM072686, K.Y.S.

Author Contributions

W.L. and M.N. carried out and interpreted all fluorescence and kinetic experiments and prepared ribosomes, elongation and initiation factors, and tRNAs; D.S. prepared thA-labeled mRNAs; K.Y.S. calculated and interpreted solvent volumes; W.L., Y.T. and B.S.C. wrote the manuscript; Y.T. and B.S.C. conceived and designed the experiments.

Conflicts of Interest

The authors declare no competing financial interests.

References

  1. Voorhees, R.M.; Ramakrishnan, V. Structural basis of the translational elongation cycle. Annu. Rev. Biochem. 2013, 82, 203–236. [Google Scholar] [CrossRef] [PubMed]
  2. Yusupova, G.; Yusupov, M. Crystal structure of eukaryotic ribosome and its complexes with inhibitors. Philos. Trans. R. Soc. Lond. B Biol. Sci. 2017, 372, 20160184. [Google Scholar] [CrossRef] [PubMed]
  3. Dever, T.E.; Green, R. The elongation, termination, and recycling phases of translation in eukaryotes. Cold Spring Harb. Perspect. Biol. 2012, 4, a013706. [Google Scholar] [CrossRef] [PubMed]
  4. Frank, J.; Gonzalez, R.L., Jr. Structure and dynamics of a processive Brownian motor: The translating ribosome. Annu. Rev. Biochem. 2010, 79, 381–412. [Google Scholar] [CrossRef] [PubMed]
  5. Wohlgemuth, I.; Pohl, C.; Mittelstaet, J.; Konevega, A.L.; Rodnina, M.V. Evolutionary optimization of speed and accuracy of decoding on the ribosome. Philos. Trans. R. Soc. Lond. B Biol. Sci. 2011, 366, 2979–2986. [Google Scholar] [CrossRef] [PubMed]
  6. Rodnina, M.V.; Wintermeyer, W. The ribosome as a molecular machine: The mechanism of tRNA-mRNA movement in translocation. Biochem. Soc. Trans. 2011, 39, 658–662. [Google Scholar] [CrossRef] [PubMed]
  7. Sanbonmatsu, K.Y. Computational studies of molecular machines: The ribosome. Curr. Opin. Struct. Biol. 2012, 22, 168–174. [Google Scholar] [CrossRef] [PubMed]
  8. Chen, J.; Choi, J.; O’Leary, S.E.; Prabhakar, A.; Petrov, A.; Grosely, R.; Puglisi, E.V.; Puglisi, J.D. The molecular choreography of protein synthesis: Translational control, regulation, and pathways. Q. Rev. Biophys. 2016, 49, e11. [Google Scholar] [CrossRef] [PubMed]
  9. Loveland, A.B.; Demo, G.; Grigorieff, N.; Korostelev, A.A. Ensemble cryo-EM elucidates the mechanism of translation fidelity. Nature 2017, 546, 113–117. [Google Scholar] [CrossRef] [PubMed]
  10. Rodnina, M.V.; Fischer, N.; Maracci, C.; Stark, H. Ribosome dynamics during decoding. Philos. Trans. R. Soc. Lond. B Biol. Sci. 2017, 372, 20160182. [Google Scholar] [CrossRef] [PubMed]
  11. Chen, C.; Zhang, H.; Broitman, S.L.; Reiche, M.; Farrell, I.; Cooperman, B.S.; Goldman, Y.E. Dynamics of translation by single ribosomes through mRNA secondary structures. Nat. Struct. Mol. Biol. 2013, 20, 582–588. [Google Scholar] [CrossRef] [PubMed]
  12. Pape, T.; Wintermeyer, W.; Rodnina, M.V. Complete kinetic mechanism of elongation factor Tu-dependent binding of aminoacyl-tRNA to the A site of the E. coli ribosome. EMBO J. 1998, 17, 7490–7497. [Google Scholar] [CrossRef] [PubMed]
  13. Blanchard, S.C.; Gonzalez, R.L.; Kim, H.D.; Chu, S.; Puglisi, J.D. tRNA selection and kinetic proofreading in translation. Nat. Struct. Mol. Biol. 2004, 11, 1008–1014. [Google Scholar] [CrossRef] [PubMed]
  14. Geggier, P.; Dace, R.; Feldman, M.B.; Terry, D.S.; Altman, R.B.; Munro, J.B.; Blanchard, S.C. Conformational sampling of aminoacyl-tRNA during selection on the bacterial ribosome. J. Biol. Chem. 2010, 399, 576–595. [Google Scholar] [CrossRef] [PubMed]
  15. Schmeing, T.M.; Voorhees, R.M.; Kelley, A.C.; Ramakrishnan, V. How mutations in tRNA distant from the anticodon effect the fidelity of decoding. Nat. Struct. Mol. Biol. 2011, 18, 432–436. [Google Scholar] [CrossRef] [PubMed]
  16. Liu, W.; Chen, C.; Kavaliauskas, D.; Knudsen, C.R.; Goldman, Y.E.; Cooperman, B.S. EF-Tu dynamics during pre-translocation complex formation: EF-Tu GDP exits the ribosome via two different pathways. Nucleic Acids Res. 2015, 43, 9519–9528. [Google Scholar] [CrossRef] [PubMed]
  17. Savelsbergh, A.; Katunin, V.I.; Mohr, D.; Peske, F.; Rodnina, M.V.; Wintermeyer, W. An elongation factor G-induced ribosome rearrangement precedes tRNA-mRNA translocation. Mol. Cell 2003, 11, 1517–1523. [Google Scholar] [CrossRef]
  18. Pan, D.; Kirillov, S.V.; Cooperman, B.S. Kinetically competent intermediates in the translocation step of protein synthesis. Mol. Cell 2007, 25, 519–525. [Google Scholar] [CrossRef] [PubMed]
  19. Munro, J.B.; Wasserman, M.R.; Altman, R.B.; Wang, L.; Blanchard, S.C. Correlated conformational events in EF-G and the ribosome regulate translocation. Nat. Struct. Mol. Biol. 2010, 17, 1470–1477. [Google Scholar] [CrossRef] [PubMed]
  20. Chen, C.; Stevens, B.; Kaur, J.; Cabral, D.; Liu, H.; Wang, Y.; Zhang, H.; Rosenblum, G.; Smilansky, Z.; Goldman, Y.E.; et al. Single-molecule fluorescence measurements of ribosomal translocation dynamics. Mol. Cell 2011, 42, 367–377. [Google Scholar] [CrossRef] [PubMed]
  21. Chen, C.; Cui, X.; Beausang, J.F.; Zhang, H.; Farrell, I.; Cooperman, B.S.; Goldman, Y.E. Elongation factor G initiates translocation through a power stroke. Proc. Natl. Acad. Sci. USA 2016, 113, 7515–7520. [Google Scholar] [CrossRef] [PubMed]
  22. Belardinelli, R.; Sharma, H.; Peske, F.; Wintermeyer, W.; Rodnina, M.V. Translocation as continuous movement through the ribosome. RNA Biol. 2016, 13, 1197–1203. [Google Scholar] [CrossRef] [PubMed]
  23. Wasserman, M.R.; Alejo, J.L.; Altman, R.B.; Blanchard, S.C. Multiperspective smFRET reveals rate-determining late intermediates of ribosomal translocation. Nat. Struct. Mol. Biol. 2016, 23, 333–341. [Google Scholar] [CrossRef] [PubMed]
  24. Pan, D.; Qin, H.; Cooperman, B.S. Synthesis and functional activity of tRNAs labeled with fluorescent hydrazides in the D-loop. RNA 2009, 15, 346–354. [Google Scholar] [CrossRef] [PubMed]
  25. Shin, D.; Sinkeldam, R.W.; Tor, Y. Emissive RNA alphabet. J. Am. Chem. Soc. 2011, 133, 14912–14915. [Google Scholar] [CrossRef] [PubMed]
  26. Samanta, P.K.; Manna, A.K.; Pati, S.K. Thieno analogues of RNA nucleosides: A detailed theoretical study. J. Phys. Chem. B 2012, 116, 7618–7626. [Google Scholar] [CrossRef] [PubMed]
  27. Gedik, M.; Brown, A. Computational study of the excited state properties of modified RNA nucleobases. J. Photochem. Photobiol. A 2013, 259, 25–32. [Google Scholar] [CrossRef]
  28. Liu, W.; Shin, D.; Tor, Y.; Cooperman, B.S. Monitoring translation with modified mRNAs strategically labeled with isomorphic fluorescent guanosine mimetic. ACS Chem. Biol. 2013, 8, 2017–2023. [Google Scholar] [CrossRef] [PubMed]
  29. Erlacher, M.D.; Chirkova, A.; Voegele, P.; Polacek, N. Generation of chemically engineered ribosomes for atomic mutagenesis studies on protein biosynthesis. Nat. Protoc. 2011, 6, 580–592. [Google Scholar] [CrossRef] [PubMed]
  30. Forconi, M.; Benz-Moy, T.; Gleitsman, K.R.; Ruben, E.; Metz, C.; Herschlag, D. Exploring purine N7 interactions via atomic mutagenesis: The group I ribozyme as a case study. RNA 2012, 18, 1222–1229. [Google Scholar] [CrossRef] [PubMed]
  31. Pan, D.; Zhang, C.M.; Kirillov, S.; Hou, Y.M.; Cooperman, B.S. Perturbation of the tRNA tertiary core differentially affects specific steps of the elongation cycle. J. Biol. Chem. 2008, 283, 18431–18440. [Google Scholar] [CrossRef] [PubMed]
  32. Daviter, T.; Gromadski, K.B.; Rodnina, M.V. The ribosome’s response to codon-anticodon mismatches. Biochimie 2006, 88, 1001–1011. [Google Scholar] [CrossRef] [PubMed]
  33. Gromadski, K.B.; Rodnina, M.V. Kinetic determinants of high-fidelity tRNA discrimination on the ribosome. Mol. Cell 2004, 13, 191–200. [Google Scholar] [CrossRef]
  34. Demeshkina, N.; Jenner, L.; Westhof, E.; Yusupov, M.; Yusupova, G. A new understanding of the decoding principle on the ribosome. Nature 2012, 484, 256–259. [Google Scholar] [CrossRef] [PubMed]
  35. Khade, P.K.; Shi, X.; Joseph, S. Steric complementarity in the decoding center isimportant for tRNA selection by the ribosome. J. Mol. Biol. 2013, 425, 3778–3789. [Google Scholar] [CrossRef] [PubMed]
  36. Ogle, J.M.; Brodersen, D.E.; Clemons, W.M.; Tarry, M.J.; Carter, A.P.; Ramakrishnan, V. Recognition of cognate transfer RNA by the 30S ribosomal subunit. Science 2001, 292, 897–902. [Google Scholar] [CrossRef] [PubMed]
  37. Demeshkina, N.; Jenner, L.; Westhof, E.; Yusupov, M.; Yusupova, G. New structural insights into the decoding mechanism: Translation infidelity via a G·U pair with Watson-Crick geometry. FEBS Lett. 2013, 587, 1848–1857. [Google Scholar] [CrossRef] [PubMed]
  38. Sanbonmatsu, K.Y.; Joseph, S. Understanding discrimination by the ribosome: Stability testing and groove measurement of codon-anticodon pairs. J. Mol. Biol. 2003, 328, 33–47. [Google Scholar] [CrossRef]
  39. Szaflarski, W.; Vesper, O.; Teraoka, Y.; Plitta, B.; Wilson, D.N.; Nierhaus, K.H. New features of the ribosome and ribosomal inhibiters: Non-enzymatic recycling, misreading and back-translocation. J. Mol. Biol. 2008, 380, 193–205. [Google Scholar] [CrossRef] [PubMed]
  40. Voss, N.R.; Gerstein, M. 3V: Cavity, channel and cleft volume calculator and extractor. Nucleic Acids Res. 2010, 38, W555–W562. [Google Scholar] [CrossRef] [PubMed]
  41. Turner, D.H. Available online: http://rna.urmc.rochester.edu/NNDB/turner04/stack.txt (accessed on 15 August 2017).
  42. Dima, R.I.; Hyeon, C.; Thirumalai, D. Extracting stacking interaction parameters for RNA from the data set of native structures. J. Mol. Biol. 2005, 347, 53–69. [Google Scholar] [CrossRef] [PubMed]
  43. Pörschke, D.; Eigen, M. Co-operative non-enzymic base recognition. 3. Kinetics of the helix-coil transition of the oligoribouridylic—Oligoriboadenylic acid system and of oligoriboadenylic acid alone at acidic pH. J. Mol. Biol. 1971, 62, 361–381. [Google Scholar] [CrossRef]
  44. Mohan, S.; Hsiao, C.; VanDeusen, H.; Gallagher, R.; Krohn, E.; Kalahar, B.; Wartell, R.M.; Williams, L.D. Mechanism of RNA double helix-propagation at atomic resolution. J. Phys. Chem. B 2009, 113, 2614–2623. [Google Scholar] [CrossRef] [PubMed]
  45. Colizzi, F.; Bussi, G. RNA unwinding from reweighted pulling simulations. J. Am. Chem. Soc. 2012, 134, 5173–5179. [Google Scholar] [CrossRef] [PubMed]
  46. Rovira, A.R.; Fin, A.; Tor, Y. Chemical mutagenesis of an emissive RNA alphabet. J. Am. Chem. Soc. 2015, 137, 14602–14605. [Google Scholar] [CrossRef] [PubMed]
  47. Pettersen, E.F.; Goddard, T.D.; Huang, C.C.; Couch, G.S.; Greenblatt, D.M.; Meng, E.C.; Ferrin, T.E. UCSF Chimera—A visualization system for exploratory research and analysis. J. Comput. Chem. 2004, 25, 1605–1612. [Google Scholar] [CrossRef] [PubMed]
Sample Availability: Not available.
Figure 1. (A) Isosteric guanosine and adenosine surrogates thG and thA; (B) mRNAs site-specifically labeled with thA (in red).
Figure 1. (A) Isosteric guanosine and adenosine surrogates thG and thA; (B) mRNAs site-specifically labeled with thA (in red).
Molecules 22 01427 g001
Figure 2. Activities of thA-containing mRNAs in complex formation. (A) 70SIC formation, measured as fMet-tRNAfMet co-sedimented per mRNA added: ribosomes, 1.0 mM; mRNA 0.5 mM; (B) Normalized PRE complex formation, with background in the absence of mRNA subtracted: ribosomes, 0.8 mM; (C) Normalized POST complex formation: ribosomes, 4.8 mM. Ranges shown are ± average deviation. Control experiments were performed, as indicated, with either no added mRNA, or with added unmodified mRNA. Values reported in (B) were normalized to the value obtained for unmodified fMKFR mRNA after subtraction of background obtained with no mRNA. Values reported in (C) were obtained following the previously published procedure for thG5 [28] and normalized to the value obtained for unmodified fMKFR mRNA.
Figure 2. Activities of thA-containing mRNAs in complex formation. (A) 70SIC formation, measured as fMet-tRNAfMet co-sedimented per mRNA added: ribosomes, 1.0 mM; mRNA 0.5 mM; (B) Normalized PRE complex formation, with background in the absence of mRNA subtracted: ribosomes, 0.8 mM; (C) Normalized POST complex formation: ribosomes, 4.8 mM. Ranges shown are ± average deviation. Control experiments were performed, as indicated, with either no added mRNA, or with added unmodified mRNA. Values reported in (B) were normalized to the value obtained for unmodified fMKFR mRNA after subtraction of background obtained with no mRNA. Values reported in (C) were obtained following the previously published procedure for thG5 [28] and normalized to the value obtained for unmodified fMKFR mRNA.
Molecules 22 01427 g002
Figure 3. Fluorescence changes of thA-containing mRNAs between ribosomal complex conversion (AD); Traces are shown for purified 70SIC (black); PRE complex (red); POST complex (green); (E) Fluorescence intensity changes were obtained by subtracting the spectrum of a purified complex (70SIC or PRE) from the spectrum of a complex prepared in situ (PRE or POST, respectively) and then normalizing the intensity differences to the peak intensity value of purified complex (70SIC or PRE, respectively).
Figure 3. Fluorescence changes of thA-containing mRNAs between ribosomal complex conversion (AD); Traces are shown for purified 70SIC (black); PRE complex (red); POST complex (green); (E) Fluorescence intensity changes were obtained by subtracting the spectrum of a purified complex (70SIC or PRE) from the spectrum of a complex prepared in situ (PRE or POST, respectively) and then normalizing the intensity differences to the peak intensity value of purified complex (70SIC or PRE, respectively).
Molecules 22 01427 g003
Figure 4. Kinetics of PRE complex formation. Ribosomes programmed with: (A) thA4; (B) thA6; (C) thA5. Fluorescence changes on rapid mixing of 70SICs (0.1 μM) With Lys-tRNALys (0.4 µM). Rates of thA fluorescence change were virtually identical using either unlabeled fMet-tRNAfMet and unlabeled aa-tRNA (green), unlabeled fMet-tRNAfMet and Lys-tRNALys(prf) (red), or fMet-tRNAfMet(prf) and unlabeled Lys-tRNALys (purple). Also shown are rates of fluorescence changes of Lys-tRNALys(prf) (dark blue) and of fMet-tRNAfMet(prf) (cyan). The inset in (C) shows an expanded time scale, making clear biphasic nature of thA5 fluorescence change. Solid black lines are fitted to Scheme 1, using rate constants provided in Table 2.
Figure 4. Kinetics of PRE complex formation. Ribosomes programmed with: (A) thA4; (B) thA6; (C) thA5. Fluorescence changes on rapid mixing of 70SICs (0.1 μM) With Lys-tRNALys (0.4 µM). Rates of thA fluorescence change were virtually identical using either unlabeled fMet-tRNAfMet and unlabeled aa-tRNA (green), unlabeled fMet-tRNAfMet and Lys-tRNALys(prf) (red), or fMet-tRNAfMet(prf) and unlabeled Lys-tRNALys (purple). Also shown are rates of fluorescence changes of Lys-tRNALys(prf) (dark blue) and of fMet-tRNAfMet(prf) (cyan). The inset in (C) shows an expanded time scale, making clear biphasic nature of thA5 fluorescence change. Solid black lines are fitted to Scheme 1, using rate constants provided in Table 2.
Molecules 22 01427 g004
Scheme 1. Kinetic scheme for thA fluorescence changes.
Scheme 1. Kinetic scheme for thA fluorescence changes.
Molecules 22 01427 sch001
Figure 5. Kinetics of POST complex formation of thA-containing mRNAs. Fluorescence changes on rapid mixing of PRE complexes (0.1 mM) programmed with thA-containing mRNAs (black) and EF-G.GTP (2 mM) compared with fluorescence changes of fMetLys-tRNALys(prf) (red), present as part of the PRE complex. (A) thA4 and thA6; (B) thA5; (C) thA1. Solid black lines are exponential fits with kapp as shown in Table 3. The inset in (B) shows an expanded time scale.
Figure 5. Kinetics of POST complex formation of thA-containing mRNAs. Fluorescence changes on rapid mixing of PRE complexes (0.1 mM) programmed with thA-containing mRNAs (black) and EF-G.GTP (2 mM) compared with fluorescence changes of fMetLys-tRNALys(prf) (red), present as part of the PRE complex. (A) thA4 and thA6; (B) thA5; (C) thA1. Solid black lines are exponential fits with kapp as shown in Table 3. The inset in (B) shows an expanded time scale.
Molecules 22 01427 g005
Table 1. Apparent rate constants (s−1) for reactions between 70SICs programmed with thA-containing mRNAs and TC (0.4 µM) 1.
Table 1. Apparent rate constants (s−1) for reactions between 70SICs programmed with thA-containing mRNAs and TC (0.4 µM) 1.
mRNAkTC 2kthAkfMet 3
thA453 ± 19 ± 15.2 ± 0.6
thA518 ± 513 ± 3, minor
0.5 ± 0.2, major
0.4 ± 0.1
thA662 ± 1111 ± 16.6 ± 0.8
fMKFR44 ± 3-8.1 ± 0.3
1 Ranges are average ± average deviation for 2 independent determinations. All fits were to a single exponential except for kthA of thA5, which was fit to a double exponential. In the latter case, the major change in fluorescence intensity was 2.2 fold greater than that of the minor change. 2 Phe-TC binding to 70SIC proceeds at ~60 s−1 at 0.4 µM TC, 20 °C [33]. 3 Phe-tRNAPhe TC accommodation proceeds at ~4 s−1 at 25 °C [31].
Table 2. Microscopic rate constants fit to Scheme 1 and relative thA fluorescence intensity changes.
Table 2. Microscopic rate constants fit to Scheme 1 and relative thA fluorescence intensity changes.
mRNAthA4thA6thA5
Rate Constants
k1 (µM−1·s−1)123 ± 13170 ± 942 ± 1
k2 (s−1)10.6 ± 0.311.3 ± 0.140 ± 16
k3 (s−1)12.0 ± 0.415.0 ± 0.50.44 ± 0.07
Relative Intensity Change from 70SIC 1
C10.140.160.05
C21.001.000.35
PRE1.001.001.00
1 Relative fluorescence intensity change on conversion of 70SIC to PRE complex is set equal to 1.00.
Table 3. Apparent rate constants for POST complex formation from PRE complexes 1.
Table 3. Apparent rate constants for POST complex formation from PRE complexes 1.
mRNAkA (s−1)kprf (s−1)
thA44.4 ± 0.34.4 ± 0.3
thA6n.d. 22.6 ± 0.2
thA11.3 ± 0.14.2 ± 0.2
thA52.3 ± 0.213 ± 1, major 3
0.94 ± 0.16, minor 3
1 Ranges are ±average deviation from 2 independent determinations. 2 not determined. 3 Fluorescence change is biphasic (top, kprf1, bottom, kprf2). The major change in fluorescence intensity was 1.7 fold greater than that of the minor change.
Table 4. Solvent volumes of A-U and G-C base pairs 1 at first elongation codon.
Table 4. Solvent volumes of A-U and G-C base pairs 1 at first elongation codon.
CodonCodon PositionVolume (Å3)
AAA5′-A4424
middle-A5271
3′-A6396
GUG5′-G4419
CGUmiddle-G5252
GUG3′-G6406
1 Ref. [28].

Share and Cite

MDPI and ACS Style

Liu, W.; Shin, D.; Ng, M.; Sanbonmatsu, K.Y.; Tor, Y.; Cooperman, B.S. Stringent Nucleotide Recognition by the Ribosome at the Middle Codon Position. Molecules 2017, 22, 1427. https://doi.org/10.3390/molecules22091427

AMA Style

Liu W, Shin D, Ng M, Sanbonmatsu KY, Tor Y, Cooperman BS. Stringent Nucleotide Recognition by the Ribosome at the Middle Codon Position. Molecules. 2017; 22(9):1427. https://doi.org/10.3390/molecules22091427

Chicago/Turabian Style

Liu, Wei, Dongwon Shin, Martin Ng, Karissa Y. Sanbonmatsu, Yitzhak Tor, and Barry S. Cooperman. 2017. "Stringent Nucleotide Recognition by the Ribosome at the Middle Codon Position" Molecules 22, no. 9: 1427. https://doi.org/10.3390/molecules22091427

APA Style

Liu, W., Shin, D., Ng, M., Sanbonmatsu, K. Y., Tor, Y., & Cooperman, B. S. (2017). Stringent Nucleotide Recognition by the Ribosome at the Middle Codon Position. Molecules, 22(9), 1427. https://doi.org/10.3390/molecules22091427

Article Metrics

Back to TopTop