# Probing the dynamic stalk region of the ribosome using
solution NMR

Xiaoling Wang, John Kirkpatrick, Helene Launay, Alfonso de Simone, Daniel

Häussinger, Christopher M. Dobson, Michele Vendruscolo, Cabrita Lisa,

Christopher A Waudby, John Christodoulou

► To cite this version:

Xiaoling Wang, John Kirkpatrick, Helene Launay, Alfonso de Simone, Daniel Häussinger, et al. Probing the dynamic stalk region of the ribosome using solution NMR. Scientific Reports, 2019, 9 (1), 10.1038/s41598-019-49190-1. hal-02466067

HAL Id: hal-02466067

https://amu.hal.science/hal-02466067

Submitted on 4 Feb 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
OPEN

# Probing the dynamic stalk region of the ribosome using solution NMR

Xiaolin Wang¹, John P. Kirkpatrick ID ¹, Hélène M. M. Launay¹, Alfonso de Simone², Daniel Häussinger ID ³, Christopher M. Dobson ID ⁴, Michele Vendruscolo ID ⁴, Lisa D. Cabrita¹, Christopher A. Waudby ID ¹ & John Christodoulou¹

We describe an NMR approach based on the measurement of residual dipolar couplings (RDCs) to probe the structural and motional properties of the dynamic regions of the ribosome. Alignment of intact 70S ribosomes in filamentous bacteriophage enabled measurement of RDCs in the mobile C-terminal domain (CTD) of the stalk protein bL12. A structural refinement of this domain using the observed RDCs did not show large changes relative to the isolated protein in the absence of the ribosome, and we also found that alignment of the CTD was almost independent of the presence of the core ribosome particle, indicating that the inter-domain linker has significant flexibility. The nature of this linker was subsequently probed in more detail using a paramagnetic alignment strategy, which revealed partial propagation of alignment between neighbouring domains, providing direct experimental validation of a structural ensemble previously derived from SAXS and NMR relaxation measurements. Our results demonstrate the prospect of better characterising dynamical and functional regions of more challenging macromolecular machines and systems, for example ribosome-nascent chain complexes.

In recent years, X-ray crystallography and cryo-electron microscopy (cryo-EM) have elucidated the details of high-resolution structures of ribosomes, revealing intricate mechanistic information about their function during the translation process¹,². In parallel, NMR-based observations of nascent polypeptide chains emerging from the ribosome are providing unique structural and mechanistic insights into co-translational folding processes³⁻⁵. In order to develop further solution-state NMR spectroscopy as a technique for structural studies of dynamic regions within large complexes, we have explored the measurement of residual dipolar couplings (RDCs) within intact ribosomes, focusing in particular on the mobile bL12 protein from the GTPase-associated region (GAR) of the prokaryotic 70S ribosome. RDCs have been used to characterise other macromolecular machines and assemblies, including HIV-1 capsid protein, bacterial Enzyme I, and the 20S proteasome⁶⁻⁸. These developments are particularly relevant as macromolecular complexes tend to exhibit a wide variety of functional motions that are challenging to characterise by methods such as X-ray crystallography or cryoelectron microscopy.

The GAR is a highly conserved region of both prokaryotic and eukaryotic ribosomes, and is so named to reflect its role in both the recruitment and the stimulation of the GTPase activity of several auxiliary factors associated with the key steps of protein synthesis: initiation (initiation factor 2, IF2), elongation (elongation factors EF-Tu and EF-G) and termination (release factor 3, RF3)⁹. The prokaryotic GAR includes helices 42–44 and 95 of the 23S rRNA, and the ribosomal proteins bL10, bL11 and bL12 (Fig. 1a). bL12, the focus of the present work, is a 120 residue dimeric protein consisting of an N-terminal dimerisation domain (NTD), which binds to the extended bL10 helix of the core ribosome particle, and a C-terminal domain (CTD), which interacts with GTPase substrates to facilitate their recruitment to the ribosome. The bL12 CTD interacts with four major translational GTPases, IF2, EF-Tu, EF-G and RF3, through a highly conserved and positively charged binding site in helices 4 and 5, identified through NMR mapping and mutagenesis¹⁰⁻¹⁴. GTPase binding has been found by mutagenesis to occur through the G4–G5 helix in the G domain of IF2 (and, by homology, likely also via the G domain in other GTPases), which presents a complementary negatively charged surface and results in extremely rapid binding driven by favourable electrostatics¹⁴. The CTD is mobile in solution, being separated from the NTD by a flexible

Received: 9 April 2019

Accepted: 31 July 2019

Published online: 19 September 2019

¹Institute of Structural and Molecular Biology, UCL and Birkbeck College London, Gower Street, London, WC1E 6BT, UK. ²Department of Life Sciences, Imperial College London, South Kensington Campus, London, SW7 2AZ, UK. ³Department of Chemistry, University of Basel, St. Johannsring 19, 4056, Basel, Switzerland. ⁴Centre for Misfolding Disease, Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge, CB2 1EW, UK. Correspondence and requests for materials should be addressed to C.A.W. (email: c.waudby@ucl.ac.uk) or J.C. (email: j.christodoulou@ucl.ac.uk)
<figure><img src="image_2.png" /><figcaption>**Figure 1.** RDC measurement and analysis of ribosome-bound and free bL12 aligned in Pf1 phage. (a) Schematic illustration of the 50S subunit of the bacterial ribosome illustrating the position and composition of the stalk region. (b) Excerpt from ¹H, ¹⁵N HSQC (red) and TROSY (blue) spectra of the ribosome-bound bL12 CTD in isotropic and aligned conditions (light and dark colouring, respectively), acquired at 700 MHz and 298 K. (c) Correlation plot of measured amide RDCs between free and ribosome-bound bL12. The error bars represent the uncertainties in the measured RDCs, as estimated from the signal-to-noise ratios and linewidths of the spectral peaks. (d) Amide RDCs measured for ribosome-bound and isolated bL12. (e) Refinement of the bL12 template structure (PDB code: 1rqu) using observed RDCs for isolated and ribosome-bound bL12. (f) Sanson-Flamsteed projection of the principal axes of the alignment tensor calculated for free and ribosome-bound bL12.</figcaption></figure>

