Neutron protein crystallography is a structural technique that permits the localization of hydrogen atoms, thereby providing important mechanistic details of protein function. We present here the workflow for mounting a protein crystal, neutron diffraction data collection, structure refinement and analysis of the neutron scattering length density maps.
Neutron crystallography is a structural technique that allows determination of hydrogen atom positions within biological macromolecules, yielding mechanistically important information about protonation and hydration states while not inducing radiation damage. X-ray diffraction, in contrast, provides only limited information on the position of light atoms and the X-ray beam rapidly induces radiation damage of photosensitive cofactors and metal centers. Presented here is the workflow employed for the IMAGINE and MaNDi beamlines at Oak Ridge National Laboratory (ORNL) to obtain a neutron diffraction structure once a protein crystal of suitable size (> 0.1 mm3) has been grown. We demonstrate mounting of hydrogenated protein crystals in quartz capillaries for neutron diffraction data collection. Also presented is the vapor exchange process of the mounted crystals with D2O-containing buffer to ensure replacement of hydrogen atoms at exchangeable sites with deuterium. The incorporation of deuterium reduces the background arising from the incoherent scattering of hydrogen atoms and prevents density cancellation caused by their negative coherent scattering length. Sample alignment and room temperature data collection strategies are illustrated using quasi-Laue data collection at IMAGINE at the High Flux Isotope Reactor (HFIR). Furthermore, crystal mounting and rapid freezing in liquid nitrogen for cryo-data collection to trap labile reaction intermediates is demonstrated at the MaNDi time-of-flight instrument at the Spallation Neutron Source (SNS). Preparation of the model coordinate and diffraction data files and visualization of the neutron scattering length density (SLD) maps will also be addressed. Structure refinement against neutron data-only or against joint X-ray/neutron data to obtain an all-atom structure of the protein of interest will finally be discussed. The process of determining a neutron structure will be demonstrated using crystals of the lytic polysaccharide monooxygenase Neurospora crassa LPMO9D, a copper-containing metalloprotein involved in the degradation of recalcitrant polysaccharides via oxidative cleavage of the glycosidic bond.
Neutron macromolecular crystallography is a technique that provides a unique window into the structure and underlying chemistry of proteins. Conceptually similar to X-ray diffraction, neutron diffraction provides atomistic details of macromolecular structure, however, the interaction of neutrons with nuclei enables localization of light atoms, often difficult to detect with X-ray diffraction1. During X-ray diffraction, X-rays scatter from the electron cloud, making light atoms such as hydrogen (H) poorly visible in electron density maps that do not have near sub-Ångström resolution2. In contrast, the scattering intensity of neutrons depends on complex interactions with the nucleus, with isotopes of the same element displaying different scattering lengths. Therefore, light atoms and their isotopes, such as hydrogen (1H) and deuterium (2H or D), have comparable visibility to the backbone carbon, nitrogen and oxygen atoms in neutron scattering length density (SLD) maps. Furthermore, since the magnitude of neutron scattering is independent of the number of electrons, scattering from light elements is not obscured by heavy elements when they are in close vicinity to each other, as is observed in X-ray scattering. The enhanced visibility of H and its isotope D when employing neutron diffraction provides valuable information about the protonation state of catalytically important residues, cofactors and ligands and aids the orientation of water molecules, revealing important information about catalytic mechanisms and protein chemistry3. Neutron diffraction also offers the advantage of being a non-destructive technique, particularly suited to biological samples sensitive to ionization such as proteins with metal centers or photosensitive redox cofactors2. The primary focus of this article is to provide an overview of the workflow to obtain a high-quality neutron protein crystal structure. We refer the interested reader to Podjarny et al.4, Blakeley5, Blakeley et al.6 and O‘Dell et al.3 for an excellent overview of neutron protein diffraction and Ashkar et al.7 for further applications of neutron scattering.
Neutrons are primarily generated during nuclear reactions employing either of two processes: nuclear fission at reactor sources or spallation at accelerator-based sources8. Reactor sources provide a continuous neutron beam by employing nuclear fission of the 235U isotope while spallation neutron sources produce a pulsed neutron beam by bombarding a target, for example a liquid metal such as mercury, with protons9. Oak Ridge National Laboratory (ORNL) in Oak Ridge, Tennessee, hosts both a steady-state neutron source at the High Flux Isotope Reactor (HFIR) and a 60 Hz pulsed source at the Spallation Neutron Source (SNS). The IMAGINE beamline, located at the HFIR, is a neutron diffractometer optimized for biological macromolecules (Supplementary Figure 1)10. IMAGINE employs a neutron image plate detector to measure quasi-Laue data using a narrow bandpass in the range of 2.8 – 4.5 Å from single crystals with unit cell edges <150 Å. The Macromolecular Neutron Diffractometer (MaNDi), located at the SNS, is a time-of-flight (TOF) Laue neutron diffractometer equipped with a spherical detector array frame (DAF) (Supplementary Figure 2)11. MaNDi measures data from single crystals with unit cell edges in the range of 10 – 300 Å by employing a tunable 2 Å-wavelength bandwidth between 2.0 – 6.0 Å12.
The process of generating neutrons is highly energy intensive, resulting in relatively weak neutron beam fluxes when contrasted to X-ray beam fluxes at synchrotron sources13. To ensure sufficient signal-to-noise ratios during data collection, it is necessary to grow crystals of suitable size and quality14. Typically, crystals with volumes > 0.1 mm3 are needed to collect data with adequate statistics15. In addition to lower fluxes, inherent properties of the interaction between neutrons and the sample nuclei must be taken into consideration16. The scattering length of neutrons differs for isotopes of the same element, a property which can be advantageously exploited in small angle neutron scattering (SANS) to mask or highlight regions of a sample – a process known as contrast matching17. In diffraction experiments, the negative coherent neutron scattering length of H (-3.741 fm for 1H) can lead to cancellation of neutron scattering density map features since the coherent neutron scattering lengths of other biologically relevant atoms, including carbon (6.6511 fm for 12C), nitrogen (9.37 fm for 14N), oxygen (5.803 fm for 16O), phosphorus (5.13 fm for 31P) and sulfur (2.804 fm for 32S), are positive (Table 1)12,14. Furthermore, the large incoherent scattering length of H (25.274 fm), increases the background during data collection, hampering the quality of the dataset and compromising data resolution7. To circumvent these limitations introduced by H it is necessary, for neutron diffraction, to exchange H for its isotope deuterium, 2H(D), which has a positive coherent neutron scattering length (6.671 fm) and significantly lower incoherent scattering length (4.04 fm)19. This can be achieved by perdeuteration, a process in which the protein is expressed by organisms grown in fully deuterated media ensuring complete incorporation of D at H sites20. It is also possible to partially deuterate protein by replacing H with D solely at the exchangeable sites (titratable groups) while the non-exchangeable carbon-bound sites remain hydrogenated21. This can be achieved by growth of hydrogenated protein crystals in deuterated mother liquor22. Most commonly however, H/D exchange of hydrogenated proteins is performed by vapor exchange following growth of suitably large crystals in H2O-based buffer23. In such cases, crystals are mounted in a quartz capillary and vapor-equilibrated with a D20-based mother liquor.
The limited neutron fluxes at neutron sources result in longer data collection times, ranging from days to several weeks24. At ORNL, both IMAGINE and MaNDi employ a narrow wavelength bandpass in the 2–6 Å range to optimize data collection25. Data can be collected at room temperature or at cryo-temperature. Cryo-data collection can potentially improve data quality and opens up the possibility for freeze-trapping catalytic intermediates. Following neutron diffraction data collection, an X-ray dataset is typically collected on the same crystal at the same temperature or on a crystal grown under identical conditions26. Data collection at the same temperature allows structure refinement to be performed against both X-ray and neutron data, preventing any potential temperature-induced artefacts such as changes in the visibility and position of waters or the occupancies of residues with alternate conformations27. Joint X-ray neutron data refinement increases the data-to-parameter ratio and provides the advantage of allowing the protein backbone coordinates to be refined against the X-ray data, while the neutron diffraction data is used to refine the position of the H/D atoms28. This is particularly useful when using partially deuterated samples, where density cancellation due to H atoms at non-exchangeable sites on the protein is present. Although the number of X-ray structures far exceeds the number of neutron structures deposited in the Protein Data Bank (PDB), software packages initially designed for the refinement of X-ray data have been expanded to encompass neutron data as well3,29,30. Following data collection, models can be refined using refinement packages such as phenix.refine, CNSsolve (nCNS) or SHELXL28,31,32,33. During the refinement process, neutron scattering density maps can be visualized for manual fitting using COOT34. Following structure solution, the coordinates and the neutron and/or X-ray diffraction data files can be submitted to the PDB, who will validate and deposit the model, making it available for public access18,29,30.
Structural analysis of proteins is a multifaceted approach in which numerous techniques are used to probe their function and mechanism35. Neutron protein crystallography provides valuable chemical insights to expand on and complement findings from additional studies such as X-ray diffraction, spectroscopy, nuclear magnetic resonance (NMR) or micro crystal electron diffraction (microED)36. Neutron protein diffraction is uniquely positioned to provide insights into enzymatic mechanisms, since H atoms are central to their chemistry. The absence of radiation damage induced by neutrons make them a probe exceptionally suited to the study of metalloproteins37. We present here a representative example of the process of neutron protein diffraction from sample preparation to data collection, refinement and analysis (Figure 1). Crystals of sufficient size for neutron diffraction experiments have been grown of the metalloprotein Neurospora crassa LPMO9D (NcLPMO9D). NcLPMO9D is a copper-containing metalloprotein involved in the degradation of recalcitrant cellulose by oxygen atom insertion at the glycosidic bond38,39. The NcLPMO9D active site contains a mononuclear copper center within a characteristic “histidine-brace” composed of the N-terminal histidine and a second conserved histidine (Supplementary Figure 3)40. The N-terminal of fungal LPMOs is methylated but the post-transitional modification does not occur during recombinant expression in yeast. In the NcLPMO9D resting state, the copper center is present in a Cu2+ oxidation state and is activated by single electron reduction to Cu1+, allowing molecular oxygen to bind and be activated by rapidly being reduced to a superoxide species41,42. The overall NcLPMO9D reaction requires further addition of one electron and two protons to form the hydroxylated polysaccharide product43. The identity of the activated oxygen species responsible for hydrogen atom abstraction (HAA) from the polysaccharide substrate has not been identified and intensive structural and computational studies are currently ongoing44,45. Given the redox chemistry at the NcLPMO9D active site, mitigation of radiation damage is particularly pertinent. We illustrate here room temperature and cryo-temperature data collection on NcLPMO9D crystals to determine the NcLPMO9D structure in the resting state and in the activated reduced form, respectively46. Emphasis will be given to protein crystal mounting, beamline instrument setup for data collection, the preparation of the data and coordinate files and the refinement steps necessary to model an all-atom neutron structure.
1. Crystal size evaluation
2. Preparation of deuterated crystallization buffer
3. Crystal harvesting
4. Crystal mounting
NOTE: Capillary mounting protocols vary with experimentalist preferences. To prevent damage to crystals, capillaries that need to be shortened should be scored with a cutting stone or sandpaper to ensure a smooth break.
5. Vapor exchange
6. Neutron protein diffraction
NOTE: Readers interested in the IMAGINE beam line specifics are encouraged to consult Meilleur et al. 2013, Meilleur et al. 201810,47.
7. Structure refinement
Neutron diffraction data on crystals of a lytic polysaccharide monooxygenase from Neurospora crassa (NcLPMO9D) were collected on IMAGINE at the HFIR at room temperature and on MaNDi at the SNS under cryo-conditions following the protocol described above. Crystals of the hydrogenated protein grown in H2O-based buffer with volume greater than 0.1 mm3 were used (illustrative example of large crystals are shown in Supplementary Figure 4 and figures thereafter). Crystals were mounted in quartz capillaries and vapor exchange with the D2O-based buffer was performed for three weeks prior to data collection (Figure 4).
Room temperature data collection was performed on the IMAGINE beamline (Figure 1). A four hour white-beam test led to high resolution diffraction suggesting that the crystal was of suitable size and quality for a full dataset to be collected. In addition to providing preliminary information on the diffraction quality of the crystal, the initial broad bandpass exposure can be used to index the diffraction pattern and determine the crystal orientation matrix. Given the P21 space group of the crystal, a data collection strategy of 18 frames with a collection time of 20 hours per frame was implemented. As with X-ray diffraction data collection, higher symmetry space groups require fewer frames (i.e. less angular coverage) to collect a complete dataset. The data were collected in quasi-Laue mode using a wavelength range of 2.8 – 4.0 Å. Following data collection, the data were indexed, integrated scaled and merged to give a neutron SLD file in MTZ format at a resolution of 2.14 Å. Data were evaluated to be of sufficient quality following similar guidelines for X-ray data analysis, although a completeness of 80 % and a CC1/2 of at least 0.3 were considered acceptable since neutron protein diffraction is a flux-limited technique.
Following room-temperature neutron diffraction data collection, the same crystal was used to collect a room temperature X-ray diffraction dataset at 1.90 Å resolution (Supplementary Figure 13). The X-ray data were used to determine the positions of the “heavier” atoms including C,N,O and S. The structure refined against the X-ray data alone was then used as the starting model to perform a joint refinement against the X-ray and neutron data. Phenix ReadySet was used to add H atoms at non-exchangeable sites, H and D atoms at exchangeable sites and D atoms to water molecules of the starting X-ray model. Following this model preparation, iterative refinements were performed against both datasets (Supplementary Figure 19 and Supplementary Figure 20). Interactive model building was performed in Coot by visually inspecting the density maps to orientate side-chains and water molecules accordingly (Supplementary Figure 22). The neutron data were primarily used to determine protonation states and water molecule orientations. Comparison of the electron density map of residues such as serine and tryptophan and the corresponding neutron SLD map illustrate the information that can be gained on protonation states at H/D exchangeable sites from neutron protein diffraction (Figure 7). A map overlay of electron and neutron SLD maps for water molecules also indicate that while hydrogen bond interactions can be inferred from X-ray data, neutrons provide clear information regarding the orientation of these hydrogen bonds (Figure 8). Neutron SLD FO-FC omit maps were generated to determine protonation states and H/D orientation of side-chains. Illustrated are the neutron SLD maps obtained for tyrosine and threonine residues, in which the neutron Fo-FC maps clearly indicate positive peaks signifying the presence of H/D (Figure 9). The collected neutron diffraction data also provided valuable information about multiple protonation states, such as the -ND3+ group of Lys (Figure 10). Refinement statistics (Rwork and Rfree) were closely monitored during model optimization to prevent over-fitting. Final statistics gave an X-ray Rwork of 12.77 % and an Rfree of 18.21%, and a neutron Rwork of 14.48% and an Rfree of 21.41% with 389 water molecules present (Supplementary Figure 28).
Cryo-temperature data were collected on NcLPMO9D following an ascorbate soak to reduce the copper active site from CuII to CuI on the MaNDi beamline (Supplementary Figure 2 and Supplementary Figure 15)45. Data were collected using TOF Laue mode following a neutron diffraction test using a 4 hour exposure to verify the quality of diffraction. Given the space group of the crystal, a data collection strategy of 18 frames with a collection dose of 80 Coulombs per frame was devised. The data were collected in TOF-Laue mode at a wavelength range of 2.0 – 4.0 Å. Following data collection, the data were indexed, integrated, scaled and merged to give a reflection file in MTZ format at a resolution of 2.40 Å51,52.
Following data collection, the 2.40 Å cryo-temperature NcLPMO9D neutron diffraction dataset was used for neutron-only data refinement. The neutron data were phased by molecular replacement using PDB 5TKH as the starting model. Phenix ReadySet was used to add H atoms at non-exchangeable sites and H/D atoms with partial occupancies at exchangeable sites. Water molecules were removed from the starting model with PDB Tools (Supplementary Figure 23). Model preparation was followed by refinement with phenix.refine using the neutron scattering table (Supplementary Figure 24). Interactive model building was performed in Coot, with water molecules being added using the positive peaks of the FO-Fc map and positioned according to potential hydrogen bond interactions (Figure 11A and Figure 11B). When analyzing neutron SLD maps, water molecules are clearly visible if they are highly ordered, however their density may be spherical or ellipsoidal if they are not well-ordered (Figure 11C-E). Neutron SLD maps were used to provide valuable information on the orientation of residues such as asparagine, in which differentiating between the carbonyl and amino groups can be challenging when using X-ray diffraction data alone (Figure 12A and Figure 12B). Peaks in FO-FC neutron SLD omit maps were also very informative in determining the protonation states of histidine residues at the Nδ– or Nε-position (Figure 12C and Figure 12D). The protonation state of residues with multiple H/D exchangeable sites can also be determined using neutron SLD maps. This was clearly illustrated with an FO-FC neutron SLD omit map of arginine, which is known to have a positive charge (Figure 12E and Figure 12F). As previously, over-fitting was prevented by monitoring Rwork and Rfree. Final statistics gave an Rwork of 22.58% and an Rfree of 30.84% (Supplementary Figure 29). Given that neutron protein diffraction is a flux limited technique in which the negative scattering length and large incoherent scattering factor of H must be taken in to account, it can be expected that a neutron data-only refinement would have poorer statistics than a joint X-ray/neutron-data refinement with fewer visible water molecules (Supplementary Figure 28 and Supplementary Figure 29).
When analyzing neutron SLD maps, it will become apparent that density cancellation due to the negative neutron scattering length of H will occur for hydrogenated proteins that were subjected to vapor exchange with D2O-containing crystallization buffer. Due to this reason, neutron SLD maps in which non-exchangeable H atoms are attached to carbon appear incomplete when compared to their electron density map counterpart (Figure 13A). The effect of cancellation is often more apparent at poorer resolutions, making it imperative to obtain protein crystals of a high quality. It is therefore preferable to perform a joint refinement of a sample with both X-ray and neutron data in which the X-ray data can be used to determine the position of the protein backbone (Figure 13B). Furthermore, sulfur atoms in cysteine and methionine may be poorly visible, requiring X-ray data for exact atom placement (Figure 13C and Figure 13D). Metals with weak neutron scattering lengths may also be challenging to model in neutron SLD maps, as is apparent in our LPMO9D maps. Collection of a low dose (free of radiation damage) X-ray dataset on the same crystal is therefore useful, since it permits metal atom positioning using electron density maps (Figure 13E and Figure 13F).
Figure 1: Flow chart of neutron protein crystallography workflow. Protein Production. In order to obtain a neutron structure, protein is first expressed. Bacterial expression in H2O- or D2O-based media is typically used to produce a high yield of hydrogenated or perdeuterated recombinant protein, respectively. The protein is purified in H2O-based buffer and then crystallized in either H2O- or D2O-based crystallization buffer to grow crystals to a minimum size of 0.1 mm3. Sample Preparation: Prior to neutron diffraction data collection, H2O-grown crystals undergo H/D exchange to exchange the protein titratable H atoms with D. H/D exchange can be done by direct soaking of the crystals in deuterated crystallization buffer, equilibration of the crystallization drop with a D2O-based reservoir, or by mounting the crystals in quartz capillaries for vapor exchange with deuterated crystallization buffer. Neutron Data Collection: Following H/D exchange, potential crystals are screened to determine the diffraction quality. Crystals with a minimum resolution of 2.5 Å are considered suitable for a full dataset to be collected. Crystals are mounted in quartz capillaries for data collection at room temperature or flash frozen in a cryo-loop for data collection at cryogenic temperature. An X-ray dataset is collected on the same (or an identical) crystal at the same temperature. Model Building: Refinement is performed using phenix.refine against both neutron and X-ray data or against the neutron data only. Manual model building of the protein structure is performed in Coot using the neutron SLD maps. Complete Structure: Following completion of the protein structure, the coordinate model is validated and deposited in the Protein Data Bank. Please click here to view a larger version of this figure.
Figure 2: Harvesting protein crystals. (A) Crystals are handled under a microscope. (B) The sealed sandwich box containing the siliconized glass plate is opened. Reservoir buffer is pipetted onto siliconized glass slides. (C) A crystal is harvested with a microloop. (D) The crystal is placed in a drop of mother liquor to wash off any debris that are often harvested along with the crystal. Please click here to view a larger version of this figure.
Figure 3: Transfer of crystal to quartz capillary. (A) The end of a quartz capillary is filled with reservoir buffer. (B) The crystal is transferred into the quartz capillary and (C) immersed in reservoir buffer. (D) The crystal is carried down capillary using reservoir buffer. Please click here to view a larger version of this figure.
Figure 4: Sealing of the quartz capillary. (A) Deuterated buffer is added at the end of the capillary to form a “plug”. (B) Wax is melted with a “wand”. (C) The capillary is placed in the melted wax to seal. (D) Wax plugs are formed on both ends to seal the capillary. (E) The crystal after mounting. (F) The sealed capillary is placed in a Petri dish and held in place with putty. Please click here to view a larger version of this figure.
Figure 5: Increased signal-to-noise of the neutron diffraction pattern. As data collection proceeds, diffracted spots become more intense. (NOTE: the live diffraction images presented here are for illustration and were taken from different crystals.) Please click here to view a larger version of this figure.
Figure 6: Interactive model building using neutron data in Coot. (A) A positive FO-FC neutron SLD density peak (green) indicating serine must be reoriented by editing chi angles. The 2FO-FC neutron SLD map is displayed in purple and 2FO-FC electron density map is displayed in blue. (B) Correctly positioned serine. (C) Positive and negative FO-FC neutron SLD density peaks (green and red, respectively) indicating that tryptophan must be rotated/translated to match difference density peak. (D) Correctly oriented tryptophan. Please click here to view a larger version of this figure.
Figure 7: Additional information from neutron SLD maps. (A) 2FO-FC electron density map (blue) displays the positions of the “heavier” atoms in serine. (B) 2FO-FC neutron SLD map (purple) clearly displays the position of the “lighter” D atom in serine. (C) 2FO-FC electron density map (blue) displays the positions of the “heavier” atoms in tryptophan. (D) 2FO-FC neutron SLD map (purple) clearly displays the position of the “lighter” D atom in tryptophan. Please click here to view a larger version of this figure.
Figure 8: Water molecule positioning. (A) The spherical shape of an 2FO-FC electron density map (blue) feature for water. (B) The 2FO-FC neutron SLD map (purple) provides information about the water orientation and hydrogen bond interaction. (C) Map overlay of electron and neutron SLD maps of water. The 2FO-FC neutron SLD map is displayed in purple and 2FO-FC electron density map is displayed in blue. Please click here to view a larger version of this figure.
Figure 9: Neutron SLD FO-FC omit maps. (A) The FO-FC neutron SLD map (green) provides clear information on the H/D orientation of tyrosine residues. The 2FO-FC neutron SLD map is displayed in purple and 2FO-FC electron density map is displayed in blue. (B) Tyrosine residue with correct H/D orientation. (C) FO-FC neutron SLD map (green) provides clear information on the H/D orientation of threonine residues. (D) Threonine residue with correct H/D orientation. Please click here to view a larger version of this figure.
Figure 10: Multiple protonation states displayed with neutron SLD maps. (A) The 2FO-FC electron density map (blue) only provides the position of the N atom of lysine ε-ammonium group. (B-E) The FO-FC neutron SLD omit map (green) clearly demonstrates the positively charged -NH3 group. The 2FO-FC neutron SLD map is displayed in purple and 2FO-FC electron density map is displayed in blue. (F) Overlay of electron density and neutron SLD maps. Please click here to view a larger version of this figure.
Figure 11: Appearance of water molecules in neutron SLD maps. (A) Water molecules are positioned according to FO-FC neutron SLD maps (green) and potential hydrogen bonds. The 2FO-FC neutron SLD map is displayed in purple. (B) Correctly positioned water molecule. (C-E) The various shapes of neutron SLD maps for water molecules depending on B-factors and hydrogen bond interactions. Please click here to view a larger version of this figure.
Figure 12: Information about amino acid orientation and protonation provided by neutron SLD maps. (A) The neutron SLD FO-FC map peaks (green) indicate incorrect orientation of an asparagine residue. The 2FO-FC neutron SLD map is displayed in purple and 2FO-FC electron density map is displayed in blue. (B) 2FO-FC neutron SLD map (purple) of the correct asparagine orientation. (C) The neutron SLD FO-FC map peak (green) indicates single protonation of the histidine at Nε. (D) 2FO-FC neutron SLD map (purple) of histidine Nε -protonation. (E) Neutron SLD FO-FC omit map peaks (green) confirm the positive charge of arginine. (F) 2FO-FC neutron SLD map (purple) of positively charged arginine. Please click here to view a larger version of this figure.
Figure 13: Discontinuous neutron SLD maps. (A) 2FO-FC neutron SLD map (purple) of a hydrogenated, vapor H/D exchanged protein. Glutamic acid displays neutron SLD map cancellation due to the negative scattering length of non-exchangeable H atoms. (B) An overlayed 2FO-FC electron density map (blue) clearly displays the density of the glutamic acid. (C) The sulfur atom in methionine is poorly visible in 2FO-FC neutron SLD maps (purple). (D) An overlayed electron density map clearly displays the density of the methionine. (E) Metal atoms, here copper, are poorly visible in neutron 2FO-FC SLD maps (purple). (F) An overlayed 2FO-FC electron density map (blue) clearly displays the density of the coordinated copper atom. Please click here to view a larger version of this figure.
Isotope | Coherent scattering length (fm) | Incoherent scattering length (fm) |
1H | -3.741 | 25.274 |
2H | 6.671 | 4.04 |
12C | 6.6511 | 0 |
14N | 9.37 | 2 |
16O | 5.803 | 0 |
23Na | 3.63 | 3.59 |
24Mg | 5.66 | 0 |
31P | 5.13 | 0.2 |
32S | 2.804 | 0 |
35Cl | 11.65 | 6.1 |
39K | 3.74 | 1.4 |
40Ca | 4.8 | 0 |
55Mn | -3.73 | 1.79 |
56Fe | 9.94 | 0 |
63Cu | 6.43 | 0.22 |
64Zn | 5.22 | 0 |
Table 1: Neutron scattering lengths and incoherent scattering values. Adapted from Sears, 199216.
Supplementary Figure 1: The IMAGINE Instrument at the High Flux Isotope Reactor. (A) The IMAGINE instrument in the cold neutron guide hall. (B) Sample in mounted in a quartz capillary attached with putty to the goniometer. The sample and detector table closes to position the crystal and the cylindrical image plate in the neutron beam. Modified with permission of the International Union of Crystallography53. Images provided with permission of Genevieve Martin, Oak Ridge National Laboratory. Please click here to download this figure.
Supplementary Figure 2: The MaNDi Instrument at the Spallation Neutron Source. (A) The MaNDi Anger camera detector array. Reproduced with the permission of International Union of Crystallography11. (B) MaNDi moveable sample stage. (C) Sample mounted in quartz capillary mounted on the goniometer at MaNDi for room temperature data collection. Images provided with permission of Genevieve Martin, Oak Ridge National Laboratory. Please click here to download this figure.
Supplementary Figure 3: Structure of the lytic polysaccharide monooxygenase NcLPMO9D. The NcLPMO9D copper-active site is located on a flat polysaccharide binding surface. The copper is coordinated by two histidine residues in a classical “histidine brace” as well as an axial tyrosine residue. Please click here to download this figure.
Supplementary Figure 4: Crystal with sufficient volume in sitting drop crystallization tray. (A) Large crystals are grown in sitting drops set up in 9-well siliconized glass plates. (B and C) Crystals are measured to identify those with volume > 0.1 mm3. Please click here to download this figure.
Supplementary Figure 5: pH meter set up for deuterated buffer readings. The pH electrode is soaked in D2O prior to use. NaOD and DCl are used to adjust the pH of deuterated buffers. Please click here to download this figure.
Supplementary Figure 6: MaNDi sample mounting guidelines. Maximum dimensions of the quartz capillary and sample position for room temperature data collection.
Reproduced from: https://neutrons.ornl.gov/mandi/sample-environment Please click here to download this figure.
Supplementary Figure 7: Removal of excess buffer. (A) Excess buffer is aspirated from the quartz capillary with microcapillary tips. (B) The remaining buffer is removed with a thin paper wick to completely dry the capillary. Please click here to download this figure.
Supplementary Figure 8: The data acquisition GUI. Input window of the “Experiment Parameters” for data collection. Please click here to download this figure.
Supplementary Figure 9: The Optics GUI. Selection of the quasi-Laue range for data collection and monitoring of the neutron count rate. Please click here to download this figure.
Supplementary Figure 10: Data collection in the data acquisition GUI. The exposure time, number of frames and angles for data collection are specified in the “Collect” tab. Data collection is then initiated using “Start Scan”. Please click here to download this figure.
Supplementary Figure 11: Diffracted neutrons detected and displayed. At the end of the exposure time, the neutron sensitive image plate detector is read out and the diffraction pattern is displayed in the data acquisition GUI. Please click here to download this figure.
Supplementary Figure 12: Data processing following neutron diffraction. Frames are indexed, integrated, wavelength normalized and scaled using Lauegen, Lscale and Scala to generate a merged reflection file following data collection. Please click here to download this figure.
Supplementary Figure 13: X-ray data collection. Home source X-ray generator set up with quartz capillary mounted crystal for room temperature data collection. Please click here to download this figure.
Supplementary Figure 14: Mounting guidelines for MaNDi cryo-data collection. Dimensions of CrystalCaps and pin height for cryo-data collection at MaNDi.
Reproduced from: https://neutrons.ornl.gov/mandi/sample-environment Please click here to download this figure.
Supplementary Figure 15: Flash freezing for cryo neutron diffraction data collection. (A) Setup for crystal soaking, harvesting with a microloop and freezing in liquid nitrogen using a cryo compatible container such as a foam Dewar. The mounted crystal is transferred directly onto the beamline cryo goniometer using precooled cryo pin tongs. (B) The wax seal is melted for crystal removal. (C) The crystal is flushed to the end of the quartz capillary for harvesting. (D) The crystal is sequentially soaked in ascorbate soak buffer and then cryoprotectant followed by flash freezing in liquid nitrogen. Please click here to download this figure.
Supplementary Figure 16: Sample alignment interface. Crystal alignment in the neutron beam, represented by the blue cross, is done by point and click centering. Please click here to download this figure.
Supplementary Figure 17: The CSS GUI for data collection. The data collection strategy, including exposure doses and angles, are uploaded in the CSS GUI. As data collection proceeds the diffracted neutrons detected on the real-time detector will be displayed in the upper panel. Please click here to download this figure.
Supplementary Figure 18: Matching R-free flags in CCP4. The R-free flags of the neutron data are matched with the R-free flags of X-ray data collected on the same or an identical crystal for joint refinement. Please click here to download this figure.
Supplementary Figure 19: Structure preparation and refinement. (A) The Phenix ReadySet tool is used to add dual H/D occupancy at exchangeable sites. (B) Both the neutron data and X-ray data are used for a joint refinement, while the initial input model was refined against the X-ray dataset collected on the same crystal or an identical crystal. Please click here to download this figure.
Supplementary Figure 20: Configuration of refinement settings. The refinement model as well as the nuclear distances are configured for joint X-ray/neutron data refinement. Please click here to download this figure.
Supplementary Figure 21: Data selection for Coot model building. The phenix MTZ file output containing X-ray and unfilled neutron data is opened in Coot to generate electron and neutron SLD maps for interactive model building. Please click here to download this figure.
Supplementary Figure 22: Interactive model building in Coot during a joint refinement. (A) A positive and negative FO-FC neutron SLD density peak (green and red, respectively) indicating that the water must be reoriented by rotation/translation. The 2FO-FC neutron SLD map is displayed in purple and 2FO-FC electron density map is displayed in blue. (B) Correctly positioned water. (C) A positive FO-FC neutron SLD map peak (green) indicate that threonine must be rotated to match difference density peak by editing chi angles. (D) Correctly oriented threonine. Please click here to download this figure.
Supplementary Figure 23: Structure preparation for neutron-only data refinement. The starting coordinate file is prepared for refinement by water atom removal in PDBTools and by addition of dual H/D occupancy at exchangeable sites. Please click here to download this figure.
Supplementary Figure 24: Neutron data-only refinement. (A) Neutron data is uploaded as well as the prepared starting model. (B) The settings for neutron data refinement use the neutron scattering table. Please click here to download this figure.
Supplementary Figure 25: Data selection for Coot model building. The unfilled neutron data are opened in Coot for interactive model building. Please click here to download this figure.
Supplementary Figure 26: Real space refinement in Coot for deuterated residues. (A) Positive and negative FO-FC neutron SLD density peaks (green and red, respectively) indicating that an arginine residue must be moved to fit the FO-FC density peak. The 2FO-FC neutron SLD map is displayed in purple and 2FO-FC electron density map is displayed in blue. (B) Utilizing Real Space Refine results in “exploding” D atoms due to missing Coot geometry restraint libraries. (C) The D atoms do not move with the rest of the residue atoms. (D) The D atom positions can be manually fixed using a text editor. Please click here to download this figure.
Supplementary Figure 27: Addition of water molecules. (A) Water molecules can be manually added to the positive FO-FC neutron SLD map density peaks (green). The inserted water molecules will initially be represented by an O atom in Coot. (B) Phenix ReadySet is used to add D atoms to the O atoms for water molecules. (C) The deuterated water molecule is successfully added. Please click here to download this figure.
Supplementary Figure 28: Refinement statistics. Final data refinement statistics following joint X-ray/neutron refinement. Please click here to download this figure.
Supplementary Figure 29: Refinement statistics. Final data refinement statistics following neutron data-only refinement. Please click here to download this figure.
Neutron protein crystallography is a highly sensitive technique to probe protonation states and water molecule orientation in proteins. This information sheds light on protein catalytic mechanisms since changes in protonation and hydrogen bonding interactions are often central to enzyme chemistry10. Neutron protein crystallography, albeit an informative technique, has a number of factors that should be taken into consideration before planning to conduct a neutron diffraction experiment, namely:
Neutron protein crystallography is a flux limited technique. In contrast to X-ray diffraction datasets, higher R-factors and lower completeness, redundancy and signal-to-noise ratios are expected for neutron datasets due to the technique inherent limitations (flux limited, quasi-Laue, longer wavelengths). Data collection of a single frame is typically 12 – 18 hours. Success of an experiment is highly dependent on sample size and quality with crystals of 0.1 mm3 often being the minimum requirement3. Neutron diffraction requires production of large amounts of protein to set up crystallization drops ranging from 10 to 800 μL. The minimum volume for growing sufficiently large crystals can be estimated using a Volume Calculator given the crystal and sample parameters (https://neutrons.ornl.gov/imagine). Growth of large crystals has most prevalently been accomplished by vapor diffusion3. Hanging drop crystallization permits growth of crystals in large drops ranging from 10-25 μL, while larger drops ranging up to ~ 50 μL can be set up using commercially available sitting drop equipment14,54. Siliconized nine-well glass plates can be used to set up very large drops, with volumes up to 800 μL. These glass plates are placed in “sandwich boxes” commercially available from Hampton Research. Further crystallization techniques include batch crystallization, in which the limit of the drop size is dictated by the vessel. Batch crystallization experiment set up can range from microliters to milliliters55. Crystallization can also be performed using the dialysis technique in which the protein is equilibrated with the precipitant via a dialysis membrane or by counter-diffusion along a precipitant concentration gradient or through a porous plug such as agarose56,57. Seeding offers another alternative to obtain crystals of the desired volume. Micro- and macroseeding have been successfully employed for large crystal growth, including large crystal of NcLPMO9D45. Some knowledge of the protein phase diagram, including the influence of temperature on solubility, aid in large crystal growth.
When planning a neutron diffraction experiment, optimization of the protein preparation to maximize signal-to-noise ratio during diffraction data collection is essential7. To circumvent density cancellation and high incoherent scattering caused by H atoms, neutron SLD maps can be improved by exchanging H atoms for its isotope D, which possesses a positive coherent scattering length and low incoherent scattering length. To accomplish this, vapor exchange of the hydrogenated protein crystal against deuterated crystallization buffer is performed. This ensures H/D exchange of solvent molecules and the labile, titratable protein H atoms23. Vapor exchange is performed by mounting the hydrogenated crystal in a quartz capillary with D2O-based, deuterated crystallization buffer “plugs” and it represents an effective, gentle technique that is most often applied14,23,35. The exchange can take several weeks and preferably requires the deuterated buffer to be frequently changed to ensure maximum H/D exchange. H/D exchange can also be performed by directly soaking the crystal in deuterated buffer. To avoid placing the crystal under stress due to D2O exposure, the soaking process should be performed gradually by incrementally increasing the D2O:H2O ratio58. In addition to this, crystallization of hydrogenated protein can also be performed in deuterated buffer for H/D exchange at labile H sites22,59. It should be noted, however, that D2O-based buffer has an effect on protein solubility requiring further adjustment of the known H2O-based conditions3,59. D2O-based buffers have also been observed to lead to smaller crystals in some cases59. Full exchange of titratable and carbon-bound H atoms to D can be achieved by expressing proteins in deuterated media to generate a perdeuterated sample20. The resulting neutron SLD maps of the perdeuterated sample will be significantly improved, no longer displaying the density cancellation of the hydrogenated sample counterpart. This is beneficial when characterizing H/D bound at non-exchangeable sites in a protein or cofactor. However, expression of perdeuterated protein is both high in cost and low in yield60. The Oak Ridge National Laboratory (ORNL) Center for Structural Molecular Biology (CSMB) offers a deuteration facility for users seeking to generate a perdeuterated sample (https://www.ornl.gov/facility/csmb). Perdeuterated expression is typically performed in a bioreactor on the 1 L scale yielding ~50 mg of purified protein61.
Following the collection of neutron diffraction data, refinement and interactive model building is performed. Refinement can be run using multiple software suites including phenix.refine, nCNS or SHELXL28,31,32,33. The Phenix suite is the most commonly utilized software for refinement of neutron diffraction data in conjunction with Coot which is used to manually build the model from the neutron SLD maps34. Although both Phenix and Coot allow for the processing of neutron diffraction data, they may lack certain features necessary to process the idiosyncrasies associated with neutron data and deuterated samples. For example, Coot does not contain geometry optimization for deuterated residues, which can lead to complications during model building since the “Real Space Refine” feature results in “exploding” residues (Supplementary Figure 26)62. This can be resolved by generating restraint files for all the deuterated residues. However, this is an intensive process and such libraries are not currently publicly available. When performing refinements in Phenix, exchangeable H/D sites will initially be set to 0.50 occupancy for H and D. As refinements are performed, the occupancy of H and D will be refined according to the neutron SLD maps. During interactive model building, difference density Fo-Fc maps are very informative in assessing H/D occupancies. Maps can be used to determine which sites possess high D occupancy, which is particularly informative at the active site where protonation states are catalytically relevant63. Ambiguous situations do arise, however, when the H:D occupancy is close to 0.70:0.30 which results in complete signal cancellation in neutron SLD maps64. It should also be taken into account that quasi-Laue neutron data sets often have completeness around 80%, which is lower than the routinely observed ≥ 98% for X-ray diffraction data. When refining neutron diffraction data in Phenix, the missing observed amplitudes (Fo) are therefore calculated from the model to complete the reflection list, thus introducing model bias. To account for this potential bias “no_fill” maps should be examined during interactive model building as opposed to “filled” maps.
Users can choose to perform a joint X-ray/neutron data refinement of their structure, or a neutron-data only refinement. Visualizing neutron SLD maps, particularly at lower resolution, may initially be disconcerting especially for a hydrogenated protein in which H is still present at non-exchangeable sites despite H/D vapor exchange. This results in neutron density map cancellations, giving the impression of discontinuous maps65,66. Collecting a corresponding X-ray dataset advantageously complements these cancellations in a joint refinement (Figure 13A and Figure 13B). A joint-refinement strategy typically involves refining the protein backbone coordinates against the X-ray data, while the neutron diffraction data is used to refine the position and occupancy of the H/D atoms at exchangeable sites28. Since introduction of joint H/D occupancy at exchangeable sites increases the number of parameters being refined, a joint refinement with X-ray data also increases the data-to-parameter ratio. A joint refinement requires a corresponding X-ray dataset to be collected at the same temperature on the same crystal or a crystal grown under the same conditions. For neutron diffraction data collected at room temperature (300K), the corresponding X-ray dataset should be collected at room temperature using a low-dose data collection strategy to limit radiation damage. Perdeuterated samples, in contrast, provide improved and continuous neutron SLD maps since they do not possess the same magnitude of H/D signal cancellation. However, the neutron scattering length of certain elements including metals and sulfur make them poorly visible in neutron SLD maps, even if the protein has been perdeuterated (Figure 13C-F)18. If a metal needs to be characterized, it is best to utilize X-ray diffraction in a joint refinement or apply spectroscopic techniques to complement diffraction experiments. Neutron-only data refinements are often performed when the neutron dataset has high resolution or if a perdeuterated protein was used. In addition, neutron-only data refinement is particularly useful if a protein highly sensitive to radiation damage is being studied, since an X-ray derived structure may possess radiation-induced artefacts. If a neutron-data-only refinement is to be performed, it must be ascertained whether the corresponding neutron dataset has sufficient completeness and resolution.
ORNL offers two facilities for collection of neutron diffraction data: the IMAGINE beamline at the HFIR as well as the MaNDi beamline at the SNS36,67. While both instruments provide effective means for collecting a neutron diffraction dataset employing similar principles, each instrument has unique specifications that should be taken into account when applying for beam time. IMAGINE collects quasi-Laue data and is optimized for room temperature data collection on crystals with unit cells up to ~100 Å. MaNDi can be used for the collection of room temperature and cryo-temperature data employing TOF-Laue collection on crystals with unit cells up to ~ 300 Å. Prior to collecting a complete dataset, a test is performed on the crystal to evaluate the quality of the obtained diffraction pattern in which the crystal is exposed to the neutron beam for a single frame. If the crystal is of sufficient quality, a full neutron diffraction dataset will be collected, indexed, integrated, scaled and merged in a process analogous to X-ray data processing. IMAGINE makes use of Lauegen and Lscale and MaNDi utilizes the Mantid package and employs three-dimensional profile fitting48,50,51,68,69,70. Scientists who become users at either of these facilities will be provided with a dataset in MTZ or HKL format for further analysis.
Neutron diffraction is a non-destructive, highly sensitive technique for probing the protonation state and hydrogen bond interactions of biological macromolecules. It is particularly useful for photo-sensitive proteins and metalloproteins. Several considerations regarding the technique as well as the processing of the data must be taken into account before conducting an experiment, however the outcome yields results which may give valuable insight into the catalytic mechanism of the protein of interest. Neutron protein crystallography complements computational, structural, biochemical and spectroscopic studies, making it a valuable tool in the biologist’s toolbox of techniques used to characterize biological macromolecules.
The authors have nothing to disclose.
Protein expression, purification and crystallization experiments were conducted at the Center for Structural Molecular Biology (CSMB), a U.S. Department of Energy Biological and Environmental Research User Facility at Oak Ridge National Laboratory. Neutron diffraction data was collected at BL-11B MaNDi at the Spallation Neutron Source (SNS) at ORNL which is sponsored by the Scientific User Facilities Division, Office of Basic Energy Sciences, U.S. Department of Energy. The authors thank Brendan Sullivan for assistance with data reduction. X-ray diffraction data was collected at the Molecular Education, Technology, and Research Innovation Center (METRIC) facilities at North Carolina State University, which is supported by the State of North Carolina. GCS acknowledges support in part from the National Research Foundation (NRF), South Africa and the Graduate Opportunities (GO!) program at ORNL. FM acknowledges support from USDA NIFA Hatch 211001.
Absorbent Paper Points Size #30-#40, 60 mm length | DiaDent/DiaVet | 218-292 | |
Capillary wax | Hampton | HR4-328 | |
CCP4 | Version 7.0.077 | ||
Conical Centrifuge Tubes (15 mL) | Corning | CLS430790 | |
Conical Centrifuge Tubes (50mL) | Corning | CLS430828 | |
Coot | Version 0.8.9.2 | ||
CrystalCap ALS | Hampton | HR4-779 | |
Curved-Tip Forceps | Mitegen | TW-CTF-1 | |
Deuterium chloride solution, 35 wt. % in D2O, ≥99 atom % D | Sigma-Aldrich | 543047 | |
Deuterium oxide 99.9 atom % D | Sigma-Aldrich | 151882 | |
Dual Thickness MicroLoops 1000 µm | Mitegen | M5-L18SP-1000 | |
FiveEasy pH meter F20-Std-Kit | Mettler Toledo | 30266626 | |
Foam Dewars Standard Vessel 800 ml | Spearlab | M-FD-800 | |
Four Color Mounting Clay | Hampton | HR4-326 | |
HEPES, BioUltra, for molecular biology, ≥99.5% (T), | Sigma-Aldrich | 54457 | |
High flux rotating anode X-ray diffractomemeter with EIGER 4M detector | Rigaku, Oxford Cryostream and Dectris | XtaLAB Synergy-R | Home source X-ray diffractometer |
Magnetic Wand Straight | Mitegen | M-R-1013198 | |
Microloader, tip for filling Femtotips and other glass microcapillaries (for research use only), 0.5 – 20 µL, 100 mm, light gray, 192 pcs. (2 racks × 96 pcs.) | Eppendorf | 930001007 | |
Microtubes volume 1.5 mL | Eppendorf | Z606340 | |
Petri Dishes with Clear Lid 100 mm diameter | Fischerbrand | FB0875713 | |
Phenix | Version 1.14-3260 | ||
Pin Tong 18 mm | Mitegen | M-R-1013196 | |
Pipette Volume 0.1-2.5 μL | Eppendorf Research | Z683779 | |
Pipette Volume 100-1000 μL | Eppendorf Research | Z683825 | |
Pipette Volume 10-100 μL | Eppendorf Research | Z683809 | |
Pipette Volume 20-200 μL | Eppendorf Research | Z683817 | |
Poly(ethylene glycol) BioXtra, average mol wt 3,350, powder | Sigma-Aldrich | P4338 | |
Quartz Capillary , 1.00 mm inner diameter, 80 mm length | Hampton | HR6-146 | Thin-walled capillary |
Research Stereomicroscope System | Olympus | SZX16 | |
Reusable B3 (SSRL/SAM Style) Goniometer Bases | Mitegen | GB-B3-R | |
Round – Miniature Hollow Glass Tubing (VitroTubes) Clear Fused Quartz / 1.00 mm inner diameter, 100 mm length | VitroCom | CV1012 | Thick-walled capillary |
Sandwich Box with cover | Hampton | HR3-132 | |
Siliconized 9 Well Glass Plate | Hampton | HR3-134 | |
Sitting Drop Crystallization Plate (24 Big Well) | Mitegen | XQ-P-24S-A | |
Sodium deuteroxide solution, 40 wt. % in D2O, 99 atom % D | Sigma-Aldrich | 176788 | |
Thick Siliconized circle cover slides (22 mm x 0.96 mm) | Hampton | HR3-247 | |
Universal Pipet Tips, 0.1 – 10 µL | VWR | 76322-528 | |
Universal Pipet Tips, 1 – 100 µL | VWR | 76322-136 | |
Universal Pipet Tips, 100 – 1000 µL | VWR | 76322-154 | |
Universal Pipet Tips, 20 – 200 µL | VWR | 76322-150 | |
Universal Pipet Tips, 1 – 20 µL | VWR | 76322-134 | |
Wax pen | Hampton | HR4-342 |