Steric interactions lead to collective tilting motion in the ribosome during mRNA-tRNA translocation

Abstract

Translocation of mRNA and tRNA through the ribosome is associated with large-scale rearrangements of the head domain in the 30S ribosomal subunit. To elucidate the relationship between 30S head dynamics and mRNA-tRNA displacement, we apply molecular dynamics simulations using an all-atom structure-based model. Here we provide a statistical analysis of 250 spontaneous transitions between the A/P-P/E and P/P-E/E ensembles. Consistent with structural studies, the ribosome samples a chimeric ap/P-pe/E intermediate, where the 30S head is rotated ∼18°. It then transiently populates a previously unreported intermediate ensemble, which is characterized by a ∼10° tilt of the head. To identify the origins of head tilting, we analyse 781 additional simulations in which specific steric features are perturbed. These calculations show that head tilting may be attributed to specific steric interactions between tRNA and the 30S subunit (PE loop and protein S13). Taken together, this study demonstrates how molecular structure can give rise to large-scale collective rearrangements.

ARTICLE

Received 1 Aug 2015 | Accepted 31 Dec 2015 | Published 3 Feb 2016

Translocation of mRNA and tRNA through the ribosome is associated with large-scale rearrangements of the head domain in the 30S ribosomal subunit. To elucidate the relationship between 30S head dynamics and mRNAtRNA displacement, we apply molecular dynamics simulations using an all-atom structure-based model. Here we provide a statistical analysis of 250 spontaneous transitions between the A/PP/E and P/PE/E ensembles. Consistent with structural studies, the ribosome samples a chimeric ap/Ppe/E intermediate, where the 30S head is rotated B18. It then transiently populates a previously unreported intermediate ensemble, which is characterized by a B10 tilt of the head. To identify the origins of head tilting, we analyse 781 additional simulations in which specic steric features are perturbed. These calculations show that head tilting may be attributed to specic steric interactions between tRNA and the 30S subunit (PE loop and protein S13). Taken together, this study demonstrates how molecular structure can give rise to large-scale collective rearrangements.

DOI: 10.1038/ncomms10586 OPEN

Steric interactions lead to collective tilting motion in the ribosome during mRNAtRNA translocation

Kien Nguyen1 & Paul C. Whitford1

1 Department of Physics, Northeastern University, Dana Research Center 111, 360 Huntington Avenue, Boston, Massachusetts 02115, USA. Correspondence and requests for materials should be addressed to P.C.W. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 7:10586 | DOI: 10.1038/ncomms10586 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10586

Ribosomes are complex multicomponent molecular machines that are responsible for the synthesis of proteins in all organisms16. In bacteria, the 70S ribosome

(2.4 MDa) is composed of two subunits: a large 50S subunit and a small 30S subunit (Fig. 1a). The 50S subunit consists of two ribosomal RNA (rRNA) molecules (23S and 5S) and over 30 proteins. The 30S subunit contains one rRNA molecule (16S) and roughly 20 proteins. Both subunits have three binding sites (A, P and E) to which transfer RNA (tRNA) molecules sequentially bind during protein synthesis. Following accommodation of an aminoacyl-tRNA (aa-tRNA) molecule into the A site, the nascent peptide chain is transferred from the P-site tRNA to the aa-tRNA in the A site, thereby extending the peptide chain by one amino acid. Subsequently, the A- and P-site tRNAs move (along with the associated mRNA) by one binding site, into the P and E sites. The displacement of tRNAs between binding sites is called translocation, and it results in a vacant A site. By vacating the A site, the next mRNA codon may be read by another incoming aa-tRNA molecule.

The process of translocation is generally described in terms of two steps: (1) tRNA hybrid-state formation and (2) translocation of the mRNAtRNA codon-anticodon pair1,4,7,8. During the rst step, the acceptor ends of the A and P site tRNAs move relative to the large subunit, from the classical A/A and P/P binding states into hybrid A/P and P/E states (Fig. 1b). In the second step, the mRNA and the tRNA anticodon stem-loops (ASLs) move relative to the small subunit, which allows the tRNAs to adopt classical P/P and E/E conformations (Fig. 1c). While hybrid-state formation occurs spontaneously9, ASL movement is signicantly inuenced

by the presence of elongation factor G (EF-G)10. Furthermore, tRNA movement between the binding sites is associated with internal motions of the ribosome, including rotations of the 30S body and head domains (Fig. 1d). Specically, the 30S body rotates relative to the 50S subunit (intersubunit rotation) during hybrid-state formation7,11, whereas rotation of the 30S head (also known as head swivel) is associated with the movement of mRNA and tRNA on the small subunit1214.

It has been accepted that Brownian noise is a major contributor to conformational changes in the ribosome during elongation15. As a result, the study of the ribosome is now expanding to questions about the ribosomes energy landscape. The energy landscape perspective has the broad utility of allowing one to quantitatively explore the structural and energetic factors that govern biomolecular dynamics. While these tools are well established for describing biomolecular folding and function1621, the development of methods to quantitatively analyse the ribosomes landscape2226 is still in its infancy. One avenue for exploring the landscape is through molecular dynamics simulations. At the highest level of detail, quantum mechanical models can provide insights into catalytic steps2729. Similarly, classical mechanical models with explicit-solvent representations are often used to predict the enthalpic contributions of intermolecular interactions26,30,31. While free energies for the full ribosome are currently not accessible with explicit-solvent models, the free energy can be calculated for smaller RNAs32, or for isolated regions of the ribosome, such as the decoding centre33,34. It is important to note that, while these explicit-solvent models are continuously being rened for RNA35, it is not yet clear which model provides the most accurate description of RNA dynamics, in general36. To complement the detailed pictures obtained from explicit-solvent calculations, coarse-grained elastic network models can describe collective, correlated motions of the ribosome. A key nding from coarse-grained models has been that the architecture of the ribosome intrinsically favours global rearrangements along specic directions3739. Since such a vast range of dynamical processes are associated with the ribosome, it is apparent that an equally wide range of models will be needed to fully characterize the numerous contributors to function.

Here we use a model that employs an intermediate level of resolution, in order to probe the physical relationship between mRNAtRNA movement and 30S head motions during the second step of translocation (30S translocation). Specically, we apply molecular dynamics simulations with an all-atom structure-based model40, and identify the roles of tRNA/ribosome sterics and exibility during this large-scale rearrangement. In this model, the unrotated conformation of the ribosome corresponds to the potential energy minimum. In addition, each tRNA molecule has afnity for two adjacent binding sites on the ribosome, where classical tRNA binding conformations are potential energy minima. With this model, we simulated 250 translocation events, which provide a description of sterically accessible pathways and intermediate ensembles between the A/PP/E and P/PE/E ensembles. While the potential energy minima are dened entirely by the unrotated, classical conformations of the ribosome and tRNA, these simulations capture the known chimeric ap/Ppe/E intermediate, where the 30S head is highly rotated. This head-rotated intermediate is in excellent agreement with numerous structures of the ribosome trapped in mid-translocational conformations1214,41. Our simulations also predict the presence of a previously unreported intermediate ensemble that is characterized by a distinct tilt of the 30S head. To identify the structural features that contribute to the head-tilted ensemble, we perform hundreds of additional simulations in which specic residues in the 30S subunit are

