Next Article in Journal
Selection of Unique Molecules for Cancer Treatment by Distance-Based Method: Hypericin Effect on Respiratory Chain
Previous Article in Journal
SARS-CoV-2 Spike-Derived Peptides Presented by HLA Molecules
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Dynamics of Subunit Rotation in a Eukaryotic Ribosome

by
Frederico Campos Freitas
1,
Gabriele Fuchs
2,
Ronaldo Junio de Oliveira
1 and
Paul Charles Whitford
3,4,*
1
Laboratório de Biofísica Teórica, Departamento de Física, Instituto de Ciências Exatas, Naturais e Educação, Universidade Federal do Triângulo Mineiro, Uberaba, MG 38064-200, Brazil
2
Department of Biological Sciences, The RNA Institute, University at Albany, 1400 Washington Ave, Albany, NY 12222, USA
3
Department of Physics, Northeastern University, 360 Huntington Ave, Boston, MA 02115, USA
4
Center for Theoretical Biological Physics, Northeastern University, 360 Huntington Ave, Boston, MA 02115, USA
*
Author to whom correspondence should be addressed.
Submission received: 19 April 2021 / Revised: 16 May 2021 / Accepted: 18 May 2021 / Published: 24 May 2021

Abstract

:
Protein synthesis by the ribosome is coordinated by an intricate series of large-scale conformational rearrangements. Structural studies can provide information about long-lived states, however biological kinetics are controlled by the intervening free-energy barriers. While there has been progress describing the energy landscapes of bacterial ribosomes, very little is known about the energetics of large-scale rearrangements in eukaryotic systems. To address this topic, we constructed an all-atom model with simplified energetics and performed simulations of subunit rotation in the yeast ribosome. In these simulations, the small subunit (SSU; ∼1 MDa) undergoes spontaneous and reversible rotation events (∼8 ). By enabling the simulation of this rearrangement under equilibrium conditions, these calculations provide initial insights into the molecular factors that control dynamics in eukaryotic ribosomes. Through this, we are able to identify specific inter-subunit interactions that have a pronounced influence on the rate-limiting free-energy barrier. We also show that, as a result of changes in molecular flexibility, the thermodynamic balance between the rotated and unrotated states is temperature-dependent. This effect may be interpreted in terms of differential molecular flexibility within the rotated and unrotated states. Together, these calculations provide a foundation, upon which the field may begin to dissect the energetics of these complex molecular machines.

1. Introduction

The ribosome is a massive molecular assembly that undergoes a wide range conformational rearrangements in order to accurately synthesize proteins [1,2,3,4]. The precise composition of a ribosome is organism-specific, though it generally contains two large RNA (rRNA) chains (thousands of nucleotides in length, each), in addition to a variable number of smaller RNA and protein molecules (∼50–100, in total). The overall architecture of the ribosome is commonly described in terms of the large subunit (LSU) and small subunit (SSU), where each has three distinct tRNA binding sites (Figure 1). During protein synthesis, aminoacyl-transfer RNA (aa-tRNA) molecules must decode messenger RNA (mRNA). Upon recognition of an mRNA codon, the incoming aa-tRNA molecule binds the ribosomal A site. The nascent protein chain that is attached to the P-site tRNA is then passed to the A-site tRNA through the formation of a peptide bond. After peptide bond formation, the A-site and P-site tRNA molecules are displaced to the P and E sites, respectively. This process, called translocation, leads to a vacant A site, which allows the ribosome to decode the next mRNA codon.
The process of tRNA translocation involves large-scale displacements of tRNA molecules (∼40 Å), which are facilitated by an elaborate sequence of collective rearrangements within the ribosome. In bacteria, countless cryo-EM and crystallographic structures have been resolved in which the SSU is rotated relative to the LSU (i.e., body rotation [5,6]), or the SSU head is rotated relative to the SSU body (head swivel [7,8] and tilting [9]). To complement these structural snapshots, biochemical measurements [10] and single-molecule studies [11,12,13,14,15] have identified coupling between SSU rearrangements and tRNA dynamics. Single-molecule measurements have also suggested aspects of hybrid-state formation are similar in eukaryotic ribosomes [16,17]. Similar to bacterial ribosomes, structural models of eukaryotic cytoplasmic [18,19] and mitochondrial [20] ribosomes have illustrated the broad range of orientations that are accessible to the SSU. In addition to rotary-like rearrangements, cryo-EM structures have also revealed that the SSU body may undergo tilt-like rearrangements in eukaryotic systems (called “rolling” [19]). These more recent insights into eukaryotic structure raise new questions into the relationship between SSU motion and tRNA dynamics. For example, how do differences in eukaryotic and bacterial ribosome structure give rise to differential dynamics? What is the relationship between SSU rotation and eukaryotic-specific tilting/rolling? While there is significant interest in understanding these dynamic properties, the biophysical features that govern eukaryotic translation are largely unexplored.
For nearly two decades, advances in structure determination have fueled the development and application of theoretical models to study subunit rotation in bacteria. In the earliest theoretical efforts, coarse-grained models were utilized to perform normal mode analysis [21] and principal component analysis [22], which illustrated how the architecture of the ribosome predisposes it to rotation-like fluctuations of the SSU. Later, highly-detailed explicit-solvent simulations (100 ns–1 μs) were applied to study small-scale structural fluctuations [23,24], which reinforced the predicted energetic accessibility of rotary-like fluctuations. Explicit-solvent simulations have also been applied to characterize how SSU-LSU bridge interactions may facilitate large-scale rotation [25]. Complementary to explicit-solvent methods, recently-developed coarse-grained models have allowed for spontaneous rotation events to be simulated [26]. With those models, it has been shown how specific subunit bridges can “hand off” the SSU during rotation. While each effort has provided insights into distinct aspects of rotation in bacterial ribosomes, simulations of spontaneous and reversible SSU rotation in a eukaryotic ribosome have not been reported previously.
In the current study, we provide a physicochemical foundation for understanding the dynamics of SSU rotation in eukaryotic ribosomes. Specifically, we developed and applied an all-atom structure-based (SMOG) model [27,28] to study the dynamics of subunit rotation in the yeast ribosome. This model includes all non-hydrogen atoms (206 k atoms, in total), and the energetics are defined to explicitly stabilize the rotated and unrotated conformations. Using this simplified model, we were able to simulate 25 reversible rotation events, where the SSU spontaneously rotated/back-rotated by ∼8 . With this data set, we provide an initial description of the free-energy landscape associated with SSU rotation. This analysis reveals a distinct sequence of rearrangements at the SSU-LSU interface, as well as a pronounced temperature dependence of the free-energy landscape. Together, these calculations establish a technical and conceptual framework for rigorously characterizing eukaryotic ribosome dynamics, both through theoretical and experimental approaches.

2. Materials and Methods

2.1. Multi-Basin Structure-Based Model

