CryoSPARC Guide
  • About CryoSPARC
  • Current Version
  • Licensing
    • Non-commercial license agreement
  • Setup, Configuration and Management
    • CryoSPARC Architecture and System Requirements
    • CryoSPARC Installation Prerequisites
    • How to Download, Install and Configure
      • Obtaining A License ID
      • Downloading and Installing CryoSPARC
      • CryoSPARC Cluster Integration Script Examples
      • Accessing the CryoSPARC User Interface
    • Deploying CryoSPARC on AWS
      • Performance Benchmarks
    • Using CryoSPARC with Cluster Management Software
    • Software Updates and Patches
    • Management and Monitoring
      • Environment variables
      • (Optional) Hosting CryoSPARC Through a Reverse Proxy
      • cryosparcm reference
      • cryosparcm cli reference
      • cryosparcw reference
    • Software System Guides
      • Guide: Updating to CryoSPARC v4
      • Guide: Installation Testing with cryosparcm test
      • Guide: Verify CryoSPARC Installation with the Extensive Validation Job (v4.3+)
      • Guide: Verify CryoSPARC Installation with the Extensive Workflow (≤v4.2)
      • Guide: Performance Benchmarking (v4.3+)
      • Guide: Download Error Reports
      • Guide: Maintenance Mode and Configurable User Facing Messages
      • Guide: User Management
      • Guide: Multi-user Unix Permissions and Data Access Control
      • Guide: Lane Assignments and Restrictions
      • Guide: Queuing Directly to a GPU
      • Guide: Priority Job Queuing
      • Guide: Configuring Custom Variables for Cluster Job Submission Scripts
      • Guide: SSD Particle Caching in CryoSPARC
      • Guide: Data Management in CryoSPARC (v4.0+)
      • Guide: Data Cleanup (v4.3+)
      • Guide: Reduce Database Size (v4.3+)
      • Guide: Data Management in CryoSPARC (≤v3.3)
      • Guide: CryoSPARC Live Session Data Management
      • Guide: Manipulating .cs Files Created By CryoSPARC
      • Guide: Migrating your CryoSPARC Instance
      • Guide: EMDB-friendly XML file for FSC plots
    • Troubleshooting
  • Application Guide (v4.0+)
    • A Tour of the CryoSPARC Interface
    • Browsing the CryoSPARC Instance
    • Projects, Workspaces and Live Sessions
    • Jobs
    • Job Views: Cards, Tree, and Table
    • Creating and Running Jobs
    • Low Level Results Interface
    • Filters and Sorting
    • View Options
    • Tags
    • Flat vs Hierarchical Navigation
    • File Browser
    • Blueprints
    • Workflows
    • Inspecting Data
    • Managing Jobs
    • Interactive Jobs
    • Upload Local Files
    • Managing Data
    • Downloading and Exporting Data
    • Instance Management
    • Admin Panel
  • Cryo-EM Foundations
    • Image Formation
      • Contrast in Cryo-EM
      • Waves as Vectors
      • Aliasing
  • Expectation Maximization in Cryo-EM
  • Processing Data in cryoSPARC
    • Get Started with CryoSPARC: Introductory Tutorial (v4.0+)
    • Tutorial Videos
    • All Job Types in CryoSPARC
      • Import
        • Job: Import Movies
        • Job: Import Micrographs
        • Job: Import Particle Stack
        • Job: Import 3D Volumes
        • Job: Import Templates
        • Job: Import Result Group
        • Job: Import Beam Shift
      • Motion Correction
        • Job: Patch Motion Correction
        • Job: Full-Frame Motion Correction
        • Job: Local Motion Correction
        • Job: MotionCor2 (Wrapper) (BETA)
        • Job: Reference Based Motion Correction (BETA)
      • CTF Estimation
        • Job: Patch CTF Estimation
        • Job: Patch CTF Extraction
        • Job: CTFFIND4 (Wrapper)
        • Job: Gctf (Wrapper) (Legacy)
      • Exposure Curation
        • Job: Micrograph Denoiser (BETA)
        • Job: Micrograph Junk Detector (BETA)
        • Interactive Job: Manually Curate Exposures
      • Particle Picking
        • Interactive Job: Manual Picker
        • Job: Blob Picker
        • Job: Template Picker
        • Job: Filament Tracer
        • Job: Blob Picker Tuner
        • Interactive Job: Inspect Particle Picks
        • Job: Create Templates
      • Extraction
        • Job: Extract from Micrographs
        • Job: Downsample Particles
        • Job: Restack Particles
      • Deep Picking
        • Guideline for Supervised Particle Picking using Deep Learning Models
        • Deep Network Particle Picker
          • T20S Proteasome: Deep Particle Picking Tutorial
          • Job: Deep Picker Train and Job: Deep Picker Inference
        • Topaz (Bepler, et al)
          • T20S Proteasome: Topaz Particle Picking Tutorial
          • T20S Proteasome: Topaz Micrograph Denoising Tutorial
          • Job: Topaz Train and Job: Topaz Cross Validation
          • Job: Topaz Extract
          • Job: Topaz Denoise
      • Particle Curation
        • Job: 2D Classification
        • Interactive Job: Select 2D Classes
        • Job: Reference Based Auto Select 2D (BETA)
        • Job: Reconstruct 2D Classes
        • Job: Rebalance 2D Classes
        • Job: Class Probability Filter (Legacy)
        • Job: Rebalance Orientations
        • Job: Subset Particles by Statistic
      • 3D Reconstruction
        • Job: Ab-Initio Reconstruction
      • 3D Refinement
        • Job: Homogeneous Refinement
        • Job: Heterogeneous Refinement
        • Job: Non-Uniform Refinement
        • Job: Homogeneous Reconstruction Only
        • Job: Heterogeneous Reconstruction Only
        • Job: Homogeneous Refinement (Legacy)
        • Job: Non-uniform Refinement (Legacy)
      • CTF Refinement
        • Job: Global CTF Refinement
        • Job: Local CTF Refinement
        • Job: Exposure Group Utilities
      • Conformational Variability
        • Job: 3D Variability
        • Job: 3D Variability Display
        • Job: 3D Classification
        • Job: Regroup 3D Classes
        • Job: Reference Based Auto Select 3D (BETA)
        • Job: 3D Flexible Refinement (3DFlex) (BETA)
      • Postprocessing
        • Job: Sharpening Tools
        • Job: DeepEMhancer (Wrapper)
        • Job: Validation (FSC)
        • Job: Local Resolution Estimation
        • Job: Local Filtering
        • Job: ResLog Analysis
        • Job: ThreeDFSC (Wrapper) (Legacy)
      • Local Refinement
        • Job: Local Refinement
        • Job: Particle Subtraction
        • Job: Local Refinement (Legacy)
      • Helical Reconstruction
        • Helical symmetry in CryoSPARC
        • Job: Helical Refinement
        • Job: Symmetry search utility
        • Job: Average Power Spectra
      • Utilities
        • Job: Exposure Sets Tool
        • Job: Exposure Tools
        • Job: Generate Micrograph Thumbnails
        • Job: Cache Particles on SSD
        • Job: Check for Corrupt Particles
        • Job: Particle Sets Tool
        • Job: Reassign Particles to Micrographs
        • Job: Remove Duplicate Particles
        • Job: Symmetry Expansion
        • Job: Volume Tools
        • Job: Volume Alignment Tools
        • Job: Align 3D maps
        • Job: Split Volumes Group
        • Job: Orientation Diagnostics
      • Simulations
        • Job: Simulate Data (GPU)
        • Job: Simulate Data (Legacy)
    • CryoSPARC Tools
    • Data Processing Tutorials
      • Case study: End-to-end processing of a ligand-bound GPCR (EMPIAR-10853)
      • Case Study: DkTx-bound TRPV1 (EMPIAR-10059)
      • Case Study: Pseudosymmetry in TRPV5 and Calmodulin (EMPIAR-10256)
      • Case Study: End-to-end processing of an inactive GPCR (EMPIAR-10668)
      • Case Study: End-to-end processing of encapsulated ferritin (EMPIAR-10716)
      • Case Study: Exploratory data processing by Oliver Clarke
      • Tutorial: Tips for Membrane Protein Structures
      • Tutorial: Common CryoSPARC Plots
      • Tutorial: Negative Stain Data
      • Tutorial: Phase Plate Data
      • Tutorial: EER File Support
      • Tutorial: EPU AFIS Beam Shift Import
      • Tutorial: Patch Motion and Patch CTF
      • Tutorial: Float16 Support
      • Tutorial: Particle Picking Calibration
      • Tutorial: Blob Picker Tuner
      • Tutorial: Helical Processing using EMPIAR-10031 (MAVS)
      • Tutorial: Maximum Box Sizes for Refinement
      • Tutorial: CTF Refinement
      • Tutorial: Ewald Sphere Correction
      • Tutorial: Symmetry Relaxation
      • Tutorial: Orientation Diagnostics
      • Tutorial: BILD files in CryoSPARC v4.4+
      • Tutorial: Mask Creation
      • Case Study: Yeast U4/U6.U5 tri-snRNP
      • Tutorial: 3D Classification
      • Tutorial: 3D Variability Analysis (Part One)
      • Tutorial: 3D Variability Analysis (Part Two)
      • Tutorial: 3D Flexible Refinement
        • Installing 3DFlex Dependencies (v4.1–v4.3)
      • Tutorial: 3D Flex Mesh Preparation
    • Webinar Recordings
  • Real-time processing in cryoSPARC Live
    • About CryoSPARC Live
    • Prerequisites and Compute Resources Setup
    • How to Access cryoSPARC Live
    • UI Overview
    • New Live Session: Start to Finish Guide
    • CryoSPARC Live Tutorial Videos
    • Live Jobs and Session-Level Functions
    • Performance Metrics
    • Managing a CryoSPARC Live Session from the CLI
    • FAQs and Troubleshooting
  • Guides for v3
    • v3 User Interface Guide
      • Dashboard
      • Project and Workspace Management
      • Create and Build Jobs
      • Queue Job, Inspect Job and Other Job Actions
      • View and Download Results
      • Job Relationships
      • Resource Manager
      • User Management
    • Tutorial: Job Builder
    • Get Started with CryoSPARC: Introductory Tutorial (v3)
    • Tutorial: Manually Curate Exposures (v3)
  • Resources
    • Questions and Support
