Contrast-Matching Detergent in Small-Angle Neutron Scattering Experiments for Membrane Protein Structural Analysis and <em>Ab Initio</em> Modeling

Ryan C. Oliver; Swe-Htet Naing; Kevin L. Weiss; Sai Venkatesh Pingali; Raquel L. Lieberman; Volker S. Urban

doi:10.3791/57901

JoVE Journal > Biochemistry

Biochimie

Contrast-Matching Detergent in Small-Angle Neutron Scattering Experiments for Membrane Protein Structural Analysis and Ab Initio Modeling

Published: October 21, 2018

doi:

10.3791/57901

Ryan C. Oliver¹, Swe-Htet Naing², Kevin L. Weiss¹, Sai Venkatesh Pingali¹, Raquel L. Lieberman², Volker S. Urban¹

¹Neutron Scattering Division,Oak Ridge National Laboratory, ²School of Chemistry and Biochemistry,Georgia Institute of Technology

Summary

This protocol demonstrates how to obtain a low-resolution ab initio model and structural details of a detergent-solubilized membrane protein in solution using small-angle neutron scattering with contrast-matching of the detergent.

Abstract

The biological small-angle neutron scattering instrument at the High-Flux Isotope Reactor of Oak Ridge National Laboratory is dedicated to the investigation of biological materials, biofuel processing, and bio-inspired materials covering nanometer to micrometer length scales. The methods presented here for investigating physical properties (i.e., size and shape) of membrane proteins (here, MmIAP, an intramembrane aspartyl protease from Methanoculleus marisnigri) in solutions of micelle-forming detergents are well-suited for this small-angle neutron scattering instrument, among others. Other biophysical characterization techniques are hindered by their inability to address the detergent contributions in a protein-detergent complex structure. Additionally, access to the Bio-Deuteration Lab provides unique capabilities for preparing large-scale cultivations and expressing deuterium-labeled proteins for enhanced scattering signal from the protein. While this technique does not provide structural details at high-resolution, the structural knowledge gap for membrane proteins contains many addressable areas of research without requiring near-atomic resolution. For example, these areas include determination of oligomeric states, complex formation, conformational changes during perturbation, and folding/unfolding events. These investigations can be readily accomplished through applications of this method.

Introduction

Membrane proteins are encoded by an estimated 30% of all genes¹ and represent a strong majority of targets for modern medicinal drugs.² These proteins perform a wide array of vital cellular functions,³ but despite their abundance and importance — only represent about 1% of total structures deposited in the Research Collaboratory for Structural Bioinformatics (RCSB) Protein Data Bank.⁴ Due to their partially hydrophobic nature, structural determination of membrane-bound proteins has been exceedingly challenging.⁵^,⁶^,⁷

As many biophysical techniques require monodisperse particles in solution for measurement, isolating membrane proteins from native membranes and stabilizing these proteins in a soluble mimic of the native membranes has been an active area of research in recent decades.⁸^,⁹^,¹⁰ These investigations have led to the development of many novel amphiphilic assemblies to solubilize membrane proteins, such as nanodiscs,¹¹^,¹²^,¹³ bicelles,¹⁴^,¹⁵ and amphipols.¹⁶^,¹⁷ However, the use of detergent micelles remains one of the most common and straightforward approaches for satisfying the solubility requirements of a given protein.¹⁸^,¹⁹^,²⁰^,²¹^,²²^,²³^,²⁴^,²⁵ Unfortunately, no single detergent or magic mixture of detergents currently exists that satisfies all membrane proteins; thus, these conditions must be empirically screened for the unique requirements of each protein.²⁶^,²⁷

Detergents self-assemble in solution above their critical micelle concentration to form aggregate structures called micelles. Micelles are composed of many detergent monomers (typically ranging from 20-200) with hydrophobic alkyl chains forming a micelle core and hydrophilic head groups arranged in a micelle shell layer facing the aqueous solvent. The behavior of detergents and micelle formation has been classically described by Charles Tanford in The Hydrophobic Effect,²⁸ and sizes and shapes of micelles from commonly used detergents in membrane protein studies have been characterized using small-angle scattering.²⁹^,³⁰ Detergent organization about membrane proteins has also been studied, and the formation of protein-detergent complexes (PDCs) is expected with detergent molecules surrounding the protein in an arrangement that resembles the neat detergent micelles.³¹

One added advantage in using detergents is that the resulting micelle properties can be manipulated by incorporating other detergents. Many detergents exhibit ideal mixing, and select properties of mixed micelles may even be predicted from the components and ratio of mixing.²² However, the presence of detergent can still present challenges for biophysical characterizations by contributing to the overall signal. For example, with X-ray and light scattering techniques, signal from detergent in the PDC is practically indistinguishable from protein.³² Investigations with single-particle cryo-electron microscopy (cryo-EM) typically rely on trapped (frozen) particles; structural details of the protein are still obscured by certain detergents or a high concentration of detergent which adds to the background.³³ Alternate approaches toward interpreting the full PDC structure (including the detergent) have been made through computational methods which seek to reconstruct the detergent around a given membrane protein.³⁴

For the case of neutron scattering, the core-shell arrangement of detergent in the micelle produces a form factor which contributes to the observed scattering. Fortunately, solution components can be altered such that they do not contribute to the net observed scattering. This "contrast matching" process is achieved by substituting deuterium for hydrogen to achieve a scattering length density that matches that of the background (buffer). A judicious choice of detergent (with available deuterated counterparts) and their ratio of mixing must be considered. For detergent micelles, this substitution can be performed using a detergent with the same head group but having a deuterated alkyl chain (d-tail instead of h-tail). Since the detergents are well-mixed,³⁵ their aggregates will have a scattering length density that is the mole-fraction weighted average of the two components (h-tails and d-tails). When this average contrast is consistent with that of the head group, the uniform aggregate structures can be fully matched to remove all contributions to observed scattering.

We present here a protocol to manipulate the neutron contrast of detergent micelles by incorporating chemically identical detergent molecules with deuterium-labeled alkyl chains.¹⁹^,³⁶^,³⁷ This permits complete simultaneous contrast matching of micelle core and shell, which is a unique capability of neutron scattering.³⁵^,³⁸ With this significantly refined level of detail, contrast matching can enable otherwise unfeasible studies of membrane protein structures. Additionally, this contrast-matching approach could be extended to other systems involving detergent, such as polymer exchange reactions³⁹ and oil-water dispersants,⁴⁰ or even other solubilizing agents, such as bicelles,⁴¹ nanodiscs,⁴² or block copolymers.⁴³ A similar approach as outlined in this manuscript, but employing a single detergent species with partial deuterium substitutions on the alkyl chain and/or head group, was recently published.³⁷ While this can be expected to improve the random distribution of hydrogen and deuterium throughout the detergent compared to the approach presented here, the limited number of available positions on the detergent for substitution and two-step detergent synthesis required poses additional challenges for consideration.

Steps 1 and 2 of the protocol detailed below often overlap since initial experiment planning must be done to submit a quality proposal. However, proposal submission is considered here as the first step to emphasize that this process should be started well in advance of a neutron experiment. It should also be noted that a prerequisite step, which should be demonstrated by the proposal, is to have biochemical and physical characterization (including purity and stability) of the sample supporting the need for neutron studies. A general discussion of small-angle neutron scattering (SANS) is beyond the scope of this article. A brief but thorough introduction is available in the reference work Characterization of Materials by Kaufmann,⁴⁴ and a comprehensive textbook focused on biological small-angle solution scattering has recently been published.⁴⁵ Further recommended reading is given in the Discussion section. Small angle scattering uses the so called scattering vector Q as the central quantity that describes the scattering process. This article uses the widely accepted definition Q = 4π sin(θ)/λ, where θ is half the angle between incoming and scattered beam and λ is the wavelength of the neutron radiation in Angstroms. Other definitions exist that use different symbols such as 's' for the scattering vector, and that may differ by a factor 2π or by using nanometers in place of Angstrom (see also discussion of Figure 10).

Protocol

1. Prepare and Submit a Neutron Facility Beam Time and Instrument Proposal

Consult online resources for identifying neutron scattering facilities that provide general user neutron beam time access, such as Oak Ridge National Laboratory (ORNL). A map of neutron facilities and information about neutron research worldwide is available online.⁴⁶ Be aware that these facilities typically have regular calls for proposals; this determines when the next beam time will be available. Neutrons are used for a variety of applications; search for small-angle neutron scattering (SANS) instruments, particularly those with capabilities for biological samples.
Use resources provided by the neutron facilities to help navigate the neutron beam time proposal submission process. Consult with a Neutron Scattering Scientist (NSS) that is associated with the instrument for which you will apply. Review the neutron facility’s current information regarding proposal calls, deadlines, and advice for submission; for Oak Ridge this is available online.⁴⁷ Submit the beam time proposal following all guidelines for a successful submission.
If deuterated protein is required (see 2.2), neutron scattering facilities are often supported by laboratories and expertise dedicated to the production of deuterated materials. To request access to the Bio-Deuteration Lab (BDL) at ORNL, select the BDL as the second instrument in the proposal system. While the SANS proposal is reviewed for feasibility and by a scientific review committee, BDL requests are only screened for feasibility. To help with the BDL feasibility review, submit a completed information request form⁴⁸ that describes the protein expression protocol and estimated yield (available online).
After successful peer-review of the beam time proposal, confirm the awarded date of beam time in advance of the experiment. Ensure that all sample and buffer details and formulas are listed correctly in the proposal for proper review. Any changes to the proposed experimental plan, including changes in sample and buffer composition, require notification of the facility.
NOTE: Changes should be discussed as soon as possible with the assigned NSS.