hinge region, and is absent in structures of ribosomes determined by X-ray crystallography and cryoelectron microscopy.

bL12 is unique among ribosomal proteins in being present in multiple copies: typically four or six monomers are present, although some cyanobacteria may contain eight¹⁵. In the particular case of the *E. coli* ribosome studied here, four copies are present¹⁶. However, while bL12 is essential for protein synthesis¹⁷, not all copies are needed. While variant ribosomes containing only one bL12 dimer show less than 50% of the protein synthesis activity of wild-type ribosomes¹⁸, engineered variants in which only one CTD is present per bL12 dimer have approximately 80% activity related to wild-type ribosomes¹⁹. These observations suggest therefore that only one CTD may be functioning within each dimer at any one time.

Aside from the CTD, the flexible hinge region of bL12 is also essential for ribosome activity²⁰. However, this flexibility, which has been proposed to facilitate the recruitment of translation factors from the cellular space to the 30S/50S interface⁹,²¹,²², has precluded the observation and structure determination of the hinge region and CTD by cryo-EM and X-ray crystallography. Solution-state NMR spectroscopy is therefore uniquely suited to the structural characterisation of this region. Indeed, ¹H, ¹⁵N-correlation spectra of 70S ribosomes were found to contain resonances arising almost exclusively from the bL12 CTD²³⁻²⁵. No chemical shift perturbations were observed between the CTD in isolation and when bound to the ribosome, implying that the domain structure
is not perturbed by ribosome binding. However, analyses of resonance intensities showed that only two of the four CTDs have sufficient mobility to be detectable by solution-state NMR, indicating that significant dynamical heterogeneity exists within the stalk region<sup>24,25</sup>. Furthermore, <sup>15</sup>N nuclear spin relaxation measurements were used to characterise the anisotropic rotational diffusion of the 50S ribosome-bound CTD, from which it was observed that the diffusion tensor of the CTD was rotated 59.7° relative to the inertia tensor of the domain, indicating that the motion of the domain is perturbed by the presence of the hinge region and potentially also interactions between the CTD and the ribosome surface<sup>24</sup>. Within isolated bL12 dimers, an NMR analysis of secondary <sup>1</sup>H<sup>10</sup> chemical shifts within the hinge region indicated a weak propensity for beta-sheet structure, indicating that the hinge may not be completely flexible<sup>26</sup>. Further measurements of NMR relaxation and small-angle x-ray scattering (SAXS) were combined with molecular modelling to create an optimized structural ensemble of the bL12 dimer, which also indicated that the hinge region was not as flexible as would be expected from a random coil, and that the orientations of the NTD and CTD were instead partially coupled, with an order parameter of S<sup>2</sup> = 0.17<sup>27</sup>.

In this work, we have continued our NMR investigations of bL12 by measuring residual dipolar couplings (RDCs) for this protein in the spectrum of intact 70S ribosomes, in order to characterize both the structure of the CTD and the dynamical properties of the hinge region. In particular, we have sought to characterise the extent to which the orientation of the bL12 CTD is coupled to the (ribosome-bound) NTD, and so test the accuracy of the optimized structural ensemble described above<sup>27</sup>. RDCs can be measured under weakly anisotropic solution conditions and report on both the local structure of the protein and the overall orientation of domains with respect to a laboratory frame of reference<sup>28,29</sup>. As such, they are useful sources of distance-independent structural information that are highly complementary to chemical shifts, scalar couplings and NOEs. Anisotropy can be induced by a variety of steric and electrostatic methods including liquid crystals, bicelles, bacteriophage and stretched polyacrylamide gels<sup>30</sup>.

## Results

**RDCs measurement of the intact 70S ribosome.** We tested the alignment of samples of ribosomes in Pf1 bacteriophage and C₁₂E₅ PEG/hexanol mixtures<sup>31</sup>. The latter induced inhomogeneous alignment (Methods), but strong alignment (12.3 Hz <sup>2</sup>H quadrupolar splitting) was obtained in Pf1 bacteriophage and this system was therefore used for further measurements (see Supplementary Fig. 1). Amide RDCs (D<sub>NH</sub>) were determined for ribosome-bound bL12 by measurement of <sup>15</sup>N frequency differences between amide resonances in <sup>1</sup>H, <sup>15</sup>N HSQC and TROSY spectra (Fig. 1b)<sup>32</sup>. We found that this approach provided the greatest precision in values of RDCs when compared with alternative methods such as the in-phase/anti-phase approach (see Supplementary Fig. 2). Measurements of additional couplings to C' and Cα atoms were also explored, but the sensitivity was not sufficient to obtain precise results. The integrity of ribosomal samples was monitored with <sup>1</sup>H stimulated-echo diffusion measurements, acquired in an interleaved fashion with HSQC and TROSY experiments (see Supplementary Fig. 3)<sup>33</sup>.