Powered by GitBook
On this page
  • At a glance
  • Description
  • Inputs
  • Particles
  • Volume
  • Mask
  • Outputs
  • Volume
  • Particles
  • Mask
  • Commonly Adjusted Parameters
  • Alignment parameters
  • Refinement parameters
  • Common Problems
  • Masking artifacts
  • Masked subvolume too small
  • Improper search extents and/or fulcrum position
  • Image Recentering
  • Common Next Steps
  • Example Workflow
  1. Processing Data in cryoSPARC
  2. All Job Types in CryoSPARC
  3. Local Refinement

Job: Local Refinement

PreviousLocal RefinementNextJob: Particle Subtraction

Last updated 1 year ago

At a glance

Refine a volume or a part of a volume to high resolution by incorporating existing pose information to prevent significant deviation from a known good alignment.

  • Use particles from an alignment with existing gold-standard half sets

  • Create a mask with a soft edge to focus the refinement on a subvolume

  • Generally a user-selected fulcrum position produces best results

  • Impose a Gaussian prior to regularize the problem when the subvolume is small or the input alignment is already of medium- to high-quality

Description

Local Refinement is a technique for generating a single 3D volume from a stack of particle images, like or . However, where Homogeneous and Non-Uniform refinements discard pre-existing alignment information, Local Refinement incorporates information from previous alignments from upstream jobs. This presents several advantages:

  • Sub-regions which are too small to align on their own may be improved by Local Refinement

  • Information about the rotation and translation needed to align the consensus volume compared to the desired sub-region provides information about flexibility

  • Combined with other job types, Local Refinement enables advanced workflows for investigation of pseudosymmetry or symmetry with flexibility

