Introduction

Ribosomes are universal cellular complexes made up of many different ribosomal proteins and long ribosomal ribonucleic acid (RNA) chains that enable decoding, successive peptide bond formation and the protection of nascent protein chains. They are present as a two-subunit structure in all living organisms with adaptive characteristics that make them fit to each physiological environment. Ribosomal proteins stabilize the ribosome structure or assist its activity, yet none plays a direct role in peptide bond formation or polypeptide elongation. The ribosome can be considered, therefore, as a giant ribozyme complex in which RNA controls catalytic activity1,2. The structural bases of mechanisms involved in ribosome function are studied intensively and observations of various conformations by cryo-electron microscopy (cryo-EM) and high-resolution X-ray crystallography established that during its functional cycle the ribosome undergoes conformational changes, ranging from local group reorientation to large domain movements3,4,5,6,7. During ribosome crystallization trials, metal clusters used to provide anomalous phasing power were found to extend dramatically the resolution of the diffraction, through a reduction of subunit internal flexibility, “by ‘gluing’ symmetrical related subunits”3. It has, furthermore, been shown that ribosome complexes with mRNA and tRNA gave better resolution than apo ribosome structures in X-ray crystallographic structures, suggesting that substrate binding helps stabilize subunits in a preferred conformation to enable more homogeneous samples for crystallization8,9. Comparisons between ribosomal subunit structures from mesophile and thermophile bacteria, as well as archaeal halophiles indicated that inherently flexible elements are involved in inter-subunit contacts and substrate binding10,11, leading Horan et al.10 to conclude that ‘inter-subunit movement is required for ribosome translocation, accounting for the universal two-subunit architecture of ribosomes’. In addition to the extensive structure-function studies, a variety of solvent conditions has been characterized and associated with halophile ribosomal subunit structural and functional integrity, during the extraordinary scientific effort that led to their crystallization11. All these data show that a detailed description of the ribosome requires an understanding in terms of dynamics.

Dynamics refers to the forces that stabilize structure and define motions. On one hand, single particle cryo-EM and X-ray structural studies, while informing on the conformational change amplitudes, do not contain information on forces or time-scales. On the other hand, molecular dynamics (MD) simulations and normal mode calculations are particularly challenging because of the size and complexity of ribosome structures12,13,14,15,16,17. Dynamics is also very sensitive to environment, which leads to a well-structured system being inactive under certain solvent conditions18, and it is now fully accepted that further to structure being adapted to function, appropriate molecular motions are also necessary. Molecular dynamics on the picosecond to nanosecond time scale acts as the ‘lubricant’ of functional conformational changes on longer time scales19,20,21. Following neutron scattering experiments, it has been proposed that adaptation occurs via evolutionary selection of dynamics on the thermal timescale22.

Neutron beams provide a soft radiation of wavelength and energy similar in magnitude to molecular thermal fluctuation amplitudes and energies, respectively, which inflicts practically no radiation damage on a biological sample. In a neutron scattering experiment, mean square displacements (MSD) and the effective force constant anchoring atoms in a structure are calculated from measured energy and momentum changes of the neutron as a function of temperature. To use a mechanical analogy: in a collision between a moving billiard ball and one that is anchored by a spring, the force constant of the spring as well as the vibration amplitude of the anchored ball can be calculated from the changes in energy and momentum of the moving ball after it is ‘scattered’. Systems studied by neutron scattering include intrinsically disordered proteins, whole cells, membranes, nanoparticles and tissue22,23,24,25,26,27.

The current study presents neutron scattering results that enabled the determination of global thermal dynamics in halophilic ribosomal subunit samples28. Halophiles were chosen because of sample stability. They also respond to solvent effects29,30, which will be applied in future experiments to explore dynamics/structure/function relations in more detail. Previous neutron scattering experiments on live cells with different physiological temperatures, demonstrated adaption of MD forces to yield similar mean square displacements (MSD) under active conditions22; e.g. psychrophiles at low temperature are ‘softer’ in order to have similar MSD to hyperthermophiles at high temperature. We would expect, therefore, that the results on halophiles would be more generally applicable than just reflecting an extremophile property, provided the measurements were performed in high salt active conditions.