a 50S

30S

b

c

E P A

E P A

Body

E P A

Head

E P A

tRNA movement

relative to 30S

d

[p10]head

[afii9835]head

30S head

30S body

50S

90

30S

Body rotation

Head rotation Head tilting

[p10]body

30S 50S

Figure 1 | Structural description of translocation on the 30S subunit. (a) Structure of the 70S ribosome with two tRNAs (red and yellow) and EF-G (orange). The 23S mRNA and 50S proteins are shown in white and ice blue. The 16S mRNA and 30S proteins are in blue and dark blue.(b) Schematic of the ribosome with mRNA (green) bound on the 30S subunit, and the P- and A-site tRNAs in hybrid P/E (red) and A/P (yellow) conformations. The 50S subunit is depicted in grey. The head and body domains of the 30S subunit are shown in blue and cyan. The region in which mRNAtRNA movement occurs on the 30S is demarcated by a dashed box. (c) Schematic of tRNAs in classical E/E and P/P conformations, after translocation on the 30S subunit. (d) Subunit rotations are described by: fbody (body rotation), fhead (head rotation) and yhead (magnitude of head tilting). Counterclockwise rotation of the body (from the perspective shown) and counterclockwise rotation of the head (viewed from above) are dened as positive. All structural representations were prepared using visual molecular dynamics (VMD)65.

2 NATURE COMMUNICATIONS | 7:10586 | DOI: 10.1038/ncomms10586 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10586 ARTICLE

perturbed. These perturbations may be considered a computational analogue of experimental mutational analysis. Using these perturbed models, we show that the tilt motion of the head is a direct consequence of steric interactions between the ribosome and tRNA. Together, our results demonstrate how a simple energetic model can be used to partition the contributors to dynamics, thereby revealing how molecular structure gives rise to biological function.

ResultsTo elucidate the interplay between 30S head dynamics and tRNA displacement on the small subunit, we performed molecular dynamics simulations of the full 70S ribosome in complex with two tRNAPhe molecules, mRNA and EF-G. For these simulations, we used a multi-basin all-atom structure-based model40. In these models, potential energy minima are dened based on (usually 23) experimentally derived structures. Here we used crystallographic structures of unrotated ribosomes to construct a potential energy landscape for translocation (see Supplementary Methods for details), where each tRNA molecule can form short-range (B10 ) stabilizing interactions with two adjacent binding sites on the ribosome (A and P sites, or P and E sites). By construction, the classical tRNA binding conformations and the unrotated ribosome conformation correspond to potential energy minima. In addition, the intramolecular interactions in EF-G were dened such that the ribosome-bound post-translocational (POST) conformation42 corresponds to the global potential energy minimum. This overall energetic representation is consistent with X-ray crystallography and cryo-electron microscopy (cryo-EM) results. That is, the unrotated and classical conformations are minima in the modelled energy landscape, which is necessarily the case in solution since these conformations may be structurally resolved.

Using our structure-based model, we rst simulated 250 independent, unrestrained (unguided) translocation events. While our focus is on the dynamics that occur after the tRNAs adopt the A/P and P/E binding conformations, we initiated simulations from a model representing an earlier point in the translocation process. Specically, each simulation was initiated with the A- and P-site tRNAs in A/A and P/E binding positions (Fig. 2a,b; models described in ref. 43). By starting the simulations from this conformation, we considered the initial A/AP/E to A/PP/E transition as an equilibration period. This equilibration step was provided, in order to reduce the inuence of initial conditions on the dynamics that follow the adoption of the A/PP/E state. Each simulation was terminated once the ribosome reached the POST state, where the tRNAs adopt the classical P/P and E/E binding conformations. The structural characteristics of all simulated events are qualitatively similar, where the A site tRNA rst moves relative to the large subunit, resulting in the formation of the A/PP/E conformation (Fig. 2c,d and Supplementary Fig. 1). Then the ASLs of the tRNAs move relative to the small subunit, which leads to the formation of the P/PE/E conformation. As noted above, rearrangements of the A-site tRNA that occur before reaching the hybrid A/PP/E conformation were not analysed. Throughout the A/PP/E to POST transition, EF-G maintains a post-translocation-like conformation (Fig. 2e and Supplementary Fig. 2). Despite the lack of a conformational rearrangement in EF-G after the A/PP/E stage (Fig. 2fh), EF-G may still facilitate subsequent translocation substeps. That is, since domain IV of EF-G is in close proximity of the 30S A site, it may serve as a steric doorstop that prevents reverse tRNA motion, which would favour forward translocation.

In contrast to previous simulations of tRNA hybrid-state formation44, targeting protocols were not used to guide the

dynamics in the current study. Rather, in the unrestrained simulations here, each transition occurred spontaneously (that is, stochastically). To illustrate the differences between unrestrained and targeted dynamics, we also performed 100 targeted molecular dynamics (TMD) simulations of translocation (Supplementary Fig. 3 and Methods). When targeted molecular dynamics was employed, the dynamics deviate signicantly from that observed in the unrestrained simulations. As discussed below, the stochastic events from unguided simulations reveal how steric interactions between the ribosome and tRNA can lead to

a b

e

E/E P/E

P/P A/A

RA-ASL

Body

A790

A/PP/E

EF-G ASLs

RA-ELB

G1338-U1341

RP-ASL

P/E A/P

Head

c

Representative trajectory (1 of 250)

f

Head rotated

30

25

RP-ASL RA-ASL RA-ELB

20

A/PP/E

R i()

Angle ()

Equilibration Equilibration

15

pe/E

[p10]head ~ 18

ap/P

10

P/PE/E

g

5

Head tilted

0

d

20

Head rotated

16

[afii9835]head ~ 10

12

Head tilted

8

h P/PE/E

4

0

4 0 2,000 4,000 6,000 8,000

E/E P/P

Time (reduced unit)

[p10]head

[afii9835]head

[p10]body

Figure 2 | Complex 30S translocation dynamics emerge from a simple energetic model. (a,b) Multiple coordinates are used to describe tRNA movement. RP ASL and RA ASL are distances of the anticodon stem-loops

(ASLs) of the P and A site tRNAs (red and yellow) to their E/E and P/P positions. These ASL distances were calculated after 30S body alignment of the simulated trajectory to the POST conformation (see Supplementary Information for details). RA ELB is the distance of the A-tRNA elbow to its