Despite the 2.4 MDa molecular-weight of the 70S ribosome, and the low sample concentration (10 µM), RDCs could be obtained for the entire CTD of ribosomal bL12 (bL12<sub>ribo</sub>), with an RMS uncertainty (⟨σ<sup>2</sup><sub>D</sub>⟩<sup>1/2</sup>) of 1.8 Hz relative to a range of ±30 Hz (Fig. 1c,d and Supplementary Tables S1, S2). Despite the multiple copies of bL12 on the ribosome, of which two are observable by solution-state NMR<sup>24,25</sup>, for all residues only single resonances were resolved, in both the presence and absence of bacteriophage, with no discernable distortions in lineshapes. This indicates either that the conformation and alignment of domains is effectively averaged over the *ca*. 100 ms timescale (1/D<sub>NH</sub>) of the residual dipolar coupling, or that for all residues the difference in residual dipolar couplings between copies is much smaller than the <sup>15</sup>N TROSY linewidth (11 ± 2 Hz, s.d.). As only a single set of resonances is observed for bL12 (i.e. isotropic chemical shifts are also averaged between multiple copies), the former possibility seems more likely, but as in the latter case any differences that exist must be smaller than the <sup>15</sup>N linewidth, they would also be smaller than the experimental uncertainty and hence unlikely to affect the conclusions of the present study. RDCs in the CTD of isolated bL12 (bL12<sub>free</sub>) were of a similar magnitude to those of the domain on the ribosome (Fig. 1c), but alignment of the isolated NTD (whose resonances are unobservable in the ribosome-bound form) appeared to be significantly weaker, with RDC values of only a few Hz (Fig. 1d). These small RDCs may be due to the weaker electric dipole moment of the NTD relative to the CTD, which contains a highly charged binding site for GTPase substrates<sup>14</sup>, leading to weaker electrostatic interactions between the NTD and the phage<sup>34</sup> and hence weaker alignment relative to the CTD.

**Structure refinement.** The measured RDCs for the CTD were fitted by singular-value-decomposition (SVD) to a template structure (PDB code: 1rq<sup>35</sup>) to determine the alignment tensor, A<sup>36</sup>. The corresponding back-calculated RDCs showed good agreement with observed RDCs for both bL12<sub>ribo</sub> and bL12<sub>free</sub> (with Q factors of 0.31 and 0.26 respectively for secondary structure elements, Table S1 and Figure S1), although for both datasets the RMS deviation (RMSD) between measured and back-calculated RDCs was larger than the RMS uncertainty in the measurements (Supplementary Table S1). Thus, while the structures of both bL12<sub>ribo</sub> and bL12<sub>free</sub> CTD are clearly similar to the template structure, some additional information appears to be present. This finding suggests the possibility of using the RDC data to derive improved structures for both bL12<sub>ribo</sub> and bL12<sub>free</sub>.

Where RDCs have been measured for multiple internuclear vectors, or in multiple alignment media, detailed structure or ensemble structure calculations may be performed by minimization of a target energy function<sup>37</sup> or using restrained molecular dynamics simulations<sup>38</sup>. Ensemble structures are particularly effective for representing flexible or disordered regions of a domain such as loops, and their analysis may be helpful in reducing bias in fitting alignment tensors, although a single structure is also generally regarded as sufficient for this purpose provided that flexible regions are excluded<sup>30</sup>. As in the present situation N-H RDCs could only be measured in a single alignment medium, to avoid overfitting we limited ourselves to refinement of secondary structure elements in a single structure (Methods). This procedure resulted in small changes for both bL12<sub>ribo</sub> and bL12<sub>free</sub> structures
(Fig. 1e), with Cα RMSDs (of all residues, relative to the template structure) of 0.35 Å and 0.52 Å respectively, while Q factors for secondary structure elements improved slightly, to 0.24 and 0.18 respectively (Supplementary Table S1). The RMSD between the measured and back-calculated RDCs remained larger than the experimental error, indicating that some structural noise relative to the true structure is still present even in the refined structure. However, further structural refinement is not feasible with N-H RDCs alone, due to the risk of overfitting and potentially distortion of the peptide bond geometry. The differences between the template and refined structures are not large, and neither is the difference between free and ribosome-bound forms, but this is perhaps to be expected given the previously observed absence of chemical shift perturbations in the CTD between free and ribosome-bound forms²⁴,²⁵. We also note that the original isolated CTD structure³⁵ was determined from a more comprehensive experimental dataset than the single set of N-H RDCs we have measured here – in analogy with data available from ribosome-bound measurements – and so is still likely to represent the more accurate structure. Nevertheless, the calculations demonstrated here are an important step towards the use of RDCs for structural refinement in more complex systems, such as disordered ribosome–nascent chain complexes, for which much less is known of the ribosome-associated structure *a priori*³,⁴.

**Alignment tensor analysis of bL12 CTD.** Following refinement, we used the resulting structures to re-fit alignment tensors to the observed RDCs and so compare the preferred orientations of free and ribosome-bound bL12. The orientations of the principal axes of the alignment tensors can be visualised in a Sanson-Flamsteed projection (Fig. 1f). Numerical uncertainties were calculated using a Monte Carlo approach³⁹ that accounts for both experimental errors in the measured RDCs and residual structural noise in the refined template structure (which is of particular importance here as a single refined structure rather than an ensemble is used for this analysis). We note that while individual differences between free and ribosome-bound RDCs may be comparable to the experimental uncertainty, the global analysis of such measurements across multiple residues can nevertheless identify significant differences, and indeed we observe a small difference in the orientations of the Aₓₓ and Aᵧᵧ axes (Fig. 1f), with a corresponding difference between the 5D angles⁴⁰ of the two tensors of 13.8 ± 7.6°. This result indicates that the orientation of the bL12 ribosome CTD is influenced by the presence of the core ribosome particle, demonstrating that the bL12 CTD is not completely independent.

