Unbend, correction of local beam-induced sample motion in cryo-EM images using a 3D spline model

  1. Lingli Kong  Is a corresponding author
  2. Ximena Zottig
  3. Johannes Elferich
  4. Nikolaus Grigorieff  Is a corresponding author
  1. RNA Therapeutics Institute, University of Massachusetts Chan Medical School, United States
  2. Howard Hughes Medical Institute, University of Massachusetts Chan Medical School, United States
8 figures, 2 tables and 1 additional file

Figures

Uniform cubic spline curve.

(a) Diagram showing the core components of a cubic B-spline curve, such as knots, control points, and the uniform knot positions ti. (b) Diagram showing the free-end boundary condition for a cubic B-spline.

3D spline model knot grid construction and shift field generation.

(a) Diagram showing the position of knots along the exposure accumulation direction and the shift amount along x for each movie frame. (b) A 3D diagram showing the cubic B-spline along the exposure accumulation direction, with red arrows representing the shifts where the knots are located (black dots on the gray line), and the cyan curve showing the spline curve for interpolating the shifts for each frame. (c, d) Shift field along x and y showing the shift amount of each pixel obtained from the bicubic B-splines on 30×30 pixels frame (for demonstration only). (e) Schematic image grid before (gray grids) and after (cyan grid) applying the shift fields in (c, d). The red arrows represent the shift of patches. Black dots represent knot positions.

Figure 2—source data 1

Unbend movie alignment and local motion correction algorithm.

https://cdn.elifesciences.org/articles/109119/elife-109119-fig2-data1-v1.xlsx
A micrograph of whole-cell M. pneumoniae.

(a) Motion trajectory panel in the cisTEM graphical user interface (GUI). (b) The cropped area shown in (a) without local distortion correction, experiencing blurring due to the local distortion (upper) and the power spectrum (lower). (c) The cropped area (upper) and the power spectrum of the local distortion-corrected micrograph (lower). (d) Detected particles (white) using 2D template matching (2DTM) on the full-frame aligned micrograph. (e) Detected particles (white) using 2DTM on the local-distortion corrected micrograph.

Figure 4 with 5 supplements
Patch shift in different types of samples.

(a) Box plots of the standard deviation of patch motions calculated for micrographs of different types of samples. (b) The mean patch shifts for each micrograph, colored by sample type.

Figure 4—source data 1

Summary statistics of per-micrograph mean patch shifts by sample type.

https://cdn.elifesciences.org/articles/109119/elife-109119-fig4-data1-v1.xlsx
Figure 4—source data 2

Motion and equivalent strain summary table for micrographs shown in Figure 4—figure supplements 1–5.

https://cdn.elifesciences.org/articles/109119/elife-109119-fig4-data2-v1.xlsx
Figure 4—figure supplement 1
M. pneumoniae sample.

(a) Screenshot of the cisTEM graphical user interface (GUI) displaying the motion-corrected micrograph and patch trajectories. (b) Micrograph based on full-frame alignment. (c) Von Mises equivalent strain field. (d) Total equivalent strain field.

Figure 4—figure supplement 2
BS-C-1 cell edge sample.

(a) Screenshot of the cisTEM graphical user interface (GUI) displaying the motion-corrected micrograph and patch trajectories. (b) Micrograph based on full-frame alignment. (c) Von Mises equivalent strain field. (d) Total equivalent strain field.

Figure 4—figure supplement 3
M. musculus ER-HoxB8 cell lamella sample.

(a) Screenshot of the cisTEM graphical user interface (GUI) displaying the motion-corrected micrograph and patch trajectories. (b) Micrograph based on full-frame alignment. (c) Von Mises equivalent strain field. (d) Total equivalent strain field.

Figure 4—figure supplement 4
S. cerevisiae lamella sample.

(a) Screenshot of the cisTEM graphical user interface (GUI) displaying the motion-corrected micrograph and patch trajectories. (b) Micrograph based on full-frame alignment. (c) Von Mises equivalent strain field. (d) Total equivalent strain field.

