Summary

Neutron Spin Echo Spectroscopy as a Unique Probe for Lipid Membrane Dynamics and Membrane-Protein Interactions

Published: May 27, 2021
doi:

Summary

This paper describes the protocols for sample preparation, data reduction, and data analysis in neutron spin echo (NSE) studies of lipid membranes. Judicious deuterium labeling of lipids enables access to different membrane dynamics on mesoscopic length and time scales, over which vital biological processes occur.

Abstract

Lipid bilayers form the main matrix of cell membranes and are the primary platform for nutrient exchange, protein-membrane interactions, and viral budding, among other vital cellular processes. For efficient biological activity, cell membranes should be rigid enough to maintain the integrity of the cell and its compartments yet fluid enough to allow membrane components, such as proteins and functional domains, to diffuse and interact. This delicate balance of elastic and fluid membrane properties, and their impact on biological function, necessitate a better understanding of collective membrane dynamics over mesoscopic length and time scales of key biological processes, e.g., membrane deformations and protein binding events. Among the techniques that can effectively probe this dynamic range is neutron spin echo (NSE) spectroscopy. Combined with deuterium labeling, NSE can be used to directly access bending and thickness fluctuations as well as mesoscopic dynamics of select membrane features. This paper provides a brief description of the NSE technique and outlines the procedures for performing NSE experiments on liposomal membranes, including details of sample preparation and deuteration schemes, along with instructions for data collection and reduction. The paper also introduces data analysis methods used to extract key membrane parameters, such as the bending rigidity modulus, area compressibility modulus, and in-plane viscosity. To illustrate the biological importance of NSE studies, select examples of membrane phenomena probed by NSE are discussed, namely, the effect of additives on membrane bending rigidity, the impact of domain formation on membrane fluctuations, and the dynamic signature of membrane-protein interactions.

Introduction

The understanding of cell membranes and their function has remarkably evolved over the last few decades. The former view of cell membranes as passive lipid bilayers that define cell boundaries and house membrane proteins1 has gradually transformed into a dynamic model in which lipid bilayers play an important role in regulating vital biological processes, including cellular signaling, molecular exchange, and protein function  to name a few2,3,4,5,6. This realization that cell membranes are highly dynamic, constantly undergoing remodeling and molecular redistribution, has urged scientific explorations beyond equilibrium structures of membranes7,8,9. Accordingly, multiple approaches have been developed to study the various dynamic modes in biological and bioinspired lipid membranes. To date, the majority of these studies have primarily focused on diffusive molecular motions10,11,12,13 and macroscopic shape fluctuations14,15,16, leaving a significant gap in understanding intermediate membrane dynamics, i.e., collective fluctuations of lipid assemblies consisting of few 10-100s of lipid molecules. These dynamics occur over length scales of few tens to few 100 Å and over time scales of sub-ns to few hundred ns (see Figure 1), referred to here as mesoscopic scales. It is indeed on these scales that key biological activity takes place at the membrane level17. This includes viral budding18, channel gating19, and membrane-protein interactions20. It is also important to point out that the energy landscape of membrane proteins21,22 shows that conformational changes in proteins  necessary for their regulatory role  happen over the ns time scales23 of collective membrane fluctuations, further emphasizing the importance of mesoscopic dynamics in the biological function of cell membranes and their bioinspired analogs20. This paper focuses on the two primary mesoscopic dynamic modes in lipid membranes, namely, bending fluctuations and thickness fluctuations.

The main challenge in directly probing these fluctuation modes is the difficulty in simultaneously accessing their spatial and temporal scales using standard spectroscopy methods. The other challenge is that direct contact techniques could impact the same fluctuations they are meant to measure16. This is further exacerbated by the compositional and structural complexity of biological membranes24,25, which results in non-homogeneous membrane features, including lipid domain formation26,27,28,29,30 and membrane asymmetry31,32,33 demanding selective probes to understand the dynamics of different membrane features. Fortunately, these challenges can be overcome with non-invasive neutron spectroscopy methods, such as neutron spin echo (NSE), which inherently access the required length and time scales, and further enable studies of selective membrane features without changing their physicochemical environment34. Indeed, over the last few years NSE spectroscopy has evolved into a unique and powerful probe of collective membrane dynamics35. Results from NSE studies on lipid membranes have produced new insights into mechanical36,37 and viscoelastic38,39 properties of lipid membranes and have shed new light on their potential role in biological function40,41.

The NSE spectroscopy technique is based on an interferometric instrument design, first proposed by Mezei42, using a series of spin-flippers and magnetic coils to control the precession of the neutron spin as neutrons traverse the instrument. The design rests on magnetic mirroring of the magnetic field elements with respect to the sample position (Figure 1A). This implies that in the absence of energy exchange between the neutron and the sample, the neutron performs the same number of spin precessions, in opposite directions, in the first and second half of the instrument (notice the π-flipper between the two precession coils). As a result, the final spin state of the neutron remains unchanged relative to the initial state – a phenomenon referred to as spin-echo (see transparent neutron in Figure 1A). However, when the neutron energetically interacts with the sample, the energy exchange modifies the number of spin precessions in the second half of the instrument, leading to a different final spin state (see Figure 1A). This is experimentally detected as a loss in polarization, as will be shown later in this paper. For more details on the NSE technique, the reader is referred to dedicated technical papers42,43,44,45.

Here, a simplified description is presented to provide a rough estimate of the length and time scales accessible with NSE. The length scales are determined by the range of achievable wavevector transfers, Q = 4π sin θ/λ, where 2θ is the scattering angle and λ is the neutron wavelength. One can see that Q is set by the wavelength range and the extent of rotation of the second arm of the spectrometer (see Figure 1A). A typical Q-range on NSE spectrometers is ~0.02-2 Å-1 46,47, and up to 0.01-4 Å-1 with recent upgrades48,49, corresponding to spatial scales of ~1-600 Å. On the other hand, the accessible time scale is calculated from the total precession angle (or phase) acquired by the neutron within the magnetic precession coils, and is found to be50: Equation 12. In this expression, t is the Fourier time defined as Equation 13, where Equation 50 is the neutron gyromagnetic ratio, Equation 51 is the coil length, and Equation 52 is the strength of the coil's magnetic field. It is worth pointing out that the Fourier time is a quantity that is strictly dependent on the instrument geometry, magnetic field strength, and neutron wavelength. For instance, using neutrons of wavelength Equation 70 = 8 Å and instrument settings of Equation 51 = 1.2 m and Equation 52 = 0.4 T, the Fourier time is calculated to be t ~ 50 ns. Experimentally, the Fourier time is tuned by changing the current in the precession coils (i.e., magnetic field strength) or using different neutron wavelengths, resulting in typical NSE time scales of ~ 1 ps to 100 ns. However, recent upgrades in NSE spectrometers have enabled access to longer Fourier times, up to ~400 ns on the J-NSE-Phoenix spectrometer at the Heinz Maier-Leibnitz Zentrum51 and the SNS-NSE spectrometer at Oak Ridge National Lab48, and up to ~1,000 ns at the IN15 NSE spectrometer at the Institut Laue-Langevin (ILL)49.

Besides direct access to the length and time scale of membrane dynamics, NSE has the inherent capabilities of neutron isotope sensitivity52. Specifically, the ability of neutrons to interact differently with the isotopes of hydrogen, the most abundant element in biological systems, results in a different neutron scattering length density,34 or NSLD (the equivalent of the optical index of refraction50), when protium is substituted by deuterium. This enables an approach known as contrast variation, which is commonly used to highlight specific membrane features or conceal others  the latter scenario is referred to as contrast matching. A frequent application of contrast variation/matching is the substitution of water (NSLD = -0.56 × 10-6 Å-2) by heavy water or D2O (NSLD = 6.4 × 10-6 Å-2) to amplify the neutron signal from protiated lipid membranes (NSLD ~ 0 × 10-6 Å-2). This approach is highly effective in studies of membrane structure because the penetration of D2O into the headgroup region of the membrane allows accurate determination of the membrane thicknesses (see Figure 2A, left panel) and of the location of different lipid subgroups when more sophisticated models are applied53,54. This paper highlights some examples on the use of contrast variation for studies of collective dynamics in biomimetic membranes and select membrane features.

Here, the effectiveness of NSE in providing unique insights into dynamical and functional membrane properties is illustrated through tangible examples of NSE studies on model and biologically relevant lipid membrane systems with emphasis on mesoscale dynamics in free-standing membranes, in the form of liposomal suspensions. For NSE measurements of in-plane membrane dynamics, the reader is referred to dedicated publications on grazing-incidence neutron spin-echo spectroscopy (GINSES)55,56 and other studies of aligned multilamellar membrane stacks57,58,59,60.