The influence of the ribosome on the orientation of the bL12 CTD could be due to incomplete flexibility of the linker region of bL12, or could be the result of transient interactions between the bL12 CTD and the ribosome core particle. We have previously shown that measurements of transferred transverse relaxation provide a sensitive probe of even weak interactions between nascent chains and the surface of the core ribosome particle⁴. However, in the case of the bL12 CTD, given that the observed ¹⁵N transverse relaxation rates of the ribosome-associated CTD (17 ± 7 Hz, s.d.) are only marginally larger than those of the free CTD (13 ± 2 Hz, s.d.), any transiently bound states must be extremely sparsely populated, and so unable to affect the observed alignment of the domain. Therefore, we focused our attention on the nature of the linker region, and sought to develop measurements of the extent to which it might couple the alignment of the CTD and the core ribosome particle.

**Measurement of alignment propagation with a lanthanide-binding tag.** To investigate further the properties of the linker region, and in particular to quantify the extent to which alignment of the CTD and core ribosome particle (including the bound NTD) might be coupled, we investigated the propagation of alignment between domains in isolated bL12 dimers following the attachment of a paramagnetic lanthanide-binding tag (LBT), Tm-DOTA-M8-4R4S⁴¹, to bL12 domains (Fig. 2a). We initially prepared three cysteine variants of bL12, containing the mutations S15C and S24C (on the surface of the NTD), and S89C (on the surface of the CTD). ¹H,¹⁵N HSQC spectra of all three variants showed dispersed resonances with only small chemical shift perturbations relative to the WT protein (Δδ<sub>NH</sub> < 0.1 ppm), and similar translational diffusion coefficients were measured with ¹H STE experiments for the WT and all variants (7.1 ± 0.6 × 10<sup>-11</sup> m² s<sup>-1</sup> (s.d.)), indicating that the stability of the domains and dimer formation had not been disrupted (Supplementary Figs 5, 6). However, upon attachment of the LBT to S15C and S24C variants, new resonances were observed with ¹H chemical shifts between 8 and 8.5 ppm, and these were found to overlay with a spectrum of the oxidized WT protein, which has previously been shown to be monomeric⁴² (Supplementary Fig. 6). ¹H STE diffusion measurements also showed an increase in translational diffusion coefficients for S15C and S24C (to 10.0 ± 0.1 × 10<sup>-11</sup> m² s<sup>-1</sup> and 11.6 ± 0.3 × 10<sup>-11</sup> m² s<sup>-1</sup> respectively). No such changes were observed in S89C upon attachment of the LBT. We therefore concluded that attachment of the LBT in both S15C and S24C variants disrupted dimer formation, and so only the S89C variant, containing the tag in the CTD, was used for further measurements.

Two samples were prepared using mixed ¹⁵N/LBT labelling schemes to distinguish between direct and indirect domain alignment: untagged domains within the bL12 dimer will only exhibit alignment if there is propagation of the alignment from the tagged CTD via the linker(s) (Fig. 2a). The directly tagged CTD showed strong alignment (maximum observed RDCs of ca. 20 Hz, Fig. 2b and Supplementary Table S2) with RDCs that fitted well to the template structure (Q=0.15). Alignments of the NTD and of the untagged CTD (separated by two linkers from the LBT-CTD), were weaker but non-zero (Fig. 2b and Supplementary Table S2), indicating that some alignment is propagated through the linker region.

**Analysis of bL12 CTD motion.** The alignment strengths of the different domains can be quantified from the magnitude of alignment tensors determined by SVD fitting of the observed RDCs (Fig. 2c). We note that while the relative uncertainties in individual RDC measurements may be large, this global fitting approach, combined with Monte Carlo error analysis, clearly finds that alignment tensor magnitudes are non-zero (Fig. 2c). The ratio of alignment tensor magnitudes between the directly aligned CTD domain, (2.5 ± 0.2) × 10⁻³, and the indirectly aligned NTD, (3.5 ± 1.3) × 10⁻⁴, then defines an order parameter, S² = 0.14 ± 0.05, that describes the degree of coupling taking place through the linker region. This measurement can be compared with previous
<figure><img src="image_3.png" /><figcaption>**Figure 2.** RDC measurement and analysis of bL12 aligned with a paramagnetic lanthanide-binding-tag. (a) Mixed ¹⁵N/lanthanide-binding-tag labelling schemes for measurement of the paramagnetic alignment of domains. (b) Amide RDCs measured for directly and indirectly paramagnetically aligned bL12 domains. (c) Magnitude of fitted alignment tensors for directly and indirectly aligned bL12 domains. Error bars represents uncertainties derived from the Monte-Carlo "mcDc" routine available in PALES³⁹. (d) bL12 linker order parameters determined experimentally as a ratio of the alignment magnitude between the CTD (direct) and the NTD (red, with standard error) and calculated from random and SAXS/NMR-optimised ensembles (blue)²⁷.</figcaption></figure>

predictions based upon: (i) a random ensemble with a fully flexible linker, and (ii) an optimised ensemble refined using SAXS and ¹⁵N nuclear spin relaxation measurements²⁷ (Fig. 2d). Our measurements are consistent with the latter ensemble, an interesting feature of which is the anti-correlated extension and retraction of the two CTDs. The linker extension is beyond that expected for a random coil linker, and may facilitate the efficient recruitment of elongation factors from the surrounding cytosol to the ribosome, while the retracted form may promote the activity of recruited CTPases by stabilizing the GTPase on the ribosome²² and remaining bound during GTP hydrolysis⁴³. In the context of ribosome-bound bL12, these anticorrelated motions may also be associated with the differential observability by solution NMR of two of the four copies of bL12 present on the bL10 stalk²⁴,²⁵,²⁷.