2. Determine Neutron Contrast Match Points and Necessary Contrast for Protein Measurement

Gather information (from supplier’s product information, online databases, published values, etc.) relating to atomic compositions and volumes for the detergent system components to be contrast-matched. This information is used to determine neutron scattering length densities (SLDs) and thus contrast match points (CMPs) in solution, which is critical to obtaining quality SANS data and structural information for the membrane protein of interest in the context of a PDC. A table summarizing the physical properties and contrast match points for detergents commonly used in biophysical studies of membrane proteins is included (Table 1).
Determine neutron SLDs using the web application MULCh: ModULes For The Analysis Of Contrast Variation Data⁴⁹ (available online⁵⁰ free-of-charge). A link to the user manual is located on the home page. Access the website and read the accompanying documentation. Figure 1 and Figure 2 provide an overview of Contrast module input and output for two related examples.
NOTE: Other websites for calculating SLDs are available, such as the one maintained by the National Institute of Standards and Technology (NIST) Center for Neutron Research.⁵¹
1. Open the Contrast module by clicking ‘Contrast’ located in the left navigation pane. Provide text for a project title and enter details below for two components (e.g., detergent head group and tail) as ‘subunit 1’ and ‘subunit 2’.
2. Enter the molecular formula for each component in the ‘Formula’ text box.
  NOTE: The radio button ‘M’ must be selected under ‘Substance Type’ to enter a molecular formula. Additionally, use the letter ‘X’ in the formula to indicate readily exchangeable hydrogen atoms. Protein, RNA, or DNA sequences can be entered if the corresponding ‘P’, ‘R’, or ‘D’ radio buttons are selected. If multiple copies of a component exist within the subunit, the value for ‘Nmolecules’ can be changed accordingly. A practical example here relates to lipid-like detergents containing two identical tails, allowing the formula for one tail to be entered with Nmolecules equal to 2.
3. Enter the volume in cubic Angstroms for each component in the box below ‘Volume (Å³)’. Use Tanford’s formula²⁸ for alkyl chains of length n, V_tail = 27.4 + n * 26.9, to calculate approximate tail volume, or obtain molecular volumes from product information provided or published values in the literature.
4. Enter details for buffer components. Change the ‘Number dissolved species in the solvent’ using the dropdown menu to add or remove rows for each buffer component. In the case of micelle-forming detergents, include the free detergent monomers as separate buffer components at the detergent’s critical micelle concentration (CMC).
5. Click ‘Submit’ to perform the neutron contrast calculations and generate a results page that provides a table of scattering length density and contrast match points as well as formulae for determining these parameters at any given percentage of D₂O in the buffer. Review the scattering parameters with particular attention to the component CMPs. If the match points are similar (within 10% D₂O) and the free micelle concentration is low, then the average CMP can be used for contrast-matching the detergent contributions. In most cases, the CMP of the detergent head group and tail will differ by more than 10-15%, producing a core-shell form factor in the SANS data which obfuscates direct collection of the scattering signal from the protein alone.
Assess the contrast between detergent head group and tail. When the detergent head group and alkyl chain tail CMPs are not well-matched, the availability and incorporation of a deuterium-substituted detergent can permit scattering length densities of select components to be manipulated. In most cases, this will require forming a mixed micelle through the incorporation of an identical detergent with deuterium-substituted alkyl chains. The addition of deuterium raises the scattering length density of the core formed by the alkyl chain tails. The desired endpoint, in this case, is such that the head group shell CMP and alkyl chain core CMP have approximately equal values.
1. Raise the CMP of the detergent micelle core to be approximately equal to the CMP of the head group shell by determining the appropriate ratio of mixing between h-tails and d-tails that achieves complete contrast matching with the head groups. A prerequisite for this determination is knowledge of the CMP for a deuterium-substituted alkyl chain tail. A commercially-available counterpart to n-Dodecyl β-D-maltoside (DDM) is available with all alkyl chain hydrogen replaced by deuterium (d25-DDM). Repeat the contrast calculations outlined in step 2.1 for a ‘d-tail’ in place of the ‘h-tail’.
2. Calculate the mole ratio of h-tail to d-tail required for contrast-matching with the head group CMP. The following formula can assist with this determination:
  CMP_head = CMP_h-tail * χ_h-tail + CMP_d-tail * χ_d-tail
  where CMP is the scattering length density and χ is the volume fraction for the detergent head, h-tail, or d-tail component. For buffered solutions of DDM, incorporation of 44% by weight (43% by mole) of the total detergent as d25-DDM with deuterium-substituted alkyl chains produces a single CMP in 48.5% D₂O for both head and tail components.
Consider the degree of deuterium substitution for the membrane protein. To measure sufficient scattering signal from the protein, the CMP for the protein should be at least 15% D₂O in solution away from the detergent CMP and measurement conditions. The signal increases with the square of this% difference. Since the CMP of an unlabeled protein typically occurs with ~42% D₂O in solution (a CMP similar to many detergent head groups), deuterium labeling of the protein is almost always a necessity.⁵²

3. Express and Purify the Membrane Protein of Interest

Prepare LB (lysogeny broth) medium and grow Escherichia coli (E. coli) cells harboring an inducible expression vector for the target protein sequence.
1. Prepare minimal media. In either H₂O or D₂O, dissolve 7.0 g/L (NH₄)₂SO₄, 5.25 g/L Na₂HPO₄, 1.6 g/L KH₂PO₄, 0.50 g/L diammonium hydrogen citrate, 5.0 g/L glycerol, 1.0 mL/L of 20% w/v MgSO₄·7H₂O, and 1.0 mL/L of Holme trace metals (0.50 g/L CaCl₂·2H₂O, 0.098 g/L CoCl₂, 0.102 g/L CuSO₄, 16.7 g/L FeCl₃·6H₂O, 0.114 g/L MnSO₄·H₂O, 22.3 g/L Na₂EDTA·2H₂O, and 0.112 g/L ZnSO₄·H₂O).⁵³^,⁵⁴
2. Prepare appropriate antibiotic stock solutions using H₂O or D₂O.
3. Prepare an isopropyl-β-D-1-thiogalactopyranoside stock solution using H₂O or D₂O.
4. Sterile-filter all solutions into dry, sterile containers using bottle-top vacuum and syringe filters with a 0.22 micron pore size.
Adapt E. coli cells to deuterium-labeled medium prior to scale-up for overexpression of a deuterated protein.
1. Inoculate 3 mL of LB medium with an isolated colony from a Lysogeny Broth (LB) agar plate or cells from a frozen glycerol stock. Grow cells at 37 °C in a shaking incubator at 250 RPM or reduce the temperature to 30 °C or below to avoid overgrowth of the culture when incubating overnight.
  NOTE: The adaptation procedure can be varied.⁵⁵^,⁵⁶^,⁵⁷
2. Once the LB culture has grown to an optical density at 600 nm (OD₆₀₀) of ~1, dilute it 1:20 in 3 mL of H₂O minimal medium and grow to an OD₆₀₀ of ~1. Repeat the 1:20 dilutions using minimal medium containing 50, 75 and 100% D₂O (or up to the desired D₂O percentage).
  NOTE: As the D₂O content is increased, growth rates will decrease. The relationship between the D₂O percentage in the growth media and deuterium substitution in the protein has been reported in the literature.⁵⁸^,⁵⁹^,⁶⁰
3. Continue growing the culture in a bioreactor to increase the yield of deuterated cell mass (Figure 3).
  NOTE: Step-by-step procedures for bioreactor operation have been reviewed elsewhere.⁶¹^,⁶²^,⁶³^,⁶⁴
Harvest and lyse E. coli cells for membrane extraction and protein purification
1. Pellet the cells by centrifugation at ~6,000 x g for ~30-45 min at 4 °C.
2. Soften the cell pellet by soaking in buffer on ice for ~30 min. Then, gently resuspend the cell pellet in buffer by gently pipetting with a serological pipet, or gentle rocking on a 2D rocker, to a final concentration of ~1 g wet cell mass/10 mL buffer.
  NOTE: An example buffer is 50 mM HEPES (4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid), pH 7.5, 200 mM NaCl supplemented with a protease inhibitor tablet.
3. Lyse resuspended cells with three passes through a high-pressure homogenizer at 10,000 psi while chilling the outflow.
  NOTE: Depending on the scale required, other lysis methods may be used (e.g., by French press or sonication).