For simplicity, this paper highlights three different schemes of membrane deuteration illustrated on a well-studied domain-forming, or phase separating, lipid bilayer system of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and 1,2-distearoyl-sn-glycero-3-phosphocholine (DSPC) mixtures61,62. The two lipids are characterized by a mismatch in their hydrocarbon chain length (14 carbons/tail in DMPC vs 18 carbons/tail in DSPC) and their gel-fluid transition temperature (Tm, DMPC = 23 °C vs Tm, DSPC = 55 °C). This results in lateral phase-separation in DMPC:DSPC membranes at temperatures between the upper and lower transition temperatures of the mixture63. The deuteration schemes considered here are chosen to demonstrate the different dynamic modes accessible in NSE measurements on liposomal membranes, namely, bending fluctuations, thickness fluctuations, and selective bending/thickness fluctuations of lateral domains. All lipid compositions are reported for DMPC:DSPC bilayers prepared at a mole fraction of 70:30, using commercially available protiated and perdeuterated variants of DMPC and DSPC. All sample preparation steps are based on 4 mL of liposomal suspension, in D2O, with a lipid concentration of 50 mg/mL, for a total lipid mass of Mtot = 200 mg per sample.

Protocol

1. Deuteration scheme required for the experiment

  1. For bending fluctuation measurements, make fully protiated liposomes in D2O (D 99.9%) or D2O-buffer (e.g., phosphate buffer prepared with D2O instead of H2O). Use fully protiated DMPC (C36H72NO8P) and DSPC (C44H88NO8P) with Equation 1 133.4 mg, where XDMPC and XDSPC are the mole fractions of DMPC and DSPC, here set to 0.7 and 0.3, respectively, and MwDMPC and MwDSPC are the molar weights given by 677.9 g/mol and 790.1 g/mol, respectively. Similarly, mDSPC = 66.6 mg. This deuteration scheme increases the scattering contrast between the membrane (NSLD ~ 0 × 10-6 Å-2) and the deuterated buffer (NSLD ~ 6.4 × 10-6 Å-2) and amplifies the signal from membrane undulations (see Figure 2A left panel).
  2. To measure the bending dynamics of select lateral membrane features, e.g., matrix dynamics in phase-separating DMPC:DSPC membranes, use protiated DMPC (C36H72NO8P) and deuterated, DSPC-d83 (C44H5NO8PD83, Mw 873.7 g/mol), such that mDMPC = 128.8 mg and mDSPC-d83 = 71.2 mg. This deuteration scheme minimizes the scattering from the undesired DSPC-rich domains, enabling selective measurements of bending fluctuations from the DMPC-rich matrix (see Figure 2B middle).
    NOTE: To find the optimal lipid deuteration required for a specific contrast matching scheme, utilize available web-based scattering length density (SLD) calculators, such as the one developed by the NIST Center for Neutron Research64. These web-based interfaces are equipped with user-friendly tools for easy calculation of the SLD of lipids with various degrees of deuteration, as well as that of lipid mixtures.
  3. For NSE measurements of average membrane thickness fluctuations (with no lateral contrast), use tail-deuterated variants of the constituent lipids, i.e., DMPC-d54 (C36H18NO8PD54, 732.3 g/mol) and DSPC-d70 (C44H18NO8PD70, 860.1 g/mol)35,38, such that mDMPC-d54 = 133.0 mg and mDSPC-d70 = 67.0 mg. This contrast scheme (Figure 2A, right panel) amplifies the scattering signal from the lipid headgroups (NSLD ~ 4.5 × 10-6 Å-2) by contrast-matching the tail-group (NSLD ~ 6.4 × 10-6 Å-2) to the deuterated buffer enabling the detection of fluctuations in membrane thickness.
  4. For thickness fluctuation studies of select membrane compartments, e.g., DMPC-rich matrix, use the same strategy described in step 1.2 by substituting protiated DMPC lipids with their tail-deuterated analogs, i.e., DMPC-d54, such that the DSPC-rich domains are contrast-matched to the deuterated buffer and the primary scattering signal is from the headgroup region of the tail-deuterated DMPC-rich matrix.

2. Preparation of lipid suspension for extrusion

  1. Calculate the mass of each constituent in the sample, depending on sample composition. As a rule of thumb, for samples with multiple molecular components, the mass of a component is given by its molar mass, Mwi, weighted by its mole fraction, Xi, and normalized over all components such that: Equation 2  where Mtot is the total mass, set here to 200 mg. See the example above for DMPC-DSPC lipid bilayers with different deuteration schemes.
  2. Using a digital semi-microbalance, weigh the calculated masses of lipids (and other sample constituents, e.g., proteins, nanoparticles, etc.) and add them to a vial or round-bottom flask – remember to weigh the vial or flask beforehand. Add 1 mL of solvent to dissolve the weighed components by manually mixing inside a hood. For pure lipid samples, use chloroform or ethanol. For samples with additional, non-lipid components (e.g., nanoparticles), choose a common solvent that disperses all components.
    1. For small lipid amounts (<10 mg), prepare a stock solution and pipette the required volume into the mixture.
      NOTE: Do not add excessive amounts of solvent as it will significantly slow down the solvent drying step described below.
  3. Dry the lipid solution, inside a hood, by gently streaming an inert gas (e.g., nitrogen, argon) in the vial while slowly rotating the vial at an angle. Keep the vials in tilted position to create a thin film of dried lipid on the vial walls, which will allow for even drying. Intermittently place the vial in a water bath at 35 °C to circumvent evaporation-mediated cooling, which will slow down the solvent evaporation.
  4. Place the vials overnight in a vacuum oven at ~35 °C to fully remove the residual solvent. For unsaturated lipids, purge vacuum with an inert gas to minimize oxidation.
  5. To ensure full solvent removal, weigh the vial after lipid drying and confirm that there is no excess mass beyond the measured amounts of materials. Do this by subtracting the mass of the vial from the measured mass after drying. If there is excess mass, dry the sample under vacuum for another 6 h. Repeat this process as needed.
  6. Hydrate the lipid film with 4 mL of D2O or D2O-buffer to obtain a lipid concentration of 50 mg/mL. For lipids with high transition temperatures, such as DMPC-DSPC mixtures, heat the buffer to above the transition temperature (60 °C) to ensure even mixing.
    NOTE: Since NSE experiments require relatively large sample volumes (~4 mL), consider hydrating the sample using half of the required buffer, i.e., 2 mL, to minimize the number of extrusions per sample (see section 3). In this case, add the remaining half of the buffer post extrusion. Notice that the capacity of syringes used in extrusion is limited to 1 mL. Thus, hydrating with 4 mL of buffer would require four sets of extrusion.
  7. Vortex-mix the hydrated lipid solution until the lipid film is fully dissolved and is no longer visible on the walls of the vial. At this stage, the hydrated lipids form multilamellar vesicles and micron sized multilamellar stacks and the suspension appears milky white.
  8. To facilitate the breaking of the lipid stacks and to reduce multilamellarity, perform five freeze/thaw cycles by placing the vial of hydrated lipid solution in a lab grade freezer (preferably -80 °C freezer) until fully frozen and then transferring the vial to a 35 °C water bath until the lipid solution is fully thawed. Vortex the thawed solution until homogenous. Repeat four more times.
    ​NOTE: Alternatively, a dry ice bath can be prepared for rapid freezing by combining acetone and dry ice.

