Saturday, October 5, 2024
No menu items!
HomeNatureThe ribosome lowers the entropic penalty of protein folding

The ribosome lowers the entropic penalty of protein folding

Protein expression and purification

DNA constructs of FLN5 were previously described7,11. Coding sequences for titin I27 and HRAS were introduced into the pLDC-17 vector using standard procedures. Further mutations were introduced using site-directed mutagenesis; for 19F labelling, amber stop codons were introduced6 in position 32 in HRAS, and residue 14 with an additional K87H point mutation in I27. FLN5 variants were expressed as His-tagged proteins and isotopically labelled in Escherichia coli BL21 DE3-Gold cells as previously described6,7; an identical protocol was used to produce purified samples of I27 and HRAS. RNC constructs comprised an arrest-enhanced variant of the SecM stalling sequence, FSTPVWIWWWPRIRGPP, as previously described6. Purification of isolated FLN5 A3A3 was performed by affinity chromatography followed by size-exclusion chromatography in the presence of 6 M urea prior to buffer exchange into Tico buffer (10 mM Hepes, 30 mM NH4Cl, 12 mM MgCl2, 1 mM EDTA). The full protein sequence of the FLN5 A3A3 is deposited together with its chemical shift assignment on the BMRB (entry 51023). For the RDC, pulse-field gradient NMR (PFG-NMR) and PRE-NMR experiments, the additional mutation C747V (referred to FLN5 A3A3 V747) was introduced to yield a cysteine-less construct for site-specific spin labelling. The protein concentration was determined using the BCA assay according to the manufacturer’s instructions. RNCs were expressed, isotopically labelled uniformly with 15N, or site-specifically with 19F, and purified as previously described6,7. For samples for intermolecular PRE-NMR experiments involving ribosome labelling, we generated modified E. coli BL21 strains with cysteine mutations in uL23 and uL24 using CRISPR as previously described2. RNC samples were prepared in Tico buffer for experiments. Western blot analyses were undertaken with an anti-hexahistidine horseradish peroxidase-linked antibody (Invitrogen, 1:5,000 dilution).

Fluorescent and PEG-maleimide labelling of 70S and RNC samples

Ribosomes and RNCs were first reduced using 2 mM TCEP overnight at 277 K, then buffer exchanged into labelling buffer. For fluorescein-5-maleimide and PEG-maleimide, labelling was performed in Tico at pH 7.5. ABD-MTS labelling was performed in labelling buffer (50 mM HEPES, 12 mM MgCl2, 20 mM NH4Cl, 1 mM EDTA, pH 8.0). Samples were labelled using a 10x molar excess of ABD-MTS, or fluorescein-5-maleimide. Cysteine mass-tagging by PEGylation was performed as previously described with 10,000-fold molar excess of PEG over sample7. ABD-MTS and PEGylation reactions were analysed using 12% Bis-Tris SDS–PAGE gels61. The fluorescein-labelled reactions were run on a 20% Tricine SDS–PAGE gel, modified from ref. 62.

NMR spectroscopy

All NMR experiments were recorded with Topspin 3.5pl2. NMR experiments of FLN5 A3A3 were performed in Tico buffer at pH 7.5 and 283 K. Chemical shifts were previously assigned7 and obtained from data recorded on a Bruker Avance III operating at 700 and 800 MHz equipped with TCI cryoprobes. All samples contained 10% (v/v) D2O and 0.001% (w/v) DSS as a reference. Data were processed analysed using NMRPipe63 (v11.7), CCPN64 (v2.4) and MATLAB (R2017a, Mathworks).

Amide 1H and 15N chemical shifts were obtained from two-dimensional 1H–15N SOFAST-HMQC experiments65 using an acquisition time of 50 ms in the direct dimension. The inter-scan delay was 50 ms. Cα chemical shifts were obtained from 3D BEST-HNCA experiments recorded at 800 MHz with acquisition times of ~50 ms and inter-scan delays of 150 ms. C’ chemical shifts were obtained from BEST HNCO experiments recorded at 700 MHz using acquisition times of ~50 ms and inter-scan delays of 200 ms. RNC samples were doped with 20 mM NiDO2A (Ni(ii) 1,4,7,10-tetraazacyclododecane-1,7-bis(acetic acid)) to enhance sensitivity66. Cosine-squared window functions were used in processing the spectra.

For PRE-NMR experiments, we used a cysteine-less construct with the C747V mutation and introduced six and eight labelling sites in the isolated and ribosome-bound protein, respectively. Samples were reduced overnight at 277 K in Tico supplemented with 2 mM TCEP. TCEP was then removed by buffer exchange into labelling buffer (50 mM HEPES, 12 mM MgCl2, 20 mM NH4Cl, 1 mM EDTA, pH 8.0) and subsequently labelled overnight at 277 K with 10× molar excess of MTSL. Following labelling, excess MTSL was removed by buffer exchanging the sample back into Tico buffer for NMR. The same labelling protocol was used for isolated protein and RNC samples. To measure the PREs, we recorded the signal intensities with MTSL in the paramagnetic and diamagnetic state. Direct measurements of relaxation rates proved not feasible for RNC samples due sensitivity limitations. 2D 1H–15N SOFAST-HMQC experiments65 were recorded at 800 MHz and 283 K using ~100 μM of protein or ~10 μM of RNC. Experiments were recorded with an acquisition time of 100 ms and 35 ms, in the direct and indirect dimension, respectively. The inter-scan delay was 450 ms to allow for complete relaxation. To acquire the diamagnetic data, the sample was reduced with 2.5 mM (RNC) or 100× molar excess (isolated) sodium ascorbate. Following complete reduction, the same HMQC experiment was recorded. To extract the PREs, spectral peaks were first fitted to a Lorentzian shape in both the direct and indirect dimension using NMRPipe63. Errors were obtained from the spectral noise (RMSE). From the fitted peaks, intensity ratios of Ipara/Idia were calculated and converted to PRE rates for Bayesian ensemble reweighting by numerically solving equation S34 (see Supplementary Notes 3–4) for Γ2. Sample integrity was monitored using interleaved 1H,15N SORDID diffusion measurements as previously described7.

PFG-NMR experiments were used to measure the diffusion coefficients and the Rh of FLN5 variants. 1D 1H,15N-XTSE diffusion measurements were recorded at 700 (FLN5, FLN5 Δ6) and 800 MHz (FLN5 A3A3). Eight to sixteen gradient strengths ranging linearly from 5% to 95% of the maximum gradient strength of 0.556 T m−1 were used. By measuring the signal intensity at each gradient strength, diffusion coefficients could be obtained by fitting the data to the Stejskal–Tanner equation67, which were converted to Rh using the Stokes–Einstein equation.