In the current study, we developed a multi-basin structure-based model, where knowledge of the rotated and unrotated conformations (PDB: 3J77 and 3J78 [18]) were used to define the potential energy function. For this, we applied a SMOG-AMBER variant [29] of the SMOG class of structure-based models [28], where the potential energy is given by:
V = b o n d s ϵ r 2 r i j r i j , 0 2 + a n g l e s ϵ θ 2 θ i j k θ i j k , 0 2 + + i m p r o p e r s ϵ χ i m p 2 χ i j k l χ i j k l , 0 2 + p l a n a r ϵ χ p l a n a r 2 F P φ i j k l + + b a c k b o n e ϵ b b F D ϕ i j k l ϕ i j k l , 0 + s i d e c h a i n s ϵ s c F D ϕ i j k l ϕ i j k l , 0 + + c o n t a c t s ϵ C σ i j r i j 12 2 σ i j r i j 6 + n o n c o n t a c t s ϵ n c σ n c r i j 12 .
Here, F D ( ϕ ) = [ 1 cos ( ϕ ) ] + 1 2 [ 1 cos ( 3 ϕ ) ] and F P φ = [ 1 cos ( 2 φ ) ] . The bonded parameters ( r i j , 0 and θ i j k , 0 ) were obtained from the AMBER03 force field [30]. The planar ring dihedrals were maintained by cosine potentials of periodicity 2. The position of the minimum of each dihedral angle ϕ i j k l , 0 was defined as the mean value adopted in the rotated and unrotated structures. This ensures that the dihedral energies in the two structures are isoenergetic. Combined with the AMBER-based bonded geometry, these terms ensure there is no intrasubunit bias toward either endpoint configuration. Dihedral energies were assigned as described previously [28]. For completeness, we will summarize the details, here. To define the dihedral interaction weights ( ϵ b b and ϵ s c ), dihedrals were first grouped based on the composition of the middle bond. Each dihedral group was given a summed weight of ϵ b b or ϵ s c . The ratio R b b / s c = ϵ b b ϵ s c was set to 1 (nucleic acids) or 2 (proteins). ϵ b b was defined to be equal for protein and nucleic acids. As a note, it was previously shown that folding mechanisms and free-energy barriers are largely robust for 1 R b b / s c 4 [27]. Contact and dihedral strengths were scaled, such that
R C / D = ϵ C ϵ b b + ϵ s c = 2
and
ϵ C + ϵ b b + ϵ s c = N ϵ ,
where N is the number of atoms in the system and ϵ is the reduced energy unit, which is equal to 2 k B T . Contact pairs were defined using the Shadow Contact Map algorithm with default values [31]. σ i j was set to 0.96 σ i j where σ i j is the interatomic distance between atoms that are in contact in the rotated (or unrotated) configuration. This scaling of contacts was introduced to avoid the artificial expansion of the ribosome that can arise from configurational entropy [29,32].
To construct the multi-basin force field, any inter-subunit contact pair that is common to both structures was assigned a single 6–12 interaction. The position of the minimum was defined for each atom pair, such that each common contact is isoenergetic (with value ϵ iso ), with respect to the rotated and unrotated conformations (Figure S1). If ϵ iso / ϵ C > 1 / 2 , the contacting pair of atoms was classified as a “common” contact. Applying this criterion, all inter-subunit pairs whose distances are similar in both structures were assigned an isoenergetic distance. All intra-subunit contact pairs were also considered common contacts. Consistent with the dihedral parameters, this isoenergetic assignment of contacts ensures that common interactions do not favor either rotation state. The weights of common contacts (inter-subunit and intra-subunit) followed the scaling rules described above [28]. For the remaining interface contacts (i.e., unique contacts), σ i j was assigned the value found in the conformation for which the contact is defined: (un)rotated contacts are given the distances found in the (un)rotated conformation. Contacts unique to the rotated configuration were then given an energetic weight of 0.21 and contacts unique to the unrotated configuration were given weights 0.19. An initial parameter sweep was performed to identify weights for which the rotated and unrotated conformations represent pronounced free-energy minima of comparable depths.
Since the available cryo-EM reconstructions resolved different numbers of atoms, minor components of the structure had to be structurally modeled. Consistent with previous efforts [32], position-restraint-based techniques were applied to model the missing regions, where the completed region in the alternate structure was aligned to the missing region. With this approach, residues 65–135 of eL24 were reconstructed in the unrotated structure (3J78), and residues 2061–2075 of the 25S rRNA were reconstructed in the rotated structure. After these steps, both models contained 206,389 non-hydrogen atoms. See Supplementary Materials for details.

2.2. Simulation Details

Simulations were performed using the GROMACS software package (v5.1.4) [33,34]. Force field files were generated using SMOG v2.3 [28] with the SBM_AA-amber-bonds force field templates [29], which are available through the smog-server.org force field repository. The contacts and dihedrals were then subsequently modified by custom scripts. Two simulations were initiated from the unrotated structure (PDB: 3J78 [18]). To ensure robustness of the results, an additional simulation was initiated from the rotated configuration (PDB: 3J77 [18]). The temperature was maintained via Langevin Dynamics protocols, with a value of 60 in GROMACS units, or a reduced value of 0.49887, for all simulations. Each simulation was continued for at least 7 · 10 9 time steps. Using the estimate of Yang et al. [35], each simulation represents an effective timescale of approximately 15 ms. To allow for equilibration, the first 10 7 time steps were excluded from analysis. The time step was 0.002 reduced units and configurations were saved every 5000 time steps.

2.3. SSU Rotation Measures

To describe the orientation of the SSU, we extended the definition of rotation angles that were originally introduced for bacterial ribosomes [32,36]. Consistent with previous descriptions [32], the overall strategy is to identify sets of residues that are structurally conserved between the rRNA of the LSU and SSU body in E. coli and yeast. Using these structurally-conserved residues, a reference E. coli structure is separately aligned to the LSU and SSU body of the yeast ribosome. This allows one to approximate the position of the SSU body and head in terms of rigid-body rotations, relative to E. coli. The reference E. coli model is PDB entry 4V6D [6], where the axis of pure rotation is defined by the structures of the classical and rotated body. With this definition, we decompose the rotation in terms of Euler angles ( ϕ , θ , ψ ). Here, we report the net body rotation as ϕ body = ϕ + ψ , body tilting/rolling as θ body = θ and the tilt direction as ψ body = ϕ + C . C is an arbitrary constant that ensures ψ body = 0 describes tilting motions that are roughly about the long axis of helix 44 in the SSU body.
In the current study, we extended our protocol for defining angles, in order to apply Euler Angle decomposition to the yeast ribosome, while still allowing for direct comparisons with SSU orientations in other organisms. First, we performed STAMP alignment [37] to determine the corresponding E. coli numbering of the yeast rRNA residues. STAMP was applied separately to the LSU and SSU body. Then, for the LSU and SSU body, we applied a second, more stringent, condition to define which residues in the rotated yeast structure are structurally conserved with the reference E. coli structure. Specifically, we performed least-squares alignment of the STAMP-aligned residues and calculated the spatial deviation of each P atom. Any P atom that deviated by more than 1 Å was then excluded in a second round of fitting. After the second round of alignment, all P atoms (including those that were not used for fitting) that are within (above) 1 Å are included in (excluded from) the fitting group. This process was repeated until the set of included residues converged. These sets of residues, which we will call the “core” groups, were then used for all subsequent angle calculations in the study. All angle calculations were performed using in-house scripts written for use with VMD [38].

3. Results

3.1. Simulating Spontaneous Subunit Rotation Events in a Complete Ribosome

To study the dynamics of subunit rotation for a eukaryotic ribosome, we applied molecular dynamics simulations with an all-atom model (206,389 atoms) that employs a simplified energetic representation. Specifically, we used a multi-basin structure-based model, which is inspired by similar models for the study of multi-domain proteins [39]. In this model, all non-hydrogen atoms are represented, and the stabilizing energetic interactions are explicitly defined to stabilize the rotated and unrotated configurations of the ribosome. Here, we define the interactions based on cryo-EM structures of a yeast ribosome in rotated and unrotated conformations (Figure 1). With this representation, we were able to simulate spontaneous (i.e., without targeting techniques) and reversible transitions between rotated and unrotated conformations (See supplementary movie).
When interpreting the physical significance of the simulated dynamics, one should recognize that the modeled interactions are intended to reflect the effective energetics of the system [40,41,42]. That is, structures that have been resolved necessarily represent free-energy minima. In the structure-based model, we directly encode these free-energy minima by defining contacts and dihedrals to stabilize the pre-assigned structures. This general approach has been applied recently to study rotation in a bacterial ribosome with a coarse-grained model [26]. In the current study, we extend to an all-atom representation, such that the presented model may later be used to study the precise influence of sterics during tRNA rearrangements in the ribosome. This consideration is motivated by previous simulations of bacterial ribosomes, which have demonstrated the critical influence of molecular sterics on tRNA dynamics during accommodation [43], A/P hybrid formation [44], P/E hybrid formation [45] and tRNA translocation [32]. In order for our model to have future utility to address these motions in eukaryotic ribosomes, it is necessary to employ atomic resolution. Here, we focus on rotation in the absence of a bound tRNA molecule.
In order to characterize the dynamics of SSU rearrangements, one must define appropriate collective coordinates that distinguish between rotated and unrotated orientations of the small subunit (SSU). To this end, we employed Euler Angle decomposition, as used previously to describe bacterial ribosomes [32,36]. Here, this method was generalized (see Methods for details) for use with non-bacterial ribosomes. Consistent with previous efforts to quantify rotation angles in the ribosome [32], we first defined sets of residues within the LSU and SSU body that undergo minimal intra-domain rearrangements during rotation (called the “core” residues). We then aligned reference structures of the core residues to the LSU and SSU body. The orientations of the aligned structures were then quantified in terms of Euler angles, which may be expressed as a net body rotation angle ( ϕ body ), a body tilt angle ( θ body ) and a corresponding tilt direction ( ψ body . See Methods for details). Pure body rotation leads to ϕ body 0 , while θ body = 0 , where the rotation axis is parallel to that defined by reference structures of E. coli [6]. Non-zero values of the tilt angle θ body represent any level of deviation from pure rotation. In relation to other studies of eukaryotic ribosomes, the tilt angle measures the so-called “rolling” rearrangement of the SSU body [19]. For consistency with descriptions of SSU head motion [32], we use the term “tilting” to describe such motions. To complete this description of subunit rotation, the direction of tilting is given by ψ body , where ψ body 50 corresponds to the tilting/rolling rearrangement that is apparent in cryo-EM models of the yeast ribosome [18].
Using our structure-based model, we performed several independent simulations, in which a total of 25 rotation and backrotation events were observed (Figure 2C,D). While applying atomic resolution is computationally demanding, the simulated events were obtained without the use of enhanced sampling methods, or artificial targeting forces. Rotation events are apparent in the time traces, where there are sharp changes in ϕ body of approximately 8 (Figure 2C). There are also abrupt transitions in the tilting angle θ body that coincide with changes in the rotation angle. Thus, this model suggests that subunit rotation and tilting/rolling are not kinetically-separable processes. While one could expect two-state-like behavior in this type of simplified model, it is common for these models to reveal sterically-induced free-energy minima [46]. That is, while the models define specific conformations as stable, the imperfect complementarity of steric interactions can impede the motion and lead to long-lived intermediates. Despite this possibility, we find that the steric composition of the SSU-LSU interface appears to be sufficiently smooth that rotation and tilting can occur simultaneously and in a two-state manner.
With regards to tilting dynamics, we find there are large variations in the tilt direction ψ body within the unrotated ensemble (Figure 2D). These large variations are expected, since the tilt direction is undefined when the tilt angle is zero. Accordingly, when θ body is small, structural fluctuations associated with thermal energy can give rise to minimal changes in θ body and large changes in ψ body .
Overall, the presented simulations demonstrate how an all-atom structure-based model can be used to provide a first-order approximation to the dynamics of rotation in a eukaryotic ribosome. Since these simulations describe spontaneous rotation events, the data set provides an opportunity to gain initial insights into the relative timing of inter-subunit contact formation during rotation, as well as the impact of specific interactions on the free-energy landscape.