3. Extrusion of the hydrated lipid solution

  1. Assemble the extruder setup using a polycarbonate membrane between two membrane supports and adding two paper filters on each side to provide additional support. Use a polycarbonate membrane with a pore size that matches the target liposomal size (common pore sizes for NSE experiments are 50 nm and 100 nm – typically, 100 nm-diameter liposomes allow for less constrained membrane fluctuations, but smaller 50 nm liposomes could be used for curvature studies). Ensure that the polycarbonate membrane is fully stretched before completing the assembly and tightening the external extruder casing.
  2. Hydrate the polycarbonate membrane by passing ~0.3 mL of D2O or D2O-buffer a few times through the membrane assembly using airtight glass syringes. Utilize the same buffer used in sample preparation. Leave it for at least 10 min, then completely suck the buffer out before introducing the sample.
  3. Fill a 1 mL gas-tight syringe with the prepared lipid solution and insert into one end of the extruder apparatus. Then, insert an empty syringe into the opposite end. Once the syringes are connected to the extruder assembly, place it into the extruder block.
  4. If elevated temperatures are needed for extrusion, as in the case of saturated lipids with high transition temperatures (e.g., DSPC, Tm = 55 °C), preheat the extruder heating block above the lipid transition temperature (e.g., 60 °C), by placing the heating block on a hot plate or by using a circulation bath as shown in Figure 3A.
    ​NOTE: This step is crucial to ensure homogeneous mixing of lipids and to avoid exerting extreme pressure during extrusion, which could rupture the polycarbonate membrane. For lipid samples with low transition temperatures (<25 °C), perform the extrusion at room temperature.
  5. To extrude the lipid solution, attach the extruder set to a programmable syringe pump with an aluminum/steel frame as shown in Figure 3A. For temperature-controlled extrusions, add a custom-built extruder base with a fluid channel and attach to a circulating water bath.
  6. Program the syringe pump to perform 15-20 extrusion cycles following the manufacturer's manual. When extruded, the color of the lipid solution changes from milky white to transparent opal blue (Figure 3B,C), indicating a final liposomal size that is smaller than the wavelength of visible light, as expected. For the type of syringe pump shown in Figure 3A, follow the steps below.
    1. Start by adjusting the pump settings. Hold down the Rate button and enter the extrusion rate (50.99 mL/h), then press the Diameter button and enter the syringe diameter (4.606 mm). Use the up arrows under each digit on the screen to change that digit value.
    2. Place the extruder set with the sample syringe to the right (see Figure 3A). Press the Withdraw button until the withdraw light turns on. Press Start and wait for the sample to dispense into the left (empty) syringe.
    3. Hit the Stop button just before the sample (right) syringe is fully empty. Record the dispensed volume and use it to program the extrusion cycle. Hold down the Rate button until phase 1 (PH:01) appears on the screen. Press the Volume button to enter the dispensed volume recorded earlier. In this phase, make sure that the Withdraw light is off – this dispenses the sample in the right direction.
    4. Press the Rate button again and use the rightmost up arrow to access phase 2 (PH:02). Press Volume to enter the same value of the dispensed volume recorded earlier. In this phase, press the Withdraw button until the Withdraw light is on – this dispenses the sample to the left.
    5. To repeat this cycle, press the Rate button again and use the rightmost up arrow to access phase 3 (PH:03). Press the Volume button until LP:SE appears on the screen and set it to 20. This is the number of loops or repeats that the pump will perform. Finally, press the Rate button, access phase 4 (PH:04), and hit the Volume button to get to the Stop function. The pump is now set up for automated extrusion.
    6. Press Start to start the extrusion cycle.
  7. Empty the syringe containing the extruded lipid suspension in a clean vial and prepare for storage or measurements. For lipid samples with high melting temperature, store the sample above the fluid phase transition until measured. Otherwise, keep samples at room temperature.
  8. Do not freeze extruded samples as freezing will cause the vesicles to burst (the suspension will turn milky white again).

4. NSE measurements for the sample(s) and reduction of the collected data

  1. Prior to the NSE experiment, characterize the extruded liposomal sample from step 3.7 using available methods to ensure adequate sample quality. A list of potential charcaterization methods that can be used to assess the quality of liposomal suspensions for NSE experiments, e.g., size distribution, multilamellarity, lateral membrane structure, is included in the discussion section.
  2. Determine the Q-range and corresponding instrument settings required for the experiment. For bending rigidity measurements of lipid bilayers, use a Q-range of ~(0.04 – 0.2) Å-1. For studies of membrane thickness fluctuations, use a Q-range of ~(0.04 – 0.2) Å-1 corresponding to the membrane thickness35,66,67.
    NOTE: Discuss the experimental setup with the instrument scientist before the start of the experiment. As mentioned earlier, SANS characterization of the sample is necessary, especially if prior information of the scattering signal is not available, as in selectively deuterated membranes. Alternatively, run static (also known as diffraction) measurements over a limited Q-range on the NSE instrument, with the caveat that such measurements take much longer compared to SANS.
  3. Using a syringe or a transfer pipette, load the extruded liposomal suspension(s) in the designated sample cells available at NSE beamlines. Note that standard NSE sample cells come in thicknesses of 1, 2, 3, and 4 mm. Choose the cell thickness in such a way to optimize the scattering signal while keeping the incoherent background signal to a reasonable intensity.
    NOTE: As a rule of thumb, use sample cells with 1 or 2 mm pathlength for protiated liposomes in deuterated buffer – thicker cells could result in multiple scattering effects that are difficult to correct for. For liposomes with higher levels of deuterations (e.g., tail contrast-matched liposomes or asymemetric liposomes with single protiated leaflets), consider using a thicker sample cell (e.g., 3 or 4 mm pathlength) to enhance the counting statistics if the sample is available in larger quantities – sometimes this can be cost prohibitive.
  4. Prepare an identical sample cell for the buffer. Use the same buffer as in the liposomal suspension. Measurements on the buffer are necessary for intensity normalization and background (BKG) corrections.
  5. Place the sample cell(s) in the sample holder of the NSE spectrometer, program the measurement runs, and collect echo data. Consult with the instrument scientist about programming the measurements if a first-time NSE user.
  6. Perform two additional sets of measurements needed for the data reduction: Resolution (R) and transmission (T) measurements.
    1. Perform Resolution (R) measurement on an elastic scattering reference (e.g., carbon)  to be run under the same settings; i.e. same wavevector and Fourier times as the sample and buffer measurements.
    2. Perform transmission (T) measurements on the sample and buffer to calculate the intensity of the transmitted neutron beam (see step 4.9. below). The transmission is calculated as the the ratio of neutron counts from the sample or buffer divided by the neutron counts for an open beam (i.e., with an empty sample position).
  7. Use the dedicated data reduction software for the NSE spectrometer on which the measurements are performed to reduce the collected data.
    NOTE: Different spectrometers might utilize different software or user interfaces. Below is an example of NSE data reduction using the Data Analysis and Visualization Environment (DAVE)68 software specifically written for the NSE spectrometer at the NIST Center for Neutron Research.
    1. Open the DAVE software and select Reduce NSE Data from the data reduction menu. Several pop-up windows will appear.
    2. Upload the data files over different Q-values using the Open .echo Files from the file menu. These files correspond to the raw data files with the spin echo signals and have the extension .echo in the file name. Once the file upload is complete, the files will show under the available data sets.
    3. Right-click on the selected file and label it according to the measurement it corresponds to; i.e., Sample, Cell (for empty cell or buffer), or Resolution.
    4. Group the detector pixles in 2 x 2 to improve the signal-to-noise ratio using the Data Set tab. Apply the same binning to all files; i.e., Resolution, Cell, and Sample.
    5. Inspect the data over all pixel groups and mask those with poor signals (see Figure 4B) by pressing the m key on the keyboard. Press Enter to access a pop up window to apply the same mask to all Fourier times or subsequent Fourier times. This can also be applied to individual pixels at any point during data reduction. Masked pixels will turn green.
    6. Ensure that the collected data is in the form of an echo signal, i.e., cosine function in terms of the phase current, over each detector pixel (see Figure 4A).
      NOTE: The phase current is proportional to the precession angle of the neutron spin; hence, it is common to represent the phase current as a phase angle as shown in Figure 4A. For measurements on pulsed sources, additional time of flight calculations are applied to the data to obtain the echo signals as a function of incident neutron wavelength within a neutron pulse.
    7. Start by fitting the resolution file. Select a resolution file from the uploaded file list and right-click on the file. Select Fit Operations: Fit Echoes (Resolution) from the pop-up menu.
    8. Ensure that the fits of the echo signals yield a number of fitting parameters, including the parameter, A, required in step 4.8. The fits are automatically performed using following expression.
      Equation 8
      Here, ζ is the period of the echo signal (i.e., cosine function in Figure 4A), σ is the width of the Gaussian envelope determined by the mean wavelength and wavelength spread of the incident neutron beam, Φc is the phase current, and Φ0 is the echo point which depends on the field path experienced by neutrons50. Physical information about the sample is encoded in the amplitude, A, of the cosine function in equation (1).
      ​NOTE: The width of the Gaussian envelope is based on values predetermined by the instrument scientist and should not be changed. The other paramters are variables that are fitted to the specific echo signal over each pixel.
    9. Inspect the fit results by clicking on each pixel to show the resulting fitting parameters, the quality of the fit, and the mean square deviation of the fit. To inspect the error associated with each fitting parameter over the entire detector, select Image Options and then select the fitting parameter of interest. This will generate a map with the value of the fitting paramter over each pixel. Right-click on the detector image. A pop up window will appear showing an error bar map of the selected fitting parameter.
    10. If the fit over a specific pixel is unsatisfactory (e.g., fit parameeters with large error bars), refit the signal over that specific pixel. Select that pixel, press the Fitting tab, and then press Fit Pixel. Input new starting parameters for the phase (Φ0) and period (ζ) in the Fitting tab to obtain a more satisfactory fit.
      NOTE: It is useful to plot the fitted phase as a function of Fourier Time. To do so, go to the main plot window and select Fit Phase v. Fourier Time. This plot should be smooth and continuous. Inspect discontinuities in this plot and refit the pixels that they correspond to.
  8. Reduce the Sample or Cell file by selecting the corresponding file from the uploaded and labeled file list.
    1. Inspect all pixels and mask the ones with bad statistics as described in step 4.7.5.
    2. Right-click on the file and select Fit Operations: Import Phases (Sample, Cell). This imports the phases and the applied mask from the Resolution file.
    3. Fit the echo signals using the same procedure described before for the Resolution file (steps 4.7.8-4.7.10). In fitting the Sample and Cell files, do not change the values of the period and echo phase point imported from the Resolution fits. These parameters depend on the instrumental settings and should not vary with samples.
    4. Before proceeding to data reduction, input the beam center for all data files. Select the data file, go to the General tab and enter X and Y beam center values. These values are recorded during the experiment.
    5. Once the fits to the Sample, Cell, and Resolution files are complete, calculate the normalized intermediate scattering function to be used later in data analysis and interpretation. To do that, right click on the Sample file to be reduced from the list of fitted files, and select Calculate I(Q) from the pop up menu. A window will appear with entry choices for the Resolution and Cell (i.e., buffer) files, and the number of Q-arcs (see step 4.9). After entering all the required informaiton, press the OK button. The results will appear in a new window.
      NOTE: The data reduction is performed according to the following equation to obtain the normalized intermediate scattering function69.
      Equation 4
      where t is Fourier time, Nup and Ndown are the neutron counts in the non-spin-flip and spin-flip configurations (measured with the π/2-flippers off and the π-flipper off and on, respectively), and the superscripts, BKG and R, correspond to the background and resolution measurements, respectively, as defined in steps 4.4 and 4.6. Note that the beam polarization Equation 3, thus changes in the spin state due to energy exchange between the neutron and the sample is detected as a drop in the polarization (from unity).
  9. Finally, group the detector pixels into Q-arcs as shown in Figure 4B to obtain the Q-dependence of the normalized intermediate scattering function, S(Q,t) / S(Q,0). This is technically referred to as data binning and should be done judiciously, i.e., taking into account the counting statistics from the sample and the expected standard deviation of the data over the grouped pixels.
  10. For strongly scattering samples, divide the detector into more Q-arcs while maintaining reasonable error bars on the resultant intermediate scattering function, S(Q,t) / S(Q,0). This yields more Q data points and is important for the data analysis procedure described below. Be aware that for weakly scattering samples, excessive binning results in poor decay signals, i.e., large error bars on S(Q,t) / S(Q,0), which could result in large uncertainties.