location when in the P/P conformation. Elbow distance is determined after 23S alignment of the trajectory to the POST conformation. In a and b, the P/E and A/A tRNA conformations (initial conguration in the simulation) are shown as opaque. The endpoint E/E and P/P conformations are in ghost representation. (c,d) Representative time traces for a single simulated event (1 of 250) shows large-scale movements of the tRNAs (2030 , c) and transient subunit rotations in the ribosome (d). In all 250 simulations, formation of the A/PP/E conformation, where RA ELB (grey line)

decreases from B30 to B0 , precedes displacement of the ASLs (red and yellow lines). In c and d, the dashed white-faded boxes outline the initial equilibration period, where the A-site tRNA elbow relaxes to the 50S P site and adopts an A/P conformation. (eh) Structural snapshots of the mRNA binding track (from the perspective of the 50S subunit) during translocation illustrate the tRNA ASL positions, relative to the 30S body and head. During the initial equilibration period, the system adopts an A/PP/E conformation (e). During translocation, the 30S head rotates (f), which is followed by back rotation and tilting of the head (g). Translocation is completed when the classical P/PE/E conformation (h) is reached.

NATURE COMMUNICATIONS | 7:10586 | DOI: 10.1038/ncomms10586 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10586

rotational and tilt movements of the 30S head during the second

step of translocation.

Sterics and exibility give rise to a known intermediate. The dynamics of biomolecular function are governed by the free-energy landscape, which is a result of stabilizing interactions, entropy and steric repulsion. While in a structure-based model the dominant potential energy minima are dened a priori, a more complex free-energy landscape can arise from entropic contributions and steric restrictions on tRNA/ribosome movement. As an example, simulations of aa-tRNA accommodation have shown that structural effects lead to pronounced free-energy barriers, even when the potential energy surface is smooth25. Since the current simulations were terminated once the POST state was reached, the free energy is not directly accessible here. However, one may infer the structural characteristics of likely intermediates and transition states from these simulations. Specically, when the ribosome encounters a free-energy barrier, molecular movements are transiently impeded, which leads to increased sampling of the local conguration space. Hence, a poorly populated region between two highly sampled regions is the signature of a free-energy barrier (a transition state) that separates two minima (intermediates or endpoint ensembles).

To identify intermediates associated with translocation, we calculated probability distributions as functions of structural coordinates. These distributions represent the statistical properties of the dynamics, obtained from hundreds of independent simulated events. We characterize the movement of the 30S head during mRNAtRNA translocation, by calculating the probability as a function of the ASL position of the P-site tRNA (relative to the 30S body) and the rotation of the 30S head: P(RP ASL, fhead)

(Fig. 3a). See Supplementary Methods for detailed descriptions of tRNA and ribosome coordinates. To provide a complementary description, we also calculated the probability as a function of the ASL position of the P-site tRNA and the magnitude of head tilt: P(RP ASL, yhead) (Fig. 3b). In these distributions, there are two

visibly distinct populations separating the A/PP/E and P/PE/E ensembles, which are labelled ap/Ppe/E and HT (for head-tilted). Populating the ap/Ppe/E ensemble is associated with movement of the tRNA from the P site towards the E site of the 30S body (RP ASLE46 ) and a large rotation of the 30S head

(fhead 16.51.5;means.d.). For reference, RP

ASL

a

ln(P)

c Head tilt axis direction: [afii9851]head

20

ap/Ppe/E

HT

A/PP/E

P/PE/E

5

15

4

[p10] head()

[afii9835] head()

5

180

0

Parallel axis

10

3

Head

Body

5

2

0

1

0

90 90

Axis perpendicular to mRNA

b

ln(P)

d

Tilt axis direction in HT ens.

16

6

0.04

0 180 90 0 90 180[afii9851]head ()

HT

5

Parallel to mRNA

FWHM=24

12

4

0.03

8

3

Probability

0.02

2

4

0.01

1

0 0 4 8 12 16

0

RP-ASL ()

Figure 3 | A/PP/E to P/PE/E transition involves large-scale rotation and tilting of the 30S head. Probability distributions calculated from250 unrestrained (unguided) simulations highlight the relationship between tRNA and 30S head motions. (a) For the transition from theA/PP/E to the P/PE/E ensemble, the probability P(RP ASL,fhead)

indicates the presence of two intermediate ensembles: a chimeric ap/Ppe/E ensemble and a head-tilted (HT) ensemble. The ap/Ppe/E ensemble includes highly rotated head congurations (fheadE1420).

Adopting the HTensemble is associated with nearly complete back rotation of the head (Dfhead E 13). (b) P(RP ASL,yhead) shows that there

is a large degree of head tilting (yheadE10) in the HT ensemble.

(c) Representative directions (parallel or perpendicular to the mRNA) of the tilt axis of the 30S head are shown. The direction of the tilt axis is measured by whead (see Supplementary Information for details).

By construction, whead 0 corresponds to tilting about an axis that is

parallel to the mRNA, where the head is displaced away from the 30S50S interface. Similarly, whead 180 corresponds to tilting towards the 30S50S

interface. In contrast, whead 90 corresponds to tilting that is

perpendicular to the mRNA. (d) P(whead) shows that head tilting in the HT

ensemble occurs predominantly along the whead 0 direction (that is,

about the mRNA axis), with a full-width at half-maximum (FWHM; red bar) of 24. As a technical note, simulations were initiated from an A/AP/E conformation. Initial relaxation into the A/PP/E ensemble was allowed before data was included for analysis. Accordingly, probability distributions were calculated for RA ELBo4 .

0 when

the tRNA is in the 30S E site. After sampling the ap/Ppe/E enemble, the system transitions to the head-tilted ensemble, which involves nearly complete back rotation of the head (DfheadE 13). Interestingly, this reverse head rotation is

accompanied by a previously unreported tilt of the 30S head (yhead 9.52.2) (Fig. 3b), which results in displacement of the

head away from the 30S50S intersubunit interface. The translocation process is completed when the tRNA molecules reach the POST state, where the head rotation and tilt angles (fhead and yhead) return to near-zero values.

The ap/Ppe/E ensemble observed in our simulations is in excellent agreement with models derived from cryo-EM and X-ray crystallography1214,41 (Supplementary Table 1). In particular, the simulated ap/Ppe/E ensemble is most similar to the cryo-EM reconstruction reported in ref. 12, which describes the ribosome with EF-G bound and tRNAs in ap/P and pe/E binding positions. Before comparing this model with our intermediate ensemble, it is instructive to discuss the overall characteristics of the cryo-EM model. In the ap/Ppe/E ribosome, the tRNA molecules are partially displaced relative to the 30S body, and the 30S head is highly rotated. This conguration allows the tRNAs to contact the A and P sites of the head, while

simultaneously contacting the P and E sites of the body. This chimeric form of intrasubunit binding is reected in the naming convention, as introduced previously14: pe/E indicates that the P-site tRNA contacts the P site of the 30S head (p), the E site of the 30S body (e), and the E site of the 50S subunit (E).