3.2. Quantifying the Energy Landscape of Rotation

A persistent challenge in molecular biophysics is to define low-dimensional measures that can accurately capture the dynamics of complex multi-dimensional processes [47,48,49,50,51]. In addition to posing an intellectual challenge, there is also a practical utility of identifying appropriate coordinates. Specifically, knowledge of appropriate one-dimensional measures can allow one to precisely characterize the relative contributions of individual interactions to biological kinetics. To this end, we explored multiple approaches for describing the simulated SSU rotation events, which together help establish a physical-chemical foundation for the analysis of eukaryotic ribosome dynamics.
Visual inspection of individual simulated time traces suggests there is a strong correlation between SSU rotation and tilting/rolling (Figure 2C). To better understand the relationship between these motions, we calculated the two-dimensional free energy ( k B T ln ( P ) ) as a function of ϕ body and θ body (Figure 3A). Consistent with the time traces, there are two distinct minima corresponding to the rotated and unrotated ensembles. However, rotation and tilting are not perfectly correlated, and there is a notable degree of variability in θ body within each ensemble. In the unrotated ensemble (low ϕ body values), the tilt angle samples values that range from 1 to 4 . Similarly, in the rotated ensemble (high ϕ body values), θ body spans a wide range of values (2.5–6.5 ). Interestingly, as the body tilts/rolls, the direction of tilting also shifts. Specifically, the direction of tilting ( ψ body ) within the unrotated ensemble is typically around 50 , while the rotated ensemble represents configurations for which the tilt direction is towards higher ψ body values. This illustrates how, upon reaching the rotated/tilted ensemble, the accessible tilting fluctuations change in character. Based on the current simulations, it is not clear whether this change in tilting dynamics has a specific biological role, though it may help coordinate tRNA dynamics during hybrid formation and/or translocation.
Since rotation and tilting are found to be correlated, we asked whether a one-dimensional description can be sufficient to quantify the free-energy barrier and kinetics associated with SSU motion. To determine whether the rotation angle ϕ body provides a kinetically-meaningful approximation to the free-energy barrier, we employed multiple independent analysis strategies. First, we applied transition path analysis [47,52] to assess whether ϕ body can unambiguously identify the transition state ensemble (TSE). That is, we calculated the probability that the system is undergoing a transition (i.e., is on a Transition Path), as a function of ϕ body : P ( TP | ϕ body ) . If ϕ body is able to unambiguously identify the TSE, then P ( TP | ϕ body ) will adopt the diffusion-limited value of 0.5. Here, we find that P ( TP | ϕ body ) reaches a peak value of 0.5 , indicating that ϕ body is suitable for describing the underlying barrier. Further, this allows us to identify the position of the TSE: ϕ body 3 . 5 . As a point of comparison, in the study of protein folding [47,52], a coordinate ρ is often considered “good” if P ( TP | ρ ) reaches values that are greater than ∼0.4.
In addition to capturing the TSE, we find that ϕ body is able to recapitulate the long-time dynamics of SSU rotation. To illustrate this, we used a Bayesian inference approach [52] to calculate the diffusion coefficient as a function of ϕ body . This produced a nearly-uniform value of the diffusion coefficient D that was approximately 0.0025 ( degrees ) 2 / τ ru . τ ru is the reduced time unit. We used this value of D to estimate the mean first passage time via the relation:
τ = ρ initial ρ final d ρ ρ min ρ d ρ exp [ ( F ( ρ ) F ( ρ ) ) / k B T ] D ( ρ ) ,
where F ( ρ ) is the free-energy as a function of ρ = ϕ body . ρ final and ρ initial are the values of ϕ body at the free-energy minima, and ρ min is the smallest accessible value of the coordinate. According to this calculation, the inferred timescale is estimated at 10 6 τ ru . This is comparable to the apparent timescale (total simulated time divided by the number of observed transitions) of 1.8 × 10 6 τ ru . Agreement between the apparent and inferred timescales is consistent with P ( TP | ϕ body ) reaching the diffusion-limited value of 0.5. Interestingly, even when using intuitively-defined coordinates, there is no guarantee they will be able to accurately recapitulated the long-time kinetics. For example, in a recent study of SSU rotation in bacteria [26], it was shown that a similar coordinate is unable to capture the TSE, and the coordinate dramatically underestimates the height of the free-energy barrier. Fortunately, here, our analysis indicates that ϕ body may be used to probe the factors that govern SSU rotation in a eukaryotic ribosome.

3.3. Simulations Implicate Millisecond-Scale Dynamics of Rotation

Since ϕ body can be used to precisely describe the free-energy barrier for rotation, we next compared the predicted free-energy barrier with biological kinetics. For this, we use the expression:
1 / τ = C exp ( Δ F / k B T ) ,
where Δ F is the height of the free-energy barrier and C is the average barrier-crossing attempt frequency. The free-energy barrier along ϕ body is Δ F 6–7 k B T . Here, we used a barrier-crossing attempt frequency for rotation C that was previously estimated from explicit-solvent simulations [23]. This conversion between the free-energy barrier and kinetics implicates a mean first-passage time τ that is in the millisecond regime (1–5 ms). When interpreting this estimate of the timescale, it is important to recall that C was obtained for a bacterial ribosome, which is ∼70% the mass of the yeast ribosome. In accordance with Stokes-Einstein scaling for rotational diffusion, this would lead one to anticipate that rotational diffusion in yeast will be reduced by a factor of roughly 1.5. Further, since C is linearly proportional to the diffusion coefficient [23], the estimated timescale is likely underestimated by a comparable factor. Accordingly, a timescale of 1–5 ms represents a lower-bound on the mean first-passage time for SSU rotation in our model.
Even though the employed model is not intended to provide a complete description of ribosome energetics, the predicted timescale for SSU rotation is generally compatible with the rate of protein synthesis in the cell. The average rate of protein synthesis in yeast has been estimated to be approximately 10 amino acids per second [53], which imposes an upper limit of 100 ms on the timescale of any individual substep. However, there is no evidence that SSU rotation is rate limiting, which would be consistent with rotation kinetics occuring on substantially shorter timescales. In terms of energetics, this empirically-imposed upper bound on the timescale indicates that the barrier height in vivo is necessarily less than 10 k B T (estimated according to ref. [23]). In terms of modeling considerations, a difference in barrier height that is only a few k B T can easily be accounted for by minor changes in the energetic representation. For example, increasing the rotated and unrotated contact strength by ∼2% would increase the barrier by 4 k B T . To make this estimate, we assume that most contacts will primarily stabilize the endpoint conformations. Since there are approximately 500 unique contacts in each endpoint, each of strength ∼0.2 reduced units (1 reduced energy unit is equal to 2 k B T ), an increase of 2% would stabilize the endpoints by 4 k B T . Since the mechanistic aspects of ribosome dynamics are often robust to parameters changes on this scale [54], we will consider the deviations between the predicted and in vivo upper-bound barrier heights to be minimal. Taken together, the compatibility of the timescales for elongation in vivo and the simulated rotation events suggests the model represents an appropriate first-order approximation to the energetics of subunit rotation.