5. Data analysis and interpretation

  1. Fit the normalized intermediate scattering functions, S(Q,t) / S(Q,0), obtained from the data reduction above to a stretched exponential function with a stretching exponent of 2/370.
    Equation 9
    NOTE: An example of these fits is provided in Figure 5B. Fits of S(Q,t) / S(Q,0) to equation (3) yield the Q-dependent relaxation rates Γ(Q).
  2. Plot Γ(Q) as a function of Q and fit to a suitable model to extract relevant membrane parameters.

Representative Results

NSE studies accessing bending fluctuations are typically performed over a Q-range of ~ (0.04 – 0.2) Å-1. This Q-range corresponds to intermediate length scales between the membrane thickness and the liposomal radius, where bending dynamics dominate. Measurement over an extended Q-range can give access to additional dynamic modes, including liposomal diffusion and intramembrane dynamics. For more details on the cross-over in membrane dynamics accessed by NSE, check these relevant publications25,71. It is important to emphasize that NSE signals are proportional to: Equation 5, where Icoh and Iinc are, respectively, the coherent and incoherent scattering intensity from the sample. Therefore, it is advisable to prepare NSE liposomal samples in deuterated buffers (i.e., buffers prepared with D2O instead of H2O) to minimize the incoherent scattering signal, mainly contributed by the hydrogen content of the sample. However, in some cases intermediate deuteration schemes (i.e., using mixtures of D2O and H2O) might be necessary to obtain optimal contrast conditions. Typically, NSE measurements of membrane bending fluctuations are performed on fully protiated liposomes in deuterated buffer, referred to as fully contrasted liposomes in Figure 5. This deuteration scheme results in a large NSLD difference between the membrane core (~0 × 10-6 Å-2) and its deuterated fluid environment (~6.4 × 10-6 Å-2), which significantly enhances the scattering signal from the liposomal membranes and improves the measurement statistics of bending dynamics. This contrast scheme (Figure 2A left panel) is frequently utilized in studies of bending rigidity of lipid membranes with single38,72 and multiple39,66 lipid components and in studies of membrane softening/stiffening by biological inclusions (e.g., cholesterol, drug molecules, peptides/proteins)36,37,73,74,75, and synthetic additives (e.g., nanoparticles)76,77.

Measurements of bending fluctuations result in relaxation rates that follow a Q3 dependence, as predicted by Zilman and Granek for thermally undulating elastic thin sheets70. A refined form of this Q-dependence is obtained from theoretical corrections by Watson and Brown78, which take into account the effects of intermonolayer friction proposed by Seifert and Langer79. By additionally defining the neutral plane to be at the interface between the hydrophilic headgroups and the hydrophobic tails of the membrane, the bending relaxation rates can then be fitted to the following expression38.

Equation 6

where ηbuff is the buffer viscosity, kBT is the thermal energy, κ and is the bending rigidity of the measured membrane (or of the contrasted portion of the membrane in selectively deuterated systems). This type of measurement enables direct calculation of membrane elastic properties in the form of the bending rigidity modulus. Note that κ is extracted from the slope of the linear fit of Γ vs. Q3, as shown in Figure 5C.

On the other hand, NSE measurements of membrane thickness fluctuations show deviations from the Q3-dependence in Γ(Q) around Q values that correspond to the membrane thickness (see Figure 2 in ref.66). To isolate the thickness fluctuation signal, one can divide Γ(Q) by Q3, as shown in Figure 5D. The resulting data show that the excess dynamics due to thickness fluctuations follow a Lorentzian function in Q, as recently corroborated in coarse-grained molecular dynamics (MD) simulations67. To fit the observed excess dynamics, Nagao et al.38 developed an expression based on the theoretical framework of membrane fluctuations by Bingham et al.80 as follows.

Equation 7

In this expression, Q0 is the peak Q-value corresponding to the membrane thickness (which can be independently obtained from SANS measurements), μ is the in-plane membrane viscosity, AL is the area per lipid (measured with SANS/SAXS), and KA is the area compressibility modulus. Assuming that KA can be calculated from κ using the polymer brush model, this expression reduces to one fit parameter, namely, the membrane viscosity μ, presenting a new approach to measure membrane viscosity without the need for fluorescence labeling or particle tethering/tracking13. The premise is that according to deformation models of elastic thin sheets81, κ and KA are interdependent such that: Equation 10, where tm is the mechanical (or deformable) membrane thickness and β is a constant that describes interleaflet coupling. The assumption is that β = 12 for fully coupled leaflets, β = 48 for completely uncoupled leaflets, and β = 24 for intermediately coupled leaflets. The latter is referred to as the polymer brush model81 and has been shown to apply in single-component and binary fluid lipid membranes39. However, this needs to be approached with caution. For instance, recent simulations by Doktorova et al.82 showed that for the polymer brush model to hold in unsaturated lipid membranes containing cholesterol, a modified expression of the mechanical membrane thickness must be used. Ideally, if an independent measurement of KA is possible, e.g., using micropipette aspiration83, then combining KA results with NSE bending rigidity measurements would present a unique opportunity to investigate interleaflet coupling in model and biological membranes – a long standing question in membrane biophysics and structural biology. Once the values of KA are validated, they can be used in equation 5 to obtain the mesoscopic membrane viscosity.

Figure 1
Figure 1: NSE instrument design and synergistic overlap with length/times scales of mesoscopic membrane dynamics. (A) Schematic of the different magnetic elements of an NSE instrument, used to manipulate the spin of neutrons traversing the instrument from left to right. The highlighted neutron indicates change in spin orientation (or polarization loss) due to energy exchange between the neutron and the sample, whereas the transparent neutron represents spin-echo, i.e., no change in the neutron spin due to zero energy exchange. The grey arrow indicates the possibility of rotating the second arm of the spectrometer to access larger scattering angles. (B) Pictorial representation of hierarchical dynamics in lipid membranes, showing various dynamic modes that span multiple length and time scales. The shaded area represents the length and time scales accessed by NSE, which overlap with the mesoscales of collective membrane fluctuations, namely, bending and thickness fluctuations. Please click here to view a larger version of this figure.

Figure 2
Figure 2: Examples of possible deuteration schemes in NSE experiments on lipid membranes. (A) Left: Fully contrasted membranes, e.g., protiated membranes in deuterated buffer, showing the NSLD profile along the normal to the membrane surface. The difference in the NSLD between the tail region (~0 × 10-2 Å-2) and headgroup region (~4.5 × 10-6 Å-2) of the membrane is due to the headgroup hydration with deuterated buffer. Right: Tail-contrast matched membranes such that the hydrocarbon tail region of the membrane has the same NSLD as the buffer, as shown in the corresponding NSLD profile along the membrane normal. (B) Domain-forming membranes with two neutron contrast schemes where the domains (center) or the matrix (left) are contrast-matched to the buffer, enabling selective studies of matrix or domain dynamics, respectively. This figure has been modified from Nickels et al., JACS 201541. (C) Asymmetric membranes prepared by cyclodextrin exchange between protiated and deuterated lipid vesicles, resulting in the deuteration of one membrane leaflet while keeping the other leaflet protiated. This allows studies of the bending dynamics of the protiated leaflet and provides insights into the mechanical coupling between opposing leaflets in asymmetric membranes. This figure has been modified from Rickeard et al., Nanoscale 202040. Please click here to view a larger version of this figure.