H30S (Haloarcula marismortui 30 S subunit) and H50S (Haloarcula marismortui 50S subunit) dynamics was measured on the complementary ~10 picosecond and nanosecond timescales, corresponding, respectively, to water diffusion and fast local vibrations20,31, on the one hand, and side-chain conformational sampling, on the other. Free and hydration water are found with similar motion parameters in both H30S and H50S. A marginal solvent effect explored on H30S samples suggested proportionally more bound water in 3 M NaCl compared to 3 M KCl, correlated with stronger hydration interactions of the Na+ ion, as observed previously in halophilic proteins30. On the nanosecond timescale, effective force constants (resilience) from the temperature dependence of the fluctuations were calculated from the temperature dependence of the fluctuations. H30S is revealed to have softer resilience and larger mean square fluctuations than H50S. The results are discussed in terms of the role of dynamics in RNA activity and of the differential functional flexibility of 30 S versus 50 S subunits, for example, to facilitate conformational adjustments required for messenger and transfer RNA binding. The measured dynamics in the different timescales can, furthermore, serve as quantitative experimental input for coarse-grained MD simulations, to provide further understanding of the ribosome as a vital molecular machine.

Results

Thermal vibrations on the ~10 picosecond timescale are similar for both subunits

Quasi-elastic neutron scattering (QENS) spectra measured on the timescale of fast thermal dynamics are given in Fig. 1.

Figure 1
figure 1

Quasi-elastic neutron scattering intensities, at 298 K, binned over the whole wave-vector range accessible on the IN5 spectrometer (0.2–2 Å−1) and normalized to unity.

Left panel: H30S-3 M NaCl (solid line) and H50S-3 M NaCl (dotted line), Right panel: – H30S-3 M NaCl (solid line) and H30S-3 M KCl (dashed line). Insets: same spectra at reference Q = 0.3 Å−1 (fluctuation amplitudes > 20 Å) and Q = 1.2 Å−1 (amplitudes ~5 Å). A broader intensity distribution indicates scattering from faster motions (see text).

The left panel shows the comparison H30s versus H50S, both in 3 M NaCl in H2O. Visual inspection reveals a slightly broader spectrum for H30S-3 M NaCl vs H50S-3M NaCl, indicating marginally faster average dynamics in the 30S subunit. All samples were equilibrated in H2O, avoiding the use of D2O (see Methods). In H2O containing samples on the 10 ps timescale – a few Å length scale, the scattering signal is dominated by the diffusion of free and hydration H2O in the samples and fast motions within the ribosomal subunits (see ‘Theoretical background’ in Methods). Water content in the samples was 65% for H30S-3 M NaCl, and 67% for H50S-3 M NaCl (see Methods), so that the apparently ‘faster’ 30S subunit cannot be accounted for by the total water content in the sample, since it is 2% lower than for the 50 S, and is likely due to faster motions of groups within the subunit.

In further analysis given in Supplementary Information (SI), the QENS was fitted to yield two diffusive populations with parameters given in Table S1. The model used for the analysis is described in detail in SI. From the values of the diffusion coefficients, the populations are likely to correspond to free and hydration water, respectively, as previously published by Bellissent-Funel and collaborators32, with motion parameters that lie within errors for H30S and H50S. Analysis of the elastic intensity yields the proportion of H atoms seen as ‘immobile’ (i.e. moving too slowly to be resolved in the experiment, see SI). These were 6.3 ± 0.5% for H30S and 9.3 ± 0.3% for H50S, indicating a larger proportion of atoms in H50S compared to H30S that are from domains in which atomic motions are too slow to be resolved on the ~10 ps time scale. The faster average motions in 30 S could, therefore, be a reflection of fewer ‘immobile’ atoms on this timescale compared to the 50 S. Considering also that the ‘faster’ 30 S contains 2% less total water, a complementary interpretation may be that the ratio free-to-bound water in the 30S is higher than in the 50S (i.e. the 50S binds more water).