Discussion

In summary, we have shown that RDCs can be successfully measured by solution-state NMR spectroscopy for the flexible bL12 region of the 2.4 MDa 70S ribosome. The analysis of these measurements indicated that, consistent with the observed absence of chemical shift changes between free and ribosome-associated forms²⁴,²⁵, the refined structures of bL12free and bL12ribo are very similar (Cα rmsd 0.37 Å). Importantly, this study has laid the foundations for detailed characterization using RDCs of modified bL12 variants on the ribosome, and also for more detailed structural investigations into the co-translational folding of ribosome–nascent chain complexes, for which the structural consequences of ribosomal tethering are at present poorly understood³,⁴. Our analysis of
RDCs also indicated that the alignment of the bL12 CTD was weakly perturbed by the presence of the ribosome, either through direct effects of a partially structured hinge region, or through interactions between the CTD and the ribosome surface. As ¹⁵N relaxation measurements indicated that any such interactions must be extremely weak, we therefore investigated the nature of the interdomain linker in more detail in the absence of the core ribosome particle, and found that specific paramagnetic alignment of domains could be used to measure the coupling between N- and C-terminal domains due to the interdomain linker, which was characterized by an S² order parameter of 0.14 ± 0.05. These results validate a previously determined structural ensemble²⁷, indicating that the interdomain linker is not fully flexible and instead allows the tethered C-terminal domains to sample a more expanded conformational space, resulting in more efficient recruitment of elongation factors.

## Methods

**NMR sample preparation.** Uniformly ¹⁵N-labeled ribosomes were prepared by growing *E. coli* cells in MDG⁴⁴ medium containing ¹⁵NH₄Cl as the nitrogen source for 20 hours at 310 K before harvesting. In order to escape stationary phase and to ensure that the ribosomes are not hibernating, the cell pellet was then resuspended in an equal volume of fresh medium and incubated at 37 °C until a steady increase of OD₆₀₀ allowed harvesting at growth phase. The resulting cell pellet was lysed by French press in purification buffer (30 mM HEPES, 12 mM Mg(OAc)₂, 1 M KOAc, 5 mM EDTA, 2 mM beta-mercaptoethanol (BME), pH 7.6) supplemented with trace of DNase. The soluble fraction of the lysate was purified with a 35% (w/v) sucrose cushion prepared in purification buffer supplemented with 5 mM ATP. The resulting ribosomal pellet was further purified through a linear density gradient of 10–35% (w/v) sucrose. Fractions containing 70S ribosomes were identified by SDS-PAGE and pure 70S ribosomes were concentrated and buffer exchanged into Tico buffer (10 mM HEPES, 12 mM MgCl₂, 30 mM NH₄Cl, 5 mM EDTA, 2 mM BME, pH 7.6). The concentration of ribosomes was determined from the absorbance at 260 nm using the extinction coefficient ε = 4.2 × 10⁷ M⁻¹cm⁻¹³³.

Uniformly ¹⁵N-labeled bL12 proteins were prepared by transformation into *E. coli* cells and growth in M9 minimal medium containing ¹⁵NH₄Cl as the nitrogen source at 37 °C. When the OD₆₀₀ reading reached 0.6, expression was induced with 1 mM IPTG for 4 h at 37 °C. The harvested cell pellet was lysed and purified with a sucrose cushion as described above. bL12 protein was purified from the supernatant using a HiTrap Q HP column (loaded in 50 mM Tris-HCl, pH 8.0; elution with a linear gradient of 0–1 M NaCl), followed by a HiTrap Phenyl HPcolumn (loaded in 25 mM Na₂HPO₄, 1.7 M (NH₄)₂SO₄, pH 7.5; elution with a reverse gradient of 1.7–1 M (NH₄)₂SO₄) and finally size exclusion chromatography with an Superdex-200 column (50 mM Na₂HPO₄, 100 mM KCl, pH 7.0). The concentration of bL12 protein was determined using the bicinchoninic acid (BCA) protein assay reagent (Pierce). Unlabeled bL12 proteins were prepared the same apart from using unlabeled NH₄Cl in the growth media.

NMR sample concentrations were 10 µM for the ribosomal particle and 50 µM for isolated bL12. To measure RDCs on the 70S ribosome particle, aligned 70S ribosome samples were prepared by addition of 15 mg ml⁻¹ Pfl filamentous phage (ASLA). For isolated bL12, alignment was achieved by addition of 6.5 mg ml⁻¹ filamentous phage. The homogeneity of the resulting anisotropic liquid crystalline medium was confirmed by monitoring the splitting and lineshape of ²H solvent signal⁴⁵ (see Supplementary Fig. 1). The other alignment medium tested (C₁₂E₅ PEG and hexanol³¹) was found not to be suitable due to a strong interaction with isolated bL12 protein. In addition, the presence of ribosomes appeared to interfere with the formation of the homogeneously aligned lamellar phase, as judged from the irregular lineshape of the ²H solvent signal (see Supplementary Fig. 1). This latter effect is most likely to arise from the diameter of the ribosomal particle being larger than the inter-lamellar spacing at the PEG concentration used.

To induce alignment using the Tm-DOTA-M8-4R4S⁴¹ lanthanide binding tag (LBT), cysteine mutations were introduced at positions S15, S24 and S89. 2D ¹H,¹⁵N correlation spectra and ¹H STE diffusion measurements were acquired for all variants to assess possible perturbations of structure and dimer formation. For labeling, bL12 (S15C, S24C and S89C) protein was incubated in buffer containing 10 mM DTT and then buffer exchanged into a non-DTT containing buffer on a PD-10 desalting column and reacted with the LBT immediately. The reaction was followed by ¹H, ¹⁵N SOFAST-HMQC NMR at 298 K until there was no further change in the spectrum.

