
==== Front
bioRxiv
BIORXIV
bioRxiv
2692-8205
Cold Spring Harbor Laboratory

39257752
10.1101/2024.08.29.610380
preprint
1
Article
Energy filtering enables macromolecular MicroED data at sub-atomic resolution
http://orcid.org/0000-0002-5466-6508
Clabbers Max T.B. 12
http://orcid.org/0000-0002-8936-0912
Hattne Johan 12
http://orcid.org/0000-0003-0055-230X
Martynowycz Michael W. 2†
http://orcid.org/0000-0002-9254-4069
Gonen Tamir 123*
1 Howard Hughes Medical Institute, University of California, Los Angeles, CA 90095
2 Department of Biological Chemistry, University of California, Los Angeles, CA 90095
3 Department of Physiology, University of California, Los Angeles, CA 90095
† Present address: Hauptman-Woodward Medical Research Institute, Buffalo, NY 14203; Department of Structural Biology, University at Buffalo, The State University of New York, Buffalo, NY 14203

* Corresponding author: tgonen@g.ucla.edu
29 8 2024
2024.08.29.610380https://creativecommons.org/licenses/by/4.0/ This work is licensed under a Creative Commons Attribution 4.0 International License, which allows reusers to distribute, remix, adapt, and build upon the material in any medium or format, so long as attribution is given to the creator. The license allows for commercial use.
nihpp-2024.08.29.610380.pdf
High resolution information is important for accurate structure modelling. However, this level of detail is typically difficult to attain in macromolecular crystallography because the diffracted intensities rapidly fade with increasing resolution. The problem cannot be circumvented by increasing the fluence as this leads to detrimental radiation damage. Previously, we demonstrated that high quality MicroED data can be obtained at low flux conditions using electron counting with direct electron detectors. The improved sensitivity and accuracy of these detectors essentially eliminate the read-out noise, such that the measurement of faint high-resolution reflections is limited by other sources of noise. Inelastic scattering is a major contributor of such noise, increasing background counts and broadening diffraction spots. Here, we demonstrate that a substantial improvement in signal-to-noise ratio can be achieved using an energy filter to largely remove the inelastically scattered electrons. This strategy resulted in sub-atomic resolution MicroED data from proteinase K crystals, enabling accurate structure modelling and the visualization of detailed features. Interestingly, filtering out the noise revealed diffuse scattering phenomena that can hold additional structural information. Our findings suggest that combining energy filtering and electron counting can provide more accurate measurements at higher resolution, providing better insights into protein function and facilitating more precise model refinement.
==== Body
pmcAt atomic resolution, commonly defined as beyond 1.2 Å, individual atoms are fully resolved revealing an accurate model of the underlying structure1. Unfortunately, noise in the measurement often makes it difficult to obtain such high-resolution information from macromolecular crystals. The mean diffracted intensity decreases rapidly at higher resolution, and radiation damage limits the total information that can be recovered from each crystal. This necessitates the use of a low fluence. The signal-to-noise ratio is therefore poor, limiting the attainable resolution. Sufficiently strong data can still be retrieved in microcrystal electron diffraction (MicroED)2, supported by the use of fast and highly sensitive cameras that are effective even at low flux densities2–4. Recently, we reported a substantial improvement in data quality by recording data using electron counting on direct electron detectors5. The improved accuracy and resolution enabled ab initio phasing using a fragment-based approach5, and allowed identification of hydrogen bond networks and protonation state of atoms6. However, this was demonstrated for triclinic lysozyme, which is relatively small and forms crystals with low solvent content, while proteinase K, which is much larger and contains more solvent, did not reach similarly high resolution. Electron counting detectors do not eliminate all noise, and further optimization of the experimental setup and data collection strategies is needed to routinely obtain more high-resolution information.

Inelastic scattering is a major factor limiting the signal-to-noise ratio7,8. Upon an inelastic interaction, electrons lose a small fraction of their energy and are no longer coherent. This contributes to increased background noise and a broadening of the Bragg peaks, thereby affecting the kinematic diffraction signal. This poses a significant challenge, as inelastic events are 3–4 times more likely than elastic scattering in biological specimens8,9. Whereas lattice filtering or subtracting a smooth radial background can partially correct for this10,11, post-processing does not reduce the noise inherently present in the measurement and cannot recover weak high-resolution reflections that are at or beneath the noise level. Removing inelastically scattered electrons experimentally based on their energy loss is therefore preferable. Previous reports assessed the effects of energy filtering on macromolecular diffraction data, demonstrating a significant reduction in background noise and sharper spots3,12,13. Interestingly, whereas the intensity and model statistics of the filtered data improved, the resolution of the reported structure did not13.