4. Pellet cell debris by centrifugation at ~4,000 – 5,000 x g for 15 min at 4 °C. Repeat this step until the supernatant is clarified.
5. Ultracentrifuge (~150,000 x g for 30-45 min) the supernatant from the prior step, which contains the soluble and membrane fractions. The pellet after ultracentrifugation contains the membrane fraction.
6. Resuspend the membrane pellet in a large Dounce homogenizer and isolate the membrane fraction again by ultracentrifugation (step 3.3.5). Repeat this step at least once to remove loosely bound proteins.
7. Solubilize membranes by gentle rocking for ~30 min at 4 °C in a buffer containing detergent suitable for purification. Solubilization is apparent when the suspension turns from cloudy to transparent. Remove insoluble material by ultracentrifugation (~150,000 x g for 30-45 min).
  NOTE: An example buffer is 50 mM HEPES, pH 7.5, 500 mM NaCl, 20 mM imidazole and 4% (w/v) DDM for purification by Ni²⁺ affinity chromatography.
Purify the protein following a protocol previously developed for protonated forms.
NOTE: See Naing et al. for an example purification protocol for a hexahistidine tagged M. marisnigri JR1 intramembrane aspartyl protease (d-MmIAP).⁶⁵ The final yield was ~ 1 mg of deuterated from 5 L deuterated cell culture growth, all of which was used in the SANS experiment (Figure 4).
Exchange the purified protein into the “Final Exchange Buffer” for contrast matching.
1. Prepare 200 mL of “Final Exchange Buffer” containing the correct mixture for contrast matching, e.g.,20 mM HEPES, pH 7.5, 250 mM NaCl, and 0.05% total DDM (0.028% DDM and 0.022% d25-DDM) in 49% D₂O.
2. Equilibrate a size exclusion column with >2 column volumes (CV) of “Final Exchange Buffer” in 3.5.1 using a Fast Protein Liquid Chromatography (FPLC) system.
3. After injecting the sample, run the column and isolate fractions corresponding to the protein of interest.
4. Concentrate the protein to the desired final concentration (2 – 5 mg/mL).
  NOTE: A volume of ~300 μL is required to fill a 1 mm cylindrical quartz cuvette used for SANS.
5. Reserve an aliquot of matched buffer for SANS background subtraction.
  NOTE: The aliquot can be taken directly from the prepared buffer, or ideally from a clean elution fraction at the end of the FPLC run.

4. Make Final Preparations for Beam Time and Collect SANS Data

After the proposal has been accepted and site access has been approved, complete all required training in advance of assigned beam time. This includes prearrival, web-based training as well as on-site radiological, laboratory, and instrument-specific training.
NOTE: Training requirements may dictate advance arrival for first time users. More information is available online.⁶⁶
Load samples and buffers for data collection.
NOTE: An overview of the biological small-angle neutron scattering (Bio-SANS) instrument sample environment area is shown in Figure 5.
1. Load samples and buffers into quartz cells. Quartz cells are available from the instrument scientist or facility. Use gel-loading tips to access the narrow opening of most sample cells.
2. Ensure that the beam shutter is closed, then approach the sample environment area and place the quartz cells in the sample changer. Record the sample changer position for each sample cell.
3. Check the area and ensure the beam path is free from obstruction, then leave the sample environment area and open the beam shutter.
  NOTE: Carefully follow all facility rules and regulations, obey all postings, and heed advice of the instrument scientist. Samples may not be removed from the beam and manipulated after exposure to the beam without following special precautions. Consult with a Radiological Control Technician (RCT) prior to these procedures.
Execute table scan for automated data collection
1. Operate the instrument through custom control LabVIEW-based software, Spectrometer Instrument Control Environment (SpICE). Follow instructions for instrument operation provided by the Neutron Scattering Scientist. An overview of the table scan operation windows is provided in Figure 6.
2. After entering the information required in the appropriate areas, execute the table scan and an automated data collection process will begin. Monitor progress via the ‘Table Scan Status’ tab.

5. Reduce SANS Data from 2D Image to 1D Plot

Identify each sample and buffer data file from the recorded scan numbers. Each measurement will be assigned a unique scan number. This information is helpful during data reduction to identify each corresponding sample and buffer pair.
Use MantidPlot software and Python script for reduction
1. Neutron Scattering Users are provided with account details to access their data on the Remote Analysis Cluster.⁶⁷ Log in to the Remote Analysis Cluster and execute “MantidPlot” from the command line. Figure 7 and Figure 8 are provided to assist these steps.
2. Obtain a User Reduction Script from the Neutron Scattering Scientist (this will usually happen during the experiment); this script is Python-based and will contain all the necessary calibration and scaling adjustments to reduce the 2D scattering image into a 1D plot of scattered intensities as a function of scattering angle, I(Q).
3. Open the provided User Reduction Script in MantidPlot and place the corresponding scan numbers or identification for the sample and buffer pairs in the appropriate list.
4. Execute the script to generate reduced data as a four-column text file (column order is Q, I(Q), I(Q) error, Q error) in the specified location. Right click the appropriate workspace in MantidPlot and select “Plot with errors…” for an initial examination of the scattering profiles.

6. Analyze Data for Structural Parameters of the Scattering Particle

Transfer the reduced data file from the analysis server to a local computer using secure FTP file transfer. This can be performed with the use of software such as FileZilla or CyberDuck. Additional instructions and details for connecting to the analysis server file system are provided on the analysis.sns.gov login page.
Download the ATSAS software suite⁶⁵^,⁶⁸ (whole ATSAS suite).⁶⁹ Individual programs⁷⁰ are also available online for analyzing the SANS data, namely determining structural parameters and ab initio models.
NOTE: Many options are available for SANS data analysis of biomacromolecules.⁴⁵
Use PRIMUS⁷¹ for data plotting, buffer scaling and subtraction, and Guinier analysis.
1. Launch the ATSAS ‘SAS Data Analysis’ application and load the reduced data files corresponding to the sample and buffer pair.
  1. To scale the buffer properly, select a data range at high Q (Q > 0.5 Å^-1) where both profiles are similar and flat, and click the ‘Scale’ button located under the ‘Operations’ tab. A scale factor will be applied to the buffer such that these two flat regions should overlap. Use caution: there should be a plausible physical reason that justifies scaling the buffer intensity to match the sample. If the discrepancy is large, consult an instrument scientist for valid approaches to background subtraction.⁴⁴^,⁴⁵^,⁷²
  2. Increase the data range to view all points. Click ‘Subtract’ to perform this operation. Right click the buffer-subtracted datafile in the legend and save this file for subsequent analysis. Figure 9 outlines the use of the ‘Operations’ tab.
2. Perform a Guinier analysis of the buffer-subtracted sample data using the ‘Analysis’ tab of the ‘SAS Data Analysis’ application.
  1. Be sure the proper file is selected in the list and click ‘Radius of Gyration’. An automated attempt to perform a Guinier fit will be provided by clicking on the ‘Autorg’ button.
  2. Expand the range of data used to include all of the low-Q data and begin narrowing the data range by taking away high-Q points until the Qmax*R_g limit is below 1.3. Use the plot of residuals to verify that data are linear in the fit range. The Guinier fit yields an approximate R_g and I₀ for the scattering particles.
  3. Make small adjustments to the fit region and monitor the sensitivity of these values to the range of data used in the fit. Figure 10A demonstrates the use of the ‘Radius of Gyration’ tool.
Obtain the probability distribution function (P(r)) in GNOM.⁷³ The output file generated here will be used as an input for the ab initio modeling process.
1. Start the Distance Distribution Wizard within the ‘Analysis’ tab. Obtaining a good fit to the data is essential to obtaining a quality model. For more information on obtaining accurate fits, please see the review⁷⁴ by Putnam, et al.
  NOTE: The GNOM program provides additional information about the fit beyond the Distance Distribution Wizard. Figure 10B demonstrates proper use of the ‘Distance Distribution’ tool, and Figure 11 illustrates some of the common errors encountered.
2. Determine Dmax, the maximum interatomic distance within the molecule.
  1. Estimate a value for Dmax by unchecking the box to force Rmax=0 and entering a large value for Dmax (150 Å, for example). The first x-intercept in the plot of P(r) yields this estimate.
  2. Make incremental changes to this Dmax value and the range of data used or number of points used in the fit to optimize the GNOM fit to the data and the resulting P(r) curve.
3. Continue to refine the GNOM fit and save the GNOM output files with good fit parameters for subsequent ab initio modeling steps.
Simulate SANS profiles from high resolution PDB models using the program CRYSON.⁷⁵
1. Open the program and select option ‘0’. Select the PDB file name in the popup file browser. Press enter to accept the default settings, except for specifying chain deuteration fractions and the fraction of D₂O in the solvent to achieve the proper contrast parameters.
2. Fit the model to an experimental data set by entering ‘Y’ and selecting the data file in the popup file browser. Accept the remaining default values, and output files will be written in the PDB file location. The ‘.fit’ file contains information for the predicted scattering profile from the high-resolution 3D model.

7. Create Ab Initio Models from the SANS Data.

NOTE: DAMMIF⁷⁶ and DAMMIN⁷⁷ within the ATSAS software suite is used to reconstruct dummy atom models (DAMs) using a simulated annealing process from the GNOM output, which contains P(r) data, or information about the probability or frequency of interatomic distances within the scattering particle. These programs may be run in batch mode or on the ATSAS-Online web server.