Since a given experimental structure represents the average conguration within a sample, we calculated the average structure of all simulated snapshots that were in the ap/Ppe/E ensemble. The root-mean-square deviation (RMSD) between the average ap/Ppe/E conguration and the cryo-EM model12 is 1.93 (calculated for all core atoms, see Supplementary Methods). While this average structure provides an overall description of the ensemble, the averaging process yields a simplied perspective that suppresses information about uctuations. Accordingly, subsequent comparisons of structural metrics include the mean and s.d. of each coordinate, calculated for a specied ensemble. In terms of subunit rotations, our ap/Ppe/E ensemble is characterized by a large rotation of the head (fhead 16.51.5; Fig. 3a) and a modest degree of body

4 NATURE COMMUNICATIONS | 7:10586 | DOI: 10.1038/ncomms10586 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10586 ARTICLE

rotation (fbody 1.91.2; Supplementary Fig. 4). These values

agree well with those obtained from cryo-EM and X-ray models1214,41, particularly with those of ref. 41, for which fhead

16.6 and fbody 1.6 (see Supplementary Table 1 for

rotation angles calculated for other experimental structures). Consistent with the experimental structures mentioned above, we also nd that EF-G adopts a post-translocation-like conformation, where its domain IV extends towards the A site of the 30S subunit. Specically, the average RMSD of EF-G in the ap/Ppe/E ensemble from its POST conformation is 2.20.3 (Supplementary Fig. 2). Furthermore, in the simulated ap/Ppe/E ensemble, there is a compaction of the ASLs of the tRNAs. To describe this compaction, we calculated the distance between the P atoms of the A31 residues in the A- and P-site tRNAs: RA31 (Supplementary Fig. 5a). When the tRNA molecules are in classical binding positions, RA31 is 30.0 (A/AP/P) and

25.2 (P/PE/E)45. In contrast, in the ap/Ppe/E intermediate from cryo-EM12, RA31 adopts a smaller value of 21.3 . This experimentally observed compaction of the ASLs is very similar to that found in our simulated ap/Ppe/E ensemble, where RA31 19.11.6 (Supplementary Fig. 5b). As a nal point of

comparison, the ASL position of the P-site tRNA relative to the 30S head is also consistent with the experimentally obtained structures of this intermediate1214,41. Specically, the distance between the P atoms of A31 in the P-site tRNA and A1229 in the 16S rRNA (Rpe) is 12.2 in the structural model12, and Rpe 11.51.3 in the simulated ap/Ppe/E ensemble

(Supplementary Table 1). Overall, the high level of agreement between our simulated ap/Ppe/E ensemble and the models derived experimentally demonstrates that the simplied energetic model is sufcient to capture this en-route translocation intermediate.

After the ap/Ppe/E ensemble is reached, the simulated system populates another intermediate ensemble, the head-tilted ensemble (labeled HT in Fig. 3a,b). The most striking structural feature of this ensemble is that the 30S head is highly tilted (yhead 9.52.2, Fig. 3b), which has not been reported in the

context of translocation. However, a similar degree of head tilting has been observed in the context of transfer-messenger RNA rescue46 (yhead 11.5). Below, we discuss the structural and

energetic features that give rise to this large-scale collective tilt motion, as well as the degree to which tilting is correlated with head rotation and tRNA displacement.

mRNAtRNA translocation involves 30S head rotation and tilt. During the transition from the ap/Ppe/E to the head-tilted ensemble, there is an apparent correlation between back rotation and tilting of the 30S head (DfheadE 13 in Fig. 3a, and

DyheadE6 in Fig. 3b). To quantify the magnitude and direction

of head tilting, we used the angles yhead and whead (shown

schematically in Figs 1d and 3c; for further details see Supplementary Figs 6 and 7, and Supplementary Methods). The tilt direction (whead) is the orientation of the axis about which tilting occurs. We dened whead 0 such that it corresponds to

tilting that is approximately about the axis of the mRNA backbone, where the 30S head is displaced away from the 30S50S interface (Fig. 3c, and Supplementary Figs 6 and 7). In addition to the direction of tilting, the magnitude of head tilting is described by yhead, as noted above. To characterize the direction

of head tilting in the head-tilted ensemble (for example, parallel or perpendicular to the mRNA), we calculated the probability P(whead) for head-tilted conformations (Fig. 3d). Since there were no explicit restrictions on the direction of tilting, it would be conceivable that tilting could occur in directions that are parallel or perpendicular to the axis of the mRNA. However, we nd that

P(whead) peakes at whead 0, where the full-width at

half-maximum is 24. This narrow distribution of tilt-axis direction demonstrates that tilting in the head-tilted ensemble occurs predominantly about the mRNA axis.

The apparent balance between rotation and tilt in the head-tilted ensemble can be interpreted as enabling opening of the mRNA binding track on the 30S subunit. This is consistent with previous structural studies that have found substantial changes in the steric features of the mRNA binding track upon head rotation47. As discussed in ref. 47, residues G1338 to U1341 in the head (that is, the PE loop) sterically separates the P and E sites of the 30S subunit. This steric obstacle has been described as a gate that opens when the PE loop (in the head) is displaced relative to residue A790 (in the body)47. Here to measure gate opening, we used the distance between the P atoms of residues U1340 and A790: Rgate (Fig. 4a). In a crystallographic model of the classical conformation, where the 30S head is unrotated and untilted, Rgate 16.7 (ref. 45). In contrast, in the ap/Ppe/E structural

model12, where the head is highly rotated, Rgate has a signicantly larger value of 25.1 . This substantial change in the gate region has led to the argument that opening may allow tRNA molecules to transition to the POST conformation4,6. While gate dynamics likely contributes to translocation, the fact that the head-rotated ap/Ppe/E conformation is stable when the gate is open12,13 suggests that head rotation alone is not sufcient to confer tRNA passage.

To describe the gate dynamics in the simulations, we evaluated the probability as a function of the ASL position of the P-site tRNA and the opening of the gate: P(RP ASL, Rgate) (Fig. 4b).

Consistent with the notion that gate-opening facilitates translocation, we observe a large degree of gate opening (Rgate 27.71.1

) in the ap/Ppe/E ensemble. Similarly, we nd that gate separation remains large in the head-tilted ensemble (Rgate 24.51.5 ). This is particularly interesting because, if

the head were unable to tilt, one would expect an approximately linear relationship between head rotation (fhead) and gate

opening (Rgate). However, during the transition from the ap/P pe/E to the head-tilted ensemble, we nd a marginal change in gate separation (DRgateE 3 ) over a large range of head

rotation angles (DfheadE 13). In contrast, during the

subsequent transition from the head-tilted to the P/PE/E

16 0 4 8 12 16

a

b

ln(P)

32

Body

ap/Ppe/E

A/PP/E

5

A790

E P

Rgate

28

4

R gate()

24

3

HT

2

20

1