We hypothesized that combining energy filtering with the improved accuracy and sensitivity of direct electron detection would enhance data quality and resolution in MicroED. Two important considerations to this approach are the total fluence that can be withstood by the crystals, and the flux tolerated by the camera before coincidence loss prevents an accurate representation of the counts. The diffraction signal therefore is proportionally weaker. To mitigate this, we devised the following protocol that continuously rotates the crystals slower, uses the same total fluence spread over a smaller wedge, and utilizes the energy filter to further improve the signal-to-noise ratio. We applied this strategy to collect electron-counted and energy-filtered MicroED data from crystalline lamellae of proteinase K that were machined to an optimal thickness of ~300 nm14. Previously, a maximum resolution of 1.4 Å was achieved from similar lamellae without the use of energy filtration15.

Initial diffraction images were taken to assess lamellae quality and achievable resolution (Fig. 1). Diffraction spots extended to the edge of the detector beyond 1.0 Å resolution, showing less background noise and sharper Bragg peaks compared to unfiltered data (Extended Data Fig. 1)5,15. Interestingly, energy filtering revealed diffuse scattering from the bulk solvent that is largely obscured in unfiltered data (Fig. 1)5,15. Continuous rotation MicroED data were collected from the best crystals, showing high quality information and spots extending out ~1.0 Å resolution (Extended Data Fig. 2-3). Individual datasets were processed and showed significant information up to 1.06 Å resolution (Extended Data Fig. 4). Not all lamellae diffracted to the same resolution depending on crystallinity, quality of the lamellae and lattice orientation. Data from 17 crystals were merged and truncated at 1.09 Å resolution to ensure high multiplicity and complete sampling of reciprocal space (Extended Data Table 1). Intensity statistics show a clear improvement in attainable resolution and a more gradual decrease in crystallographic quality indicators compared to unfiltered data (Fig. 2A)15. Importantly, improving the resolution more than doubled the number of unique reflections used for structure refinement (Extended Data Table 1). At atomic resolution, the electrostatic potential map showed highly resolved features, enabling accurate interpretation of the structural model and visualization of hydrogen atoms (Fig. 2B).

We demonstrate that optimizing MicroED data collection by integrating an energy filter with direct electron detection cooperatively increases data quality and resolution by mitigating the effects of inelastic scattering. The filtered data exhibited reduced background noise, sharper spots, and an improved signal-to-noise ratio. Better high-resolution data provides more detailed insight into protein structure and function, and may prove useful in visualization of hydrogen atoms and hydrogen bonding networks6, as well as charge distribution3,10. With this improved data collection approach, the attainable resolution for proteinase K is similar to the best reported proteinase K structure by X-ray crystallography16, indicating that as technologies improve MicroED can rival more established structural biology methods. Diffraction spots extended beyond 1.0 Å resolution on individual frames, suggesting that there is still room for improvement. A data collection strategy that combines a larger number of datasets3, or a serial approach using individual snapshots11, may further improve data quality and resolution. Furthermore, the possibility to separate inelastic scattering from the signal is unique for electrons. It uncovers diffuse scattering which can originate from protein disorder or dynamics, providing additional information that improves accurate modeling and refinement17. As energy filtering removes a major obstacle to better data quality, other sources of noise that remain can be addressed more effectively and potentially lead to further improved accuracy in macromolecular electron crystallography.

Methods

Crystallization

Proteinase K powder from Engyodontium/tritirachium album was purchased from Sigma Aldrich and used without further purification. Crystals were grown in batch by dissolving 40 mg/ml proteinase K in 20 mM MES-NaOH pH 6.5. The protein solution was mixed at a 1:1 ratio with a precipitant solution of 0.5 M NaNO3, 0.1 M CaCl2, 0.1 M MES-NaOH pH 6.5. The mixture was incubated at 4 °C. Microcrystals with dimensions ranging between 8–12 μm grew within 24h.

Sample preparation