Start the PRIMUS Shape Wizard from the ‘Dammif’ button on the ‘Analysis’ tab of ‘SAS Data Analysis’, and use a manual selection of parameters. The input for the wizard is illustrated in Figure 12.
1. Define the Guinier range from the fit (step 6.3.2) and proceed to the next step using the navigation buttons. Define fit values from the P(r) plot (step 6.4) and proceed to the next step.
2. Provide parameters for the ab initio modeling process. Enter a prefix for the model output filenames (e.g., SANSEnvelope), choose either ‘fast’ or ‘slow’ mode modeling, enter a value for the number of models (17 recommended for statistical significance), select from available options for particle symmetry, anisometry, and angular scale (if known). Ensure that the boxes are checked to average models with DAMAVER and refine the final model with DAMMIN.
3. Initiate the process using the ‘commit’ button. Once the process is complete, view and save the working directory.
  NOTE: The typical modeling process for a single, small protein can take up to a few hours with a typical personal computer. Files are stored in temporary directories until saved.
Overlay and superimpose the final SANS envelope with a related high-resolution model using SUPCOMB within ATSAS. Place a copy of the high-resolution PDB model to be fit to the SANS envelope in the working directory. Execute SUPCOMB from the command line using the PDB filenames as two arguments with the template structure listed first: $ supcomb HiRes.pdb SANSEnvelope.pdb
NOTE: The second filename listed will be rewritten as a new file with an “r” appended to the end and having new coordinates for superimposition on the template structure.
Visualize the ab initio model results using PyMOL (pymol.org), a 3D molecular graphics viewing program. Figure 13 provides an overview of the PyMOL operations.
NOTE: There are publication guidelines for structural modeling of small-angle scattering data from biomolecules in solution.⁷⁸
1. Launch the PyMOL application and open the .pdb files corresponding to the 3D SANS envelope and the SUPCOMB-aligned high-resolution structure to be viewed.
2. Visualize the SANS envelope by representing this model as a surface. Click the ‘S’ button next to the model and select ‘Show: As: Surface’.
3. Visualize the high-resolution protein backbone structure by representing this model as a cartoon. Select the ‘S’ button next to this model and select ‘Show: As: Cartoon’. Display the chain as a rainbow from N- to C- terminus by selecting the ‘C’ button and ‘By Chain: Chainbows’.
4. Modify the transparency of the surface representation for clarity. From the ‘Setting’ menu bar, select ‘Transparency » Surface » 40%’.
5. Make the background white and non-opaque. From the ‘Display’ menu bar, select ‘Background »’ and select ‘Opaque’ to uncheck this option and ‘White’ to mark this color for background.
  NOTE: Additional operations and uses for PyMOL can be found at PyMolWiki.org.

Representative Results

A beam time and instrument proposal should clearly convey all information needed to the review committee so that a valid assessment of the proposed experiment can be made. Communication with an NSS is highly suggested for inexperienced users. The NSS can assess initial feasibility and guide proposal submission to emphasize feasibility, safety, and the potential for high-impact science. The information provided in the proposal should include background information and context for the significance of the research; knowledge that is expected to be gained and how this impacts the current understanding in the related field of science; a description of the work, samples, methods, and procedures that will be employed; and, if applicable, previous productivity of the team at the facility, including relevant publications and results. Useful resources, such as proposal templates and tips for preparing the proposal, are available to Users via the Neutron Science User Portal (neutrons.ornl.gov/users).

Experiment planning is a dynamic process that often begins during the initial stages of proposal submission, but may not be fully conceived until just prior to the experiment. However, keep in mind that any changes that deviate significantly from the description in the proposal (including changes to buffer conditions or sample composition) must be approved by the facility prior to the start of the experiment.

This protocol assumes that a method for expressing and purifying the membrane protein of interest into detergent micelles in solution has been successfully demonstrated. In this case, the membrane protein of interest is an intramembrane aspartyl protease (IAP), which has an established purification protocol and has been previously determined to be soluble and active in buffer containing DDM micelles.

Here, we demonstrate neutron contrast calculations using the MULCh: Contrast module for a solution of IAP protein in contrast-matched mixed DDM / d25-DDM micelles. The strategy outlined herein was to first determine the degree of mixing between the two detergents necessary to achieve a complete contrast match of the micelle. The end result was to determine the relative contrast of each component and the background (H₂O/D₂O ratio) necessary to contrast-match the scattering from detergent in order to observe only the protein scattering.

Figure 1 demonstrates a proper input for the MULCh Contrast module and the resulting output for contrast calculations of the DDM detergent head group and tail as subunits 1 and 2, respectively. The formula used for the head group is C₁₂H₁₄X₇O₁₁ with a volume of 348 Å³, and the alkyl chain tail formula is given as C₁₂H₂₅ with a volume of 350 Å³.

It is apparent from the Contrast module output that the DDM head group and tail have different scattering length densities, and thus contrast between the two components will be observed. For DDM, the head groups have a CMP in 49% D₂O while the alkyl chain tails have a CMP in 2% D₂O, with their average CMP occurring in 22% D₂O. Therefore, the aim of the next step will be to design a mixed micelle that incorporates deuterium-substituted alkyl chains to increase the average CMP of the detergent tails to match the CMP of the head groups. The contrast calculation was repeated for a substituted detergent, DDM with a fully-deuterated tail (d25-DDM), which similarly results in contrast between the head group and tail. However, contrast values of the h-tails and d-tails are significantly different. Recall that detergent mixtures are known to produce well-mixed micelles in solution,²² thus a blend of these h- and d-tails in the micelle core should produce an average SLD equal to that of the head group shell, yielding a micelle with single CMP. Since the head group is common to both DDM and d25-DDM, the strategy is to find a mixture of these detergents that produces an average contrast value from the mixed tails that matches the contrast of the head groups.

The target average CMP for the detergent tails is that of the maltoside head groups, or 49% D₂O. To estimate the ratio of mixing the average CMP of each h-tail and d-tail component must be known. These values for some common detergents and commercially-available deuterium labelled counterparts are provided in Table 1. Using these CMP values for DDM and d25-DDM, the mole fraction of h-tails and d-tails that satisfies the equation provided in step 2.2 is χ_d-tails = 0.43.

Figure 2 demonstrates a more advanced input to the MULCh Contrast module that can be used to confirm the final mixed-micelle contrast match conditions and determine the degree of deuteration necessary on the protein. Here, subunit 1 refers to the contrast-matched mixed micelle with 2 components, DDM and d25-DDM, while subunit 2 refers to the M. marisnigri JR1 intramembrane aspartyl protease (MmIAP) given by its amino acid sequence. The SLD of the protein has been elevated by expression in deuterated growth media to yield a CMP for the protein in buffer containing ~92% D₂O. The degree of separation between the protein's match point in 92% D₂O and the measurement condition (48.5% D₂O) suggests that sufficient scattering signal will be obtained from the protein.

Production of d-MmIAP was carried out with Rosetta 2 E. coli cells harboring the pET22b-MmIAP vector. Minimal medium with unlabeled glycerol in 90% D₂O was selected as the growth medium. After adaptation to 90% D₂O minimal medium, the culture volume was scaled up to 400 mL and used to inoculate 3.6 L of fresh 90% D₂O minimal medium in a 5.5 L bioreactor vessel.

Figure 3 provides a trace of the process values during the fed-batch cultivation. At the time of inoculation, the temperature set point was 30 °C, the dissolved oxygen (DO) set point was 100%, the agitation set point was 200 rpm, and the flow rate of compressed air was 4 L/min. The pH was held above a set point of 6.9 using a 10% w/v solution of sodium hydroxide in 90% D₂O. Once the DO spike signaled depletion of the initial 5 g/L glycerol, feeding(30% w/v glycerol, 0.2% MgSO₄ in 90% D2O) was initiated, which continued throughout the cultivation. After around 7 hours of feeding, the culture temperature was reduced to 18 °C and isopropyl-β-D-1-thiogalactopyranoside was added to a final concentration of 1 mM. Upon harvesting the culture, approximately 145 g wet weight of cell paste was collected via centrifugation (6,000 x g for 45 min)

Purification of d-MmIAP proceeded as for the protonated enzyme except that enzyme yield was significantly lower and final purity was somewhat lower than typical. From 5 L, ~20 g of wet membranes was isolated, all of which was solubilized in DDM for purification and loaded onto the first Ni²⁺ affinity column. To maximize purification yield, the flow-through fraction was diluted and re-rerun on Ni²⁺ affinity column. The procedure was repeated a third time before polishing the concentrated d-MmIAP sample by size exclusion chromatography (Figure 4).

Figure 5 depicts a quartz sample cell and its related sample changer setup. Sample environments are managed by the Bio-SANS operations team and NSS. A variety of different environments can be arranged to perform measurements with temperature control, humidity control, mechanical tumbling, high temperature, and sustained pressure, among others.