**NMR data acquisition.** All NMR experiments were recorded at 298 K on a Bruker Avance III 700 MHz spectrometer equipped with a TCI cryoprobe. Data were acquired by recording a repeated series of relatively short (time) experiments in an interleaved manner. This approach allows the use of 1D and diffusion experiments to monitor the ribosome integrity at regular intervals. The final spectra for measurement of coupling constants were obtained by summing the individual experiments over the time period for which the ribosome was deemed to be intact, as assessed by ¹H STE diffusion experiments⁴⁶. These were recorded using an encoding/decoding gradient length, δ, of 4.8 ms, a diffusion delay, Δ, of 400 ms, and gradient strengths that were 5, 35, 65 and 95% of the maximum, 55.57 G cm⁻¹. To exclude interfering low molecular weight components, the integrity of ribosome samples was monitored using the intensity ratio between 65% and 95% gradient strengths (see Supplementary Fig. 3). An initial decrease in the diffusion coefficient for the isotropic sample can be attributed to a poor baseline due to the residual water signal (an issue which in the future could be ameliorated by the use of isotope-edited XSTE experiments⁴⁷), as the diffusion coefficients calculated from experiments later in time drops back to ca. 2 × 10⁻¹¹ m² s⁻¹, which corresponds to intact 70S ribosomes³³. The slower diffusion observed for the aligned sample is a result of elevated viscosity due to the presence of phage (see Supplementary Fig. 3).

Both in-phase/anti-phase (IPAP) ¹H-coupled ¹⁵N-HSQC spectra and the pairwise combination of ¹⁵N-HSQC/¹⁵N-TROSY spectra were recorded for the ribosome, and the combination of ¹⁵N-HSQC/¹⁵N-TROSY spectra was found to yield slightly smaller relative uncertainties in the calculated RDCs than with the approach where two sub-spectra containing the downfield and upfield lines of the ¹⁵N doublets are generated by taking the sum and difference between in-phase and anti-phase ¹H-coupled ¹⁵N-HSQC spectra (see Supplementary Fig. 2).
Backbone ¹⁵N R₂ relaxation rates were measured on ribosome-bound and isolated bL12 using established methods⁴⁸,⁴⁹.

**RDC data analysis.** Residual dipolar couplings, *D*, were measured as the difference between the isotropic and anisotropic splittings in the ¹⁵N dimension of HSQC and TROSY spectra for the ribosome, and IPAP ¹H-coupled ¹⁵N-HSQC spectra for isolated bL12 (Supplementary Table S2). Uncertainties, $\sigma_D$, were propagated using standard methods from the uncertainty in peak positions, determined from the signal-to-noise ratio (SN) and ¹⁵N linewidth (LW), $\sigma_v$ = LW/2SN. Initial data analysis was performed by fitting the measured RDCs to the available NMR-derived structure of the bL12 protein (PDB code 1rq³), as implemented by the SVD (singular-value-decomposition) routine in the program PALES³⁶, in order to derive the alignment tensor (A) and the corresponding set of back-calculated RDCs. Orientations of alignment tensors are relative to the PDB frame of the template structure. To be able to restrict the analysis of RDCs to the structure and orientational probability distribution of the domain, but not the local dynamics within the domain, residues with local flexibility (order parameter $S^2$ < 0.85) or undergoing conformational exchange are excluded from the structure and alignment tensor analysis. In addition, residues whose resonances are overlapped in the spectra and residues with experimental uncertainties $\sigma_D$ > 5 Hz are excluded.

**Structure refinement.** The template structure for the CTD of bL12 was refined using Xplor-NIH⁵⁰. The refinement was implemented in a simulating- annealing molecular-dynamics simulation by introducing an additional term to the force-field whose energy depends on the difference between the experimental RDC and that back-calculated from the structure. In the refinement, only the orientation of the secondary-structure elements (individual α-helices and β-strands) is optimized⁵¹ and the local structure of the secondary structure elements was fixed using non-crystallographic symmetry (NCS) restraints, and only those RDCs corresponding to rigid N-H vectors within the secondary structure elements were included as restraints. The values for the magnitude and rhombicity of the alignment tensor for back-calculation of the RDCs were taken as those of the best-fit alignment tensor derived using the template structure. The orientation of the alignment tensor is allowed to float during the refinement. The RDC restraints were implemented as soft-square-well susceptibility anisotropy (SANI) restraints, with the width of the according to the experimental uncertainty ($\sigma_D$) of each measured RDC. Euler angles for rotation of the alignment tensor from the PDB frame into the principal axis frame before and after structural refinement of bL12 free and bL12 ribo CTD are listed in Supplementary Table S3. The five components of the irreducible tensor representations of the alignment tensors and their general magnitudes are shown in Supplementary Table S4.

**References**

1. Fischer, N. et al. Structure of the E. coli ribosome-EF-Tu complex at <3 Å resolution by Cs-corrected cryo-EM. *Nature* **520**, 567–570 (2015).

2. Selmer, M. et al. Structure of the 70S ribosome complexed with mRNA and tRNA. *Science* **313**, 1935–1942 (2006).

3. Cabrita, L. D. et al. A structural ensemble of a ribosome-nascent chain complex during cotranslational protein folding. *Nat. Struct. Mol. Biol.* **23**, 278–285 (2016).

4. Deckert, A. et al. Structural characterization of the interaction of α-synuclein nascent chains with the ribosomal surface and trigger factor. *Proc. Natl. Acad. Sci. USA* **113**, 5012–5017 (2016).