Figure 3
Figure 3: Illustration of the setup for automated extrusion of liposomes. (A) Custom-built automated extruder using a syringe pump, a mini extruder set, and an aluminum/steel frame to enable cyclic extrusions. (B) and (C) show the difference in visual appearance of lipid suspensions before (milky white) and after (transparent opal blue) extrusion. This is due to the initial formation of micron-sized lipid stacks or giant vesicles which are on the order of, or larger than, the wavelength of visible light. After extrusion, the suspension will comprise nanoscopic vesicles (~100 nm), which are smaller than the wavelength of visible light, yielding a transparent suspension. Please click here to view a larger version of this figure.

Figure 4
Figure 4: Representative data from NSE experiments on liposomal suspensions. (A) Example of an echo-signal over a single detector pixel (marked pixel in panel B), showing the fits of the echo signal using equation (1), with an illustration of the different parameters required in the echo fit. Note that the echo signal is plotted as a function of the phase angle instead of the phase current as discussed in step 4.7 of the protocol. (B) NSE detector image showing the variation in neutron counts per pixel. The image also shows eliminated detector pixels (green) due to poor echo signals. The binning of the detector pixels in Q-arcs (also known as Debye-Scherrer rings) yields the Q-dependence of the intermediate scattering function, needed for analyzing and interpreting NSE data. This figure was modified from Ashkar, J. Appl. Phys. 202050. Please click here to view a larger version of this figure.

Figure 5
Figure 5: Representative results from NSE experiments on liposomal suspensions with different deuteration schemes. (A) Scattering geometry of a neutron interacting with a liposome, showing the scattering angle, 2θ, and the wavevector transfer, Equation 11. (B) Intermediate scattering functions, S(Q,t) / S(Q,0), exhibit decays as a function of the Fourier time. Fit of the measured decays to a stretched exponential function given by equation 3 yields the relaxation rates, Γ(Q). (C) For fully protiated liposomes in deuterated buffer, Γ(Q) follows a Q3 dependence, typical of bending dynamics. The linear fit of the obtained data to a Zilman-Granek model yields the bending rigidity modulus of the membrane. (D) For tail deuterated liposomes, excess dynamics are observed in addition to bending fluctuations and are most pronounced at Q-values that correspond to the membrane thickness. Fitting the excess dynamics to a Lorentzian function (equation 5) allows extraction of the membrane viscosity. Data sets were collected on the NSE spectrometer at NIST. Please click here to view a larger version of this figure.

Discussion

NSE is a powerful and unique technique in measuring mesoscopic dynamics of lipid membranes under various conditions. The effective utilization of NSE depends on sample quality, neutron contrast, and the range of accessible dynamics that can be probed for a given sample. Thus, several critical steps are required for performing successful NSE experiments and collecting high-quality data. A key step in ensuring the effective use of neutron beam time during an NSE experiment is to characterize the liposomal suspensions with lab-based methods prior to the NSE experiment. For exmaple, the size distribution (or diffusion constant) of extruded liposomes can be determined by dynamic light scattering (DLS), readily available in individual labs or in shared facilities84. Cryo-electron microscopy is another charcaterization method recently validated on liposomal samples, where high-resolution microscopy images on cryomicrotomed sections of liposomal suspension can be effectively used to examine liposomal unilamellarity65, domain formation85,86, or the incorporation of additives such as nanoparticles76 and proteins87. Alternatively, small-angle x-ray scattering (SAXS) can be used to characterize the membrane structure88, assess liposomal multilamellarity65, or evaluate the effects of additives on membrane structural properties89. Besides these lab-based techniques, it is highly advisable that NSE measurements on liposomal samples are paired with structural studies using small-angle neutron scattering (SANS)54,90. SANS is an excellent complement to NSE, not only for acquiring structural membrane information but also for examining the intensity of the neutron scattering signal from the sample, confirming the contrast scheme, and making an informed choice about the Q-range over which the NSE measurements should be performed. Therefore, it is recommended that NSE users request SANS beamtime when applying for NSE experiments.

However, NSE suffers from sample limitations in studies of biological membranes. One of the major limiting factors of such experiments is the standard amount of sample required for NSE measurements (2-4 mL) and the high sample concentrations amounting to 100-200 mg of membrane material (lipids and proteins) to obtain high quality data. In many cases, the production of such amounts of biological material is not feasible or is cost prohibitive. In such scenarios, it is possible to reduce the concentration to 20-25 mg/mL, but this would require at least a 4-fold increase in the acquisition time in order to obtain comparable statistics to samples with 50 mg/mL concentrations. These stringent requirements on sample volume and concentration could be alleviated with the next generation of NSE spectrometers on higher-flux neutron sources, such as the second target station at Oak Ridge National Lab and the European Spallation Source. Another critical limitation in performing NSE experiments on lipid membranes requiring selective deuteration schemes is the lack of commercial availabilty of some deuterated variants of lipid molecules or their exorbitant prices, if available. In some cases, these limitations can be circumvented by requesting the synthesis of deuterated lipids (or cholesterol, proteins) through user deuteration facilities, such as the bio-deuteration lab at Oak Ridge National Lab, the national deuteration facility at ANSTO, or the deuteration facility at the ISIS Neutron and Muon Source. Access to these facilities and their synthesis capabilities is available through submitted user proposals which are peer-reviewed based on the scientific merit of the proposed material synthesis and its intended use in isotope-sensitive studies.

Despite these limitations, the application of NSE spectroscopy in studies of membrane mechanics has led to the determination of the bending rigidity moduli of membranes with various degrees of complexity, from single-component lipid membranes35,38 to multicomponent biomimetic membranes41,66,91, all of which have advanced our understanding of the dynamic nature of lipid membranes. For instance, NSE bending rigidity measurements of lipid membranes with different molecular units, e.g., lipids of different acyl chain lengths and chain saturation38,72,92, have provided essential information about the role of molecular chemistry in membrane mechanics. When paired with structural information, such as membrane thickness or molecular packing93, these measurements start to provide new perspectives on the interdependence between membrane structure and dynamics and how they influence membrane function. The mesoscopic scales of NSE uniquely position it for such fundamental investigations of structure-property relations, most relevant on the length scale of molecular assemblies. This topic was recently explored in two NSE studies on cholesterol-rich lipid membranes36 and in binary lipid membranes with hydrophobic mismatch between the two lipid components39. Both studies found strong evidence that membrane mechanics scale with the area per lipid, corroborating the conclusions from a recent all-atom MD simulations by Doktorova et al.82. These findings emphasize the self-assembled nature of lipid membranes and provide a unifying picture of molecular packing as a key parameter in defining membrane dynamical and functional properties.

Other applications of NSE involve studies of the mechanical response of membranes to small additives, including biological molecules such as cholesterol36,37, trehalose92, and melittin73,94, or inorganic additives such as nanoparticles for drug delivery applications76. NSE has also been used to understand how membrane mechanics respond to changes in their environment, including temperature92, pH74, and the presence of crowding macromolecules96. Such studies are contributing to a better understanding of the factors that influence the softening or stiffening of lipid membranes, under biological conditions related to health and disease, and in controlled settings for therapeutic applications. Notably, NSE measurements have also been used to probe the effect of antimicrobial peptides on membrane dynamics73,94,95. Further examples of NSE applications on biomembranes include studies of the dynamics of flattened membrane structures, called thylakoids, which house the photosynthetic machinery in cyanobacterial cells97,98.

One can also utilize selective lipid deuteration in NSE studies to investigate the dynamics of specific membrane features that are relevant to biological function. For example, Nickels et al. used selective lipid deuteration in domain-forming lipid membranes to generate lateral contrast within the membrane, as previously illustrated by Heberle et al.28. This deuteration scheme enabled independent measurements of the bending rigidity of the lipid domains and the host lipid matrix41 (see Figure 2B). The findings confirmed that the two membrane compartments have distinct bending rigidity moduli, which could be a driving mechanism for domain formation in cellular membranes. In a more recent study, Rickeard et al. used cyclodextrin exchange between protiated and deuterated liposomes to obtain asymmetric liposomes with isotopically labeled leaflets40 (Figure 2C). Their end liposomes had a protiated leaflet and a deuterated leaflet that is contrast matched to the buffer, enabling studies of individual leaflet dynamics and providing a first direct experimental account of the effect of asymmetry and leaflet coupling on membrane bending fluctuations.