G1338-U1341 (PE loop)

P/PE/E

0

Head

RP-ASL ()

Figure 4 | Changes in steric features of the mRNA binding track are associated with tRNA movement. (a) The PE loop (residues G1338-U1341) in the head and residue A790 in the body (pink) form a steric gate separating the E and P sites. Rgate is the distance between the P atoms of

U1340 and A790 (black beads), and is used to probe opening and closing of the gate. (b) The probability distribution P(RP ASL,Rgate) implicates a large

degree of gate opening (Rgate 2630 ) in the ap/Ppe/E ensemble,

which is associated with head rotation (fheadE18, cf. Fig. 3a,b). In the

subsequently populated HT ensemble, gate opening is only slightly decreased. While there is a reduced degree of head rotation in the HT ensemble, large gate opening (RgateE2326 ) is maintained by a

compensatory tilt of the head (yheadE10, cf. Fig. 3a,b).

NATURE COMMUNICATIONS | 7:10586 | DOI: 10.1038/ncomms10586 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10586

ensemble, there is a signicant reduction in gate separation (DRgateE 7 ) over a smaller range of rotation values

(DfheadE 5). While initial opening of the gate is associated

with head rotation, as implicated by structural data1214,41, our simulations predict that gate closure only occurs when the head relaxes from a tilted orientation. This suggests an extended description of gate opening, where there are compensatory rotation and tilting uctuations that together provide a sufciently wide corridor for tRNA to translocate.

Isolating the origins of head tilting. From visual inspection of the simulated trajectories, it is apparent that the PE loop and protein S13 transiently interact with the tRNA molecules during translocation (Fig. 5ac and Supplementary Movie 1). Consistent with this observation, we veried that in all simulations the tRNAs transiently contact (are within 3 of) the C-terminal tail of S13 and the PE loop. To evaluate the extent to which these tRNAribosome interactions contribute to head tilting, we performed additional simulations in which specic steric interactions were removed. Specically, we tested the relative inuence of protein S13 and the PE loop by excluding each of their steric contributions in separate sets of simulations. For these tests, we used the following variations of our multi-basin structure-based model: Model (1) Same model as described above; Model (2) Identical to Model 1, except steric interactions between S13 and tRNA were not included; Model (3) Identical to Model 1, except steric interactions between the PE loop and tRNA were not included; Model (4) Identical to the above models, except both S13 and PE loop sterics were not included. See Supplementary Methods for technical descriptions of the models. These simulations mimic an ideal mutation experiment. That is, mutational studies aim to modulate a single variable at a time, in order to elucidate the relative contributions of specic interactions. Here we removed specic steric interactions between residues without

introducing indirect perturbations. Thus, changes in the dynamics can be directly attributed to the steric interactions that are varied.

For each of the four models, we simulated B200 independent, unrestrained translocation events. To compare the degree of head titling in each model, we calculated the probability as a function of head rotation and head tilting: P(fhead, yhead) (Fig. 5dg). Since

the focus of this exercise was to explore the role of specic interactions during formation of the head-tilted ensemble, these additional simulations were initiated from an ap/Ppe/E conformation and terminated when the system adopted the P/PE/E conformation (fheadE0 and yheadE0). We nd that

steric interactions have signicant effects on the degree of head tilting, as highlighted by the average tilt angle as a function of head rotation: yheadfhead (dashed lines in Fig. 5dg). When all

steric interactions are included (Model 1), yhead reaches a maximum value of B10. When the steric interactions between protein S13 and tRNA are excluded (Model 2), the maximum value of yhead is reduced to B8. In contrast, there is a more signicant reduction in head tilting when the sterics of the PE loop are not included (Model 3). Under those conditions, yhead only reaches a value of B6, corresponding to a 40% reduction in the scale of head tilting. Finally, when the sterics of both S13 and the PE loop are not included (Model 4), head tilting is further attenuated with a maximum yhead value of only 34. This residual head tilting indicates that, while the dominant contributions stem from PE loop and S13 interactions, additional structural factors also inuence tilting. Together, comparison of the models demonstrates that both the PE loop and protein S13 contribute to head tilting during back rotation of the head, though the relative contribution of the PE loop is larger.

Previous biochemical and genetic assays have shown that protein S13 is involved in stabilizing the pre-translocational state48, and that ribosomes lacking S13 exhibit increased rates of factor-free translation49. Consistent with these observations, our

a b c

ap/Ppe/E HT P/PE/E

Body

Head

[afii9835]head~ 10

PE loop

[p10]head~ 18 [p10]head~ 6 [p10]head = 0, [afii9835]head = 0

S13

d e f g

P/PE/E

ap/Ppe/E

ln(P)

15

S13 sterics removed

PE loop sterics removed S13 and PE loop sterics removed

4

10

3

[afii9835] head ()

2

5

1

0 5 0 5 10 15 20

5 0 5 10 15 20 5 0 5 10 15 20 5 0 5 10 15 20

0

[p10]head () [p10]head () [p10]head () [p10]head ()

Figure 5 | Head tilting results from tRNA interactions with the PE loop and protein S13. Representative structural snapshots of (a) ap/Ppe/E,(b) HT and (c) P/PE/E congurations. The PE loop (pink) and protein S13 (grey) are in close proximity to the P- and A-site tRNAs (red and yellow) during translocation. To elucidate the extent to which the PE loop and protein S13 contribute to tilting, perturbations were introduced to the model. (dg) Probability distributions P(fhead, yhead) illustrate the relative inuence of steric composition on tilting. When all sterics are included (d), yheadfhead

(the average tilt as a function of rotation; dashed grey line) reaches a maximum value of B10. When steric interactions between protein S13 and tRNA are not included (e), the maximum value of yhead is reduced to B8. When the sterics of S13 are present, while steric interactions between the PE loop and tRNA are not included (f), the maximum tilt is reduced to B6. When both S13 and PE loop sterics are not included (g), head tilting is most strongly attenuated ( yheado5 ).

6 NATURE COMMUNICATIONS | 7:10586 | DOI: 10.1038/ncomms10586 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10586 ARTICLE

results suggest that protein S13 (in addition to the PE loop) poses a steric barrier that impedes the movement of the mRNAtRNA complex towards the POST state. That is, before the reverse rotation of the 30S head (before, and including the ap/Ppe/E state), the C terminus of S13 is positioned between the A- and P-site tRNAs (Fig. 5a). As discussed above, this positioning of S13 then leads to head tilting during reverse rotation. It is conceivable that the steric barrier introduced by S13 may also impede forward movement of the A-site tRNA relative to the 30S head. This suggests that, in addition to potentially introducing stabilizing interactions between S13 and tRNA, the steric inuence of S13 may increase the observed dwell time of the pre-translocational state.