Figure 6 demonstrates Bio-SANS operations and execution of the automated table scans. The table scan is configured to be intuitive and user-friendly. On the first tab (Load Samples), information is provided about the sample changer and setup, such as identification of samples at each position in the sample changer, cell thicknesses, and sample type (open beam, sample, background, etc.). Additional metadata tags can be applied here also. On the next tab (Plan Experiment), the instrument configuration and any associated parameters (such as temperature control) are defined. This step is arranged by the NSS at the start of the experiment. The user must simply input the measurement times for each sample (in seconds). The final tab (Execute Scans) provides a summary of the table scan data and allows the automated scan to be executed.

After the 2D scattering images are recorded, this data must be reduced to a 1D plot of the intensity versus Q – a function of the scattering angle. During data reduction, image corrections such as pixel sensitivity masks and dark background subtractions are also applied. MantidPlot software was designed to interpret the raw Bio-SANS instrument data. The NSS will generally assist in reducing and retrieving the final data at the end of the experiment.

Figure 7 provides an overview of the remote analysis cluster access to MantidPlot and the Script Window used for data reduction. Raw data are accessible at the remote analysis cluster (analysis.sns.gov) via UCAMS/XCAMS user authentication. After logging into the remote desktop, the MantidPlot software can be launched from a terminal window by simply executing 'MantidPlot' under any path. Open the Script Window in MantidPlot and edit the User Reduction Script as directed by the NSS. Executing this script will generate the reduced data files, which can then be transferred to a local machine for analysis using secure FTP.

Figure 8 demonstrates plotting and viewing the data after executing the User Reduction Script. Reduced data will appear in the Workspaces window. By default, final merged data will have the suffix _f appended. Right-click will allow Plot Spectrum with Errors to be selected and data to be plotted. Display preferences are accessed by selecting Preferences from the 'View' menu bar. Formatting of axes and plotting range is easily accessible by double-clicking the axes. Scattering profiles of buffers, which contain no scattering particles of significant size should appear as a relatively flat line (constant intensity). For samples, intensity should be observed at lower scattering angles (Q) with exponential decay to a flat incoherent background.

Figure 9 provides an overview of proper buffer subtraction using the SAS Data Analysis (or primusqt) software package following data reduction. Quite often the incoherent background (intensity at high-Q) is slightly mismatched between the sample and buffer. For proper subtraction, data in this region should overlap. The scale correction applies an intensity scale factor to the data which results in overlapping data, and the subtraction can be performed. If the scale factor results in buffer data points with larger intensity than corresponding points in the sample, then these points will be negative and not rendered in a log plot. If negative values in the buffer-subtracted file are abundant, make a minor reduction in the buffer scale factor and perform the subtraction again.

Figure 10 provides example usages of the Primus Guinier and Distance Distribution Wizards. A Guinier analysis provides a preliminary estimate of the particles radius of gyration from the low-angle scattering data. As data approach zero, a linear region should be present in the data. Fitting of a line in this region allows the R_g to be determined from the slope and an I₀ to be extrapolated from the intensity intercept at Q=0. An upward trend in the data as Q approaches zero indicates aggregation in the sample, while downward trends are usually indicative of interparticle repulsions. These trends are most apparent when the distribution of residuals is non-stochastic. The upper limit of the Guinier fit for globular particles is constrained by Q*R_g < 1.3 ('s' is used in place of Q in the ATSAS software).

The R_g determined for d-MmIAP from this Guinier fit was 15.4 ± 1.1 Å with an estimated I₀ given as 0.580 ± 0.001.

The distance distribution wizard allows systematic changes to be made to parameters of the P(r) fit. The P(r) curve is an indirect Fourier transform of the curve fit to the experimental data. Therefore, it is essential to verify that the fit in the left panel accurately reflects the experimental data. For the fit shown in this example, a Dmax of 46 Å was obtained.

Figure 11 illustrates common errors in Dmax selection. A Dmax that is artificially large often causes P(r) values to become negative as the P(r) function oscillates about the x-axis. A Dmax that is artificially small leads to a truncated P(r) curve with an abrupt transition to Rmax=0.

The first steps in the Primus Shape Wizard are identical to the Guinier and Distance Distribution Wizards. Once this information has been provided to obtain good fits to the experimental data, the ab initio model setup is configured.

Figure 12 illustrates a proper input for the Primus Shape Wizard. Here, experimental SANS data for the IAP was provided with a scale in Angstroms, and the expected file prefix of "SANSEnvelope" was given. A set of 17 initial dummy atom models were requested to be generated using "fast" model mode with no symmetry (P1) applied and no anisometry selected. The set of 17 models will be aligned, averaged, and filtered with DAMAVER, and the average model refined using DAMMIN. Inspection of the prefix-damsel.log text document provides a summary of the selection criterion and any models discarded from the set. The final refined model was written as SANSEnvelope-1.pdb.

The program SUPCOMB is used to perform a superimposition of the SANS envelope model with the high-resolution model, with only the two PDB structure files as input. By default, "r" is appended to the reoriented and superimposed file name.

Once the SANS envelope and high-resolution structure have been superimposed, these models can be visualized using any molecular graphics viewing program. The results from PyMOL are suitable for publication quality images and can be used to emphasize biophysical results relating to the structural investigation. Here, we demonstrate a few basic PyMOL operations used to provide 3D representations and views of the resulting superimposed structures.

Figure 13 provides an overview of the PyMOL visualization process. Once the PDB structure files have been opened in PyMOL, the models should be visible in the PyMOL Viewer window. Change the representations of each model using the 'S' button next to each file name. Use a surface representation for the SANS envelope and a cartoon representation of the protein backbone from the high-resolution model. Select a suitable color scheme from options available under the 'C' button. For protein chains, a "chainbow" coloring provides a color gradient from N- to C- terminus to aid interpretation. Transparency is applied to the surface representation to allow better visualization of the protein structure within the envelope. For publication images, a white background is recommended. Once these visualization queues have been applied, the 3D structures can be examined. Manipulations of the perspective and molecule rotation/translation can be performed by clicking and dragging the structure.

Figure 1. Usage of the MULCh Contrast module to determine contrast parameters for components of dodecyl maltoside (DDM). Input values are shown in the screen on the left with the subsequent output to the right. The Contrast module input page is accessed by clicking 'Contrast' (green circle) from the navigation menu on the left side of each screen. Input areas for the project title, buffer components, and subunits 1 and 2 are labelled as such. Output pages contain important particle information and contrast properties for the described system. Relevant tables and calculated match points have been boxed in orange for clarity. Please click here to view a larger version of this figure.

Figure 2. Usage of the MULCh Contrast module to determine contrast parameters for components of the membrane protein-detergent complex. In this instance, input has been given to describe the overall PDC system in buffered solution. Subunit 1 describes the contrast-matched mixed micelle composed of DDM with 43% by mole as d25-DDM. Subunit 2 describes the membrane protein, with the amino acid sequence for MmIAP input as the formula. The deuteration level (green circle) describes the protein's degree of deuterium substitution, and has a direct impact on the SLD and calculated match-point that will be estimated for the protein. Results (orange box) indicate that match-point for the mixed detergent will occur in 48.5% D₂O, while the protein's match-point occurs at ~90% D₂O. Please click here to view a larger version of this figure.

Figure 3. Sample preparation – Bioreactor trace. The process values displayed are temperature set point (orange line), dissolved oxygen (DO) set point (blue line), agitation set point (red line), and compressed air flow rate (pink line). Pump 1 (green line) was used for pH control. The DO spike at 18:40 signaled depletion of the intial glycerol and was used to initiate feeding of glycerol solution via pump 2 (black line). X-axis denotes hours. Please click here to view a larger version of this figure.

Figure 4. Sample preparation – FPLC and SDS-PAGE characterization of the sample. Final size exclusion chromatogram for d-MmIAP equilibrated with 20 mM HEPES, pH 7.5, 250 mM NaCl, 48.5% D₂O and 0.05% total DDM, of which 44% (w/v) is tail-deuterated d25-DDM. Inset: SDS-PAGE analysis with Coomassie staining. Pooled fractions are labeled A, B, and C. Region "B" was used in the SANS experiment. Annotated molecular weight marker is in kDa. Image reproduced with permission from Naing et al.⁶⁵ Please click here to view a larger version of this figure.

Figure 5. Data collection – Quartz banjo sample cell and autosampler setup. (A) An empty quartz banjo cell is shown resting against the autosampler block. Cells are filled with sample and inserted into one of the 15 available positions. (B) The sample block is mounted behind an aperture selector (not shown) between the Beam Tube Extender and Silicon window at the entrance to the detector tank. Hoses connecting to a temperature-controlled water bath and channels throughout the sample block provide temperature regulation at the sample positions. The arrow indicates the direction of the neutron beam. Please click here to view a larger version of this figure.