A small solvent effect is observed in Fig. 1 (right panel), where a marginally broader spectrum (faster dynamics) is observed for H30S-3 M KCl sample compared to H30S-3 M NaCl. The role of high K+ ion intracellular concentrations in halophiles has been discussed previously in terms of ion binding and hydration33,34. Water content in the H30S-3 M KCl and H30S-3 M NaCl was 62% and 65%, respectively (see Methods). Here again, the effect is small and the sample with less water appears to have faster dynamics. We tentatively propose an interpretation in terms of the different water binding properties of K+ and Na+ ions, respectively. Because of the stronger hydration interactions of Na+ compared to K+, we would expect a relatively higher ratio of free to bound water in the KCl sample, which would account for the ‘faster’ dynamics. Stronger hydration interactions in NaCl compared to KCl have been observed in a previous neutron scattering dynamics study of malate dehydrogenase from Haloarcula marismortui35. And different hydration effects between tRNA solutions in NaCl and KCl solvents have been reported from a small angle neutron scattering study36.

H30S is significantly softer with larger MSD than H50S (both in 3 M salt) on the nanosecond timescale of conformational sampling

Temperature dependence elastic scans on H30S and H50S were recorded on the IN16 spectrometer at ILL (see Methods). MSD versus temperature plots for the three samples are shown in Fig. 2. The observed dynamics is dominated by chemical group motions within the ribosomal subunits because free H2O and hydration water diffusion, observed by QENS on IN5, lie outside the elastic window ~1 ns - ~5 Å time- length-scale of the IN16 spectrometer32,37. Very slow water molecules that are strongly bound to ionic groups within the ribosomal subunits are considered as part of the biological complex. All samples in Fig. 2 indicate a flatter MSD below ~273 K, which may be due to a ‘stiffening’ of the ribosome subunit complexes due to freezing of free water in the samples. We, therefore, performed straight-line fits to the MSD beyond the kink at ~273 K to provide absolute scale values for the effective force constants (resilience, in N/m) in this timescale (see Methods). The values calculated from the MSD temperature dependence are in Table 1. NaCl and KCl sample values are within errors, indicating there is no solvent effect in this timescale except for slightly lower MSD for the KCl condition, which may be due to lower macromolecular plasticity in the drier sample (62% hydration compared to 65% for the sample in NaCl, see Methods). The effective force constant () is a factor of two smaller (softer) and the MSD at 37 °C significantly higher for H30S compared to the H50S both in 3 M NaCl. This difference cannot be due to water content because the ‘softer’ 30 S sample is drier (65% water content) than the 50 S (67%).

Table 1 <k> values calculated from the MSD temperature dependence (equation 3) and MSD at reference temperature value of 37 °C, in the nanosecond time scale.
Figure 2
figure 2

Mean square displacements (MSD) on the nanosecond timescale of H30S-3 M NaCl (circles), H30S-3 M KCl (squares) and H50S-3 M NaCl (triangles) extracted from elastic intensities in the wave-vector range 0.34–0.85 A−1.

The solid lines represent the best fits for data above 270 K (see text), while the dashed line corresponds to the MSD at the reference value of 37 °C. A steeper slope with temperature indicates a softer effective force constant for the motions (see text).

Compared to previous studies on the same timescale for a membrane38 and an intrinsically disordered protein39, the MSD values for the ribosomal subunits are significantly higher and the resilience values significantly lower. This is certainly due to the higher H2O content in the ribosome samples allowing the observation of large amplitude, softer diffusive motions. Indeed, an IN16 study of HDL (high density lipoprotein) in solution yielded MSD and resilience values of the same order as calculated from the ribosome data (J. Peters, private communication).

Discussion

In their 2015 review, Herschlag et al.40 write that “To understand RNA, it is necessary to move beyond a descriptive categorization towards quantitative predictions of its molecular conformations and functional behavior.” The plasticity of RNA structure has been demonstrated in the first decade after the first crystallographic structure of tRNA, including by small angle neutron scattering experiments41. Ribozyme motifs constitute evolved RNA molecules that carry out identical chemical functionality and represent ideal systems to study underlying structure-function relationships, illustrating the diversity of RNA’s functional role in biology (review by Fürtig et al.42). In a time-resolved NMR spectroscopy study of the adenine-dependent hairpin ribozyme, it has been shown that high activation barriers have to be overcome to populate active states with short life times43 and concluded that conformational dynamics, not chemistry, constitute the catalytic rate-limiting step. Fischer et al.6 have followed tRNA movement through the ribosome during translocation by time-resolved single-particle cryo-EM, and images of intermediate states revealed dynamic interactions between tRNA and ribosomal residues in a flat free energy landscape at physiological temperature. The authors concluded that “The ribosome functions as a Brownian machine that couples spontaneous conformational changes driven by thermal energy to directed movement”. The results from the various studies triggered further studies on the molecular dynamics of ribosomes in terms of forces and thermal energies by MD simulations (see references in the Introduction) and the present neutron scattering study, which delivers quantitative absolute scale values that could be used as input and/or experimental checks in computational studies of the ribosome, as was done, for example, in a previous study by Stadler et al.44 comparing dynamics of hemoglobin from different species, in which force constants measured by neutron scattering guided a coarse-grained elastic network approach to provide an understanding of the functional dynamics differences at the amino acid level44.