RDCs for isolated FLN5 A3A3 were measured in Tico buffer at 283 K and pH 7.5 in a PEG/octanol mixture68. RDCs are reported as the splitting of the isotropic splitting subtracted from the aligned splitting, corrected for the negative gyromagnetic ratio of 15N. RDCs were measured by preparing a solution containing 4.6% (w/w) pentaethylene glycol monooctyl ether (C8E5), 1-octonal (molar ratio 1-octanol:C8E5 = 0.94) and 110 μM of protein. Alignment was confirmed by measuring the D2O deuterium splitting at 283 K (17.6 Hz). All RDC NMR experiments were acquired on a Bruker Avance III HD 800 MHz spectrometer equipped with a TCI cryoprobe. A set of four different RDCs (1DNH, 1DCαCO, 1DCαHα and 2DHNCO) was measured per sample (isotropic and anisotropic) using the 3D BEST HNCO (JCOH and JCC) or BEST-HNCOCA (JCAHA) experiments69,70,71. The one-bond 1H−15N coupling was determined by recording two 15N-HSQC sub-spectra, in-phase (IP) and anti-phase (AP). For the measurement of the 1H-13CO coupling constants a BEST HNCO-JCOH experiment was used with an introduced DIPSAP filter. Such J-mismatch compensated DIPSAP spin-state filter offers an attractive approach for accurate measurement of small spin–spin coupling constants72. For that, three separate experiments were recorded with different filter lengths (2τ = 1/J) for each anisotropic and isotropic media, where the sub-spectra associated to the separated spin states (two in phase and one anti-phase) are combined using a linear relation k (IP) + (k − 1) (IP) ± (AP) with k = 0.73, the theoretical optimized scaling factor. The spectra were recorded with 144 × 104 × 1,536 complex points in the 13C(t1)/15N (t2)/1H (t3) dimensions, respectively, and with the spectral widths set to 15,244 Hz (1H), 2,070 Hz (15N) and 1,510 Hz (13C) for the HNCO-JCOH. For the HNCO-JCC and HNCOCA-JCAHA 256 × 200 × 1,536 complex points were acquired in the 13C(t1)/15N (t2)/1H (t3) dimensions, with spectral widths of 15,244 Hz (1H), 1,900 Hz (15N) and 1,214 Hz/5050 Hz (13C). The recycle delay was set to 200 ms, the acquisition time to 100 ms with 16 scans per increment, and the data was acquired in the non-uniform sampling format (2246 points for HNCO-JCOH and 7680 for the HNCO-JCC/HNCOCA-JCAHA experiments were sampled using the schedule generator from the web portal nus@HMS (http://gwagner.med.harvard.edu/intranet/hmsIST/). The time domain data was converted into the NMRPipe63 format and reconstructed using the sparse multidimensional iterative lineshape-enhanced method (SMILE)73. Coupling constants were obtained from line splitting in the 13C or 15N dimension obtained with CCPN analysis software64.

19F NMR experiments were recorded on a 500 MHz Bruker Avance III spectrometer equipped with a TCI cryoprobe at 298 K (unless otherwise indicated) using a 350 ms acquisition time and 1.5–3 s recycle delay as previously described6. We used an amber-suppression strategy to incorporate the unnatural amino acid tfmF, as previously described6. Multiple experiments were recorded in succession to monitor sample integrity over time also as previously described6. Data were processed using NMRPipe63. Spectra were baseline corrected, peaks were fit to Lorentzian functions and errors of the linewidths and integrals (that is, populations) were estimated using bootstrapping (200 iterations, calculating the standard error of the mean), or from the spectral noise for states whose resonance was not detectable, in MATLAB6. 19F-translational diffusion experiments were performed as previously described6.

Thermodynamic parameters of folding (ΔH, ΔS and ΔCp) were obtained from a nonlinear fit to a modified Gibbs–Helmholtz equation, assuming ΔCp remains constant across the experimental temperature range:

$${\rm{ln}}\left({K}_{{\rm{eq}},T}\right)=-\left(\frac{\Delta {H}_{{T}_{0}}+\Delta {C}_{{\rm{p}}}(T-{T}_{0})}{R}\right)\left(\frac{1}{T}\right)+\left(\frac{\Delta {S}_{{T}_{0}}+\Delta {C}_{{\rm{p}}}{\rm{ln}}\left(\frac{T}{{T}_{0}}\right)}{R}\right)$$

(S1)

Keq is the equilibrium constant, T is the temperature in Kelvin and T0 is the standard temperature (298 K). We also fitted the data to the linear van’t Hoff equation (assuming ΔCp = 0).

$${\rm{ln}}\left({K}_{{\rm{eq}}}\right)=-\left(\frac{\Delta H}{R}\right)\left(\frac{1}{T}\right)+\frac{\Delta S}{R}$$

(S2)

The Scipy package with optimize.curve_fit function was used to perform the fits74 and errors were estimated as one s.d. from the diagonal elements of the parameter covariance matrix. All parameters (ΔH, ΔS, ΔCp) generally showed strong correlations with each other (r ≥ 0.8), and thus, their uncertainties correlate also. These parameter correlations are expected75. The magnitudes of ΔH and −TΔS are also expected to correlate because we study the temperature dependence of folding in a range where ΔG of folding is close to 0.

Folding free energies were calculated from the experimental populations using ΔG = −RTln(K). The folding free energy of the FLN5+67 wild-type RNC, where no unfolded state is observable, was estimated on the basis of two destabilizing mutants FLN5(V664A/F665A) and FLN5(V707A). The stability of these mutants was measured using 19F NMR on and off the ribosome. The FLN5+67 wild-type folding energy (ΔGN-U) was then calculated as the average from the V664A/F665A and V707A mutants using ΔGWT,+67 = ΔGmut,+67 − ΔΔGmut-WT,iso, where ΔΔGmut-WT,iso is the experimentally measured destabilization in isolation. Given that at FLN5+34, both mutants show a weaker destabilization than in isolation, we reasoned that this estimate of the FLN5+67 wild-type folding free energy is its lower bound (most negative).

19F transverse relaxation rate (R2) measurements were recorded using a Hahn-echo sequence and acquired as pseudo-2D experiments with relaxation delays of 0.1 to 200 ms. Data were processed using NMRPipe and analysed using MATLAB. Data were fit to lineshapes and R2 was obtained by fitting the integrals to single exponential functions. We also orthogonally determined R2 from linewidth measurements of spectra acquired by 1D 19F pulse-acquire experiments, which showed excellent correlations. The lineshape-derived R2 values also showed good correlation with previously determined rotational correlation times25 (τC). We additionally determined the S2τC of FLN5 in 60% glycerol at 278 K by measurements of triple quantum build-up and single quantum relaxation as previously described76,77. Thus, our R2 values can be used to determine rotational correlation times (τC,exp). The obtained τC,exp was used to estimate the bound population as \(\frac{{\tau }_{{\rm{C}},\exp }-{\tau }_{{\rm{C,iso}}}}{{\tau }_{{\rm{C,bound}}}-\,{\tau }_{{\rm{C,iso}}}}\), where τC,iso is the rotational correlation time of the isolated protein25 (7.7 ns at 298 K) and τC,bound is the expected rotational correlation time of the bound state. τC,bound is taken as the rotational correlation time of the ribosome itself (~3,000 ns at 298 K) for a fully rigid bound state. From the bound populations (pB), the resulting change in the folding free energies of the intermediates was calculated as ΔΔGI-U,RNC-iso = RT(ln(1 − pB)). We report the estimate for a fully rigid bound state in the main text (\({S}_{{\rm{bound}}}^{2}=1.0\)) but note that even one order of magnitude more flexibility in the bound state (\({S}_{{\rm{bound}}}^{2}=0.1\)) only accounts for up to 1.1 ± 0.6 and 0.4 ± 0.1 kcal mol−1 of stabilization for I1 and I2 on the ribosome at FLN5+47, respectively. These estimates still cannot account for the >4 kcal mol−1 of intermediate stabilization observed on the ribosome6.

All NMR experiments of RNCs are recorded and continuously interleaved with a series of 1D 1H/19 F spectra and 1H,15N/19F diffusion measurements6,7,61,78. These provide the most sensitive means to assess changes in the sample, and when alterations in signal intensities or linewidths (that is, transverse relaxation rates), chemical shifts or translational diffusion measurements of the nascent chain are observed, data acquisition is halted. Only data corresponding to intact RNCs are summed together and subjected to a final round of analysis. Where signal-to-noise remains low, datasets from multiple samples are compared to ensure identical spectra, before summation together into a single NMR spectrum. Biochemical assays provide an orthogonal means to assess nascent chain attachment to the ribosome. Identical samples incubated in parallel with NMR samples are analysed by SDS–PAGE (under low pH conditions61) and detected with nascent-chain-specific antibodies. Ribosome-bound species migrate with an addition ~17-kDa band-shift relative to released nascent chains due to the presence of the tRNA covalently linked to the nascent chain. Combined with time-resolved NMR measurements, these analyses confirm that the reported NMR resonances originate exclusively from intact RNCs.

Mass spectrometry

FLN5 A3A3 was buffer exchanged into 100 mM (NH4)2CO3 at pH 6.8 (using formic acid for pH adjustment). Analyses were run on the Agilent 6510 QTOF LC–MS system at the UCL Chemistry Mass Spectrometry Facility. Samples contained ~10–20 μM of protein and 10 μl were injected onto a liquid chromatography column (PLRP-S, 1,000 Å, 8 μm, 150 mm × 2.1 mm, maintained at 60 °C). The liquid chromatography was run using water with 0.1% formic acid as mobile phase A and acetonitrile with 0.1% formic acid as phase B with a gradient elution and a flow rate of 0.3 ml min−1. ESI mass spectra were continuously acquired. The data were processed to zero charge mass spectra with the MassHunter software, utilizing the maximum entropy deconvolution algorithm.

Small-angle X-ray scattering

We measured SAXS of an isolated FLN5 A3A3 C747V sample in Tico buffer supplemented with 1% (w/v) glycerol. Data collection was performed at the DIAMOND B21 beamline (UK)79 with a beam wavelength of 0.9408 Å, flux of 4 × 1012 photons s−1 and an EigerX 4 M (Dectris) detector distanced at 3.712 m from the sample. A capillary with a 1.5 mm diameter kept at 283 K was used for data acquisition. We acquired SAXS data at multiple protein concentrations (5.5, 2.75, 1.38, 0.69, 0.34 and 0.17 mg ml−1) to assess whether the sample exhibited signs of aggregation or interparticle interference. At 5 mg ml−1, we observed weak signs of interparticle interference in the low q region of the scattering profile, which is also reflected in the Rg obtained by Guinier analysis (using the autorg tool from ATSAS80; Supplementary Table 1). Data were recorded as a series of frames, non-defective frames were averaged, and buffer subtracted with PRIMUS80. Size-exclusion chromatography–SAXS (SEC–SAXS) experiments were additionally performed in Tico buffer with 1% (w/v) glycerol using a KW402.5 (Shodex) column to confirm the monodispersity of the sample. We chose the 2.75 mg ml−1 dataset as the final dataset to compare with our molecular dynamics simulations. This dataset exhibited the highest signal to noise ratio and did not show signs of interparticle interference, and accordingly, the Rg obtained from the 2.75 mg ml−1 dataset is consistent with the value obtained from lower concentrations and the main SEC–SAXS peak (Supplementary Table 1).

Circular dichroism spectroscopy

The circular dichroism (CD) spectrum of isolated FLN5 A3A3 V747 was acquired in 10 mM Na2HPO4 pH 7.5 at 283 K. A Chirascan-plus CD spectrometer (Applied Photophysics), a protein concentration of 44 μM and a cuvette with a 0.5 cm pathlength were used.

HRAS refolding experiments

HRAS refolding experiments were performed with the HRAS G-domain (residues 1–166). The protein was unfolded overnight at 298 K in Tico buffer with 2 mM β-mercaptoethanol, 8 M urea and protein concentration of 15 μM. The protein was then refolded by rapidly diluting into Tico buffer (supplemented with 2 mM β-mercaptoethanol and 50 μM GDP) to reach final urea and protein concentrations of 0.94 M and 1.76 μM, respectively, and allowed to incubate at 298 K for 24 h. For NMR analyses of refolded samples, we prepared 18 μM of refolded protein with the same urea concentrations and dilutions.

We assayed the functional/activity state of HRAS using GDP/GTP nucleotide exchange (‘activity’) assay81 with fluorescently labelled GTP that exhibits higher fluorescence when bound to HRAS than free in solution (BODIPY FL GTP, ThermoFisher). 0.4 μM of HRAS, 0.01 μM of BODIPY GTP and 1 μM of SOScat (the catalytic domain of Son of sevenless) were incubated at room temperature and the maximum (plateau) fluorescence recorded and normalized by the signal of the buffer (signal/noise ratio). SOScat was produced as previously described82. Fluorescence measurements were performed using the CLARIOstar microplate reader (BMG Labtech) with excitation and emission wavelengths set to 488 and 514 nm, respectively.

The proteolytic stability of HRAS was assayed with thermolysin at a concentration of 0.05 mg ml−1 incubated with HRAS samples over the course of 5 h in vitro and 9 h in rabbit reticulocyte lysate (RRL, TNT coupled reticulocyte lysate, Promega). Reactions were quenched with 23 mM EDTA. Timepoints were analysed by western blot analysis using a pan-RAS polyclonal antibody (ThermoFisher, 1:1,000 dilution), utilizing an anti-rabbit IgG horseradish peroxidase-linked secondary antibody (Cell Signaling Technology, 1:1,000 dilution). Densitometry analyses were performed with ImageJ83. For the RRL experiments, refolding reactions were performed in RRL for 24 h at 298 K and a final HRAS concentration of 1.6 μM followed by pulse proteolysis and we quantified the relative band intensities (refolded/control) for each time point to account for increased background on the western blot during the proteolysis reaction.

Molecular dynamics simulations

We used the FLN5 A3A3 C747V sequence for all simulations. A reliability and reproducibility checklist is provided in Supplementary Table 8. GROMACS (version 2021)84 was used for all all-atom molecular dynamics simulations in explicit solvent. We employed the Charmm36m force field in combination with the CHARMM TIP3P water model (C36m) and the CHARMM TIP3P water model with an increased water hydrogen LJ well-depth (denoted here as C36m+W)85. We also used the a99sb-disp force field together the a99sb-disp TIP4P-D water model86. Default protonation states were used in all cases. Starting from a random extended conformation, for all force field combinations the system was solvated in a dodecahedron box with 151,135 water molecules and 12 mM MgCl2 (resulting in 455,116 atoms and an initial box volume of 4,688 nm3). Systems were then energy minimized using the steepest-decent algorithm. For the following dynamics simulations, we used the LINCS algorithm87 to constrain all bonds connected to hydrogen and a timestep of 2 fs using the leap-frog algorithm for integration. Nonbonded interactions were calculated with a cut-off at 1.2 nm (including a switching function at 1.0 nm for van der Waals interactions) and the particle mesh Ewald (PME) method88 was used for long-range electrostatic calculations. We then equilibrated the systems in two phases. First, we performed a 500 ps equilibration simulation in the NVT ensemble with position restraints on all protein heavy atoms. The temperature was kept at 283 K using the velocity rescaling algorithm89 and a time constant of 0.1 ps. Next, we further equilibrated the systems for 500 ps in the NPT ensemble at 283 K and a pressure of 1 bar with a compressibility of 4.5 × 10−5 bar−1 using the Berendsen barostat90. Following equilibration, we relaxed our initial structure for 100 ns at 283 K without any position restraints using the Parrinello–Rahman algorithm91 and then picked five structures from this simulation for production simulations. We ran a total of 5× 2 μs (with different initial coordinates and velocities) yielding a total of 10 μs of sampling per force field. For the C36m+W combination we ran an additional 5× 2 μs starting from 5 new starting structures yielding 20 μs in total.

We also generated a prior ensemble with a physics-based coarse-grained (C-alpha) model. We generated the C-alpha model template from the FLN5 crystal structure using SMOG 2.392, where all bonded terms have a global energy minimum at the values taken in the crystal structure93. Nonbonded van der Waals interactions were modelled using a 10–12 Lennard-Jones potential with σ and λ parameters described in the M1 parameter determined by Tesei et al.94 (equation (S3)). We used the arithmetic mean of two residues to determine σ and λ. Electrostatic interactions were modelled using the Debye–Hückel theory with parameters described previously7. Interactions between Cα beads separated by less than four residues were excluded. We ran initial simulations at a range of reduced temperatures to determine the effect on the average compactness and ran final simulations at a reduced temperature of 1.247 (150 K in GROMACS) as we did not observe a significant increase in average Rg beyond this temperature. Simulations were run for a total of 3 × 109 steps with GROMACS (v2018.3).

$${u}_{{\rm{LJ}}}=\mathop{\sum }\limits_{i}^{N}\lambda \left[{5\left(\frac{\sigma }{r}\right)}^{12}-6{\left(\frac{\sigma }{r}\right)}^{10}\right]$$

(S3)

After simulations, the coarse-grained ensemble was backmapped to an all-atom structure using PULCHRA (v3.06)95.

RNC simulations were parameterized using the C36m+W force field/water model combination85,96. We modelled the ribosome using the structure PDB 4YBB97 as a template, which we previously refined against a cryo-EM map containing an FLN5 RNC98. As in our previous work, we only retained ribosome atoms around the nascent chain exit tunnel and accessible surface outside the vestibule7. The FLN6 linker and SecM sequence were initially modelled using a cryo-EM map of a FLN5+47 RNC (Mitropoulou et al., manuscript in preparation). The rest of the nascent chain (MHHHHHAS N-terminal tag and FLN5) was then built using PyMol version 2.3 (The PyMol Molecular Graphics System, Schrödinger) and we generated a random initial starting structure with a short simulation using a structure-based force field, SMOG2.392, without native contacts. The FLN5+31 A3A3 RNC (containing the C747V mutation) complex was then centred in a dodecahedral box, solvated using 1,030,527 water molecules and neutralized with 706 Mg2+ ions, resulting in a final system size of 3,163,127 atoms. The initial box volume was 32,117 nm3. The large box size was necessary to accommodate the highly expanded unfolded state. We then used the same cut-offs and simulation methods as for the isolated protein. We initially also ran a 500 ps equilibration simulation in the NVT ensemble using position restraints on all heavy atoms using a force constant of 1,000 kJ mol−1 nm2 in along the x, y and z axes. We used a temperature of 283 K, which was held constant using the velocity rescaling algorithm89 and a time constant of 0.1 ps. Then, we ran a 500 ps equilibration simulation in the NPT ensemble at 283 K using the same position restraints. The pressure was kept at 1 bar with a compressibility of 4.5 × 10−5 bar−1 using the Berendsen barostat90. The position restraints for all nascent chain atoms (except the terminal residue at the PTC in the ribosome) were then removed, while the ribosome atoms kept being position restrained. In this setup, we ran a 1 ns equilibration simulation at 283 K and 1 bar, using the Parrinello–Rahman algorithm91. All production simulations were performed using position restraints for the ribosome atoms and C-terminal nascent chain residue at the PTC. Using the equilibrated configuration, we then ran two simulations of ~100 ns to picked ten starting structures for production simulations. Then, ten production simulations of 1.5 μs each (15 μs) were initiated from these different starting structures using random initial velocities. Before the production simulation, each structure was re-equilibrated at 283 K and 1 bar with a 500 ps NVT and 500 ps NPT simulation.

Lastly, to compare our C36m+W simulations with a model that only considers steric exclusion as a nonbonded interaction, we also ran simulations of a simple all-atom model, based on a structure-based model template92. We used the FLN5 crystal structure to define the energy minima of all bond and dihedral angles and removed all native contacts. Simulations of isolated and ribosome-bound FLN5 A3A3 were run for 1 × 109 steps and 100,000 frames were harvested for analysis. This ensemble was used to compare the expansion of the ensemble, ribosome interactions and conformational entropy with the C36m+W simulations.

Calculation of PREs

The transverse PRE rates of backbone amide groups, Γ2, were back-calculated from the ensembles using the Solomon–Bloembergen equation99,100

$${\varGamma }_{2}=\frac{1}{15}{\left(\frac{{\mu }_{0}}{4\pi }\right)}^{2}{\gamma }_{{\rm{H}}}^{2}\,{g}_{{\rm{e}}}^{2}{\mu }_{{\rm{B}}}^{2}S(S-1)[4J(0)+3J({\omega }_{{\rm{H}}})]$$

(S4)

where μ0 is the permeability of space, γH is the gyromagnetic ratio of the proton, ge is the electron g-factor, γB is the Bohr magneton, S is the proton nuclear spin and J(ω0) is the generalized spectral density function. For flexible spin labels attached via rotatable bonds the spectral density can be expressed as in equation (S5)101.

$$J({\omega }_{{\rm{H}}})=\langle {r}^{-6}\rangle \left[\frac{{S}^{2}{\tau }_{{\rm{c}}}}{1+{({\omega }_{{\rm{H}}}{\tau }_{{\rm{c}}})}^{2}}+\frac{(1-{S}^{2}){\tau }_{{\rm{t}}}}{1+{({\omega }_{{\rm{H}}}{\tau }_{{\rm{t}}})}^{2}}\right]$$

(S5)

where \(\langle {r}^{-6}\rangle \) is the average of the electron–hydrogen distance (r) distribution, S2 is the generalized order parameter for the electron–hydrogen interaction vector, τC is the correlation time defined in terms of the rotational correlation time of the protein (τr) and the electron spin relaxation time (τs):

$${\tau }_{c}={\left({\tau }_{r}^{-1}+{\tau }_{s}^{-1}\right)}^{-1}$$

(S6)

τt is the total correlation time defined as:

$${\tau }_{{\rm{t}}}={\left({\tau }_{r}^{-1}+{\tau }_{s}^{-1}+{\tau }_{{\rm{i}}}^{-1}\right)}^{-1}$$

(S7)

τi is the internal correlation time of the spin label. Since for nitroxide labels electron spin relaxation occurs on a much slower timescale than rotational tumbling101,102, τC can be approximated to τr such that expression for τt simplifies to

$${\tau }_{{\rm{t}}}={\left({\tau }_{{\rm{C}}}^{-1}+{\tau }_{i}^{-1}\right)}^{-1}$$

(S8)

Given that τC is not known a priori, we iteratively scanned τC values in the range of 1 to 15 ns to find a value for which optimal agreement with the experimental data is achieved (as judged by the reduced χ2)94,103. The spin label correlation time104, τi was set to 500 ps, in agreement with molecular dynamics simulations105 and electron spin resonance measurement106.

The generalized order parameter S2 for the electron–hydrogen interaction vector can be decomposed into its radial and angular components107:

$${S}_{{\rm{PRE}}}^{2}\approx {S}_{{\rm{PRE}},{\rm{angular}}}^{2}{S}_{{\rm{PRE}},{\rm{radial}}}^{2}$$

(S9)

where the individual components are defined as

$${S}_{{\rm{PRE,angular}}}^{2}=\frac{4\pi }{5}\mathop{\sum }\limits_{m=-2}^{2}{\left|\langle {Y}_{2}^{m}({\varOmega }^{{\rm{mol}}})\rangle \right|}^{2}$$

(S10)

$${S}_{{\rm{PRE}},{\rm{radial}}}^{2}=\langle {r}^{-6}{\rangle }^{-1}\langle {r}^{-3}{\rangle }^{2}$$

(S11)

and \({Y}_{2}^{m}\) are the second order spherical harmonics and Ωmol are the Euler angles in the frame. A weighted ensemble average of S2 can be calculated by taking a weighted ensemble average of the individual radial and angular components.

A previously published rotamer library containing 216 MTSL rotamers108 was used to explicitly model the flexibility of the spin label, similar to other existing methods109,110. The rotamer library was aligned to all employed labelling sites for each conformer using the backbone atoms of the labelling site and Cys-MTSL moiety. Clashing rotamers were discarded, where a steric clash between the rotamer and the protein was defined using a 2.5 Å cut-off distance. Only backbone and Cβ atoms were considered for the protein, assuming sidechains can rearrange to accommodate the MTSL rotamer111. For MTSL, only the sidechain was included (heavy atoms beyond the Cβ atom). Protein frames for which at least one labelling position cannot sterically allow any MTSL rotamers were discarded. The rotamer library was used to calculate a weighted ensemble-averaged Γ2 over the rotamer ensemble for each protein conformer in the protein ensemble using equations (S3−S11). The protein ensemble average can then be calculated by averaging Γ2 over the ensemble.

PRE intensity ratios were then calculated from the ensemble-averaged PRE rate, \(\langle {\varGamma }_{2}\rangle \), using

$$\frac{{{\rm{I}}}_{{\rm{p}}{\rm{a}}{\rm{r}}{\rm{a}}}}{{{\rm{I}}}_{{\rm{d}}{\rm{i}}{\rm{a}}}}=\frac{{R}_{2}{{\rm{e}}}^{-2\Delta \langle {\varGamma }_{2}\rangle }}{{R}_{2}+\langle {\varGamma }_{2}\rangle }\times \frac{{R}_{2,{\rm{M}}{\rm{Q}}}}{{R}_{2,{\rm{M}}{\rm{Q}}}+\langle {\varGamma }_{2}\rangle }$$

(S12)

where R2 is the linewidth in the proton dimension (residue-specific), R2,MQ is the linewidth in the nitrogen dimension (multiple-quantum term) and Δ is the delay time in the HMQC experiment (5.43 ms). See Supplementary Note 3 for additional details.

For RNCs, we considered that that ribosome tethering may increase the correlation time of the electron–amide interaction vector due to restricted molecular tumbling near the exit tunnel. We therefore calculated an order parameter, \({S}_{{\rm{NC}}}^{2}\), which quantifies the motion of the electron-interaction vector over the entire nascent chain conformer ensemble (\({S}_{{\rm{NC}}}^{2}\) is distinct from the order parameter S2 that quantifies the motion of the MTSL rotamer library attached to a labelling site for a specific protein conformer; equation (S9)). S2 is given by

$${S}_{{\rm{NC}}}^{2}\approx {S}_{{\rm{NC,angular}}}^{2}{S}_{{\rm{NC,radial}}}^{2}$$

(S13)

where \({S}_{{\rm{NC}},{\rm{angular}}}^{2}\) and \({S}_{{\rm{NC}},{\rm{radial}}}^{2}\) are given by

$${S}_{{\rm{PRE}},{\rm{angular}}}^{2}=\frac{4\pi }{5}\mathop{\sum }\limits_{m=-2}^{2}{| \langle {Y}_{2}^{m}({\varOmega }^{{\rm{mol}}})\rangle | }^{2}$$

(S14)

$${S}_{{\rm{PRE}},{\rm{radial}}}^{2}=\langle {r}^{-6}{\rangle }^{-1}\langle {r}^{-3}{\rangle }^{2}$$

(S15)

and \({Y}_{2}^{{\rm{m}}}\) are the second order spherical harmonics and Ωmol are the Euler angles in the frame. We approximated the position of the free electron with the Cα atom of the labelling site in this case. A \({S}_{{\rm{NC}}}^{2}\) value of 0 indicates that the vector tumbles completely independent of the ribosome and that the correlation time of the electron–amide vector is the same as for the isolated protein, τC,iso. A \({S}_{{\rm{NC}}}^{2}\) value of 1 means that the vector tumbles with the same rotational correlation time as the ribosome (τr,70S = 3.3 μs per cP, as determined by fluorescence depolarization112, and τr,70S = 4.3 μs at 283 K in H2O). The effective correlation time, τC,eff, of each amide-electron vector is given by

$${\tau }_{{\rm{C,eff}}}\,={S}_{{\rm{NC}}}^{2}{\tau }_{{\rm{r,70S}}}+(1-{S}_{{\rm{NC}}}^{2}){\tau }_{{\rm{C,iso}}}$$

(S16)

We used a value of 3 ns for τC,iso, which was the optimal value determined for isolated FLN5 A3A3. Generally, τC (equation (S6)) is approximated as τC ≈ τr because the electron spin relaxation time, τs, occurs on a much slower timescale. In fact, measurements of the spin relaxation time of nitroxides have been measured to be on a timescale from hundreds of nanoseconds to several microseconds113,114,115. The calculated values of τC,eff are predominantly below 100 ns except for labelling sites C744, uL23 G90C and uL24 N53C, where values of up to ~250 ns are observed (Supplementary Tables 5 and 6). Thus, we still expect τC to be dominated by τr and make use of the τC ≈ τ approximation.

Finally, reference PRE profiles for a fully extended peptide were calculated from a linear polyalanine chain using a τC of 5 ns and R2,H/R2,MQ of 100 Hz.

Bayesian inference reweighting

We performed ensemble refinement by reweighting the molecular dynamics-derived ensembles against the experimentally deduced Γ2 rates using the Bayesian Inference of Ensembles (BioEn) software and method described in the corresponding paper116,117. These calculations were performed using in-house scripts of the software with the modification to incorporate upper and lower bound restraints in addition to regular restraints with gaussian errors. To this end, these inequality restraints were treated as normal gaussian restraints but subjected to a conditional statement. Lower bound restraints (Γ2 > 64.5 s−1 for isolated FLN5 A3A3; Γ2 > 96.0 s−1 for the RNCs) were applied only if the back-calculated Γ2 was below the lower bound value. Similarly, upper bound restraints (Γ2 < 2.2 s−1 for isolated FLN5 A3A3; Γ2 < 3.7 s−1 for the RNCs) were applied only if the back-calculated average was above the upper bound. This effectively allows the back-calculated value to vary freely above the lower bound and below the upper bound but imposes a penalty if the inequality condition is not met. The errors of the lower and upper bound values were taken as the combined relative error of that datapoint (that is, the intensity ratio).

As described by Köfinger et al.117, the reweighting optimization problem can be efficiently solved by minimizing the negative log-posterior function (L).

$$L=\theta {S}_{{\rm{KL}}}+\mathop{\sum }\limits_{i=1}^{M}\frac{{\left({\sum }_{\alpha =1}^{N}{w}_{\alpha }{y}_{i}^{\alpha }-{Y}_{i}\right)}^{2}}{2{\sigma }_{i}^{2}}$$

(S17)

θ expresses the confidence in the initial ensemble, N is the ensemble size, M is the number of experimental restraints, wα is the vector of weights for the conformers in the ensemble, \({y}_{i}^{\alpha }\) is the back-calculated experimental value i, Yi is the experimental restraint i, σi is the uncertainty of experimental restraint i, and SKL is the Kullback–Leibler divergence defined as

$${S}_{{\rm{KL}}}=\mathop{\sum }\limits_{\alpha =1}^{N}{w}_{\alpha }{\rm{ln}}\frac{{w}_{\alpha }}{{w}_{\alpha }^{0}}$$

(S18)

\({w}_{\alpha }^{0}\) is the vector of initial weights (which were uniform). We used the log-weights method to minimize the negative log-posterior117 and performed reweighting calculations for a range of θ values, as the optimal value of θ cannot be known a priori. Therefore, we conduct L-curve analysis117,118 by plotting SKL (entropy) on the x axis and the goodness of fit, quantified by the reduced χ2 value, on the y axis. The reduced χ2 was calculated against the experimental intensity ratios (Ipara/Idia). This is an effective method to prevent overfitting and introducing a minimal amount of bias into the prior ensemble117,119. After reweighting, we also calculated the effective fraction of frames contributing to the ensemble average119 as an indication of the extent of fitting.

$${N}_{{\rm{eff}}}=\exp (-{S}_{{\rm{KL}}})$$

(S19)

For RNCs, we used the same approach with an additional modification. Since the PRE depends on τC,eff and \({S}_{{\rm{NC}}}^{2}\) which are a function of the weights of individual structures in the ensemble, this consequently leads to changes in τC,eff and \({S}_{{\rm{NC}}}^{2}\) when reweighting is performed. Therefore, the conformer-specific PRE values that were used for reweighting are not the same anymore after reweighting. To account for this, we performed 20 iterative rounds of reweighting where each additional round receives input weights and τC,eff from the previous round. We found that this leads to convergence of the weights and conformer-specific PREs.

We found that for the ribosomal labelling sites, uL23 G90C and uL24 N53C, the reweighting results are sensitive to the specific ribosome structure used to fit the MTSL rotamer library to, since small variations in the local structure of the labelling site can lead to different rotamer distributions. We tested two different rotamer distributions for the ribosomal labelling sites (Extended Data Fig. 5a), finding that one of them (referred to as R2) gives better agreement with the intermolecular PRE data after reweighting and fits better into the expected density or MTSL rotamers when rotamers are fitted to ten high-resolution ribosome structures (Extended Data Fig. 5a). The R2 rotamer distribution is more representative of the expected variation from structural changes in the labelling sites and was therefore used for our final reweighting calculations.

Calculation of RDCs, R
h and chemical shifts from molecular dynamics simulations

To back-calculate the Rh from static structures we used an approximate relationship between Rg and Rh values120, the latter being calculated from the programme HYDROPRO121. Thus, we calculated the Rg from Cα atoms using MDAnalysis122 and then converted it to Rh using

$${R}_{{\rm{h}}}=\frac{{R}_{{\rm{g}}}}{\frac{{\alpha }_{1}\left({R}_{{\rm{g}}}-{\alpha }_{2}{N}^{0.33}\right)}{{N}^{0.60}-{N}^{0.33}}+{\alpha }_{3}}$$

(S20)

N is the number of amino acids, α1 takes a value of 0.216 Å−1, α2 takes a value of 4.06 Å, and α3 has a value of 0.821. The estimated value of Rh (relative to the HYDROPRO calculation) has an average relative uncertainty120 of 3%. HYDROPRO itself has a relative uncertainty of ±4% with respect to experimental values121. Therefore, we treat the back-calculated ensemble-average Rh with a total relative uncertainty of ±5%. The ensemble average was calculated as previously described by Ahmed et al. for back-calculation of PFG-NMR derived values123 of Rh

$$\langle {R}_{{\rm{h}}}\rangle ={\rm{ln}}{\left(\langle \exp \left(-{R}_{{\rm{h}}}^{-1}\right)\rangle \right)}^{-1}.$$

(S21)

Chemical shifts were calculated using the SHIFTX2 software124 and RDCs were calculated using the global alignment prediction method implemented in PALES125. We then scaled the magnitude (that is, the extent of alignment) of the calculated RDCs by a global factor to optimize the Q-factor for each ensemble.

Calculation of SAXS profiles from molecular dynamics simulations

We used Pepsi-SAXS126 to compute the theoretical scattering profiles of each conformer in the molecular dynamics ensembles. We treated the contrast of the hydration layer (δp) and the effective atomic radius (r0) as global parameters and used values of 3.34 e− nm−3 and 1.025 × rm (rm = average atomic radius of the protein) in line with previous work that showed these parameters to well suited for flexible proteins127. The constant background and scale factor were also fitted globally using least-squares regression103,127. The goodness of fit was assessed using the reduced χ2 metric, where n is the number of datapoints, q is the scattering angle, \({I}_{q}^{{\rm{calc}}}\) and \({I}_{{q}}^{\text{exp}}\) are the calculated and experimental scattering intensities, respectively, and σq is the experimental error:

$${\chi }_{r}^{2}=\frac{1}{n}\mathop{\sum }\limits_{{q}}^{n}\frac{{({I}_{{q}}^{\text{calc}}-{I}_{{q}}^{\text{exp}})}^{2}}{{\sigma }_{{q}}^{2}}$$

(S22)

Structural analysis

The Python package MDAnalysis122 and MDTraj128 were used for general analysis of the ensembles involving atomic coordinates. For native contact analysis, we calculated the fraction of native contacts (relative to the native FLN5 crystal structure) as129

$$Q(X)=\frac{1}{N}\sum _{ij}\frac{1}{1+{{\rm{e}}}^{\left(\beta \left({r}_{i,j}-\lambda {r}_{i,j}^{0}\right)\right)}}$$

(S23)

where ri,j and \({r}_{i,j}^{0}\) are the distances between atoms i and j in frame X and the template structure, respectively, β modulates the smoothness of the switching function (default value 5 Å−1 used) and λ is a factor allowing for fluctuations of the contact distance (default value 1.8 used).

Asphericity was calculated using MDAnalysis122 as defined by Dima and Thirumalai130:

$$\Delta =\frac{3}{2}\frac{{\sum }_{n=1}^{3}{\left({\lambda }_{i}-\bar{\lambda }\right)}^{2}}{tr{T}^{2}}$$

(S24)

\(\bar{\lambda }\) represents the mean eigenvalue obtained from the inertia tensor \(\bar{\lambda }=\,\frac{{\rm{trT}}}{3}\).

For the intrachain contact analysis, we defined contacts between Cα–Cα distances of less than 10 Å. The contact features qualitatively were unchanged when using lower cut-off values or when calculating contacts between all heavy atoms. Secondary structure populations were calculated using DSSP131 implemented in MDTraj. The SASA was calculated using GROMACS84. Clustering was also performed in GROMACS using the GROMOS algorithm132 and Cα RMSD cut-offs in the range of 1.2–1.8 nm.

Error analysis from ensembles

Errors from the molecular dynamics ensembles were estimated using a block analysis of the full concatenated ensembles (composed of multiple statistically independent simulations). We performed block analysis for the concatenated ensembles to verify that the estimate of the standard error of the mean (s.e.m.) plateaus/fluctuates at block sizes larger than blocks corresponding to the individual trajectories. The final block size was chosen either in the plateau region of block analysis plots or corresponding to the blocks of the statistically independent simulation (10 independent simulations were run for the isolated and RNC systems with the C36m+W force field, and thus 10 blocks were chosen for block analysis and error estimation). The error after reweighting with PRE-NMR data was calculated the same way using a weighted standard error, where blocks are weighted according to the weights obtained from reweighting with PRE-NMR data. Exemplar block analysis plots are shown in Supplementary Fig. 8.

Energetic analyses from structural ensembles

The conformational entropy was calculated as defined by Baxa et al.133. Proline, glycine and alanine entropies were calculated from the backbone probability distribution Pi(Φ,Ψ). Residues with a maximum of two sidechain torsion angles, Xn, the entropy was calculated from the probability distribution Pi(Φ,Ψ,X1,X2), while residues with more sidechain torsion angles was calculated from the sum of entropies obtained using the Pi(Φ,Ψ,X1), and Pi(Xn), after subtraction of the entropy obtained from Pi(X1). Entropies were calculated from probability distributions using \(S={-k}_{{\rm{B}}}{\sum }_{i=1}^{n}{P}_{i}{\rm{ln}}({P}_{i})\). We used a block analysis from the pooled ensembles (i.e., all individual trajectories concatenated together) to check that the entropy difference between on and off the ribosome is robust with respect to sampling by calculating entropy changes with increasing amounts of total sampling (from the 15 and 20μs of concatenated sampling for the RNC and isolated protein, respectively). The errors were then also estimated from the same sampling/block sizes up to 7.5 μs of molecular dynamics sampling. This is because the estimate of entropy differences trend increases up to total sampling times of 7.5 μs (Extended Data Fig. 6g).

The energetic contributions due to changes in solvation were estimated based on empirical relationships between changes in the polar and apolar accessible surface area75,134 (ΔASApolar and ΔASAapolar). The apolar and polar surface area of the protein were defined based on the atomic partial charges in the C36m force field85. Atoms with an absolute charge of less than or equal to 0.3 were defined as apolar. The change in heat capacity of hydration is related to these quantities by

$$\Delta C=\Delta {C}_{{\rm{apolar}}}+\Delta {C}_{{\rm{polar}}}=\alpha \times \Delta {{\rm{ASA}}}_{{\rm{apolar}}}+\beta \times \Delta {{\rm{ASA}}}_{{\rm{polar}}}$$

(S25)

where α and β are 0.34 ± 0.11 and −0.12 ± 0.12 cal mol−1 K−1 Å−2, respectively. We obtained these values as an average and standard deviation of parameters previously reported in the literature as summarized in ref. 135 to account for the uncertainty of the parameters in addition to the uncertainty coming from conformational sampling in our simulations. The enthalpy change due to solvation is then obtained from75

$$\Delta {H}_{{\rm{solv}}}\left(333\,{\rm{K}}\right)=\gamma \times \Delta {{\rm{ASA}}}_{apolar}+\delta \times \Delta {{\rm{ASA}}}_{polar}$$

(S26)

$$\Delta {H}_{{\rm{solv}}}\left(T\right)=\Delta {H}_{{\rm{solv}}}\left(333\,{\rm{K}}\right)+\Delta C(T-333\,{\rm{K}})$$

(S27)

T is the temperature and γ and δ constants taking on values of −8.44 and 31.4 cal mol−1 Å−2, respectively. While we are not aware of alternative parameter sets for the solvation enthalpy (equation (S26)) in the literature, we treated these parameters with a relative uncertainty of 50% to show that even with such high levels of uncertainty our conclusions are not affected. Finally, the solvation entropy and change in free energy are then calculated using

$$\Delta {S}_{333{\rm{K}},{\rm{solv}}}=\Delta {C}_{{\rm{apolar}}}\,{\rm{ln}}\left(\frac{T}{{T}_{{\rm{apolar}}}}\right)-\Delta {C}_{{\rm{polar}}}\,{\rm{ln}}\left(\frac{T}{{T}_{{\rm{polar}}}}\right)$$

(S28)

$$\Delta {G}_{{\rm{solv}}}=\Delta {H}_{{\rm{solv}}}-T\Delta {S}_{{\rm{solv}}}$$

(S29)

where Tapolar and Tpolar are the temperatures at which ΔSsolv,apolar and ΔSsolv,polar are 0 (385 K and 335 K, respectively). Our previous work indicated ribosome solvation changes during coTF is not a major factor in coTF thermodynamics (see Supplementary Note 9), we estimated the above quantities using surface areas calculated excluding the ribosome. We regard these absolute quantities as an estimated upper bound for ΔGsolv because it is likely that folding intermediates and the native state also interact with the ribosome6, thus effectively cancelling out any reduction in SASA of the unfolded state due to ribosome interactions. However, in the following section we describe an alternative, more direct approach for the solvation entropy that does not rely on this assumption.

Calculation of solvation entropy changes using the 2PT method

The water and solvation entropy changes were also assessed more directly from molecular dynamics simulations using the two-phase thermodynamic (2PT) method136 implemented in the DoSPT code (https://dospt.org/index.php/DoSPT)137. For these calculations, we chose five snapshots from our isolated FLN5 A3A3 V747 simulations detailed above (that is, with different initial protein conformations and solvent configuration) and use these to initiate short molecular dynamics simulations for entropy calculations. We first re-equilibrated the boxes for 10 ns at the target temperature in the NPT ensemble at 1 bar using the velocity rescaling algorithm89 and the Parrinello–Rahman algorithm91 as detailed above and the velocity Verlet integration algorithm (md-vv in GROMACS84). Production simulations were then run in the NVT ensemble at 283 K and 298 K (to assess the effect of temperature on the water entropy calculations) for 20 ps using the md-vv integrator and saving coordinates and velocities for analysis every 4 fs. Control simulations of pure TIP3P (CHARMM TIP3P) water in a cubic box with a box vector length of 5 nm, resulting in 4,055 water molecules. Five independent simulations were performed by first energy minimizing the system using the steepest-decent algorithm. Then, using a 2 fs timestep and thermostat/barostat settings as for the protein and the md-vv integrator we equilibrated the water box first in the NVT ensemble for 1 ns, followed by 1 ns in the NPT ensemble using the Berendson barostat90. The water box was then further equilibrated in the NVT ensemble for 1 ns prior to the production simulation in the NVT ensemble for 20 ps, saving coordinates and velocities every 4 fs. These production simulations were also performed at 283 K and 298 K and then used to calculate the molar entropies of pure water at these temperatures with DoSPT.

For water entropy calculations in the protein system, we first analysed the radial distribution function water surrounding the protein molecule using our 15 μs and 20 μs molecular dynamics ensembles of the isolated protein and RNC and the GROMACS rdf functionality84 to identify the region of the first two hydration shells that show significantly reduced water dynamics. Using this analysis, we chose a distance cut-off of 3.5 Å between the protein and water centre of mass to define the hydration layer around the protein. With this criterion we then calculated the probability distribution and average number of water molecules in the hydration layer to assess the difference in solvation on and off the ribosome. Water molecules that remain within a defined distance range from the protein during the entire 20 ps production simulation were then selected to calculate the average molar entropy per molecule of water in different environments with DoSPT. The accessible volume for this subsystem was estimated by using the average volume occupied per water molecule in a pure water box under identical conditions multiplied by the number of molecules. To obtain the change in solvation entropy (difference between the RNC and isolated system, ΔSsolv,RNC-iso), we used

$$\Delta {S}_{{\rm{solv,RNC}}-{\rm{iso}}}={N}_{{\rm{diff}}}\Delta {S}_{{\rm{solv,water}}}$$

(S30)

where Ndiff is the average difference in the number of water molecules in the hydration layer (RNC-iso) and ΔSsolv,water is the entropy difference between water molecules in the hydration layer (0–3.5 Å from the protein) and water molecules in bulk solution (defined here as 36–46 Å from the protein).

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

RELATED ARTICLES

Most Popular

Recent Comments