Figure 6. Data Collection – Overview of SPICE operations and table scan. The table scan operations are accessible from the 'Table Scan' button on the far left menu of the SPICE instrument control software. Green arrows denote the active tab at the top of the screen and orange boxes highlight areas in the table for active user inputs. (A) The first tab of the table scan is used to provide information about the sample changer and sample cell positions. Provide labels, sample thicknesses, and sample types for each position used in the sample changer. Checking the 'Do Scan' box will add the corresponding row to the table scan queue. (B) On the second tab, details about the instrument configuration and measurement time for each sample (in seconds) are recorded. Instrument configurations are preconfigured by the Neutron Scattering Scientist and provided to users based on details provided in the beam time and instrument proposal. Additional parameters can be added to the instrument control, such as the ability to define temperature setpoints and hold times, throughout the measurement queue. (C) The final tab provides an overview of the measurements to be made for each row in the table. If no changes are needed, the automated scan is initiated by clicking the 'Execute Scans' button. Please click here to view a larger version of this figure.

Figure 7. Data Reduction – Overview of reduction script and MantidPlot operations. MantidPlot software is accessible on the analysis cluster by typing "mantidplot" in a terminal command prompt at any active directory (the yellow circle and arrow are added for emphasis). Once the software is open, access the script window (toggled by button marked in green) and open the script provided by the Neutron Scattering Scientist. Follow the instructions provided in the script, which should only require the scan numbers for each sample and buffer to be entered in a list for automated reduction, as well as scan numbers for the empty cell for subtraction and empty beam for transmission and beam center determination. Please click here to view a larger version of this figure.

Figure 8. Data Analysis – Plotting of data using MantidPlot. Data are referenced as 'Workspaces' in MantidPlot. The workspace area will be populated with filenames as produced by the user reduction script. Workspace data for 1D data sets can be plotted by right-clicking the workspace and selecting "Plot Spectrum with Errors…". Additional data can be added to the current plot by simply dragging and dropping workspaces. Formatting of the plot window (such as for log-log plots) can be performed by selecting 'Preferences' from the 'View' menu bar, or double-clicking on the axis labels. Please click here to view a larger version of this figure.

Figure 9. Data Analysis – ATSAS software: basic operations and buffer subtraction. The primusqt application (SAS Data Analysis) provides visualization for proper background (buffer) subtraction. The buffer file should be scaled for subtraction such that data overlap at high-Q (where scattering is flat as a result of remaining signal from incoherent background). When this high-Q range of data is selected, the Scale operation will apply a scale factor to the lower data file in the table that produces overlap in the defined region. After scaling, the buffer file can be subtracted to yield the net scattering profile of the particle. Please click here to view a larger version of this figure.

Figure 10. Data Analysis – Guinier and P(r) (distance distribution) determination. (A) The Primus Guinier Wizard is used to perform a Guinier analysis, providing an initial estimate of I₀ and R_g. The particle type and range of data to fit are used to define the fit. A plot of the fit and corresponding residuals are shown to the left, while Guinier terms (I₀ and R_g), Q*R_g limits, and quality of linear fit are provided as output in the area marked by the orange box. (B) The Primus Distance Distribution Wizard is used to perform the distance distribution analysis, providing the P(r) curve which defines the probability of interatomic distances within the scattering particle and includes the Dmax or maximum interatomic distance. Parameters indicated by the orange box can be systematically investigated to determine a proper distance distribution function. Please click here to view a larger version of this figure.

Figure 11. Improper results from the distance distribution analysis. A properly selected Dmax should produce a P(r) curve that peaks with a gradual decay to zero. Dmax values that are too large will likely produce negative P(r) values or oscillations near the x-axis at high values of r. Dmax values that are artificially small often result in P(r) curves where the upper bound appears truncated. Please click here to view a larger version of this figure.

Figure 12. Modeling – Ab initio model setup using the Primus Shape Wizard. Ab initio modeling parameters are provided in a single input window. A prefix is supplied for output file names with other options for processing available. Annealing procedure can be fast (bigger beads with faster cooling) or slow (smaller beads with slower cooling). Number of repetitions should be of sufficient size to examine reproducibility of model features. Particle symmetry and anisometry can be provided, if known. Angular scale should correspond to the units of the measured data. DAMAVER averaging will align and average all dummy atom models, and then apply a selection criterion which excludes any outlier models. A core of fixed atoms from the averaged model will be further refined using DAMMIN if this option is selected. This refined model should represent the 3D low-resolution envelope of the experimental SANS profile. Please click here to view a larger version of this figure.

Figure 13. Visualization – Experimental SANS envelope and overlay with high-resolution model using PyMOL. After opening the PDB structure files in PyMOL, models appear in the PyMOL Viewer. Active models are listed in the table to the right, along with action buttons which can be used to manipulate the model and its representation. Basic operations are provided for visualizing these models which allow regions of agreement (or mismatch) between the SANS envelope and high-resolution structure to be identified. Please click here to view a larger version of this figure.

Table 1. Physical properties of detergents commonly used in membrane protein investigations. Please click here to download this file.

Discussion

Structural biology researchers take advantage of complementary structural techniques like solution scattering to obtain biochemical and structural details (such as overall size and shape) from biomolecules in solution. SANS is a particularly attractive technique for determining low resolution structures of membrane proteins, a core focus of modern structural biology and biochemistry. SANS requires quantities of purified proteins comparable to those of crystallographic trials (1 mg/sample). The recently expanding commercial availability of high-purity deuterated detergents relevant to membrane protein studies makes accessible a means to manipulate this hydrogen/deuterium content for SANS experiments of membrane protein-detergent complexes, allowing the protein signal to be recorded directly. When successful, an ab initio model molecular envelope can be calculated from data collected over just several hours. Thus, membrane protein biochemists, biophysicists, and structural biologists can readily take advantage of SANS to obtain coveted initial 3D models of membrane proteins in solution, structures in complex with binding partners or substrates, and models obtained by SANS can be used for phasing diffraction data. Routine usage of SANS to characterize membrane proteins would be transformative for the membrane protein biochemistry field and have ripple effects from basic structure function to drug discovery and development. The added advantage of neutron contrast matching with deuterium substitution makes SANS a valuable technique for studying protein-protein, protein-DNA, or other biomolecular complexes, which can be readily manipulated in their hydrogen/deuterium content.

For additional details about small-angle scattering instrument design and theory, the following reviews are recommended: The Bio-SANS instrument at the High Flux Isotope Reactor of Oak Ridge National Laboratory,⁷⁹ Small-angle scattering for structural biology – expanding the frontier while avoiding the pitfalls,⁷² and Small-angle scattering studies of biological macromolecules in solution.⁸⁰ The Neutron Scattering Scientist can also discuss current instrument configurations and the optimal parameters for a given system. For instance, data at lower Q (which is a function of the scattering angle and neutron wavelength) is generally desired for larger systems. The most common instrument configuration at Bio-SANS provides a Q_min of 0.003 Å^-1 and is suitable for larger protein complexes up to hundreds of kDa. A solution containing ~3 mg mL^-1 of membrane protein of approximate size of 30 kDa (monomer) in contrast-matched micelles with 48.5% D₂O in solution requires a typical measurement time around 8 hours each for sample, buffer, and empty cell (24 hours = 1day total). Instrument measurement times at a given contrast condition can be approximately scaled according to the product of the scattering particle's molecular mass and concentration. However, scattering particles must be measured under non-interacting conditions, including any non-specific protein aggregation, which places a practical upper limit on the concentration of the sample. Hour-long measurements place limits on the stability of certain proteins. Fortunately, SANS measurements can be done routinely at refrigerated temperatures to improve protein life time, and the employed beam of low-energy neutrons does not cause any radiation damage during the measurement.

An overview of this process and determination of many key factors toward the successful recording of the neutron scattering from a sample measured at the contrast match point of the detergent are presented in this manuscript. This includes the critical step toward obtaining a complete match of the aggregate detergent by designing a detergent mixture with a uniform neutron contrast between the detergent head group and alkyl chain components. Measurements at this single detergent match point provide a scattering profile attributed to only the membrane protein of interest and allowed an ab initio envelope representing the membrane protein to be reconstructed from the data. The data analysis and ab initio modeling protocols also demonstrate the potential information obtained from such investigations, which could aim to address overall structure, conformational changes, oligomeric states, among others. One limitation is that to date only very few detergents are commercially available with deuterium substitutions.¹⁹

While perdeuteration may not be necessary, the aim in this case should be to achieve a protein with >65% deuterium labeling in the protein. Alternatively, if detergent with selectively-deuterated head groups and tails is available, and the CMP of the detergent micelle occurs near 100% D₂O, then sufficient contrast from the protein can be achieved without the need for deuterium labeling. Determine the appropriate deuteration level of the protein to achieve sufficient scattering when measured at the contrast-match point of the detergent.There are numerous challenges associated with membrane protein expression and purification,⁸¹ which remain beyond the scope of this article. Unfortunately, high levels of deuteration are currently not (yet) possible for eukaryotic systems. While we realize this is a limitation for some eukaryotic proteins, most membrane proteins have bacterial orthologs for which this method is suitable.

Divulgations

The authors have nothing to disclose.

Acknowledgements

The Office of Biological and Environmental Research supported research at ORNL's Center for Structural Molecular Biology (CSMB) and Bio-SANS using facilities supported by the Scientific User Facilities Division, Office of Basic Energy Sciences, US Department of Energy. Structural work on membrane proteins in the Lieberman lab has been supported by NIH (DK091357, GM095638) and NSF (0845445).

Materials