Selective membrane deuteration has also been used in NSE studies of membrane thickness fluctuations, a long predicted dynamic mode in lipid membranes99 that was only recently observed with the advent of NSE spectroscopy35,100. These measurements utilize tail-deuterated membranes to amplify the signal from the membrane headgroup regions and resolve the thickness fluctuation signal. This type of NSE experiments is relatively recent, but it has been effectively used to understand the interdependence of membrane elastic and viscous properties38, to explore the scaling of bending rigidity and viscosity with molecular packing in mixed lipid membranes39, and to probe the local effects of cholesterol on membrane viscosity36. Another area of biological significance in which this dynamic mode could have far-reaching implications is mesoscopic membrane-protein interactions95. It is known that the function of membrane proteins is tightly linked to hydrophobic matching between the protein and the host membrane. Thus, variations in membrane thickness, due to thickness fluctuations, could act as a regulatory mechanism for the function of membrane proteins. NSE is extremely well suited for such studies as it can directly probe the effects of protein binding and insertion on membrane thickness fluctuations. Recent NSE measurements from our group (unpublished) suggest that transmembrane protein insertion could significantly suppress membrane thickness fluctuations and could present a potential mechanism for regulating signaling events. This is a pressing, yet underdeveloped, area of research where NSE can have significant impact in understanding the dynamic responses of membranes to protein binding and insertion on the length and time scales of key biological functions imparted by the interactions of proteins with cell membranes.

In summary, NSE has evolved over the last few years as a powerful tool for interrogating membrane dynamics over spatial and temporal scales of vital biological functions. The technique is rapidly gaining widespread interest and its potential in answering key questions in membrane function is becoming well recognized. The contrast variation capabilities within NSE have positioned it as a unique approach to measure mesoscopic membrane properties that would otherwise be challenging to obtain. Another significant advantage of NSE over traditional spectroscopy methods in studies of membrane dynamics is its overlap with the length and time scales accessible with MD simulations, allowing for synergistic experimental/computation studies to gain a molecular level understanding of the different molecular components making up membranes. Despite its promise, there are still some limitations in the use of NSE in biological membrane studies, including the requirement for large sample volumes, the difficulty in selective deuteration in biological systems, and the relatively low neutron flux on NSE spectrometers, which results in longer measurement times and limited beamtime availability. However, these shortcomings could be overcome in the near future with constant developments in neutron sources and instrumentation along with advances in deuteration facilities.

Disclosures

The authors have nothing to disclose.

Acknowledgements

R. Ashkar thanks M. Nagao, L.-R. Stingaciu, and P. Zolnierczuk for many useful discussions and for their frequent assistance with NSE experiments on their respective beamlines. The authors acknowledge the use of neutron spin echo spectrometers at NIST and ORNL. The NSE spectrometer at NIST is supported by the Center for High Resolution Neutron Scattering, a partnership between the National Institute of Standards and Technology and the National Science Foundation under agreement no. DMR-1508249. The NSE spectrometer at ORNL's Spallation Neutron Source is supported by the Scientific User Facilities Division, Office of Basic Energy Sciences, US Department of Energy. Oak Ridge National Laboratory is managed by UT-Battelle, LLC under US DOE Contract No. DE-AC05-00OR22725.

Materials

Chloroform (biotech grade) Sigma Aldrich 496189 Biotech. grade, ≥99.8%, contains 0.5-1.0% ethanol as stabilizer
Circulating water bath Julabo SE-12 Heating Circulator with smart pump, programmable temperature settings, and external sensor connection for measurement and control
Deuterium Oxide Cambridge Isotopes Laboratories DLM-4 Deuterated water; Heavy water (D2O) (D, 99.9%)
Digital Semi-Microbalance Mettler Toledo MS105 Semi-micro balance with 120 g capacity, 0.01 mg readability, high resolution weighing cell, ergonomic doors, and pipette-check application
Ethanol (molecular biology grade) Sigma Aldrich E7023 200 proof ethanol for molecular biology applications
Glass Pipets VWR 36360-536 Disposable Soda Lime glass Pasteur pipets
Glass Vials Thermo Scientific B7990-1 Borosilicate glass vials with PTFE/Silione septum caps
Lab grade freezer Fisher Scientific IU2886D Ultra-low temprature freezer (-86 to -50 C) for long-term storage of lipids and proteins
Lipids (protaited or perdeuterated) Avanti Polar Lipids varies by lipid Lipids can be purchased from Avanti in powder form or in a chloroform solution with the required amounts and deuteration schemes.
Millipore water purifier Millipore Sigma ZRQSVP3US Direct-Q® 3 UV Water Purification System which deliver both pure and ultrapure water with a built-in UV lamp to reduce the levels of organics for biological  applications
Mini Extruder Set Avanti Polar Lipids 610020 Mini-extruder set includes mini-extruder, heating block, 2 GasTight Syringes, and 2 O-rings, Polycarbonate Membranes, and Filter Supports
Quick Connect Fittings Grainger 2YDA1 and 2YDA7 Push-button tube fittings for QuickConnect water circulation applications, e.g. high temperature vesicle extrusion
Syringe Pump SyringePump.com New Era-1000 Fully programmable syringe pump for infusion and withdrawal; programs up to 41 pumping phases with adjustable pumping rates, dispensed volumes, and extrusion cycles
Ultrasonic bath Fisher Scientific CPX2800 Temperature controlled ultra sonic bath with programmable functionality for degassing and ultrasonic applications
Vacuum Oven Thermo Scientific 3608 0.7 cu ft vaccum oven with built-in-high-limit thermostat guards against overheating
Vortex Mixer Fisher Scientific 02-215-414 Variable speed, analog control that allows low rpm start-up for gentle shaking or high-speed mixing for vigorous vortexing of samples