Figure 4—figure supplement 5
BS-C-1 cell lysate sample.

(a) Screenshot of the cisTEM graphical user interface (GUI) displaying the motion-corrected micrograph and patch trajectories. (b) Micrograph based on full-frame alignment. (c) Von Mises equivalent strain field. (d) Total equivalent strain field.

Magnitude of deformation correction of the micrograph shown in Figure 3a.

(a) Von Mises equivalent strain and (b) total equivalent strain.

Evaluation of local motion correction using Unbend across different sample types.

(a) Box plots displaying the distribution of the number of 2D template matching (2DTM) detections per micrograph. (b) Histograms of the 2DTM signal-to-noise ratio (SNR) of detected 2DTM targets. Unbend improves the overall distribution of SNR values for all sample types. (c) 2DTM SNR difference between Unbend and Unblur. Across all sample types, Unbend led to improvements in either the number of detected targets or 2DTM SNR values, or both.

Impact of patch-based motion correction on the detection of 60S subunits in focused ion beam-milled lamellae analyzed in a large dataset.

(a) Histogram showing the percentage of increase in 60S detection in montages processed using Unbend, compared to Unblur. (b) Scatterplot of the number of detected 60S per exposure in a representative montage that did not show a substantial increase of detection after processing with Unbend, compared to Unblur. (c) Similar to panel (b), but for a representative montage with a more substantial increase of detected targets after processing with Unbend, compared to Unblur.

Comparison of local motion correction methods.

(a) Box plots showing the distribution of the number of 2D template matching (2DTM) detections per micrograph for each motion correction method. (b) Distribution of 2DTM signal-to-noise ratio (SNR) differences between Unbend and the other tested methods for particles detected in common (ΔSNR=SNRUnbendSNRother). Vertical lines mark the zero-crossing in the sorted ΔSNR distribution for each method. Percentages in the legend indicate the fraction of commonly detected particles with higher SNR in Unbend-processed micrographs.

Figure 8—source data 1

Statistics table for detected particles from micrographs processed by different software.

https://cdn.elifesciences.org/articles/109119/elife-109119-fig8-data1-v1.xlsx

Tables

Table 1
Unbend runtimes on movies analyzed in Figure 6.
M. pneumoniaeAethiops: BS-C-1(cell edge)S. cerevisiaeMusculus:ER-HoxB8Aethiops: BS-C-1(cell lysate)
Input movie size
(pixels)
5760×40925760×409211520×818411520×818411520×8184
No. of frames per movie2430507530
Exposure per frame
(e2)
1.2961.0210.6000.6001.004
No. of movies6465305976
Input pixel size
(Å/pixel)
1.0531.330.530.4150.415
Defocus
(µm) (mean ± standard deviation)
1.3±0.51.5±0.30.4±0.11.0±0.51.2±0.4
Patch number8×610×88×67×57×5
Output image size
(pixels)
4044×28735107×36284070×28923187×22643187×2264
Output pixel size
(Å/pixel)
1.51.51.51.51.5
Runtime per movie
(mean ± standard deviation)
1m34s±5 s3m21s±35 s7m58s±26 s12m30s±43 s4m49s±9 s
Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Cell line (M. pneumonia)M129O’Reilly et al., 2020;
Lucas et al., 2021
ATCC 29342
Cell line (C. aethiops)BS-C-1ATCCCCL-26;
RRIDs:CVCL_0607
Cell line (M. musculus)ER-HoxB8Elferich et al., 2022
Cell line (S. cerevisiae)BY4741Lucas et al., 2022ATCC S288C
Cell line (C. albicans)PY6413Serrano et al., 2026
Software, algorithmcisTEMLucas et al., 2021;
Grant et al., 2018;
Elferich et al., 2024
https://cistem.org/

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Lingli Kong
  2. Ximena Zottig
  3. Johannes Elferich
  4. Nikolaus Grigorieff
(2026)
Unbend, correction of local beam-induced sample motion in cryo-EM images using a 3D spline model
eLife 14:RP109119.
https://doi.org/10.7554/eLife.109119.3