Amicon Ultra MWCO 50KDa concentrator	EMD Millipore	UFC905096	labware
Ammonium citrate dibasic	Fisher Scientific	A663	medium component
Ammonium sulfate	EMD Millipore	2150	medium component
Bioflo 310 Bioreactor System	Eppendorf	M1287-2110	equipment
Calcium chloride dihydrate	Acros	423525000	medium component
Carbenicillin	IBI Scientific	IB02025	antibiotic
Chloramphenicol	EMD Millipore	3130	antibiotic
Cobalt (II) chloride	Acros	AC21413-0050	medium component
Copper (II) sulfate	Acros	AC19771-1000	medium component
Deuterium oxide	Sigma-Aldrich	756822	medium component
Drierite Gas Purifier	W.A. Hammond Drierite Co. Ltd.	27068
EDTA, disodium, dihydrate	EMD Millipore	4010	medium component
Emulsiflex-C3	Avestin	EF-C3	equipment
Äkta Purifier UPC100	GE Healthcare		equipment
Glycerol	Sigma-Aldrich	G5516	medium component
HEPES	Sigma-Aldrich	H4034
HiPrep 16/60 Sephacryl S-300 HR column	GE Healthcare	17116701
Imidazole	VWR	97064-622
IPTG	Teknova	I3325
Iron(III) chloride hexahydrate	MP Biochemicals	ICN19404590	medium component
LB Agar Miller	Fisher Scientific	BP1425-2
Magnesium sulfate heptahydrate	VWR	97062-134	medium component
Manganese(II) sulfate monohydrate	Acros	AC20590-5000	medium component
MaxQ 6000 Incubated/Refrigerated Shaker	Thermo Scientific	SHKE6000-7	equipment
n-Dodecyl-d25-β-D-maltopyranoside	Anatrace	D310T
n-Dodecyl-β-D-maltopyranoside	Anatrace	D310A
Potassium phosphate monobasic	VWR	97062-346	medium component
RC 6 Plus Centrifuge	Thermo Scientific Sorvall	46910	equipment
SIGMAFAST protease inhibitor cocktail tablets, EDTA-free	Sigma-Aldrich	S8830
Sodium chloride	Sigma-Aldrich	S3014
Sodium hydroxide	Sigma-Aldrich	795429
Sodium phosphate dibasic	Sigma-Aldrich	S7907	medium component
Sterile 25mm syringe filter with 0.2µm PES membrane	VWR	28145-501	labware
Sterile disposable bottle top filter with 0.2µm PES membrane	Thermo Scientific	596-4520	labware
Superdex 200 10/300 GL	GE Healthcare	17517501
Superose-12 10/300 GL column	GE Healthcare	17517301
Ultrospec 10 Cell Density Meter	GE Healthcare	80211630	equipment
Zinc sulfate monohydrate	Acros	AC38980-2500	medium component

References