References

  1. Singer, S. J., Nicolson, G. L. The fluid mosaic model of the structure of cell membranes. Science. 175 (4023), 720-731 (1972).
  2. Andersen, O. S., Koeppe, R. E. Bilayer thickness and membrane protein function: an energetic perspective. Annual Review of Biophysics and Biomolecular Structure. 36, 107-130 (2007).
  3. Lundbæk, J. A., Collingwood, S. A., Ingólfsson, H. I., Kapoor, R., Andersen, O. S. Lipid bilayer regulation of membrane protein function: gramicidin channels as molecular force probes. Journal of The Royal Society Interface. 7 (44), 373-395 (2010).
  4. Bradley, R. P., Radhakrishnan, R. Curvature-undulation coupling as a basis for curvature sensing and generation in bilayer membranes. Proceedings of the National Academy of Sciences of the United States of America. 113 (35), 117-124 (2016).
  5. Perozo, E., Cortes, D. M., Sompornpisut, P., Kloda, A., Martinac, B. Open channel structure of MscL and the gating mechanism of mechanosensitive channels. Nature. 418 (6901), 942-948 (2002).
  6. Jensen, M. &. #. 2. 1. 6. ;., Mouritsen, O. G. Lipids do influence protein function-the hydrophobic matching hypothesis revisited. Biochimica et Biophysica Acta (BBA) – Biomembranes. 1666 (1-2), 205-226 (2004).
  7. Rajendran, L., Simons, K. Lipid rafts and membrane dynamics. Journal of Cell Science. 118 (6), 1099-1102 (2005).
  8. Katchalsky, A., Spangler, R. Dynamics of membrane processes. Quarterly Reviews of Biophysics. 1 (2), 127-175 (1968).
  9. Rheinstädter, M. C. Collective molecular dynamics in proteins and membranes (Review). Biointerphases. 3 (2), 83-90 (2008).
  10. Fujiwara, T., Ritchie, K., Murakoshi, H., Jacobson, K., Kusumi, A. Phospholipids undergo hop diffusion in compartmentalized cell membrane. The Journal of Cell Biology. 157 (6), 1071-1082 (2002).
  11. Hac, A. E., Seeger, H. M., Fidorra, M., Heimburg, T. Diffusion in two-component lipid membranes–a fluorescence correlation spectroscopy and monte carlo simulation study. Biophysical Journal. 88 (1), 317-333 (2005).
  12. Heinrich, M., Tian, A., Esposito, C., Baumgart, T. Dynamic sorting of lipids and proteins in membrane tubes with a moving phase boundary. Proceedings of the National Academy of Sciences of the United States of America. 107 (16), 7208-7213 (2010).
  13. Hormel, T. T., Kurihara, S. Q., Brennan, M. K., Wozniak, M. C., Parthasarathy, R. Measuring lipid membrane viscosity using rotational and translational probe diffusion. Physical Review Letters. 112 (18), 188101 (2014).
  14. Dimova, R. Recent developments in the field of bending rigidity measurements on membranes. Advances in Colloid and Interface Science. 208, 225-234 (2014).
  15. Bassereau, P., Sorre, B., Lévy, A. Bending lipid membranes: Experiments after W. Helfrich’s model. Advances in Colloid and Interface Science. 208, 47-57 (2014).
  16. Monzel, C., Sengupta, K. Measuring shape fluctuations in biological membranes. Journal of Physics D: Applied Physics. 49 (24), 243002 (2016).
  17. Deserno, M. Mesoscopic membrane physics: concepts, simulations, and selected applications. Macromolecular Rapid Communications. 30 (9-10), 752-771 (2009).
  18. Reynwar, B. J., et al. Aggregation and vesiculation of membrane proteins by curvature-mediated interactions. Nature. 447 (7143), 461-464 (2007).
  19. Haswell, E. S., Phillips, R., Rees, D. C. Mechanosensitive channels: what can they do and how do they do it. Structure. 19 (10), 1356-1369 (2011).
  20. Phillips, R., Ursell, T., Wiggins, P., Sens, P. Emerging roles for lipids in shaping membrane-protein function. Nature. 459 (7245), 379-385 (2009).
  21. Dill, K. A., Chan, H. S. From Levinthal to pathways to funnels. Nature Structural Biology. 4 (1), 10-19 (1997).
  22. Henzler-Wildman, K., Kern, D. Dynamic personalities of proteins. Nature. 450 (7172), 964-972 (2007).
  23. Grimaldo, M., Roosen-Runge, F., Zhang, F., Schreiber, F., Seydel, T. Dynamics of proteins in solution. Quarterly Reviews of Biophysics. 52, 7 (2019).
  24. Lyman, E., Hsieh, C. -. L., Eggeling, C. From dynamics to membrane organization: experimental breakthroughs occasion a “modeling manifesto”. Biophysical Journal. 115 (4), 595-604 (2018).
  25. Arriaga, L. R., et al. Dissipative curvature fluctuations in bilayer vesicles: Coexistence of pure-bending and hybrid curvature-compression modes. The European Physical Journal. E, Soft Matter. 31 (1), 105-113 (2010).
  26. Honerkamp-Smith, A. R., Veatch, S. L., Keller, S. L. An introduction to critical points for biophysicists; observations of compositional heterogeneity in lipid membranes. Biochimica et Biophysica Acta (BBA) – Biomembranes. 1788 (1), 53-63 (2009).
  27. Veatch, S. L., Keller, S. L. Organization in lipid membranes containing cholesterol. Physical Review Letters. 89 (26), 268101 (2002).
  28. Heberle, F. A., et al. Bilayer thickness mismatch controls domain size in model membranes. Journal of the American Chemical Society. 135 (18), 6853-6859 (2013).
  29. Nickels, J. D., et al. The in vivo structure of biological membranes and evidence for lipid domains. PLOS Biology. 15 (5), 2002214 (2017).
  30. Simons, K., Ikonen, E. Functional rafts in cell membranes. Nature. 387 (6633), 569-572 (1997).
  31. van Meer, G., Voelker, D. R., Feigenson, G. W. Membrane lipids: where they are and how they behave. Nature Reviews. Molecular Cell Biology. 9 (2), 112-124 (2008).
  32. Liu, S. -. L., et al. Orthogonal lipid sensors identify transbilayer asymmetry of plasma membrane cholesterol. Nature Chemical Biology. 13, 268 (2016).
  33. Rothman, J., Lenard, J. Membrane asymmetry. Science. 195 (4280), 743-753 (1977).
  34. Ashkar, R., et al. Neutron scattering in the biological sciences: progress and prospects. Acta Crystallographica Section D. 74 (12), 1129-1168 (2018).
  35. Woodka, A. C., Butler, P. D., Porcar, L., Farago, B., Nagao, M. Lipid bilayers and membrane dynamics: insight into thickness fluctuations. Physical Review Letters. 109 (5), 058102 (2012).
  36. Chakraborty, S., et al. How cholesterol stiffens unsaturated lipid membranes. Proceedings of the National Academy of Sciences of the United States of America. 117 (36), 21896-21905 (2020).
  37. Arriaga, L. R., et al. Stiffening effect of cholesterol on disordered lipid phases: a combined neutron spin echo + dynamic light scattering analysis of the bending elasticity of large unilamellar vesicles. Biophysical Journal. 96 (9), 3629-3637 (2009).
  38. Nagao, M., Kelley, E. G., Ashkar, R., Bradbury, R., Butler, P. D. Probing elastic and viscous properties of phospholipid bilayers using neutron spin echo spectroscopy. The Journal of Physical Chemistry Letters. 8 (19), 4679-4684 (2017).
  39. Kelley, E. G., Butler, P. D., Ashkar, R., Bradbury, R., Nagao, M. Scaling relationships for the elastic moduli and viscosity of mixed lipid membranes. Proceedings of the National Academy of Sciences of the United States of America. 117 (38), 23365-23373 (2020).
  40. Rickeard, B. W., et al. Transverse lipid organization dictates bending fluctuations in model plasma membranes. Nanoscale. 12 (3), 1438-1447 (2020).
  41. Nickels, J. D., et al. Mechanical properties of nanoscopic lipid domains. Journal of the American Chemical Society. 137 (50), 15772-15780 (2015).
  42. Mezei, F. Neutron spin echo: A new concept in polarized thermal neutron techniques. Zeitschrift für Physik A Hadrons and Nuclei. 255 (2), 146-160 (1972).
  43. Hayter, J. B., Penfold, J. Neutron spin-echo integral transform spectroscopy. Zeitschrift für Physik B Condensed Matter. 35 (2), 199-205 (1979).
  44. Monkenbusch, M., Richter, D., Imae, T., Kanaya, T., Furusaka, M., Torikai, N. . Neutrons in Soft Matter. , 147-182 (2011).
  45. Pynn, R., Mezei, F., Pappas, C., Gutberlet, T. . Neutron Spin Echo. , 159-177 (2003).
  46. Holderer, O., et al. The JCNS neutron spin-echo spectrometer J-NSE at the FRM II. Measurement Science and Technology. 19 (3), 034022 (2008).
  47. Schleger, P., et al. The long-wavelength neutron spin-echo spectrometer IN15 at the Institut Laue-Langevin. Physica B: Condensed Matter. 241-243, 164-165 (1997).
  48. Holderer, O., Zolnierczuk, P., Pasini, S., Stingaciu, L., Monkenbusch, M. A better view through new glasses: Developments at the Jülich neutron spin echo spectrometers. Physica B: Condensed Matter. 562, 9-12 (2019).
  49. Farago, B., et al. The IN15 upgrade. Neutron News. 26 (3), 15-17 (2015).
  50. Ashkar, R. Selective dynamics in polymeric materials: Insights from quasi-elastic neutron scattering spectroscopy. Journal of Applied Physics. 127 (15), 151101 (2020).
  51. Pasini, S., Holderer, O., Kozielewski, T., Richter, D., Phoenix Monkenbusch, M. J-NSE- Phoenix, a neutron spin-echo spectrometer with optimized superconducting precession coils at the MLZ in Garching. Review of Scientific Instruments. 90 (4), 043107 (2019).
  52. Svergun, D. I., Koch, M. H. J., Timmins, P. A., May, R. P. . Small Angle X-Ray and Neutron Scattering from Solutions of Biological Macromolecules. , (2013).
  53. Eicher, B., et al. Joint small-angle X-ray and neutron scattering data analysis of asymmetric lipid vesicles. Journal of Applied Crystallography. 50 (2), 419-429 (2017).
  54. Heberle, F. A., et al. Model-based approaches for the determination of lipid bilayer structure from small-angle neutron and X-ray scattering data. European Biophysics Journal. 41 (10), 875-890 (2012).
  55. Jaksch, S., Koutsioubas, A., Mattauch, S., Holderer, O., Frielinghaus, H. Long-range excitations in phospholipid membranes. Chemistry and Physics of Lipids. 225, 104788 (2019).
  56. Jaksch, S., et al. Influence of ibuprofen on phospholipid membranes. Physical Review E. 91 (2), 022716 (2015).
  57. Armstrong, C. L., et al. Effect of cholesterol on the lateral nanoscale dynamics of fluid membranes. European Biophysics Journal. 41 (10), 901-913 (2012).
  58. Rheinstädter, M. C., Häußler, W., Salditt, T. Dispersion relation of lipid membrane shape fluctuations by neutron spin-echo spectrometry. Physical Review Letters. 97 (4), 048103 (2006).
  59. Armstrong, C. L., Häußler, W., Seydel, T., Katsaras, J., Rheinstädter, M. C. Nanosecond lipid dynamics in membranes containing cholesterol. Soft Matter. 10 (15), 2600-2611 (2014).
  60. Nickels, J. D., et al. Lipid rafts: buffers of cell membrane physical properties. The Journal of Physical Chemistry B. 123 (9), 2050-2056 (2019).
  61. Michonova-Alexova, E. I., Sugár, I. P. Component and state separation in DMPC/DSPC lipid bilayers: a Monte Carlo simulation study. Biophysical Journal. 83 (4), 1820-1833 (2002).
  62. Sugár, I. P., Thompson, T. E., Biltonen, R. L. Monte Carlo simulation of two-component bilayers: DMPC/DSPC mixtures. Biophysical Journal. 76 (4), 2099-2110 (1999).
  63. Mabrey, S., Sturtevant, J. M. Investigation of phase transitions of lipids and lipid mixtures by sensitivity differential scanning calorimetry. Proceedings of the National Academy of Sciences. 73 (11), 3862-3866 (1976).
  64. . Neutron activation and scattering calculator Available from: https://www.ncnr.nist.gov/resources/activation/ (2021)
  65. Scott, H. L., et al. On the mechanism of bilayer separation by extrusion, or why your LUVs are not really unilamellar. Biophysical Journal. 117 (8), 1381-1386 (2019).
  66. Ashkar, R., et al. Tuning membrane thickness fluctuations in model lipid bilayers. Biophysical Journal. 109 (1), 106-112 (2015).
  67. Carrillo, J. -. M. Y., Katsaras, J., Sumpter, B. G., Ashkar, R. A computational approach for modeling neutron scattering data from lipid bilayers. Journal of Chemical Theory and Computation. 13 (2), 916-925 (2017).
  68. Azuah, R. T. DAVE: a comprehensive software suite for the reduction, visualization, and analysis of low energy neutron spectroscopic data. Journal of Research of the National Institute of Standards and Technology. 114 (6), 341-358 (2009).
  69. Van Hove, L. Correlations in space and time and born approximation scattering in systems of interacting particles. Physical Review. 95 (1), 249-262 (1954).
  70. Zilman, A. G., Granek, R. Undulations and dynamic structure factor of membranes. Physical Review Letters. 77 (23), 4788-4791 (1996).
  71. Kelley, E. G., Butler, P. D., Nagao, M. . Collective dynamics in model biological membranes measured by neutron spin echo spectroscopy. , 131-176 (2019).
  72. Zheng, Y., Michihiro, N., Dobrin, P. B. Bending elasticity of saturated and monounsaturated phospholipid membranes studied by the neutron spin echo technique. Journal of Physics: Condensed Matter. 21 (15), 155104 (2009).
  73. Sharma, V. K., Qian, S. Effect of an antimicrobial peptide on lateral segregation of lipids: a structure and dynamics study by neutron scattering. Langmuir. 35 (11), 4152-4160 (2019).
  74. Boggara, M. B., Faraone, A., Krishnamoorti, R. Effect of pH and Ibuprofen on the Phospholipid Bilayer Bending Modulus. The Journal of Physical Chemistry B. 114 (24), 8061-8066 (2010).
  75. Lee, J. -. H., et al. Thermal fluctuation and elasticity of lipid vesicles interacting with pore-forming peptides. Physical Review Letters. 105 (3), 038101 (2010).
  76. Chakraborty, S., Abbasi, A., Bothun, G. D., Nagao, M., Kitchens, C. L. Phospholipid bilayer softening due to hydrophobic gold nanoparticle inclusions. Langmuir. 34 (44), 13416-13425 (2018).
  77. Hoffmann, I., et al. Softening of phospholipid membranes by the adhesion of silica nanoparticles – as seen by neutron spin-echo (NSE). Nanoscale. 6 (12), 6945-6952 (2014).
  78. Watson, M. C., Brown, F. L. H. Interpreting membrane scattering experiments at the mesoscale: the contribution of dissipation within the bilayer. Biophysical Journal. 98 (6), 9-11 (2010).
  79. Seifert, U., Langer, S. A. Viscous modes of fluid bilayer membranes. Europhysics Letters (EPL). 23 (1), 71-76 (1993).
  80. Bingham, R. J., Smye, S. W., Olmsted, P. D. Dynamics of an asymmetric bilayer lipid membrane in a viscous solvent. EPL (Europhysics Letters). 111 (1), 18004 (2015).
  81. Rawicz, W., Olbrich, K. C., McIntosh, T., Needham, D., Evans, E. Effect of chain length and unsaturation on elasticity of lipid bilayers. Biophysical Journal. 79 (1), 328-339 (2000).
  82. Doktorova, M., LeVine, M. V., Khelashvili, G., Weinstein, H. A new computational method for membrane compressibility: bilayer mechanical thickness revisited. Biophysical Journal. 116 (3), 487-502 (2019).
  83. Evans, E., Needham, D. Physical properties of surfactant bilayer membranes: thermal transitions, elasticity, rigidity, cohesion and colloidal interactions. The Journal of Physical Chemistry. 91 (16), 4219-4228 (1987).
  84. Lesieur, S., Grabielle-Madelmont, C., Paternostre, M. T., Ollivon, M. Size analysis and stability study of lipid vesicles by high-performance gel exclusion chromatography, turbidity, and dynamic light scattering. Analytical Biochemistry. 192 (2), 334-343 (1991).
  85. Heberle, F. A., et al. Direct label-free imaging of nanodomains in biomimetic and biological membranes by cryogenic electron microscopy. Proceedings of the National Academy of Sciences of the United States of America. 117 (33), 19943-19952 (2020).
  86. Cornell, C. E., Mileant, A., Thakkar, N., Lee, K. K., Keller, S. L. Direct imaging of liquid domains in membranes by cryo-electron tomography. Proceedings of the National Academy of Sciences of the United States of America. 117 (33), 19713-19719 (2020).
  87. Yao, X., Fan, X., Yan, N. Cryo-EM analysis of a membrane protein embedded in the liposome. Proceedings of the National Academy of Sciences of the United States of America. 117 (31), 18497-18503 (2020).
  88. Kučerka, N., Nieh, M. -. P., Katsaras, J. Fluid phase lipid areas and bilayer thicknesses of commonly used phosphatidylcholines as a function of temperature. Biochimica et Biophysica Acta (BBA) – Biomembranes. 1808 (11), 2761-2771 (2011).
  89. Nielsen, J. E., Bjørnestad, V. A., Lund, R. Resolving the structural interactions between antimicrobial peptides and lipid membranes using small-angle scattering methods: the case of indolicidin. Soft Matter. 14 (43), 8750-8763 (2018).
  90. Kučerka, N., et al. Lipid bilayer structure determined by the simultaneous analysis of neutron and X-ray scattering data. Biophysical Journal. 95 (5), 2356-2367 (2008).
  91. Kelley, E. G., Butler, P. D., Nagao, M. Scaling of lipid membrane rigidity with domain area fraction. Soft Matter. 15 (13), 2762-2767 (2019).
  92. Brüning, B. -. A., et al. Bilayer undulation dynamics in unilamellar phospholipid vesicles: Effect of temperature, cholesterol and trehalose. Biochimica et Biophysica Acta (BBA) – Biomembranes. 1838 (10), 2412-2419 (2014).
  93. Kučerka, N., et al. Areas of monounsaturated diacylphosphatidylcholines. Biophysical Journal. 97 (7), 1926-1932 (2009).
  94. Sharma, V. K., Mamontov, E., Anunciado, D. B., O’Neill, H., Urban, V. S. Effect of antimicrobial peptide on the dynamics of phosphocholine membrane: role of cholesterol and physical state of bilayer. Soft Matter. 11 (34), 6755-6767 (2015).
  95. Kelley, E. G., Butler, P. D., Nagao, M. Collective dynamics in lipid membranes containing transmembrane peptides. Soft Matter. , (2021).
  96. Yu, J., et al. Structure and dynamics of lipid membranes interacting with antivirulence end-phosphorylated polyethylene glycol block copolymers. Soft Matter. 16 (4), 983-989 (2020).
  97. Stingaciu, L. -. R., et al. Revealing the dynamics of thylakoid membranes in living cyanobacterial cells. Scientific Reports. 6 (1), 19627 (2016).
  98. Stingaciu, L. -. R., O’Neill, H. M., Liberton, M., Pakrasi, H. B., Urban, V. S. Influence of chemically disrupted photosynthesis on cyanobacterial thylakoid dynamics in synechocystis sp. PCC 6803. Scientific Reports. 9 (1), 5711 (2019).
  99. Miller, I. R. Energetics of fluctuation in lipid bilayer thickness. Biophysical Journal. 45 (3), 643-644 (1984).
  100. Nagao, M. Observation of local thickness fluctuations in surfactant membranes using neutron spin echo. Physical Review E. 80 (3), 031606 (2009).

Play Video

Cite This Article
Kumarage, T., Nguyen, J., Ashkar, R. Neutron Spin Echo Spectroscopy as a Unique Probe for Lipid Membrane Dynamics and Membrane-Protein Interactions. J. Vis. Exp. (171), e62396, doi:10.3791/62396 (2021).

View Video