DiscussionAs our physicalchemical understanding of the ribosome continues to be developed, an emerging theme is that elongation dynamics result from an interplay between energy accumulation and release in the ribosome and tRNA. In the words of Frank et al.50, a tRNA molecule may be regarded as a (non-linear) molecular spring, where the classical tRNA conformations correspond to low-energy states. During elongation, tRNA molecules then interconvert between nonclassical states, which is facilitated by tRNAribosome interactions and factor binding. Consistent with this perspective, our model describes the classical conformations of tRNA as energetic minima. That is, when a tRNA adopts hybrid conformations, the internal energy of the molecule is higher (B5 reduced units, Supplementary Methods), though these states are also stabilized by ribosometRNA interactions. In terms of concepts developed to describe protein function51, nonclassical binding states are associated with an accumulation of strain energy, which is released when the system reaches the low-energy P/PE/E conformation.

The strain energy framework may also be applied to the ribosome. Specically, when the tRNAribosome assembly reaches the chimeric ap/Ppe/E ensemble, the 30S head becomes highly rotated, implying that the small subunit is in a strained state. In contrast, in the ap/Ppe/E ensemble, the tRNAs adopt conformations that are near their classical P/P and E/E conformations. That is, the RMSD values of the average simulated ap/P and pe/E tRNA congurations, relative to the P/P and E/E conformations, are only 0.42 and 0.96 . Experiments have shown that thermal energy enables tRNA to transiently sample energetically strained hybrid conformations9,11,52. Our analysis extends this description, and suggests that strain energy accumulated in the tRNAs may be subsequently transferred to the ribosome as the ap/Ppe/E ensemble is reached. This strain in the ribosome is then released as the system relaxes to the unrotated POST state.

With regards to strain in EF-G, structural studies have implicated a range of different conformations of the factor bound to ribosomes in numerous states, including the A/AP/P53, A/P* P/E54, ap/appe/E13 and ap/Ppe/E12 states. These structural considerations suggest that strain accumulation and release in EF-G may likely occur during earlier stages of translocation, particularly before the system reaches the head-tilted ensemble. In the preceding ap/Ppe/E ensemble, EF-G adopts a post-translocation-like conformation, where its domain IV is adjacent to the ap/PtRNA (Supplementary Fig. 2), implying that EF-G is energetically relaxed. Therefore, as the ribosome transitions from the ap/Ppe/E via the head-tilted to the P/PE/E ensemble, strain energy is not expected to be signicant. Rather, during these rearrangements, the presence of close interactions between domain IV and tRNA suggests that EF-G may act in a more passive capacity. That is, it appears to simply serve as a steric

doorstop that prevents reverse movement of the tRNAs. By preventing reverse motion, EF-G effectively stabilizes head-rotated conformations, relative to the pre-translocational state, as suggested by ref. 13. In addition, compaction of the tRNA ASLs and domain IV may position the tRNA molecules such that the free-energy barrier for forward translocation is reduced. This raises the possibility that EF-G may accelerate translocation kinetics, even in the absence of direct energy transfer between GTP hydrolysis and the ribosome.

The picture of strain accumulation and release provides a framework for exploring the balance between molecular exibility and sterics, which may be extended to account for more detailed energetic contributions. While our models elucidate steric effects and implicate correlated uctuations, the complete energy landscape of the ribosome is complex, and it results from many factors, including chargecharge interactions and solvation. In the current model, the stabilizing contribution of these factors are implicitly described. However, when the system is far from the endpoints (for example, in transition states), many transient interactions are formed. Here we have shown that steric interactions frequently occur during translocation, primarily in the form of RNARNA backbone interactions. Since the RNA backbone is negatively charged, one could expect that the observed steric effects will be amplied by electrostatic effects. For example, the steric repulsion between the P-site tRNA and the PE loop will likely become stronger when explicit electro-statics are included, though counterion effects may partially mitigate this amplication. Similarly, the degree of tRNA compaction in the ap/Ppe/E ensemble may also be modulated by electrostatic interactions. While the current model predicts a compaction of the tRNA ASLs that is similar to cryo-EM observations, the tRNA molecules are roughly 1 closer to each other in the simulations. Repulsive chargecharge interactions between the tRNAs may decrease the degree of compaction, which would further improve the level of agreement between cryo-EM and the simulations. In addition to electrostatics, the models may be extended to explore the role of desolvation. In the context of protein folding, many theoretical and experimental ndings have implicated an essential role of desolvation20,5557, where models that include desolvation effects predict more cooperative protein folding dynamics. In terms of RNA base stacking, desolvation effects have been extensively quantied computationally58. Since the conformational rearrangements probed in the current study do not involve breaking or formation of basestacking interactions, we do not expect the current results to be sensitive to desolvation effects. However, for other steps of elongation, such as displacement of the tRNA CCA end inside the peptidyl transferase centre, the effect of desolvation may be pronounced.

In summary, this work demonstrates how simple energetic considerations can elucidate the relationship between molecular structure and functional dynamics in the ribosome. Our ndings provide a structural and energetic framework for describing and interpreting the complex dynamics of mRNAtRNA translocation. With this framework, we may now explore the physical relationship between the many molecular components of the ribosome, which can help guide the design of more precise experimental measurements. Through the continued integration of theoretical and experimental insights, it will be possible to establish a cohesive physicalchemical description that bridges detailed biochemical interactions and large-scale dynamics.

Methods

The forceeld. To simulate translocation, we used a multi-basin all-atom structure-based model. In the simplest form, a single-basin structure-based model denes an experimentally derived structure as the global potential energy

NATURE COMMUNICATIONS | 7:10586 | DOI: 10.1038/ncomms10586 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10586

minimum40,59. To construct a potential energy function for translocation, we rst generated a structure-based model for the P/PE/E conguration. In addition, contacts between the ribosome and mRNAtRNA, which are found in the A/AP/P conguration were included as stabilizing interactions. The energetic weights of the contacts between the ribosome and mRNAtRNA were rescaled by a factor of 0.3, to account for the transient nature of mRNAtRNA binding. That is, intrasubunit interactions maintain the long-timescale structural integrity of the ribosome, whereas mRNAtRNA associates and dissociates from the ribosome on much shorter timescales. Interface contacts between the 30S and 50S subunits, as well as between the 30S head and 30S body, were also weakened (see Supplementary Methods for details). These contacts were reduced in strength to mimic the observation that subunit rotations spontaneously occur11 on timescales that are far shorter than the lifetime of a fully formed ribosome. For complete details, see Supplementary Information.

Simulation details. Simulations were performed with the Gromacs (v4.6.1) software package60,61 using forceeld les generated by the SMOG-model web server (smog-server.org)62. Reduced units were used for all calculations. Each simulation was performed for a minimum of 4 106 time steps of size 0.002, and