Standard holey carbon electron microscopy grids (Quantifoil, Cu 200 mesh, R2/2) were glow discharged for 30 s at 15 mA on the negative setting. Samples were prepared using a Leica GP2 vitrification device set at 4 °C and 90% humidity. For each sample, 3 μl of crystal solution was deposited onto the grid, incubated for 10s, and any excess liquid was blotted away from the back side. Next, the sample was soaked for 30s with 3 μl cryoprotectant solution containing 30% glycerol, 250 mM NaNO3, 50 mM CaCl2, 60 mM MES-NaOH pH 6.5. After incubation, any excess solution was blotted away using filter paper for a second time. Immediately after, the grid was rapidly vitrified by plunging it into liquid ethane. Grids were stored in a liquid nitrogen dewar prior to further use.

Focused ion beam milling

Grids were loaded onto a Helios Hydra 5 CX dual-beam plasma FIB/SEM (Thermo Fisher Scientific). Prior to milling, grids were coated with a thin protective layer of platinum for 45s using the gas injection system. Microcrystals of proteinase K were machined using a stepwise protocol to an optimal thickness of approximately 300 nm using a 30 kV Argon plasma ion beam. Coarse milling steps were performed using a 2.0 nA current to a thickness of approximately 3 μm. Finer milling steps at 0.2 nA were used to thin the lamellae to 600 nm. Final polishing steps were performed at 60 pA down to 300 nm thickness. Grids were cryo-transferred immediately after to the TEM for data collection. Grids were rotated by 90° relative to the milling direction such that the rotation axis on the microscope is perpendicular to the milling direction.

Hardware setup

Data were collected on a Titan Krios G3i TEM (Thermo Fisher Scientific) equipped with a X-FEG operated at an acceleration voltage of 300 kV, a post-column Selectris energy filter, and a Falcon 4i direct electron detector. The microscope was aligned for low flux density conditions using the 50 μm C2 aperture, spot size 11, and gun lens setting 8 for a less bright but more coherent illumination. A parallel electron beam of 10 μm diameter was used for data collection. The flux density at these conditions is approximately 0.002 e-/Å2/s. The energy spread of the emitted electrons was characterized at ΔE = 0.834 eV +/− 0.006 eV at FWHM. The zero loss peak of the energy filter was first centered in imaging mode. On our system, there is an offset in the position of the zero-loss peak when switching between imaging and diffraction mode. Therefore, the energy filter slit was carefully offset manually in defocused diffraction mode to align it with the selected area (SA) aperture used for diffraction data collection.

Data collection

Electron-counted and energy-filtered MicroED data were collected using the continuous rotation method. Diffraction data were collected using a 150 μm SA aperture, corresponding to a beam diameter of ~3 μm at the sample plane as defined by the aperture. The energy filter zero-loss peak was centered in defocused diffraction at the 5 and 10 eV slit settings. The effective sample to detector distance was calibrated at 1402 mm using a standard evaporated aluminum grid (Ted Pella). Crystals were rotated at a slow angular increment of 0.0476 °/s covering a total tilt range of 20.0°. Data were collected over 420 s exposures at a total fluence of ~0.84 e-/Å2. Diffraction stills were taken using the same settings without any stage rotation from stationary crystals. No beam stop post energy filter was used. Even without a beam stop, the low resolution reflections close to the center, where most of the inelastic scattering normally accumulates, are in sharp contrast to the low background. Data were recorded on a Falcon 4i direct electron detector in electron counting mode operating at an internal frame rate of 320 Hz. The proactive dose protector was manually disabled. Raw data were written in electron event representation (EER) format with an effective readout speed of 308 frames per second.

Data processing