In each iteration of a local refinement, the volume is masked. Then, the algorithm finds the optimal pose of the masked subvolume within the given search extent of rotations around a fulcrum and translations for each particle image. The map is updated using the new poses, and the process continues.

Local Refinement should be your first choice when attempting to resolve two subdomains which are flexible relative to each other but internally rigid. For example, the head and foot regions of the bend around flexible hinges, but are otherwise relatively rigid. Local refinement significantly improved map quality for those domains.

If the domains undergo continuous deformations or conformational changes, algorithms like or might produce better results.

Inputs

Particles

Unlike other refinement jobs, particle inputs to Local Refinement must have been previously aligned in 3D, and must have a FSC half-set split. These conditions are satisfied by, for example, Homogeneous or Non-Uniform refinement.

Particle Subtraction and Local Refinement

Particle subtraction attempts to accommodate the fact that we cannot apply a mask to our particle images. It subtracts a projection of the excluded regions of a map from the particle images, leaving behind the parts of the image corresponding to what we want to align. Because this subtraction requires projection of the information in our volume by the pose of the particle, it requires a high-quality initial alignment to work properly. See the Particle Subtraction job page for more detail.

Volume

Local Refinement accepts any volume with half maps as an input.

Mask

Local Refinement requires a mask that includes the region to refine. Note that the mask does not necessarily have to exclude anything — a mask covering the entire volume is acceptable.