was extended until the P/PE/E conguration was adopted. Employing a timescale correction factor reported by Kouza et al.63, the effective timescale of each simulation may be roughly approximated as 10100 ms. For a detailed discussion on timescale estimates in all-atom structure-based models, see ref. 25 and the Supplementary Material of ref. 64. Langevin Dynamics protocols were used to ensure a constant temperature of 0.5 (reduced units). As discussed elsewhere64, the scale of structural uctuations at this temperature is consistent with that obtained using an explicit-solvent model at 300 K, and with values estimated from crystallographic B-factors. 250 simulations were initially performed for the full mRNAtRNA translocation process (results presented in Figs 24). Around200 additional simulations were performed for each of the four modied forceelds, where the sterics of the PE loop and/or protein S13 were excluded (results presented in Fig. 5). In total, 1,031 translocation events were simulated without the use of targeting or steering protocols.

References

1. Noller, H. F., Yusupov, M. M., Yusupova, G. Z., Baucom, A. & Cate, J. H. D. Translocation of tRNA during protein synthesis. FEBS Lett. 514, 1116 (2002).

2. Ramakrishnan, V. Ribosome structure and the mechanism of translation. Cell 108, 557572 (2002).

3. Bashan, A. et al. Structural basis of the ribosomal machinery for peptide bond formation, translocation, and nascent chain progression. Mol. Cell 11, 91102 (2003).

4. Frank, J., Gao, H., Sengupta, J., Gao, N. & Taylor, D. J. The process of mRNA-tRNA translocation. Proc. Natl Acad. Sci. USA 104, 1967119678 (2007).

5. Aitken, C. E., Petrov, A. & Puglisi, J. D. Single ribosome dynamics and the mechanism of translation. Annu. Rev. Biophys. 39, 491513 (2010).

6. Yamamoto, H. et al. EF-G and EF4: translocation and back-translocation on the bacterial ribosome. Nat. Rev. Microbiol. 12, 89100 (2014).

7. Valle, M. et al. Locking and unlocking of ribosomal motions. Cell 114, 123134 (2003).

8. Ermolenko, D. N. & Noller, H. F. mRNA translocation occurs during the second step of ribosomal intersubunit rotation. Nat. Struct. Mol. Biol. 18, 457462 (2011).

9. Moazed, D. & Noller, H. F. Intermediate states in the movement of transfer RNA in the ribosome. Nature 342, 142148 (1989).

10. Rodnina, M. V. & Wintermeyer, W. The ribosome as a molecular machine: the mechanism of tRNA-mRNA movement in translocation. Biochem. Soc. Trans. 39, 658662 (2011).

11. Cornish, P. V., Ermolenko, D. N., Noller, H. F. & Ha, T. Spontaneous intersubunit rotation in single ribosomes. Mol. Cell 30, 578588 (2008).

12. Ramrath, D. J. F. et al. Visualization of two transfer RNAs trapped in transit during elongation factor G-mediated translocation. Proc. Natl Acad. Sci. USA 110, 2096420969 (2013).

13. Zhou, J., Lancaster, L., Donohue, J. P. & Noller, H. F. How the ribosome hands the A-site tRNA to the P site during EF-G-catalyzed translocation. Science 345, 11881191 (2014).

14. Ratje, A. H. et al. Head swivel on the ribosome facilitates translocation by means of intra-subunit tRNA hybrid sites. Nature 468, 713716 (2010).15. Spirin, A. S. The ribosome as a conveying thermal ratchet machine. J. Biol. Chem. 284, 2110321119 (2009).

16. Bryngelson, J. D. & Wolynes, P. G. Intermediates and barrier crossing in a random energy-model (with applications to protein folding). J. Phys. Chem. 93, 69026915 (1989).

17. Onuchic, J. N., Luthey-Schulten, Z. & Wolynes, P. G. Theory of protein folding: the energy landscape perspective. Annu. Rev. Phys. Chem. 48, 545600 (1997).

18. Thirumalai, D. & Hyeon, C. RNA and protein folding: common themes and variations. Biochemistry 44, 49574970 (2005).

19. Hyeon, C. & Thirumalai, D. Capturing the essence of folding and functions of biomolecules using coarse-grained models. Nat. Commun. 2, 487 (2011).

20. Chan, H. S., Zhang, Z., Wallin, S. & Liu, Z. Cooperativity, local-nonlocal coupling, and nonnative interactions: principles of protein folding from coarse-grained models. Annu. Rev. Phys. Chem. 62, 301326 (2011).

21. Frauenfelder, H., Sligar, S. G. & Wolynes, P. G. The energy landscapes and motions of proteins. Science 254, 15981603 (1991).

22. Savir, Y. & Tlusty, T. The ribosome as an optimal decoder: a lesson in molecular recognition. Cell 153, 471479 (2013).

23. Whitford, P. C., Blanchard, S. C., Cate, J. H. D. & Sanbonmatsu, K. Y. Connecting the kinetics and energy landscape of tRNA translocation on the ribosome. PLoS Comput. Biol. 9, e1003003 (2013).

24. Dashti, A. et al. Trajectories of the ribosome as a Brownian nanomachine. Proc. Natl Acad. Sci. USA 111, 1749217497 (2014).

25. Noel, J. K., Chahine, J., Leite, V. B. P. & Whitford, P. C. Capturing transition paths and transition states for conformational rearrangements in the ribosome. Biophys. J. 107, 28812890 (2014).

26. Bock, L. V. et al. Energy barriers and driving forces in tRNA translocation through the ribosome. Nat. Struct. Mol. Biol. 20, 13901396 (2013).

27. Trobro, S. & Aqvist, J. Role of ribosomal protein L27 in peptidyl transfer. Biochemistry 47, 48984906 (2008).

28. Gindulyte, A. et al. The transition state for formation of the peptide bond in the ribosome. Proc. Natl Acad. Sci. USA 103, 1332713332 (2006).

29. Adamczyk, A. J. & Warshel, A. Converting structural information into an allosteric-energy-based picture for elongation factor Tu activation by the ribosome. Proc. Natl Acad. Sci. USA 108, 98279832 (2011).

30. Alexander, R. W., Eargle, J. & Luthey-Schulten, Z. Experimental and computational determination of tRNA dynamics. FEBS Lett. 584, 376386 (2010).

31. Trabuco, L. G. et al. The role of L1 stalk-tRNA interaction in the ribosome elongation cycle. J. Mol. Biol. 402, 741760 (2010).

32. Garcia, A. E. & Paschek, D. Simulation of the pressure and temperature folding/unfolding equilibrium of a small rna hairpin. J. Am. Chem. Soc. 130, 815817 (2008).

33. Vaiana, A. C. & Sanbonmatsu, K. Y. Stochastic Gating and Drug-Ribosome Interactions. J. Mol. Biol. 386, 648661 (2009).

34. Satpati, P., Sund, J. & Aqvist, J. Structure-based energetics of mRNA decoding on the ribosome. Biochemistry 53, 17141722 (2014).

35. Chen, A. A. & Garca, A. E. High-resolution reversible folding of hyperstable RNA tetraloops using molecular dynamics simulations. Proc. Natl Acad. Sci. USA 110, 1682016825 (2013).