Individual MicroED datasets in EER format were binned by two and converted to SMV format using the MicroED tools (available at https://cryoem.ucla.edu), after summing batches of 308 frames and applying post-counting gain corrections. Individual MicroED datasets were processed using XDS18. Data were integrated up to a cross correlation between two random half sets that was still significant at the 0.1% level19. Individual datasets were analyzed and merged using XSCALE18 and XSCALE_ISOCLUSTER20. The merged data were truncated at a mean I/σI ≥ 1.0 and a CC1/2 of 33.0 % in the highest resolution shell. Data were merged using Aimless21.

Structure solution and refinement

The structure was phased by molecular replacement using electron scattering factors in Phaser22. The model was inspected and built using Coot23. Two calcium ions and one nitrate were placed, and a total of 21 residues were modelled in alternate conformations. The structure was refined using REFMAC524 and phenix.refine25. Refinement was done using electron scattering factors and individual anisotropic B-factors for all atoms except hydrogens. Hydrogen omit maps were calculated using REFMAC5. Difference peaks higher than 2σ and 3σ in the omit map within less than 1.0 Å from any known riding positions were identified as potential hydrogen atoms.

Extended Data

Extended Data Figure 1. Unfiltered diffraction pattern of a stationary proteinase K lamella. Areas highlighted by the insets are magnified and the corresponding peak profiles are plotted.

Extended Data Figure 2. Individual frame from MicroED data collection using a 10 eV energy filter slit setting. Areas highlighted by the insets are magnified and the corresponding peak profiles are plotted.

Extended Data Figure 3. Individual frame from MicroED data collection using a 5 eV energy filter slit setting. Areas highlighted by the insets are magnified and the corresponding peak profiles are plotted.

Extended Data Figure 4. Intensity statistics for individual MicroED datasets, featuring from left to right I/σI, CC1/2, and Rpim as function of the resolution.

Extended Data Table 1. Data processing and refinement statistics

Data collection	
No. of crystals	17	
Space group	P43212	
Unit cell dimensions		
 a, b, c (Å)	66.92, 66.92, 107.56	
 α, β, γ (°)	90, 90, 90	
Resolution (Å)*	56.90–1.09 (1.12–1.09)	
Observed reflections	2,811,895 (95,185)	
Unique reflections	98,530 (6,664)	
Multiplicity	28.5 (14.3)	
Completeness (%)	97.5 (90.3)	
Rmerge	0.284 (1.154)	
<I/σI>	7.8 (1.4)	
CC1/2 (%)	99.5 (33.0)	
Refinement	
Resolution (Å)	56.82–1.09	
No. of reflections	98,228	
Rwork/Rfree	0.149/0.182	
No. of atoms		
 Protein	2193	
 Ligand/ions	6	
 Waters	357	
Hydrogen omit maps		
 No. of hydrogen atoms	1942	
 Peaks at ≥ 3σ	737 (37.9%)	
 Peaks at ≥ 2σ	1280 (69.7%)	
<B> factors (Å2)		
 Overall	11.76	
 Protein	9.17	
 Ligand/ions	11.67	
 Waters	27.67	
R.m.s. deviations		
 Bond lengths (Å)	0.016	
 Bond angles (°)	1.24	
* Values in parenthesis are for the highest resolution shell. Data were truncated at a mean I/σI ≥ 1.0 and a cross correlation between two random half sets that is statistically significant at CC1/2 ~ 33%.

Acknowledgements

This study was supported by the National Institutes of Health P41GM136508 and the Department of Defense HDTRA1-21-1-0004 and the Howard Hughes Medical Institute.

Data availability

Coordinates and structure factors have been deposited to the PDB.

Figure 1. Energy-filtered diffraction pattern of a proteinase K lamella. Areas highlighted by the insets are magnified and the corresponding peak profiles are plotted.

Figure 2. A. Comparison of intensity statics for filtered and unfiltered MicroED data, featuring from left to right: I/σI, CC1/2, and Rpim. B. Slice through the structural model showing the electrostatic potential map, the location of the slice in the structure is indicated in the inset. C. Maps and hydrogen omit maps shown for Tyr128. D. Maps and hydrogen omit maps for active site residues and waters. Electrostatic potential 2mFo-DFc maps are shown in blue and contoured at 4σ (B), and 2.2σ (C-D). Difference mFo-DFc maps are shown in green and red for positive and negative density contoured at 3σ (C-D). Hydrogen omit maps are shown in green and contoured at 2.3σ (C-D). Difference peaks marked with an asterisk indicate potential hydrogen atoms that have not been included in riding positions prior to map calculations.
==== Refs
References

1. Wlodawer A. , Minor W. , Dauter Z. & Jaskolski M. Protein crystallography for non-crystallographers, or how to get the best (but not more) from published macromolecular structures. FEBS J. 275 , 1–21 (2008).
2. Nannenga B. L. , Shi D. , Leslie A. G. W. & Gonen T. High-resolution structure determination by continuous-rotation data collection in MicroED. Nat. Methods 11 , 927–930 (2014).25086503
3. Yonekura K. , Kato K. , Ogasawara M. , Tomita M. & Toyoshima C. Electron crystallography of ultrathin 3D protein crystals: Atomic model with charges. Proc. Natl. Acad. Sci. 112 , 3368–3373 (2015).25730881
4. Clabbers M. T. B. Protein structure determination by electron diffraction using a single three-dimensional nanocrystal. Acta Crystallogr. Sect. Struct. Biol. 73 , 738–748 (2017).
5. Martynowycz M. W. , Clabbers M. T. B. , Hattne J. & Gonen T. Ab initio phasing macromolecular structures using electron-counted MicroED data. Nat. Methods 19 , 724–729 (2022).35637302
6. Clabbers M. T. B. , Martynowycz M. W. , Hattne J. & Gonen T. Hydrogens and hydrogen-bond networks in macromolecular MicroED data. J. Struct. Biol. X 6 , 100078 (2022).36507068
7. Cowley J. M. & Moodie A. F. The scattering of electrons by atoms and crystals. I. A new theoretical approach. Acta Crystallogr. 10 , 609–619 (1957).
8. Latychevskaia T. & Abrahams J. P. Inelastic scattering and solvent scattering reduce dynamical diffraction in biological crystals. Acta Crystallogr. Sect. B Struct. Sci. Cryst. Eng. Mater. 75 , 523–531 (2019).
9. Henderson R. The potential and limitations of neutrons, electrons and X-rays for atomic resolution microscopy of unstained biological molecules. Q. Rev. Biophys. 28 , 171–193 (1995).7568675
10. Mitsuoka K. The structure of bacteriorhodopsin at 3.0 Å resolution based on electron crystallography: implication of the charge distribution 1 1Edited by R. Huber. J. Mol. Biol. 286 , 861–882 (1999).10024456
11. Bücker R. Serial protein crystallography in an electron microscope. Nat. Commun. 11 , 996 (2020).32081905
12. Yonekura K. , Maki-Yonekura S. & Namba K. Quantitative Comparison of Zero-Loss and Conventional Electron Diffraction from Two-Dimensional and Thin Three-Dimensional Protein Crystals. Biophys. J. 82 , 2784–2797 (2002).11964264
13. Yonekura K. , Ishikawa T. & Maki-Yonekura S. A new cryo-EM system for electron 3D crystallography by eEFD. J. Struct. Biol. 206 , 243–253 (2019).30928615
14. Martynowycz M. W. , Clabbers M. T. B. , Unge J. , Hattne J. & Gonen T. Benchmarking the ideal sample thickness in cryo-EM. Proc. Natl. Acad. Sci. 118 , e2108884118 (2021).34873060
15. Martynowycz M. W. A robust approach for MicroED sample preparation of lipidic cubic phase embedded membrane protein crystals. Nat. Commun. 14 , 1086 (2023).36841804
16. Masuda T. Atomic resolution structure of serine protease proteinase K at ambient temperature. Sci. Rep. 7 , 45604 (2017).28361898
17. Meisburger S. P. , Case D. A. & Ando N. Diffuse X-ray scattering from correlated motions in a protein crystal. Nat. Commun. 11 , 1271 (2020).32152274
18. Kabsch W. XDS. Acta Crystallogr. D Biol. Crystallogr. 66 , 125–132 (2010).20124692
19. Karplus P. A. & Diederichs K. Linking Crystallographic Model and Data Quality. Science 336 , 1030–1033 (2012).22628654
20. Diederichs K. Dissecting random and systematic differences between noisy composite data sets. Acta Crystallogr. Sect. Struct. Biol. 73 , 286–293 (2017).
21. Evans P. R. & Murshudov G. N. How good are my data and what is the resolution? Acta Crystallogr. D Biol. Crystallogr. 69 , 1204–1214 (2013).23793146
22. McCoy A. J. Phaser crystallographic software. J. Appl. Crystallogr. 40 , 658–674 (2007).19461840
23. Emsley P. , Lohkamp B. , Scott W. G. & Cowtan K. Features and development of Coot. Acta Crystallogr. D Biol. Crystallogr. 66 , 486–501 (2010).20383002
24. Murshudov G. N. REFMAC 5 for the refinement of macromolecular crystal structures. Acta Crystallogr. D Biol. Crystallogr. 67 , 355–367 (2011).21460454
25. Afonine P. V. Towards automated crystallographic structure refinement with phenix.refine. Acta Crystallogr. D Biol. Crystallogr. 68 , 352–367 (2012).22505256