Outputs

Volume

The refined sub-volume is not masked. It is possible that new features will become visible which are outside the mask if they are rigidly associated with the refined sub-volume.

Particles

Particles are output with new poses.

Mask

If dynamic masking is enabled, the refined mask is output. Otherwise the mask is the same as the input.

Commonly Adjusted Parameters

Alignment parameters

The alignment parameters control how much the sub-volume is rotated during the search for its optimal pose. Broadly speaking, the search extent should be large enough to capture the range of motion of the sub-volume, but small enough that the prior alignment information is still useful in preventing overfitting.

Use pose/shift Gaussian prior during alignment

If true, during each iteration, poses far from the input pose are softly penalized. This increases stability while still allowing for an alignment far from the input, if the score is sufficiently high.

When this option is enabled, two additional settings are available to control the standard deviation of the prior over rotation and shifts, respectively. Wider standard deviations “weaken” the prior, while tighter ones increase the penalty on poses far from the input.

Generally, your prior should follow your expected distribution of the particles. For instance, when refining a highly flexible domain, wide shifts are appropriate. On the other hand, when improving detail in well-aligned, rigid region of the target a prior as small as 3° or 2 Å may work best.

Regardless of the width of the Gaussian, no poses outside the search extent will be considered. If priors are enabled, then the search extents are given a default value of three times the prior standard deviations, for both the rotation and shift priors.

For small masks or low SNR data, the use of these priors is encouraged to help increase stability and reduce overfitting.

Search Extents

These parameters set the maximum distance from the starting alignment that the sub-volume will be rotated (Rotation search extent) or translated (Shift search extent) in degrees and Å, respectively.

Re-centering

In each iteration, the sub-volume is assigned a new pose (rotation and shift). These settings control whether future iterations should start from the input pose (False) or the updated pose (True) for the rotations and shifts, respectively.

If either of these settings are set to True, we strongly recommend that a pose/shift Gaussian prior is applied to the alignments to prevent particles from drifting far from their initial pose in early iterations.

These settings re-center both the search grid and the Gaussian prior, if applied.

Fulcrum Location

Local refinement requires a fulcrum about which it will rotate the search volume while aligning to particles. By default, this is the center of the input mask. It can also be set to the center of the box (the fulcrum for all other refinements).

Additionally, you can specify an exact fulcrum location. This can be entered as a triplet of coordinates with the suffix px or A to indicate the triplet is being specified in pixels or Angstroms, respectively. Specify this coordinate in x, y, z order, with 0, 0, 0 px or 0, 0, 0 Acorresponding to the corner of the box.

Fulcrum Location

The choice of a fulcrum location is an important parameter to consider. An inappropriate fulcrum location can exaggerate or dampen the expected results of a given rotation by coupling it with a translation. See Improper search extents and/or fulcrum position for more information.

Refinement parameters

Marginalization

Enabling marginalization allows for insertion of particle images at multiple poses, weighted by the likelihood of that pose. Marginalization can improve results when the target is smaller, when overfitting occurs, or when the particle stack still contains some noisy outliers.

Non-Uniform Refinement Enable

Symmetry

Mask

When set to static, use the input mask. When set to dynamic, perform dynamic masking to re-generate a mask at each iteration of refinement. We generally recommend static masking to avoid mask-based overfitting.

If dynamic masking is used, the following parameters influence mask generation:

  • Regions of the map with a value greater than the Dynamic mask threshold are included in the mask.

  • Dynamic mask near and far control the distance in Å from the thresholded map edge to begin and end the soft edge, respectively.

Minimize over per-particle scale

When enabled, optimize the per-particle scale at each refinement iteration. In general, this should be kept off due to the instability of particle scale when considering a masked subregion compared to the full particle image.

Initial lowpass resolution

The input map will be low-pass filtered to this resolution before the first iteration. Subsequent iterations will be filtered to the GSFSC resolution.

The initial lowpass resolution can be an important parameter to tune. The default value of 12 Å is conservative. If the input particles are already well-aligned, lowpass filtering worsens alignments in the initial iterations by discarding good alignment information. At best, this results in slightly longer job runtimes as the initial iteration is spent getting back to the quality of the input alignments. At worst, for smaller regions of the particle which are difficult to initially align, the particles may shift too far from their input alignments and never recover the quality of the initial alignments.