36. Bergonzo, C., Henriksen, N. M., Roe, D. R. & Cheatham, T. E. Highly sampled tetranucleotide and tetraloop motifs enable evaluation of common RNA force elds. RNA 21, 114 (2015).

37. Wang, Y., Rader, A. J., Bahar, I. & Jernigan, R. L. Global ribosome motions revealed with elastic network model. J. Struct. Biol. 147, 302314 (2004).

38. Trylska, J., Tozzini, V. & McCammon, J. A. Exploring global motions and correlations in the ribosome. Biophys. J. 89, 14551463 (2005).

39. Kurkcuoglu, O., Kurkcuoglu, Z., Doruker, P. & Jernigan, R. L. Collective dynamics of the ribosomal tunnel revealed by elastic network modeling. Proteins 75, 837845 (2009).

40. Whitford, P. C. et al. An all-atom structure-based potential for proteins: Bridging minimal models with all-atom empirical forceelds. Proteins 75, 430441 (2009).

41. Zhou, J., Lancaster, L., Donohue, J. P. & Noller, H. F. Crystal structures of EF-G-ribosome complexes trapped in intermediate states of translocation. Science 340, 1236086 (2013).

42. Gao, Y.-G. et al. The structure of the ribosome with elongation factor G trapped in the posttranslocational state. Science 326, 694699 (2009).

43. Whitford, P. C. et al. Excited states of ribosome translocation revealed through integrative molecular modeling. Proc. Natl Acad. Sci. USA 108, 1894318948 (2011).

44. Whitford, P. C. & Sanbonmatsu, K. Y. Simulating movement of tRNA through the ribosome during hybrid-state formation. J. Chem. Phys. 139, 121919 (2013).

45. Jenner, L. B., Demeshkina, N., Yusupova, G. & Yusupov, M. Structural aspects of messenger RNA reading frame maintenance by the ribosome. Nat. Struct. Mol. Biol. 17, 555560 (2010).

46. Ramrath, D. J. F. et al. The complex of tmRNA-SmpB and EF-G on translocating ribosomes. Nature 485, 526529 (2012).

47. Schuwirth, B. S. et al. Structures of the bacterial ribosome at 3.5 resolution. Science 310, 827 (2005).

48. Cukras, A. R., Southworth, D., Brunelle, J., Culver, G. & Green, R. Ribosomal proteins S12 and S13 function as control elements for translocation of the mRNA: tRNA complex. Mol. Cell 12, 321328 (2003).

49. Cukras, A. R. & Green, R. Multiple effects of S13 in modulating the strength of intersubunit interactions in the ribosome during translation. J. Mol. Biol. 349, 4759 (2005).

8 NATURE COMMUNICATIONS | 7:10586 | DOI: 10.1038/ncomms10586 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10586 ARTICLE

50. Frank, J. et al. The role of tRNA as a molecular spring in decoding, accommodation, and peptidyl transfer. FEBS Lett. 579, 959962 (2005).

51. Miyashita, O., Onuchic, J. N. & Wolynes, P. G. Nonlinear elasticity, proteinquakes, and the energy landscapes of functional transitions in proteins. Proc. Natl Acad. Sci. USA 100, 1257012575 (2003).

52. Wang, L., Altman, R. B. & Blanchard, S. C. Insights into the molecular determinants of ef-g catalyzed translocation. RNA 17, 21892200 (2011).

53. Lin, J., Gagnon, M. G., Bulkley, D. & Steitz, T. A. Conformational changes of elongation factor G on the ribosome during tRNA translocation. Cell 160, 219227 (2015).

54. Brilot, A. F., Korostelev, A. A., Ermolenko, D. N. & Grigorieff, N. Structure of the ribosome with elongation factor G trapped in the pretranslocation state. Proc. Natl Acad. Sci. USA 110, 2099420999 (2013).

55. Cheung, M. S., Garca, A. E. & Onuchic, J. N. Protein folding mediated by solvation: water expulsion and formation of the hydrophobic core occur after the structural collapse. Proc. Natl Acad. Sci. USA 99, 685690 (2002).

56. Ferguson, A., Liu, Z. & Chan, H. S. Desolvation barrier effects are a likely contributor to the remarkable diversity in the folding rates of small proteins.J. Mol. Biol. 389, 619636 (2009).57. Dias, C. L. & Chan, H. S. Pressure-dependent properties of elementary hydrophobic interactions: ramications for activation properties of protein folding. J. Phys. Chem. B 118, 74887509 (2014).

58.poner, J. et al. Nature of base stacking: reference quantum-chemical stacking energies in ten unique B-DNA base-pair steps. Chem. Eur. J. 12, 28542865 (2006).

59. Clementi, C., Nymeyer, H. & Onuchic, J. N. Topological and energetic factors: What determines the structural details of the transition state ensemble and en-route intermediates for protein folding? An investigation for small globular proteins. J. Mol. Biol. 298, 937953 (2000).

60. Lindahl, E., Hess, B. & van der Spoel, D. Gromacs 3.0: A package for molecular simulation and trajectory analysis. J. Mol. Model. 7, 306317 (2001).

61. Hess, B., Kutzner, C., van der Spoel, D. & Lindahl, E. GROMACS 4: algorithms for highly efcient, load-balanced, and scalable molecular simulation. J. Chem. Theory Comput. 4, 435447 (2008).

62. Noel, J. K., Whitford, P. C., Sanbonmatsu, K. Y. & Onuchic, J. N. SMOG@ctbp: simplied deployment of structure-based models in GROMACS. Nucleic Acids Res. 38, W657W661 (2010).

63. Kouza, M., Li, M. S., OBrien, E. P., Hu, C.-K. & Thirumalai, D. Effect of nite size on cooperativity and rates of protein folding. J. Phys. Chem. A 110, 671676 (2006).

64. Whitford, P. C. et al. Accommodation of aminoacyl-tRNA into the ribosome involves reversible excursions along multiple pathways. RNA 16, 11961204 (2010).

65. Humphrey, W., Dalke, A. & Schulten, K. VMD: visual molecular dynamics.J. Mol. Graph. 14, 3338 (1996).

Acknowledgements

This work was supported by a National Science Foundation CAREER Award (grant No. MCB-1350312). We are also grateful for technical support from Nilay Roy, as well as for the generous computing resources provided by the Northeastern University Discovery Cluster.

Author contributions

K.N. performed the simulations. K.N. and P.C.W. analysed the data and prepared the manuscript. P.C.W. conceived and designed the study.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Nguyen, K. and Whitford P. C. steric interactions lead to collective tilting motion in the ribosome during mRNAtRNA translocation.

Nat. Commun. 7:10586 doi: 10.1038/ncomms10586 (2016).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

NATURE COMMUNICATIONS | 7:10586 | DOI: 10.1038/ncomms10586 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9