QENS results from IN5 and MSD results IN16 are complementary, informing on different dynamics. On the ~10 ps timescale, QENS is dominated by free and hydration water diffusion in the samples. Similar motion parameters were found for H30S and H50S, in 3 M NaCl. A comparison of H30S in 3 M NaCl and 3 M KCl suggested a higher ratio of free to bound water in the KCl sample in accordance with the stronger hydration interactions of the Na+ ion. Analysis of the elastic component on the 10ps timescale indicated a higher proportion (~9% cf ~6%) of atoms that appear to be immobile. This relative ‘immobility’ of the 50S in the fast IN5 timescale is reflected in the longer timescale of IN16, sampling large amplitude side chain motions, in which the mean effective force constant for nanosecond motions in H50S is a factor of 2 ‘stiffer’ than for H30S, while mean square displacements are a factor of 1.5 larger for H30S. These are significant differences indicating a flatter energy landscape for the small subunit that may play a functional role in initial mRNA and tRNA binding. The observations are likely be related to the larger amino acid/nucleotide ratio and solvent interaction area in H30S. H atoms in RNA groups have been shown in MD simulations as well as by neutron scattering to have stiffer dynamics45. The observation that H30S is softer than H50S is, therefore, in accordance with its higher nucleotide/(amino acid residue) ratio, 0.60 for 50S compared to 0.46 for 30 S. The larger surface to volume ratio in H30S compared to H50S would also contribute to a softer structure through a larger area of solvent interactions. They reveal the dynamic nature of conformational flexibility observed for the small subunit by crystallization trials, X-ray crystallography and cryo-EM. In the hypothesis that dynamics is similarly adapted to function under physiological conditions of temperature and solvent11,22, we expect the results obtained on the halophile ribosomal subunits in high salt to be generally applicable; i.e. 30S subunits are softer than 50 S for mesophiles under mesophile conditions, for thermophiles under thermophile conditions etc. The work could, therefore, pave the way for the design of ribosome inhibitors that would act by stiffening small subunit dynamics.

The current study opens up important perspectives for the experimental characterization of ribosome functional molecular dynamics. Since the scattering cross-section of deuterium is much weaker than that of natural abundance hydrogen, neutron scattering using specific deuterium labeling is a powerful method to study motions in parts of a complex system. This was done, for example, to measure the motions in the active core of the membrane protein bacteriorhodopsin46. Deuterium labelling of ribosomes has a long history47 and can be applied to observe separately the dynamics of different ribosome components (RNA, protein, bound tRNA…), with the aim of better characterising the dynamics interactions in a working ribosome as well as to relate the experimental observations with coarse-grained MD simulations, currently being developed by others.

Materials and Methods

Sample preparation

Ribosomes and ribosomal subunits from Haloarcula marismortui were obtained as described in Harms et al.28. Neutron scattering experiments can profit from the different scattering cross-sections of natural abundance hydrogen and deuterium in order to reduce or enhance the contribution of water to the signal. It is now well accepted, however, that the essential role of H-bonds in protein stabilization and dynamics can lead to significantly different dynamics in H2O and D2O as shown in previous thermodynamics and neutron scattering studies29,35. We were anxious to stay as close to dynamics under physiological conditions as possible so we avoided D2O, and included a slight excess of H2O in the samples. About 100 mg of each sample in either 3 M NaCl or 3 M KCl solvent, in which halophilic ribosomes are structurally stable, was equilibrated in H2O atmosphere to the following hydration levels: 65% for H30S-3 M NaCl, 62% for H30S-3 M KCl and 67% for H50S-3 M NaCl. Samples were sealed in vacuum-tight flat aluminum holders. In order to verify that no loss of material had occurred and that the hydration state was maintained, samples were weighed before and after the neutron scattering experiments.