5. Waudby, C. A. et al. Systematic mapping of free energy landscapes of a growing filamin domain during biosynthesis. *Proc. Natl. Acad. Sci. USA* **115**, 9744–9749 (2018).

6. Deshmukh, L. et al. Structure and dynamics of full-length HIV-1 capsid protein in solution. *J. Am. Chem. Soc.* **135**, 16133–16147 (2013).

7. Schwieters, C. D. et al. Solution structure of the 128 kDa enzyme I dimer from Escherichia coli and its 146 kDa complex with HPr using residual dipolar couplings and small- and wide-angle X-ray scattering. *J. Am. Chem. Soc.* **132**, 13026–13045 (2010).

8. Ruschak, A. M. & Kay, L. E. Proteasome allostery as a population shift between interchanging conformers. *Proc. Natl. Acad. Sci. USA* **109**, E3454–E62 (2012).

9. Wahl, M. C. & Möller, W. Structure and function of the acidic ribosomal stalk proteins. *Curr. Protein Pept. Sci.* **3**, 93–106 (2002).

10. Helgstrand, M. et al. The ribosomal stalk binds to translation factors IF2, EF-Tu, EF-G and RF3 via a conserved region of the L12 C-terminal domain. *J. Mol. Biol.* **365**, 468–479 (2007).

11. Kothe, U., Wieden, H.-J., Mohr, D. & Rodnina, M. V. Interaction of helix D of elongation factor Tu with helices 4 and 5 of protein L7/12 on the ribosome. *J. Mol. Biol.* **336**, 1011–1021 (2004).

12. Tourigny, D. S., Fernandez, I. S., Kelley, A. C. & Ramakrishnan, V. Elongation factor G bound to the ribosome in an intermediate state of translocation. *Science* **340**, 1235490–1235490 (2013).

13. Pallesen, J. et al. Cryo-EM visualization of the ribosome in termination complex with apo-RF3 and RF1. *Elife* **2**, e00411 (2013).

14. Ge, X., Mandava, C. S., Lind, C., Aqvist, J. & Sanyal, S. Complementary charge-based interaction between the ribosomal-stalk protein L7/12 and IF2 is the key to rapid subunit association. *Proc. Natl. Acad. Sci. USA* **115**, 4649–4654 (2018).

15. Davydov, I. I. et al. Evolution of the protein stoichiometry in the L12 stalk of bacterial and organellar ribosomes. *Nat Commun* **4**, 1387 (2013).

16. Subramanian, A. R. Copies of proteins L7 and L12 and heterogeneity of the large subunit of Escherichia coli ribosome. *J. Mol. Biol.* **95**, 1–8 (1975).

17. Traub, P. & Nomura, M. Structure and function of Escherichia coli ribosomes. I. Partial fractionation of the functionally active ribosomal proteins and reconstitution of artificial subribosomal particles. *J. Mol. Biol.* **34**, 575–593 (1968).

18. Mandava, C. S. et al. Bacterial ribosome requires multiple L12 dimers for efficient initiation and elongation of protein synthesis involving IF2 and EF-G. *Nucleic Acids Res.* **40**, 2054–2064 (2012).

19. Oleinikov, A. V., Jokhadze, G. G. & Traut, R. R. A single-headed dimer of Escherichia coli ribosomal protein L7/L12 supports protein synthesis. *Proc. Natl. Acad. Sci. USA* **95**, 4215–4218 (1998).

20. Oleinikov, A. V., Perroud, B., Wang, B. & Traut, R. R. Structural and functional domains of Escherichia coli ribosomal protein L7/L12. The hinge region is required for activity. *J. Biol. Chem.* **268**, 917–922 (1993).

21. Mohr, D., Wintermeyer, W. & Rodnina, M. V. GTPase activation of elongation factors Tu and G on the ribosome. *Biochemistry* **41**, 12520–12528 (2002).

22. Diaconu, M. et al. Structural basis for the function of the ribosomal L7/12 stalk in factor binding and GTPase activation. *Cell* **121**, 991–1004 (2005).
23. Gudkov, A. T., Gongadze, G. M., Bushuev, V. N. & Okon, M. S. Proton nuclear magnetic resonance study of the ribosomal protein L7/L12 in *situ*. *FEBS Lett.* **138**, 229–232 (1982).

24. Christodoulou, J. et al. Heteronuclear NMR investigations of dynamic regions of intact Escherichia coli ribosomes. *Proc. Natl. Acad. Sci. USA* **101**, 10949–10954 (2004).

25. Mulder, F. A. A. et al. Conformation and dynamics of ribosomal stalk protein L12 in solution and on the ribosome. *Biochemistry* **43**, 5930–5936 (2004).

26. Bocharov, E. V., Gudkov, A. T. & Arseniev, A. S. Topology of the secondary structure elements of ribosomal protein L7/L12 from E. coli in solution. *FEBS Lett.* **379**, 291–294 (1996).

27. Bernadó, P. et al. Structure and Dynamics of Ribosomal Protein L12: An Ensemble Model Based on SAXS and NMR Relaxation. *Biophys. J.* **98**, 2374–2382 (2010).

28. Tolman, J. R., Flanagan, J. M., Kennedy, M. A. & Prestegard, J. H. Nuclear magnetic dipole interactions in field-oriented proteins: information for structure determination in solution. *Proc. Natl. Acad. Sci. USA* **92**, 9279–9283 (1995).

29. Tjandra, N. & Bax, A. Direct measurement of distances and angles in biomolecules by NMR in a dilute liquid crystalline medium. *Science* **278**, 1111–1114 (1997).

30. Chen, K. & Tjandra, N. The use of residual dipolar coupling in studying proteins by NMR. *Top Curr Chem* **326**, 47–67 (2012).