3.4. Asynchronous Dynamics of Subunit Bridge Interactions during Rotation

While applying rigid-body descriptions of biomolecular dynamics is often motivated by analogies with macroscopic systems, the scale of thermal energy in molecular systems gives rise to a fundamentally distinct relationship between structure and dynamics. In the cell, solvent introduces energetic fluctuations that are of the same scale ( k B T ) as the interactions that maintain structural integrity. This leads to heterogeneous and anisotropic structural fluctuations [55] that can manifest in the form of large-scale global motions [22], as well as more localized distortions and partial unfolding [56,57]. In contrast, a corollary of rigid-body arguments would be that all intersubunit bridge interactions should simultaneously interconvert as the SSU rotates. However, such a process would likely be associated with a large free-energy barrier that would lead to prohibitively slow dynamics.
To better understand how molecular flexibility can facilitate collective rotation of the SSU, we used our simulated transitions to explore the ordering of intersubunit bridge rearrangements. Specifically, we calculated the probability that each rotated and unrotated contact is formed: P N i j . N refers to the conformation (rotated, or unrotated) and i j denotes the atom pairs involved in a contact. Here, uniquely rotated/unrotated contacts are defined as atomic interactions that are present (i.e., proximal atoms) in only one of the conformations. Further, to be classified as “unique,” the atom pair must be substantially farther apart (see Methods) in the alternate conformation. We then calculated the average probability of all unique contacts that are defined with a specific atom: P N i = Σ j P N i j N i , where N i is the number of contacts that are defined with atom i. We then calculated P N i as a function of the rotation angle ϕ body , for both the rotated (Figure 4A) and unrotated (Figure 4B) contacts. Since the contact probabilities will be proportional to the stability imparted by a specific interaction, such analysis allows one to categorize interactions in terms of their contributions to the free-energy barrier and biological kinetics.
Contact analysis reveals how specific bridge interactions can have differential effects on the global kinetics of the system. We find that most unrotated contacts break early in the rotation process (Figure 4B), where almost no interactions are still formed by the time the ribosome reaches the TSE ( ϕ body 3 . 5 ). This suggests that, while contacts that are unique to the unrotated conformation can contribute to the stability of the classical configuration, these interactions are unlikely to directly affect the free-energy barrier. In contrast, the rotated contacts form over a broad range of ϕ body values, where some contacts are likely to form prior to the system reaching the TSE (Figure 4A).
To visually depict which interface regions are likely to have the strongest influence on the free-energy barrier, we identified all atoms that have at least one contact that is formed ( P N i j > 0.8 ) in the TSE. Consistent with Figure 4A,B, there are only four unrotated contacts that are likely to be formed (yellow spheres in Figure 4C,D). This is in sharp contrast with the dynamics of the rotated contacts, for which there are clusters of formed contacts scattered across the subunit interface. In particular, there is a dense cluster of formed contacts that are centered around bridge B8, which is formed between protein L23 and SSU helix 14. This suggests that, during rotation, the flexibility of the bridge allows it to “reach out” and form rotated contacts before the SSU body has fully transitioned to a rotated configuration. This finding suggests new ways in which experiments may modulate rotation kinetics in eukaryotic ribosomes. For example, it may be possible to introduce mutations to protein L23 that will specifically impact the rotated and TS ensembles, while leaving the stability of the classical configuration unperturbed. Together, these calculations provide a physical framework that can guide the development of next-generation experimental techniques that will be able to control biological dynamics.

3.5. Molecular Flexibility Leads to Temperature-Induced Population Shift

In addition to characterizing the transition-state ensemble and the dynamics of SSU-LSU bridge interactions, we next asked whether the balance between rotated and unrotated ensembles is likely to be temperature-dependent. To explore this possibility, we calculated the relative free energy of the rotated and unrotated ensembles as a function of temperature. We find that increases in temperature are associated with increased stability of the rotated ensemble (Figure 5A). As described below, one can understand this effect as arising from an increase in mobility of a localized region within the LSU upon rotation.
To probe the temperature dependence of the energy landscape, we pursued a free-energy perturbation (FEP) approach. That is, we calculated the free-energy at arbitrary temperatures according to: F ( ϕ body ) = k B T ln ( P ( ϕ body ) ) , where P ( ϕ body ) represents a probability distribution for which each simulated frame was assigned a weight:
W = exp 1 k B T 1 k B T 0 U
U is the potential energy of the simulated snapshot, T 0 is the temperature of the simulation, and T is the temperature of interest. When T 0 = T (i.e., original simulated distribution), W reduces to 1. Before discussing the results obtained from the FEP calculations, it is important to explain the rationale and assumptions that are implicit in these types of methods. The first consideration is the scale of the simulations. That is, even though simulations with the simplified model are faster (many orders of magnitude) than simulations with explicit-solvent models, the size (>200,000 atoms) and effective timescale (milliseconds) remains very computationally demanding (∼10 10 simulated timesteps, per simulation). Thus, it is not practical to obtain equilibrium sampling for a range of temperatures. The second consideration that supports the use of FEP is that the distributions calculated from different simulations were highly reproducible (not shown). This suggests the acquired sampling provides adequate coverage of the local phase space, which is necessary for FEP techniques to provide reliable estimates. The final consideration is that we are interested in temperature effects over a relatively small temperature range. That is, we consider a temperature range of ± 2 % , where dramatic changes in the accessible conformations of the system are not likely to occur. This is an important point, since one may only extrapolate free energies using perturbation techniques if the perturbed parameter/s (i.e., temperature) would primarily redistribute the probabilities of well-sampled regions of phase space. Since the presented simulations were performed at a temperature for which only modest changes in flexibility are expected with this model [58], it is reasonable to assume that the accessible range of configurations will not be altered dramatically by small changes in temperature.
We find that the rotated ensemble is favored as temperature is increased (Figure 5A). For this, the rotated and unrotated ensembles were defined as 5.4 < ϕ body < 7 . 86 and 1.12 < ϕ body < 0 . 82 . The difference in free energy of the two ensembles was then defined as Δ F = F rot F unrot = k B T ln ( P unrot / P rot ) , where F N is the free energy of ensemble N. For the reported simulations, the free energy of the two ensembles was comparable (Figure 3D), where the relative stability of the rotated ensemble is found to increase with temperature. This suggests that differences in the configurational entropy of the rotated and unrotated ensembles can lead to a distinct temperature-dependence of the distribution between these states. Since the current model does not provide an explicit treatment of the solvent, this predicted trend is due solely to the contributions of the configurational entropy of the ribosome. Based on available data, it is not known whether solvation entropy will provide a significant contribution to this free-energy difference. If future experiments corroborate the observed trend, then one may infer that configurational entropy considerations are sufficient. If an opposite trend were found in experiments, this would suggest that solvation entropy is the dominant contributor. However, since the majority of the SSU-LSU interface is maintained in the rotated and unrotated states, it is reasonable to expect that solvation entropy changes will be minimal.
To provide a structural interpretation for the origins of the temperature-dependent distributions, we considered the flexibility of the ribosome in each ensemble. To assess the differential flexibility of the ribosome, we calculated the spatial root mean-squared fluctuations (rmsf) of each atom, where the average was calculated separately for the unrotated and rotated (Figure 5C) ensembles. In both cases, the rmsf values are highly heterogeneous, where the stalk regions are most flexible, along with some peripheral protein tails. To identify which regions of the structure contribute to the observed temperature-dependent dynamics, we calculated the difference in rmsf values Δ rmsf = rmsf rot rmsf unrot for each atom (Figure 5D). This reveals that the differences in flexibility within the ensembles are primarily centered around a small region of H63 (expansion segment 27b, ES27b: red in Figure 5D). Interestingly, there is also one region in the LSU rRNA involved in the formation of a bridge interaction (ES31a in H79; blue in Figure 5D ) that exhibits a decrease in mobility, though the attenuation of mobility is of a smaller scale than the increase in mobility of ES27b.
In summary, as the ribosome rotates, the ES27b region can dissociate from the SSU interface and adopt a wide range of configurations. This increased mobility can then entropically drive the ribosome towards the rotated ensemble. At a qualitative level, this type of behavior of reminiscent of the notion of an entropic spring in polymer physics, where extension-associated reductions in entropy lead to an effective contractile force. There have also been similar examples in the study of protein conformational rearrangements, where increased mobility has been suggested to act as an entropic counterweight [59]. Accordingly, the current study helps elucidate the interplay between flexibility and dynamics in the ribosome, while revealing a common theme with simpler molecular systems.

4. Discussion