Neutron scattering experiments

Quasi elastic neutron scattering (QENS) spectra were collected on the high-resolution time–of–flight spectrometer IN548 at the ILL in Grenoble. IN5 is a disk chopper time-of-flight spectrometer, which spans a wide range of incident wavelengths and energy resolutions. Neutrons, scattered at the sample, are detected at 4 m distance by 259 spherically ordered 3He detectors of the height of 3 m, covering an angular range between 14.5° and 132.5°. In addition, a position sensitive detector (PSD) counts neutrons in the small angle region between 2° and 8°. The instrument was used at two incident wavelengths of 5.1 and 10 Å providing the energy resolution of 75 and 12 μeV (FWHM) and the momentum transfer range 0.21 < Q < 1.98 Å–1 and 0.1 < Q < 1.01 Å–1, respectively.

Elastic (zero energy transfer, within instrumental resolution) spectra as a function of temperature were acquired on the cold backscattering spectrometer IN1649 at ILL, with an incident wavelength of λ = 6.27 Å and energy resolution of 0.9 μeV (FWHM) corresponding to a time window of about 1 ns and a momentum transfer (Q) range from 0.19 to 1.89 Å−1.

Data correction on both IN5 and IN16 was performed using LAMP50 to correct for transmission, normalize the raw data to the neutron flux, subtract the background due to the empty cell, correct for detector efficiency and energy resolution by using vanadium, a purely elastic, incoherent scatterer with a constant signal as a function of (Q)37. Transmission values were all above 0.9 so that multiple scattering effects were not taken into consideration for the data treatment.

Theoretical background

Neutron scattering data are sensitive mainly to hydrogen motions, because the incoherent neutron scattering cross-section of 1H is an order of magnitude larger than that of other nuclei in biological samples. Hydrogen atoms are homogeneously distributed in biological macromolecules, however, and are representative of global averaged molecular dynamics51. H-bond diffusive dynamics is mainly observed on the picosecond timescale, while, on longer timescales, hydrogen atoms reflect the conformational sampling motions of the groups to which they are bound51.

The scattered elastic incoherent intensity can be described in a Gaussian approximation by the dynamic structure factor at zero energy exchange52,53

where ΔE is instrumental energy resolution, related to the time window through Heisenberg’s uncertainty principle, and <u2> is the average time-dependent atomic mean square displacement (MSD)54,55 in the limit defined by ΔE, the energy resolution of the spectrometer. (Note that with respect to the radius of gyration of the motion about its center of mass (Rg), <u2> = 2Rg2).

Similarly to the Guinier approximation in small angle scattering, the Q range of validity for the Gaussian approximation depends on the geometry of the motion and could go as far as <u2> Q2 ≈ 456. The <u2> can be obtained for each temperature by the slope of the semi-logarithmic plot of the incoherent scattering function through

Note that in studies where temperature scan elastic window data are collected down to 20 K the absolute scale calibration is performed through normalization to the scattering curve at 20 K (e.g. see Wood et al.57). At this temperature the MSD is essentially zero leading to a theoretically constant ln Sel vs Q2. In the current study, it was important to remain above freezing temperature to avoid sample damage. At the lowest measured temperature, the MSD value is ~10 Å, very far from zero. The ln Sel vs Q2 slope values are, therefore, calibrated to an absolute Å2 scale through the division by the elastic incoherent vanadium scattering that is constant as a function of (Q) (see ‘Neutron scattering experiments’, above).

An effective average force constant for sample dynamics, <k>, can be calculated from the slope of <u2> as a function of temperature, by applying a quasi-harmonic approximation58:

In equation (3), the numerical constant is for <k> expressed in Newton per meter when <u2> is given in Ångstrom squared and T is the temperature in Kelvin.

Additional Information

How to cite this article: Zaccai, G. et al. The fluctuating ribosome: thermal molecular dynamics characterized by neutron scattering. Sci. Rep. 6, 37138; doi: 10.1038/srep37138 (2016).

Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.