Common Problems

Masking artifacts

Masked subvolume too small

Since local refinement uses the masked volume to search for the proper pose, if the masked subregion is too small there will not be enough information for a high-quality alignment. If the angle and shifts change by a large amount and the sub-volume becomes noisy and low-quality, this is a likely culprit.

Generally, your prior should follow your expected distribution of the particles. For instance, when refining a highly flexible domain, wide shifts are appropriate. On the other hand, when improving detail in well-aligned, rigid region of the target a prior as small as 3° or 2 Å may work best.

Pose and shift change histograms are expected to be peaked relatively close to 0, and to decay outward. A large proportion of particles near the right extreme of the histograms indicates poor alignability, which may be alleviated by the use of a gaussian prior or a larger mask. The example histograms above illustrate an overfitted, poorly-aligned sub-volume (blue) and a well-regularized alignment of the same particles (grey) using a Gaussian prior and larger subvolume.

Improper search extents and/or fulcrum position

Ideally, the distribution of angle and shift changes is smooth after the final iteration of local refinement. If either distribution has peaks, especially if those peaks are at the furthest extent of the search range, the refinement may benefit from an increased search range.

A poor choice of fulcrum position may also result in peaks in the angle and shift change plots. A fulcrum position that is far away from the true point about which the sub-volume rotates requires a large translation to correct the alignment. Therefore, especially where large shifts are observed, optimizing the fulcrum may prove beneficial.

Image Recentering

Common Next Steps

A local refinement often provides the best possible map for a given rigid sub-volume. If the volume improves, creating a mask around the new features may iteratively increase map quality.

Example Workflow

Particles can also be provided from a job, but subtraction is not required.

Ensure that your mask has a sufficiently soft edge to avoid ringing artifacts. For more information on mask creation, please see .

increased alignment quality when turning both re-centering and Gaussian priors on.

For an example of fulcrum selection, see the .

Enable cross-validation-optimal filtering, as in a . This can significantly improve results for membrane proteins or other targets with disordered regions.

Enforce point-group symmetry. Note that this parameter uses the box center, regardless of the choice of fulcrum. We recommend performing rather than using this setting, as symmetry expansion can accommodate minor breaks in symmetry.

We recommend that this value is set to the local resolution for the sub-volume to be refined (which can be determined with a job). If the sub-volume to be refined is large and well-aligned, a heuristic of 1–2 Å higher value (that is, worse resolution) than the GSFSC resolution may also give good results.

As discussed in , any time a mask “cuts through” density, a soft edge is absolutely essential. In practice, it is often best to try several widths of dilation and padding to find the optimal combination.

If a small volume must be masked and aligned, this problem can be somewhat mitigated by imposing a on the shifts and rotations. With this prior, poses receive a penalty to their score which increases the further the pose is from the input pose. The smaller the standard deviation of the Gaussian, the harsher the penalty. Penalizing large shifts in this way can reduce overfitting of small sub-volumes.

In some cases, that re-centering the images using Volume Alignment Tools significantly improves the results of Local Refinement. If the results of a Local Refinement are unexpectedly poor this may be one possible solution.

If after local refinement you suspect that the aligned sub-volume may have some heterogeneity, will quickly tease apart discrete conformations in your aligned sub-volume. Using the mask from your Local Refinement job and a solvent mask around the whole map will likely provide the best results.

Homogeneous
Non-Uniform refinement
Yeast U4/U6.U5 tri-snRNP (spliceosome)
3D Variability Analysis
3D Flexible Refinement
Particle Subtraction
the relevant guide page
Some users have reported
tri-snRNP case study
Non-Uniform Refinement
symmetry expansion
Local Resolution Estimation
Mask Selection and Creation
users have reported
3D Classification
Case Study: Yeast U4/U6.U5 tri-snRNP
Gaussian prior
A Gaussian prior penalizes poses far from the initial value.
Search extents control the maximum displacements of the sub-volume.
Recentering controls whether each iteration searches a region around the input pose or the optimal pose found in the previous iteration.
The two default fulcrum locations (left: mask_center; middle: box_center) and a user-defined fulcrum (right).
Small sub-volumes do not provide enough information for reliable alignment, resulting in noisy, over-fitted results.
Rotation about a fulcrum that is far from the true center of rotation results in an additional shift. This may result in poor alignment.