31. Rückert, M. & Otting, G. Alignment of Biological Macromolecules in Novel Nonionic Liquid Crystalline Media for NMR Experiments - Journal of the American Chemical Society (ACS Publications). *J. Am. Chem. Soc.* (2000).

32. Kontaxis, G., Clore, G. M. & Bax, A. Evaluation of cross-correlation effects and measurement of one-bond couplings in proteins with short transverse relaxation times. *J. Magn. Reson.* **143**, 184–196 (2000).

33. Cassaignau, A. M. E. et al. A strategy for co-translational folding studies of ribosome-bound nascent chain complexes using NMR spectroscopy. *Nat Protoc* **11**, 1492–1507 (2016).

34. Hansen, M. R., Hanson, P. & Pardi, A. Filamentous bacteriophage for aligning RNA, DNA, and proteins for measurement of nuclear magnetic resonance dipolar coupling interactions. *Meth. Enzymol.* **317**, 220–240 (2000).

35. Bocharov, E. V. et al. From structure and dynamics of protein L7/L12 to molecular switching in ribosome. *J. Biol. Chem.* **279**, 17697–17706 (2004).

36. Zweckstetter, M. NMR: prediction of molecular alignment from structure using the PALES software. *Nat Protoc* **3**, 679–690 (2008).

37. Schwieters, C. D., Bermejo, G. A. & Clore, G. M. Xplor-NIH for molecular structure determination from NMR and other data sources. *Protein Sci.* **27**, 26–40 (2018).

38. Camilloni, C. & Vendruscolo, M. A Tensor-Free Method for the Structural and Dynamical Refinement of Proteins using Residual Dipolar Couplings. *J Phys Chem B*, https://doi.org/10.1021/jp5021824 (2014).

39. Zweckstetter, M. & Bax, A. Evaluation of uncertainty in alignment tensors obtained from dipolar couplings. *J. Biomol. NMR* **23**, 127–137 (2002).

40. Sass, J. et al. Purple Membrane Induced Alignment of Biological Macromolecules in the Magnetic Field. *J. Am. Chem. Soc.* **121**, 2047–2055 (1999).

41. Häussinger, D., Huang, J.-R. & Grzesiek, S. DOTA-M8: An extremely rigid, high-affinity lanthanide chelating tag for PCS NMR spectroscopy. *J. Am. Chem. Soc.* **131**, 14761–14767 (2009).

42. Caldwell, P., Luk, D. C., Weissbach, H. & Brot, N. Oxidation of the methionine residues of Escherichia coli ribosomal protein L12 decreases the protein's biological activity. *Proc. Natl. Acad. Sci. USA* **75**, 5349–5352 (1978).

43. Savelsbergh, A., Mohr, D., Kothe, U., Wintermeyer, W. & Rodnina, M. V. Control of phosphate release from elongation factor G by ribosomal protein L7/12. *EMBO J.* **24**, 4316–4323 (2005).

44. Studier, F. W. Protein production by auto-induction in high density shaking cultures. *Protein Expr. Purif.* **41**, 207–234 (2005).

45. Zweckstetter, M. & Bax, A. Characterization of molecular alignment in aqueous suspensions of Pf1 bacteriophage. *J. Biomol. NMR* **20**, 365–377 (2001).

46. Tanner, J. Use of the stimulated echo in NMR diffusion studies. *J Chem Phys* (1970).

47. Ferrage, F., Zoonens, M., Warschawski, D. E., Popot, J.-L. & Bodenhausen, G. Slow diffusion of macromolecular assemblies by a new pulsed field gradient NMR method. *J. Am. Chem. Soc.* **125**, 2541–2545 (2003).

48. Kay, L. E., Torchia, D. A. & Bax, A. Backbone dynamics of proteins as studied by 15N inverse detected heteronuclear NMR spectroscopy: application to staphylococcal nuclease. *Biochemistry* **28**, 8972–8979 (1989).

49. Dayie, K. T. & Wagner, G. Relaxation-rate measurements for 15N-1H groups with pulsed-field gradients and preservation of coherence pathways. *Journal of Magnetic Resonance* **111**, 121–126 (1994).

50. Schwieters, C. D., Kuszewski, J. J., Tjandra, N. & Clore, G. M. The Xplor-NIH NMR molecular structure determination package. *J. Magn. Reson.* **160**, 65–73 (2003).

51. Kirkpatrick, J. P., Li, P. & Carlomagno, T. Probing mutation-induced structural perturbations by refinement against residual dipolar couplings: application to the U4 spliceosomal RNP complex. *Chembiochem* **10**, 1007–1014 (2009).

## Acknowledgements

We acknowledge the use of the UCL NMR Centre and the staff for their support. This work was supported by a Wellcome Trust Investigator Award to J.C. (206409/Z/17/Z), a Swiss National Science Foundation research grant (SNF 200021 130263) to D.H., and the Francis Crick Institute through provision of access to the MRC Biomedical NMR Centre. The Francis Crick Institute receives its core funding from Cancer Research UK (FC001029), the UK Medical Research Council (FC001029), and the Wellcome Trust (FC001029).

## Author Contributions

X.W., J.P.K., H.M.M.L., C.M.D., M.V., L.D.C., C.A.W. and J.C. conceived the experiments. X.W., J.P.K. and H.M.M.L. conducted the experiments. D.H. provided reagents. A.d.S. and M.V. performed structure calculations. All authors analysed and discussed the results. X.W., C.A.W. and J.C. wrote the manuscript with input from all authors.

## Additional Information

Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-019-49190-1.

**Competing Interests:** The authors declare no competing interests.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
<figure><img src="image_4.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019