While structural methods can reveal atomic details of the ribosome at various stages of function, biological kinetics are controlled by energetics. Due to the complex character of ribosomal motions, directly quantifying the energetics through experiments has proven to be extremely difficult. Thus, there is need for a quantitative physical foundation, with which precise experimental measurements may be devised and executed. In the case of bacterial ribosomes, there are various single-molecule strategies available for studying large-scale conformational changes [12,14,15,60,61,62], where each can provide a unique perspective on the dynamics. To help integrate these complementary insights, molecular simulation represents a powerful approach for reconciling apparent differences within a common framework [26,63] and guiding the development of next-generation experiments. That is, rather than attempting to predict the exact behavior of a complex assembly, simulations may be used to establish broad trends and relationships. In regards to the current study, the presented simulations could be used to identify single-molecule probes that can precisely monitor rotation in eukaryotic ribosomes. These types of analyses can suggest experimental strategies that will be able to isolate the factors that control biological dynamics in the cell.
In the presented study, we illustrate how molecular simulations of a eukaryotic ribosome can provide initial insights into the relationship between molecular flexibility and kinetics. For example, we find that rotation is best described in terms of asynchronous motions that are facilitated by heterogeneous flexibility of the ribosome. As another example, we find that the flexibility of a specific helical region (ES27b) can increase upon rotation. This increase in mobility is accompanied by an increase in configurational entropy that can drive the rotation process. In terms of experiments, these predictions suggest that it may be possible to alter the dynamics of large-scale processes by introducing localized modifications to the ribosome (e.g., mutations, small-molecular binding) that impact flexibility. In future studies, it will be interesting to see the many ways in which flexibility can help orchestrate the dynamics of these assemblies.
With the ever-increasing availability of computational facilities, the ribosome is now becoming a model system for exploring theoretical concepts in biomolecular dynamics. That is, while performing a ribosome simulation used to represent a major technical accomplishment, the field is now entering a stage where the primary challenge is to craft pointed questions and deploy suitable models that can address them. This stage of development is reminiscent of the protein folding field in the late 1990s and early 2000s. At that time, there was an endless stream of proposed theoretical models, where each could be tested by applying simulation techniques. 20 years later, similar approaches are becoming possible for large assemblies, such as the ribosome. Building on current efforts, we anticipate that the continued integration of experiments, theoretical concepts and simulation techniques will allow for the identification of the precise molecular factors that control complex biomolecular assemblies.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/biophysica1020016/s1, Figure S1: Lennard-Jones potential used to show the isoenergetic distance definition and Figure S2: Probability of contact formation for run, as a function of SSU body angle; The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/biophysica1020016/s2, Movie S1: All-atom Simulation of Subunit Rotation in a Yeast Ribosome. References [18,27,28,30,31,32,64] are cited in the Supplementary Materials.

Author Contributions

Conceptualization, P.C.W. and G.F.; methodology, F.C.F., R.J.d.O., P.C.W.; formal analysis, all authors; writing—original draft preparation, F.C.F. and P.C.W.; writing—review and editing, all authors. All authors have read and agreed to the published version of the manuscript.

Funding

PCW was supported by NSF grant MCB-1915843. Work at the Center for Theoretical Biological Physics was also supported by the NSF (Grant PHY-2019745). FCF was financed by the Coodenação de Aperfeiçoamento de Pessoal de Nível Superior—Brasil (Capes)—Finance Code 001. Financial support for RJO was provided by Fundação de Amparo à Pesquisa do Estado de Minas Gerais (FAPEMIG, APQ-00941-14) and Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq, 438316/2018-5 and 312328/2019-2). GF was supported by NSF grant MCB-2047629 and NIH R03 AI144839.

Data Availability Statement

All simulated data is available upon reasonable request.

Acknowledgments

We would like to acknowledge generous support from the Northeastern University Discovery cluster and Northeastern University Research Computing staff. We also thank Andrei Korostelev for providing the model with improved stereochemistry.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
rmsfroot mean-squared fluctuations
LSUlarge subunit
SSUsmall subunit
mfptmean first-passage time