Wallin, E., Heijne, G. V. Genome-wide analysis of integral membrane proteins from eubacterial, archaean, and eukaryotic organisms. Protein Science. 7 (4), 1029-1038 (1998).
Tautermann, C. S. GPCR structures in drug design, emerging opportunities with new structures. Bioorganic & Medicinal Chemistry Letters. 24 (17), 4073-4079 (2014).
Cournia, Z., et al. Membrane protein structure, function, and dynamics: A perspective from experiments and theory. The Journal of Membrane Biology. 248 (4), 611-640 (2015).
Doerr, A. Membrane protein structures. Nature Methods. 6, 35 (2008).
Bill, R. M., et al. Overcoming barriers to membrane protein structure determination. Nature Biotechnology. 29, 335 (2011).
Carpenter, E. P., Beis, K., Cameron, A. D., Iwata, S. Overcoming the challenges of membrane protein crystallography. Current Opinion in Structural Biology. 18 (5), 581-586 (2008).
Sanders, C. R., Sönnichsen, F. Solution NMR of membrane proteins: Practice and challenges. Magnetic Resonance in Chemistry. 44 (S1), S24-S40 (2006).
Garavito, R. M., Ferguson-Miller, S. Detergents as tools in membrane biochemistry. Journal of Biological Chemistry. 276 (35), 32403-32406 (2001).
Privé, G. G. Detergents for the stabilization and crystallization of membrane proteins. Methods. 41 (4), 388-397 (2007).
Wessel, D., Flügge, U. I. A method for the quantitative recovery of protein in dilute solution in the presence of detergents and lipids. Analytical Biochemistry. 138 (1), 141-143 (1984).
Bayburt, T. H., Grinkova, Y. V., Sligar, S. G. Self-assembly of discoidal phospholipid bilayer nanoparticles with membrane scaffold proteins. Nano Letters. 2 (8), 853-856 (2002).
Bayburt, T. H., Sligar, S. G. Membrane protein assembly into Nanodiscs. FEBS Letters. 584 (9), 1721-1727 (2010).
Skar-Gislinge, N., et al. Elliptical structure of phospholipid bilayer nanodiscs encapsulated by scaffold proteins: Casting the roles of the lipids and the protein. Journal of the American Chemical Society. 132 (39), 13713-13722 (2010).
Sanders, C. R., Schwonek, J. P. Characterization of magnetically orientable bilayers in mixtures of dihexanoylphosphatidylcholine and dimyristoylphosphatidylcholine by solid-state NMR. Biochimie. 31 (37), 8898-8905 (1992).
Vestergaard, M., Kraft, J. F., Vosegaard, T., Thøgersen, L., Schiøtt, B. Bicelles and other membrane mimics: comparison of structure, properties, and dynamics from MD simulations. The Journal of Physical Chemistry B. 119 (52), 15831-15843 (2015).
Popot, J. -. L., et al. Amphipols From A to Z. Annual Review of Biophysics. 40 (1), 379-408 (2011).
Tribet, C., Audebert, R., Popot, J. -. L. Amphipols: Polymers that keep membrane proteins soluble in aqueous solutions. Proceedings of the National Academy of Sciences. 93 (26), 15047 (1996).
Fernández, C., Wüthrich, K. NMR solution structure determination of membrane proteins reconstituted in detergent micelles. FEBS Letters. 555 (1), 144-150 (2003).
Hiruma-Shimizu, K., Shimizu, H., Thompson, G. S., Kalverda, A. P., Patching, S. G. Deuterated detergents for structural and functional studies of membrane proteins: Properties, chemical synthesis and applications. Molecular Membrane Biology. 32 (5-8), 139-155 (2015).
Krueger-Koplin, R. D., et al. An evaluation of detergents for NMR structural studies of membrane proteins. Journal of Biomolecular NMR. 28 (1), 43-57 (2004).
Linke, D., Burgess, R. R., Deutscher, M. P. . Methods in Enzymology. 463, 603-617 (2009).
Oliver, R. C., et al. Tuning micelle dimensions and properties with binary surfactant mixtures. Langmuir. 30 (44), 13353-13361 (2014).
Orwick-Rydmark, M., Arnold, T., Linke, D. The use of detergents to purify membrane proteins. Current Protocols in Protein Science. 84 (1), (2016).
Tanford, C., Reynolds, J. A. Characterization of membrane proteins in detergent solutions. Biochimica et Biophysica Acta (BBA) – Reviews on Biomembranes. 457 (2), 133-170 (1976).
Tulumello, D. V., Deber, C. M. Efficiency of detergents at maintaining membrane protein structures in their biologically relevant forms. Biochimica et Biophysica Acta (BBA) – Biomembranes. 1818 (5), 1351-1358 (2012).
Arachea, B. T., et al. Detergent selection for enhanced extraction of membrane proteins. Protein Expression and Purification. 86 (1), 12-20 (2012).
Seddon, A. M., Curnow, P., Booth, P. J. Membrane proteins, lipids and detergents: Not just a soap opera. Biochimica et Biophysica Acta (BBA) – Biomembranes. 1666 (1), 105-117 (2004).
Tanford, C. . The Hydrophobic Effect: Formation of Micelles and Biological Membranes. , (1980).
Littrell, K., Urban, V., Tiede, D., Thiyagarajan, P. Solution structure of detergent micelles at conditions relevant to membrane protein crystallization. Journal of Applied Crystallography. 33 (3 Part 1), 577-581 (2000).
Oliver, R. C., et al. Dependence of micelle size and shape on detergent alkyl chain length and head group. PLOS ONE. 8 (5), e62488 (2013).
Maire, M., Champeil, P., Møller, J. V. Interaction of membrane proteins and lipids with solubilizing detergents. Biochimica et Biophysica Acta (BBA) – Biomembranes. 1508 (1), 86-111 (2000).
Hong, X., Weng, Y. -. X., Li, M. Determination of the topological shape of integral membrane protein light-harvesting complex LH2 from photosynthetic bacteria in the detergent solution by small-angle X-ray scattering. Biophysical Journal. 86 (2), 1082-1088 (2004).
Vinothkumar, K. R. Membrane protein structures without crystals, by single particle electron cryomicroscopy. Current Opinion in Structural Biology. 33, 103-114 (2015).
Pérez, J., Koutsioubas, A. Memprot: A program to model the detergent corona around a membrane protein based on S-SAXS data. Acta Crystallographica Section D. 71 (1), 86-93 (2015).
Oliver, R. C., Pingali, S. V., Urban, V. S. Designing mixed detergent micelles for uniform neutron contrast. The Journal of Physical Chemistry Letters. 8 (20), 5041-5046 (2017).
Hiruma-Shimizu, K., Kalverda, A. P., Henderson, P. J. F., Homans, S. W., Patching, S. G. Synthesis of uniformly deuterated n-dodecyl-β-D-maltoside (d39-DDM) for solubilization of membrane proteins in TROSY NMR experiments. Journal of Labelled Compounds & Radiopharmaceuticals. 57 (14), 737-743 (2014).
Midtgaard, S. R., et al. Invisible detergents for structure determination of membrane proteins by small-angle neutron scattering. The FEBS Journal. 285 (2), 357-371 (2018).
Gabel, F., Chaudhuri, B., Muñoz, I. G., Qian, S., Urban, V. S. . Biological Small Angle Scattering: Techniques, Strategies and Tips. , 201-214 (2017).
Schantz, A. B., et al. PEE-PEO Block copolymer exchange rate between mixed micelles is detergent and temperature activated. Macromolecules. 50 (6), 2484-2494 (2017).
Liyana-Arachchi, T. P., et al. Bubble bursting as an aerosol generation mechanism during an oil spill in the deep-sea environment: Molecular dynamics simulations of oil alkanes and dispersants in atmospheric air/salt water interfaces. Environmental Science: Processes & Impacts. 16 (1), 53-64 (2014).
Dos Santos Morais, R., et al. Contrast-matched isotropic bicelles: A versatile tool to specifically probe the solution structure of peripheral membrane proteins using SANS. Langmuir. 33 (26), 6572-6580 (2017).
Maric, S., et al. Stealth carriers for low-resolution structure determination of membrane proteins in solution. Acta Crystallographica Section D. 70 (2), 317-328 (2014).
Pedersen, J. S., Svaneborg, C., Almdal, K., Hamley, I. W., Young, R. N. A small-angle neutron and x-ray contrast variation scattering study of the structure of block copolymer micelles: Corona shape and excluded volume interactions. Macromolecules. 36 (2), 416-433 (2003).
Urban, V. S., Kaufmann, E. N. . Characterization of Materials. , (2012).
Chaudhuri, B., Muñoz, I. G., Qian, S., Urban, V. S., Cohen, I. R., et al. . Advances in Experimental Medicine and Biology. 1009, 1-268 (2017).
. Submitting a Research Proposal Available from: https://neutrons.ornl.gov/users/proposals (2018)
. Accompanying Bio-Deuteration Laboratory Proposal Available from: https://www.ornl.gov/sites/default/files/BDL_info_request.docx (2014)
Whitten, A. E., Cai, S., Trewhella, J. MULCh: Modules for the analysis of small-angle neutron contrast variation data from biomolecular assemblies. Journal of Applied Crystallography. 41 (1), 222-226 (2008).
. Scattering Length Density Calculator by National Institute of Standards and Technology (NIST) Center for Neutron Research Available from: https://www.ncnr.nist.gov/resources/activation/ (2018)
Ibel, K., Stuhrmann, H. B. Comparison of neutron and X-ray scattering of dilute myoglobin solutions. Journal of Molecular Biology. 93 (2), 255-265 (1975).
Holme, T., Arvidson, S., Lindholm, B., Pavlu, B. Enzymes: Laboratory-scale production. Process Biochemistry. 5 (9), 62-66 (1970).
Larsson, G., Enfors, S. -. O. Protein release and foaming in Escherichia coli cultures grown in minimal medium. Bioprocess Engineering. 15 (5), 231-237 (1996).
Artero, J. -. B., Hartlein, M., McSweeney, S., Timmins, P. A comparison of refined X-ray structures of hydrogenated and perdeuterated rat [gamma]E-crystallin in H2O and D2O. Acta Crystallographica Section D. 61 (11), 1541-1549 (2005).
Paliy, O., Bloor, D., Brockwell, D., Gilbert, P., Barber, J. Improved methods of cultivation and production of deuteriated proteins from E. coli strains grown on fully deuteriated minimal medium. Journal of applied microbiology. 94 (4), 580-586 (2003).
Sivashanmugam, A., et al. Practical protocols for production of very high yields of recombinant proteins using Escherichia coli. Protein Science. 18 (5), 936-948 (2009).
Hoopes, J. T., Elberson, M. A., Preston, R. J., Reddy, P. T., Kelman, Z., Kelman, Z. . Methods in Enzymology. 565, 27-44 (2015).
Leiting, B., Marsilio, F., O’Connell, J. F. Predictable deuteration of recombinant proteins expressed in Escherichia coli. Analytical Biochemistry. 265 (2), 351-355 (1998).
Perkins, S. J. Estimation of deuteration levels in whole cells and cellular proteins by 1H n.m.r. spectroscopy and neutron scattering. Biochemical Journal. 199 (1), 163-170 (1981).
Obom, K. M., Magno, A., Cummings, P. J. Operation of a benchtop bioreactor. Journal of Visualized Experiments. (79), e50582 (2013).
Duff, A. P., Wilde, K. L., Rekas, A., Lake, V., Holden, P. J., Kelman, Z. . Methods in Enzymology. 565, 3-25 (2015).
Haertlein, M., Kelman, Z., et al. . Methods in Enzymology. 566, 113-157 (2016).
Meilleur, F., Weiss, K. L., Myles, D. A. A. T Micro and Nano Technologies in Bioanalysis. Methods in Molecular Biology. 544, 281-292 (2009).
Naing, S. -. H., Oliver, R. C., Weiss, K. L., Urban, V. S., Lieberman, R. L. Solution structure of an intramembrane aspartyl protease via small angle neutron scattering. Biophysical Journal. 114 (3), 602-608 (2018).
. Training Requirements for First Time and Repeat Users Available from: https://neutrons.ornl.gov/users (2018)
. Remote Analysis Cluster Available from: https://analysis.sns.gov (2018)
Franke, D., et al. ATSAS 2.8: A comprehensive data analysis suite for small-angle scattering from macromolecular solutions. Journal of Applied Crystallography. 50 (4), 1212-1225 (2017).
. Software Suite Available from: https://www.embl-hamburg.de/biosaxs/download.html (2018)
. Software Individual Programs Available from: https://www.embl-hamburg.de/biosaxs/software.html (2018)
Konarev, P. V., Volkov, V. V., Sokolova, A. V., Koch, M. H. J., Svergun, D. I. PRIMUS: A Windows PC-based system for small-angle scattering data analysis. Journal of Applied Crystallography. 36 (5), 1277-1282 (2003).
Jacques, D. A., Trewhella, J. Small-angle scattering for structural biology-Expanding the frontier while avoiding the pitfalls. Protein Science. 19 (4), 642-657 (2010).
Svergun, D. Determination of the regularization parameter in indirect-transform methods using perceptual criteria. Journal of Applied Crystallography. 25 (4), 495-503 (1992).
Putnam, C. D., Hammel, M., Hura, G. L., Tainer, J. A. X-ray solution scattering (SAXS) combined with crystallography and computation: defining accurate macromolecular structures, conformations and assemblies in solution. Quarterly Reviews of Biophysics. 40 (3), 191-285 (2007).
Svergun, D. I., et al. Protein hydration in solution: Experimental observation by x-ray and neutron scattering. Proceedings of the National Academy of Sciences of the United States of America. 95 (5), 2267-2272 (1998).
Franke, D., Svergun, D. I. DAMMIF, a program for rapid ab-initio shape determination in small-angle scattering. Journal of Applied Crystallography. 42 (2), 342-346 (2009).
Svergun, D. I. Restoring Low Resolution Structure of Biological Macromolecules from Solution Scattering Using Simulated Annealing. Biophysical Journal. 76 (6), 2879-2886 (1999).
Jacques, D. A., Guss, J. M., Svergun, D. I., Trewhella, J. Publication guidelines for structural modelling of small-angle scattering data from biomolecules in solution. Acta Crystallographica Section D. 68 (6), 620-626 (2012).
Heller, W. T., et al. The Bio-SANS instrument at the High Flux Isotope Reactor of Oak Ridge National Laboratory. Journal of Applied Crystallography. 47 (4), 1238-1246 (2014).
Svergun, D. I., Koch, M. H. J. Small-angle scattering studies of biological macromolecules in solution. Reports on Progress in Physics. 66 (10), 1735 (2003).
Johnson, J. L., Kalyoncu, S., Lieberman, R. L., Mus-Veteau, I. . Heterologous Expression of Membrane Proteins: Methods and Protocols. , 281-301 (2016).

Play Video

PDF

DOI

DOWNLOAD MATERIALS LIST

Citer Cet Article

Oliver, R. C., Naing, S., Weiss, K. L., Pingali, S. V., Lieberman, R. L., Urban, V. S. Contrast-Matching Detergent in Small-Angle Neutron Scattering Experiments for Membrane Protein Structural Analysis and Ab Initio Modeling. J. Vis. Exp. (140), e57901, doi:10.3791/57901 (2018).