References

  1. Korostelev, A.; Noller, H.F. The ribosome in focus: New structures bring new insights. Trends Biochem. Sci. 2007, 32, 434–441. [Google Scholar] [CrossRef] [PubMed]
  2. 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] [Green Version]
  3. Schmeing, T.M.; Ramakrishnan, V. What recent ribosome structures have revealed about the mechanism of translation. Nature 2009, 461, 1234–1242. [Google Scholar] [CrossRef] [PubMed]
  4. 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] [Green Version]
  5. Valle, M.; Zavialov, A.; Sengupta, J.; Rawat, U.; Ehrenberg, M.; Frank, J. Locking and unlocking of ribosomal motions. Cell 2003, 114, 123–134. [Google Scholar] [CrossRef] [Green Version]
  6. Dunkle, J.; Wang, L.; Feldman, M.B.; Pulk, A.; Chen, V.B.; Kapral, G.J.; Noeske, J.; Richardson, J.S.; Blanchard, S.C.; Cate, J.H.D. Structures of the bacterial ribosome in classical and hybrid states of tRNA binding. Science 2011, 332, 981–984. [Google Scholar] [CrossRef] [PubMed]
  7. Ratje, A.; Loerke, J.; Mikolajka, A.; Brünner, M.; Hildebrand, P.W.; Starosta, A.; Dönhöfer, A.; Connell, S.; Fucini, P.; Mielke, T.; et al. Head swivel on the ribosome facilitates translocation by means of intra-subunit tRNA hybrid sites. Nature 2010, 468, 713–716. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Mohan, S.; Donohue, J.P.; Noller, H.F. Molecular mechanics of 30S subunit head rotation. Proc. Natl. Acad. Sci. USA 2014, 111, 13325–13330. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Ramrath, D.J.F.; Yamamoto, H.; Rother, K.; Wittek, D.; Pech, M.; Mielke, T.; Loerke, J.; Scheerer, P.; Ivanov, P.; Teraoka, Y.; et al. The complex of tmRNA–SmpB and EF-G on translocating ribosomes. Nature 2013, 485, 526–529. [Google Scholar] [CrossRef] [PubMed]
  10. Horan, L.H.; Noller, H.F. Intersubunit movement is required for ribosomal translocation. Proc. Natl. Acad. Sci. USA 2007, 104, 4881–4885. [Google Scholar] [CrossRef] [Green Version]
  11. Ermolenko, D.N.; Majumdar, Z.K.; Hickerson, R.P.; Spiegel, P.C.; Clegg, R.M.; Noller, H.F. Observation of intersubunit movement of the ribosome in solution using FRET. J. Mol. Biol. 2007, 370, 530–540. [Google Scholar] [CrossRef]
  12. Cornish, P.V.; Ermolenko, D.N.; Noller, H.F.; Ha, T. Spontaneous intersubunit rotation in single ribosomes. Mol. Cell 2008, 30, 578–588. [Google Scholar] [CrossRef] [Green Version]
  13. Cornish, P.V.; Ermolenko, D.N.; Staple, D.W.; Hoang, L.; Hickerson, R.P.; Noller, H.F.; Ha, T. Following movement of the L1 stalk between three functional states in single ribosomes. Proc. Natl. Acad. Sci. USA 2009, 106, 2571–2576. [Google Scholar] [CrossRef] [Green Version]
  14. Fei, J.; Bronson, J.E.; Hofman, J.M.; Srinivas, R.L.; Wiggins, C.H.; Gonzalez, R.L. Allosteric collaboration between elongation factor G and the ribosomal L1 stalk directs tRNA movements during translation. Proc. Natl. Acad. Sci. USA 2009, 106, 15702–15707. [Google Scholar] [CrossRef] [Green Version]
  15. Marshall, R.; Dorywalska, M.; Puglisi, J. Irreversible chemical steps control intersubunit dynamics during translation. Proc. Nat. Acad. Sci. USA 2008, 105, 15364–15369. [Google Scholar] [CrossRef] [Green Version]
  16. Ferguson, A.; Wang, L.; Altman, R.B.; Terry, D.S.; Juette, M.F.; Burnett, B.J.; Alejo, J.L.; Dass, R.A.; Parks, M.M.; Vincent, C.T.; et al. Functional Dynamics within the Human Ribosome Regulate the Rate of Active Protein Synthesis. Mol. Cell 2015, 60, 475–486. [Google Scholar] [CrossRef] [Green Version]
  17. Flis, J.; Holm, M.; Rundlet, E.J.; Loerke, J.; Hilal, T.; Dabrowski, M.; Bürger, J.; Mielke, T.; Blanchard, S.C.; Spahn, C.M.T.; et al. tRNA Translocation by the Eukaryotic 80S Ribosome and the Impact of GTP Hydrolysis. Cell Rep. 2018, 25, 2676–2688.e7. [Google Scholar] [CrossRef] [Green Version]
  18. Svidritskiy, E.; Brilot, A.F.; Koh, C.S.; Grigorieff, N.; Korostelev, A.A. Structures of yeast 80S ribosome-tRNA complexes in the rotated and nonrotated conformations. Structure 2014, 22, 1210–1218. [Google Scholar] [CrossRef] [Green Version]
  19. Budkevich, T.V.; Giesebrecht, J.; Behrmann, E.; Loerke, J.; Ramrath, D.J.; Mielke, T.; Ismer, J.; Hildebrand, P.W.; Tung, C.S.; Nierhaus, K.H.; et al. Regulation of the mammalian elongation cycle by subunit rolling: A eukaryotic-specific ribosome rearrangement. Cell 2014, 158, 121–131. [Google Scholar] [CrossRef] [Green Version]
  20. Itoh, Y.; Naschberger, A.; Mortezaei, N.; Herrmann, J.M.; Amunts, A. Analysis of translating mitoribosome reveals functional characteristics of translation in mitochondria of fungi. Nat. Commun. 2020, 11, 1–10. [Google Scholar] [CrossRef]
  21. Tama, F.; Valle, M.; Frank, J.; Brooks, C.L. Dynamic reorganization of the functionally active ribosome explored by normal mode analysis and cryo-electron microscopy. Proc. Natl. Acad. Sci. USA 2003, 100, 9319–9323. [Google Scholar] [CrossRef] [Green Version]
  22. Wang, Y.; Rader, A.J.; Bahar, I.; Jernigan, R.L. Global ribosome motions revealed with elastic network model. J. Struct. Biol. 2004, 147, 302–314. [Google Scholar] [CrossRef]
  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. 2013, 9, e1003003. [Google Scholar] [CrossRef]
  24. Bock, L.V.; Blau, C.; Schröder, G.F.; Davydov, I.I.; Fischer, N.; Stark, H.; Rodnina, M.V.; Vaiana, A.C.; Grubmüller, H. Energy barriers and driving forces in tRNA translocation through the ribosome. Nat. Struct. Mol. Biol. 2013, 20, 1390–1396. [Google Scholar] [CrossRef] [PubMed]
  25. Bock, L.V.; Blau, C.; Vaiana, A.C.; Grubmüller, H. Dynamic contact network between ribosomal subunits enables rapid large-scale rotation during spontaneous translocation. Nucleic Acid Res. 2015, 43, 6747–6760. [Google Scholar] [CrossRef]
  26. Levi, M.; Whitford, P.C. Dissecting the energetics of subunit rotation in the ribosome. J. Phys. Chem. B 2019, 123, 2812–2923. [Google Scholar] [CrossRef] [PubMed]
  27. Whitford, P.C.; Noel, J.K.; Gosavi, S.; Schug, A.; Sanbonmatsu, K.Y.; Onuchic, J.N. An all-atom structure-based potential for proteins: Bridging minimal models with all-atom empirical forcefields. Prot. Struct. Func. Bioinfo. 2009, 75, 430–441. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Noel, J.K.; Levi, M.; Raghunathan, M.; Lammert, H.; Hayes, R.L.; Onuchic, J.N.; Whitford, P.C. SMOG 2: A Versatile Software Package for Generating Structure-Based Models. PLoS Comput. Biol. 2016, 12, e1004794. [Google Scholar] [CrossRef] [Green Version]
  29. Whitford, P.C.; Jiang, W.; Serwer, P. Simulations of Phage T7 Capsid Expansion Reveal the Role of Molecular Sterics on Dynamics. Viruses 2020, 12, 1273. [Google Scholar] [CrossRef]
  30. Duan, Y.; Wu, C.; Chowdhury, S.; Lee, M.C.; Xiong, G.; Zhang, W.; Yang, R.; Cieplak, P.; Luo, R.; Lee, T.; et al. A point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum mechanical calculations. J. Comp. Chem. 2003, 24, 1999–2012. [Google Scholar] [CrossRef]
  31. Noel, J.K.; Whitford, P.C.; Onuchic, J.N. The shadow map: A general contact definition for capturing the dynamics of biomolecular folding and function. J. Phys. Chem. B 2012, 116, 8692–8702. [Google Scholar] [CrossRef] [Green Version]
  32. Nguyen, K.; Whitford, P.C. Steric interactions lead to collective tilting motion in the ribosome during mRNA-tRNA translocation. Nat. Commun. 2016, 7, 10586. [Google Scholar] [CrossRef] [Green Version]
  33. Lindahl, E.; Hess, B.; van der Spoel, D. GROMACS 3.0: A package for molecular simulation and trajectory analysis. J. Mol. Mod. 2001, 7, 306–317. [Google Scholar] [CrossRef]
  34. Abraham, M.J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J.C.; Hess, B.; Lindahl, E. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 2015, 1–2, 19–25. [Google Scholar] [CrossRef] [Green Version]
  35. Yang, H.; Bandarkar, P.; Horne, R.; Leite, V.B.; Chahine, J.; Whitford, P.C. Diffusion of tRNA inside the ribosome is position-dependent. J. Chem. Phys. 2019, 151. [Google Scholar] [CrossRef]
  36. Hoffer, E.D.; Hong, S.; Sunita, S.; Maehigashi, T.; Gonzalez, R.L.J.; Whitford, P.C.; Dunham, C.M. Structural insights into mRNA reading frame regulation by tRNA modification and slippery codon–anticodon pairing. eLife 2020, 9, e51898. [Google Scholar] [CrossRef]
  37. Russell, R.B.; Barton, G.J. Multiple protein sequence alignment from tertiary structure comparison. Proteins 1992, 14, 309–323. [Google Scholar] [CrossRef]
  38. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  39. Okazaki, K.; Koga, N.; Takada, S.; Onuchic, J.N.; Wolynes, P.G. Multiple-basin energy landscapes for large-amplitude conformational motions of proteins: Structure-based molecular dynamics simulations. Proc. Nat. Acad. Sci. USA 2006, 103, 11844–11849. [Google Scholar] [CrossRef] [Green Version]
  40. Hyeon, C.; Thirumalai, D. Capturing the essence of folding and functions of biomolecules using coarse-grained models. Nat. Commun. 2011, 2, 487. [Google Scholar] [CrossRef]
  41. Pierro, M.D.; Potoyan, D.A.; Wolynes, P.G.; Onuchic, J.N. Anomalous diffusion, spatial coherence, and viscoelasticity from the energy landscape of human chromosomes. Proc. Natl. Acad. Sci. USA 2018, 115, 7753–7758. [Google Scholar] [CrossRef] [Green Version]
  42. 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. 2011, 62. [Google Scholar] [CrossRef]
  43. Noel, J.K.; Whitford, P.C. How EF-Tu can contribute to efficient proofreading of aa-tRNA by the ribosome. Nat. Commun. 2016, 7, 13314. [Google Scholar] [CrossRef] [Green Version]
  44. Nguyen, K.; Yang, H.; Whitford, P.C. How the Ribosomal A-Site Finger Can Lead to tRNA Species-Dependent Dynamics. J. Phys. Chem. B 2017, 121, 2767–2775. [Google Scholar] [CrossRef]
  45. Levi, M.; Walak, K.; Wang, A.; Mohanty, U.; Whitford, P.C. A steric gate controls P/E hybrid-state formation of tRNA on the ribosome. Nat. Commun. 2020, 11, 5706. [Google Scholar] [CrossRef]
  46. Levi, M.; Noel, J.K.; Whitford, P.C. Studying ribosome dynamics with simplified models. Methods 2019, 162-163, 128–140. [Google Scholar] [CrossRef]
  47. Best, R.B.; Hummer, G. Reaction coordinates and rates from transition paths. Proc. Natl. Acad. Sci. USA 2005, 102, 6732–6737. [Google Scholar] [CrossRef] [Green Version]
  48. Das, P.; Moll, M.; Stamati, H.; Kavraki, L.; Clementi, C. Low-dimensional, free-energy landscapes of protein-folding reactions by nonlinear dimensionality reduction. Proc. Natl. Acad. Sci. USA 2006, 103, 9885–9890. [Google Scholar] [CrossRef] [Green Version]
  49. Best, R.; Paci, E.; Hummer, G.; Dudko, O. Pulling direction as a reaction coordinate for the mechanical unfolding of single molecules. J. Phys. Chem. B 2008, 112, 5968–5976. [Google Scholar] [CrossRef]
  50. Krivov, S.V. On reaction coordinate optimality. J. Chem. Theory Comput. 2013, 9, 135–146. [Google Scholar] [CrossRef]
  51. Makarov, D.E. Barrier Crossing Dynamics from Single-Molecule Measurements. J. Phys. Chem. B 2021, 125, 2467–2476. [Google Scholar] [CrossRef] [PubMed]
  52. Hummer, G. Position-dependent diffusion coefficients and free energies from Bayesian analysis of equilibrium and replica molecular dynamics simulations. New J. Phys. 2005, 7, 34. [Google Scholar] [CrossRef]
  53. Lacroute, F. RNA and protein elongation rates in Saccharomyces cerevisiae. MGG Mol. Gen. Genet. 1973, 125, 319–327. [Google Scholar] [CrossRef] [PubMed]
  54. Noel, J.K.; Chahine, J.; Leite, V.B.; Whitford, P.C. Capturing transition paths and transition states for conformational rearrangements in the ribosome. Biophys. J. 2014, 107, 2881–2890. [Google Scholar] [CrossRef] [Green Version]
  55. Atilgan, A.R.; Durell, S.R.; Jernigan, R.L.; Demirel, M.C.; Keskin, O.; Bahar, I. Anisotropy of fluctuation dynamics of proteins with an elastic network model. Biophys. J. 2001, 80, 505–515. [Google Scholar] [CrossRef] [Green Version]
  56. Miyashita, O.; Onuchic, J.; Wolynes, P. Nonlinear elasticity, proteinquakes, and the energy landscapes of functional transitions in proteins. Proc. Natl. Acad. Sci. USA 2003, 100, 12570–12575. [Google Scholar] [CrossRef] [Green Version]
  57. Shan, Y.; Arkhipov, A.; Kim, E.T.; Pan, A.C.; Shawa, D.E. Transitions to catalytically inactive conformations in EGFR kinase. Proc. Natl. Acad. Sci. USA 2013, 110, 7270–7275. [Google Scholar] [CrossRef] [Green Version]
  58. Jackson, J.; Nguyen, K.; Whitford, P.C. Exploring the balance between folding and functional dynamics in proteins and RNA. Int. J. Mol. Sci. 2015, 16, 6868–6889. [Google Scholar] [CrossRef] [Green Version]
  59. Müller, C.W.; Schlauderer, G.J.; Reinstein, J.; Schulz, G.E. Adenylate kinase motions during catalysis: An energetic counterweight balancing substrate binding. Structure 1996, 4, 147–156. [Google Scholar] [CrossRef] [Green Version]
  60. Munro, J.B.; Altman, R.B.; O’Connor, N.; Blanchard, S.C. Identification of two distinct hybrid state intermediates on the ribosome. Mol. Cell 2007, 25, 505–517. [Google Scholar] [CrossRef] [Green Version]
  61. Chen, J.; Petrov, A.; Tsai, A.; O’Leary, S.E.; Puglisi, J.D. Coordinated conformational and compositional dynamics drive ribosome translocation. Nat. Struct. Mol. Biol. 2013, 20, 718–727. [Google Scholar] [CrossRef] [Green Version]
  62. 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]
  63. Levi, M.; Nguyen, K.; Dukaye, L.; Whitford, P.C. Quantifying the Relationship between Single-Molecule Probes and Subunit Rotation in the Ribosome. Biophys. J. 2017, 113, 2777–2786. [Google Scholar] [CrossRef] [Green Version]
  64. Whitford, P.C.; Ahmed, A.; Yu, Y.; Hennelly, S.P.; Tama, F.; Spahn, C.M.T.; Onuchic, J.N.; Sanbonmatsu, K.Y. Excited states of ribosome translocation revealed through integrative molecular modeling. Proc. Natl. Acad. Sci. USA 2011, 108, 18943–18948. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Subunit rotation in a eukaryotic ribosome. The elongation cycle of the ribosome involves numerous large-scale conformational rearrangements. (left) In the ribosome, the nascent protein chain is attached to the P-site tRNA molecule (red). A “classical” P/P configuration is shown. In this state, the small subunit (SSU; rRNA:cyan, protein:blue) is described as being in an “unrotated” conformation. Perspective shown in the bottom panel is rotated ∼90 about the horizontal axis. (right) After peptide bond formation, where the nascent chain is transferred to an incoming tRNA molecule, the P-site tRNA adopts a hybrid P/E conformation, where it is displaced toward the E site of the large subunit. P/E formation is also accompanied by a ∼8 counterclockwise rotation of the SSU, relative to the large subunit (LSU; rRNA:gray, proteins:white). White dashed lines are shown to highlight the relative rotation of the SSU. In the current study, we apply molecular dynamics simulations to probe the dynamics and energetics of SSU rotation in yeast. Structures shown are PDB entries 3J78 (unrotated, classical) and 3J77 (P/E, rotated) [18].
Figure 1. Subunit rotation in a eukaryotic ribosome. The elongation cycle of the ribosome involves numerous large-scale conformational rearrangements. (left) In the ribosome, the nascent protein chain is attached to the P-site tRNA molecule (red). A “classical” P/P configuration is shown. In this state, the small subunit (SSU; rRNA:cyan, protein:blue) is described as being in an “unrotated” conformation. Perspective shown in the bottom panel is rotated ∼90 about the horizontal axis. (right) After peptide bond formation, where the nascent chain is transferred to an incoming tRNA molecule, the P-site tRNA adopts a hybrid P/E conformation, where it is displaced toward the E site of the large subunit. P/E formation is also accompanied by a ∼8 counterclockwise rotation of the SSU, relative to the large subunit (LSU; rRNA:gray, proteins:white). White dashed lines are shown to highlight the relative rotation of the SSU. In the current study, we apply molecular dynamics simulations to probe the dynamics and energetics of SSU rotation in yeast. Structures shown are PDB entries 3J78 (unrotated, classical) and 3J77 (P/E, rotated) [18].
Biophysica 01 00016 g001
Figure 2. Simulations of subunit rotation in a eukaryotic ribosome. Using an all-atom structure-based model [28], we simulated spontaneous rotation and backrotation of the SSU in yeast. (A) Euler angle decomposition was used to describe the orientation of the SSU, where the angle ϕ body measures rotation of the SSU body, relative to the LSU. ϕ body is defined as rotation about the yellow vector. (B) The tilt angle θ body describes rotation that is orthogonal to the body rotation angle ϕ body . The tilt axis (red) may be in any direction perpendicular to the rotation axis (yellow), where the direction of the tilting axis is given by ψ body . In the figure, ψ body 60 is shown. (C) ϕ body and θ body shown for a single simulation (of 3, in total). In this model, there are distinct and sharp transitions between the unrotated ( ϕ body 1 ) and rotated ( ϕ body 7 ) orientations. There are also concommittent changes in θ body . For reference, the cryo-EM structures [18] correspond to ϕ body 1 . 8 and θ body 2 . 4 for the unrotated state and 8 . 6 and θ body 5 . 0 for the rotated state. Negative and non-zero angles for the unrotated conformation reflect the orientation of the SSU in yeast, relative to the reference bacterial system (E. coli). (D) In the simulation, the direction of tilting ( ψ body ) shifts to slightly higher values as the SSU rotates and tilts. This reveals how the direction of structural fluctuations depends on the global conformation of the ribosome. The effective simulated times are estimated based on previous comparisons with explicit-solvent simulations [35].
Figure 2. Simulations of subunit rotation in a eukaryotic ribosome. Using an all-atom structure-based model [28], we simulated spontaneous rotation and backrotation of the SSU in yeast. (A) Euler angle decomposition was used to describe the orientation of the SSU, where the angle ϕ body measures rotation of the SSU body, relative to the LSU. ϕ body is defined as rotation about the yellow vector. (B) The tilt angle θ body describes rotation that is orthogonal to the body rotation angle ϕ body . The tilt axis (red) may be in any direction perpendicular to the rotation axis (yellow), where the direction of the tilting axis is given by ψ body . In the figure, ψ body 60 is shown. (C) ϕ body and θ body shown for a single simulation (of 3, in total). In this model, there are distinct and sharp transitions between the unrotated ( ϕ body 1 ) and rotated ( ϕ body 7 ) orientations. There are also concommittent changes in θ body . For reference, the cryo-EM structures [18] correspond to ϕ body 1 . 8 and θ body 2 . 4 for the unrotated state and 8 . 6 and θ body 5 . 0 for the rotated state. Negative and non-zero angles for the unrotated conformation reflect the orientation of the SSU in yeast, relative to the reference bacterial system (E. coli). (D) In the simulation, the direction of tilting ( ψ body ) shifts to slightly higher values as the SSU rotates and tilts. This reveals how the direction of structural fluctuations depends on the global conformation of the ribosome. The effective simulated times are estimated based on previous comparisons with explicit-solvent simulations [35].
Biophysica 01 00016 g002
Figure 3. The free-energy landscape of rotation. Equilibrium simulations of the intact yeast ribosome (Figure 2) were used to characterize the energetics of subunit rotation. (A) The free energy as a function of ϕ body and θ body suggests the barrier for rotation is 7 k B T . Rotation and tilting are correlated, though the simulations reveal that tilting motions are widely accessible within the unrotated ensemble. When ϕ body 1 , the tilting angle θ body can adopt values up to 4 . Similarly, there is a range of 4 in the tilt angle for the rotated ensemble. This indicates that, while rotation and tilting are correlated, the motions are only weakly coupled. (B) Free energy as a function of the tilt angle θ body and direction of tilting ψ body . Comparison of cryo-EM structures would suggest a single direction of tilting, though the simulations indicate that as the ribosome adopts more tilted orientations, the direction of tilting shifts to slightly higher values of ψ body . As a result, it is appropriate to describe tilting in terms of a twist-like rearrangement that is concomitant with rotation. (C) While rotation/tilting involves a complex combination of motions, the kinetics of rotation is described well by the single coordinate ϕ body . Specifically, the probability of being on a transition path as a function of ϕ body , P ( TP | ϕ body ), adopts a peak value of 0.5 . Gray area indicates the identified TSE region. (D) The free energy as a function of ϕ body also yields a barrier that is comparable to that of the two-dimensional landscape shown in panel A.
Figure 3. The free-energy landscape of rotation. Equilibrium simulations of the intact yeast ribosome (Figure 2) were used to characterize the energetics of subunit rotation. (A) The free energy as a function of ϕ body and θ body suggests the barrier for rotation is 7 k B T . Rotation and tilting are correlated, though the simulations reveal that tilting motions are widely accessible within the unrotated ensemble. When ϕ body 1 , the tilting angle θ body can adopt values up to 4 . Similarly, there is a range of 4 in the tilt angle for the rotated ensemble. This indicates that, while rotation and tilting are correlated, the motions are only weakly coupled. (B) Free energy as a function of the tilt angle θ body and direction of tilting ψ body . Comparison of cryo-EM structures would suggest a single direction of tilting, though the simulations indicate that as the ribosome adopts more tilted orientations, the direction of tilting shifts to slightly higher values of ψ body . As a result, it is appropriate to describe tilting in terms of a twist-like rearrangement that is concomitant with rotation. (C) While rotation/tilting involves a complex combination of motions, the kinetics of rotation is described well by the single coordinate ϕ body . Specifically, the probability of being on a transition path as a function of ϕ body , P ( TP | ϕ body ), adopts a peak value of 0.5 . Gray area indicates the identified TSE region. (D) The free energy as a function of ϕ body also yields a barrier that is comparable to that of the two-dimensional landscape shown in panel A.
Biophysica 01 00016 g003
Figure 4. Contact analysis reveals sequential displacements at the SSU-LSU interface. (A) Average fraction of rotated-specific inter-subunit contacts formed, by atom, as a function of ϕ body . Many rotated-specific contacts are formed by the time the ribosome reaches the TSE ( ϕ body 3 . 5 ). Contact probabilities for each simulation may be found in Figure S2. (B) In contrast to the rotated-specific contacts (panel A), the unrotated-specific contacts are rarely formed in the TSE. (C) Structure of the SSU, viewed from the SSU-LSU interface. Red (yellow) spheres indicate which atoms have at least one rotated-specific (unrotated-specific) contact formed in more than 80 % of the TSE frames. (D) Structure of the LSU, viewed from the SSU-LSU interface, shown in the same representation as in panel C. While there are only four atoms involved in unrotated contacts in the TSE (yellow), there are clear clusters of rotated contacts formed in the TSE (red). In particular, contacts between h14 and protein L23 (also called uL14, found in bridge B8; circled) are frequently formed, indicating they can have a strong influence on the free-energy barrier and kinetics of rotation.
Figure 4. Contact analysis reveals sequential displacements at the SSU-LSU interface. (A) Average fraction of rotated-specific inter-subunit contacts formed, by atom, as a function of ϕ body . Many rotated-specific contacts are formed by the time the ribosome reaches the TSE ( ϕ body 3 . 5 ). Contact probabilities for each simulation may be found in Figure S2. (B) In contrast to the rotated-specific contacts (panel A), the unrotated-specific contacts are rarely formed in the TSE. (C) Structure of the SSU, viewed from the SSU-LSU interface. Red (yellow) spheres indicate which atoms have at least one rotated-specific (unrotated-specific) contact formed in more than 80 % of the TSE frames. (D) Structure of the LSU, viewed from the SSU-LSU interface, shown in the same representation as in panel C. While there are only four atoms involved in unrotated contacts in the TSE (yellow), there are clear clusters of rotated contacts formed in the TSE (red). In particular, contacts between h14 and protein L23 (also called uL14, found in bridge B8; circled) are frequently formed, indicating they can have a strong influence on the free-energy barrier and kinetics of rotation.
Biophysica 01 00016 g004
Figure 5. Temperature-dependent free-energy landscape. The presented simulations reveal how the free-energy landscape is influenced by changes in molecular flexibility upon rotation. (A) Difference in free energy of the rotated ( 5 . 40 < ϕ body < 7 . 83 ) and unrotated ( 1 . 12 < ϕ body < 0 . 82 ) ensembles, as a function of temperature. Positive values indicate the free energy of the rotated ensemble is higher. As T is increased, the free energy of the rotated ensemble monotonically decreases, relative to the unrotated ensemble. (B) Spatial root mean squared fluctuations (rmsf), by atom, for all simulated frames within the rotated ensemble. Peaks in the protein regions are typically due to marginally ordered tails at the peripheral regions of the ribosome. (C) Difference between rmsf values calculated for the rotated and unrotated ensembles: Δ rmsf . Positive values indicate elevated mobility in the rotated ensemble. Consistent with the temperature-induced shift in population towards the rotated ensemble (panel A), there is one region of the LSU rRNA that is significantly more flexible in the rotated ensemble (expansion segment 27b, ES27b). (D) Structural representation of the LSU rRNA, colored by Δ rmsf values (blue to red). There is a large increase in flexibility of ES27b, suggesting that rotation allows this region to adopt a wider range of configurations and provide an entropic drive towards the rotated ensemble. This observation is consistent with the lack of electron density obtained for this region in a cryo-EM reconstruction of the rotated state [18].
Figure 5. Temperature-dependent free-energy landscape. The presented simulations reveal how the free-energy landscape is influenced by changes in molecular flexibility upon rotation. (A) Difference in free energy of the rotated ( 5 . 40 < ϕ body < 7 . 83 ) and unrotated ( 1 . 12 < ϕ body < 0 . 82 ) ensembles, as a function of temperature. Positive values indicate the free energy of the rotated ensemble is higher. As T is increased, the free energy of the rotated ensemble monotonically decreases, relative to the unrotated ensemble. (B) Spatial root mean squared fluctuations (rmsf), by atom, for all simulated frames within the rotated ensemble. Peaks in the protein regions are typically due to marginally ordered tails at the peripheral regions of the ribosome. (C) Difference between rmsf values calculated for the rotated and unrotated ensembles: Δ rmsf . Positive values indicate elevated mobility in the rotated ensemble. Consistent with the temperature-induced shift in population towards the rotated ensemble (panel A), there is one region of the LSU rRNA that is significantly more flexible in the rotated ensemble (expansion segment 27b, ES27b). (D) Structural representation of the LSU rRNA, colored by Δ rmsf values (blue to red). There is a large increase in flexibility of ES27b, suggesting that rotation allows this region to adopt a wider range of configurations and provide an entropic drive towards the rotated ensemble. This observation is consistent with the lack of electron density obtained for this region in a cryo-EM reconstruction of the rotated state [18].
Biophysica 01 00016 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Freitas, F.C.; Fuchs, G.; de Oliveira, R.J.; Whitford, P.C. The Dynamics of Subunit Rotation in a Eukaryotic Ribosome. Biophysica 2021, 1, 204-221. https://0-doi-org.brum.beds.ac.uk/10.3390/biophysica1020016

AMA Style

Freitas FC, Fuchs G, de Oliveira RJ, Whitford PC. The Dynamics of Subunit Rotation in a Eukaryotic Ribosome. Biophysica. 2021; 1(2):204-221. https://0-doi-org.brum.beds.ac.uk/10.3390/biophysica1020016

Chicago/Turabian Style

Freitas, Frederico Campos, Gabriele Fuchs, Ronaldo Junio de Oliveira, and Paul Charles Whitford. 2021. "The Dynamics of Subunit Rotation in a Eukaryotic Ribosome" Biophysica 1, no. 2: 204-221. https://0-doi-org.brum.beds.ac.uk/10.3390/biophysica1020016

Article Metrics

Back to TopTop