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
  • Definition of helical symmetry in CryoSPARC
  • Advanced notes on helical symmetry
  • Exploration of symmetry parameters
  • Using asymmetric reconstructions and symmetry search utility job to explore symmetry parameters
  • Citations
  1. Processing Data in cryoSPARC
  2. All Job Types in CryoSPARC
  3. Helical Reconstruction

Helical symmetry in CryoSPARC

Explanation of treatment of helical symmetry in CryoSPARC.

PreviousHelical ReconstructionNextJob: Helical Refinement

Last updated 1 year ago

Many proteins that form filaments can be characterized as having some form of helical symmetry. In cryo-EM, symmetry is a powerful tool to both increase the effective signal in a dataset, and to provide a structural constraint on the reconstruction. With helical symmetry, the latter point is especially important in order to arrive at correct, high-resolution structures. However, a major challenge in reconstruction of helical assemblies is determining the parameters that define the helical symmetry. While it is often necessary, there is currently no universal method for determining or validating helical symmetry parameters. This page details how helical symmetry is defined and treated in CryoSPARC, as well as some tools developed in CryoSPARC to assist in the exploration of symmetry parameters.

Definition of helical symmetry in CryoSPARC

This section covers some of the math involved in representing helical symmetry. Similar to point group symmetry, proteins with helical symmetry are comprised of a single asymmetric unit that exists in multiple different spatial positions. With point group symmetry, the transformations that relate one copy of the subunit to all others are characterized by pure rotations around a fixed point. With helical symmetry, these transformations are characterized instead by a simultaneous rotation ("twist") and translation ("rise") along the same axis (called the helical axis, "meridian", or the "screw" axis). The image below shows an example of a simple one-start helical lattice on the right, compared to a continuous helix on the left [1].

Global helical symmetry can be defined through two sets of equivalent parameters:

Advanced notes on helical symmetry

CryoSPARC also supports additional point group symmetry associated with a helical assembly, namely cyclic symmetries (with a cyclic axis coincident with the helical axis), and dihedral symmetries, (with the dyad axis perpendicular to the helical axis). These are often referred to as "n-start" helices, where n is the cyclic order.

It should also be noted that in proteins that form a helical lattice, there are many ambiguities in the characterization of helical symmetry. For the purposes of symmetry imposition during refinement, the helical twist and rise pair should:

  • Correspond to a valid symmetry transform that translates an individual asymmetric unit to another identical asymmetric unit on the helical lattice, and

  • Have the minimum helical rise out of all such valid transformations

The first condition ensures that asymmetric units are properly averaged together, and the second condition ensures that symmetry is used to its maximum extent.

In addition to these ambiguities, there may also be ambiguities in the hand of helical symmetry. Since the hand of a cryo-EM map is not constrained (projection images carry no information of the correct hand), it's possible for maps to reach high resolution but have an incorrect hand. It's generally only possible to observe the hand of a density map once the resolution is high enough that structural motifs such as alpha helices can be seen – if these are left-handed, that likely indicates the hand of the map is inverted. If symmetry was applied, and the map has the wrong hand, the helical twist should have its sign inverted in order to produce the correct symmetry. Note that the Volume Tools job can invert the hand of a map.

Finally, more advanced cases exist, such as:

  • helical proteins with multiple starts but without cyclic symmetry (i.e. "screw symmetry", or "n-start C1 helical symmetry")

  • pseudo-helical symmetry (e.g. microtubules with a symmetry-breaking seam)

The bottom line is that is most cases of global helical symmetry, a valid twist/rise representation will exist. In the first case of screw symmetry, these are reducible into a valid representation within the global twist/rise parameterization, as described by He and Scheres [2]. On the other hand, microtubules with seams may or may not be able to fit into a global helical symmetry formulation. In cases where you suspect the helical symmetry is not global, it may be insightful to refine the structures both with and without symmetry applied, and observe the differences.

Exploration of symmetry parameters

Note: Reliable determination of helical symmetry is an open problem in cryo-EM reconstruction, and it has been observed that the error landscape over helical symmetry parameters contains many false minima [3]. Furthermore, validation of helical symmetry parameters poses an even greater challenge, as standard methods of validation (including FSC) often are not sufficient when comparing two reconstructions with differing helical symmetry. For these reasons, the tools and workflows presented on this page are not intended to immediately provide unambiguous solutions to, and/or validate, the helical symmetry parameters. Rather, these tools are primarily intended to help guide the initial exploration of candidate symmetry parameters, and to enable refining these candidate solutions during a helical refinement. Though the success of these methods is often dataset-dependent, they may be useful as part of the process of determining helical symmetry.

Where possible, validation of imposed helical symmetry parameters should be done by cross-referencing with other methods (e.g. prior knowledge, known similar structures, Fourier-Bessel analysis, etc.). In addition, density maps with imposed symmetry should always be inspected for expected structural motifs and secondary structure. Often times, imposing incorrect symmetry parameters can still produce refined maps with relatively high claimed resolutions, but structural motifs may appear highly distorted, may be "smeared" out along the azimuth, or may not exist altogether. Finally, FSC plots should also always be inspected. High amplitude oscillations in the FSC curves are typical in early iterations when the resolution is still fairly low, but if they persist into the final iterations of the refinement, that may indicate that the wrong symmetry was imposed.

Using asymmetric reconstructions and symmetry search utility job to explore symmetry parameters

Using helical refinement

  • Initial lowpass resolution (A) : This controls the resolution used in the first iteration.

  • Number of images for initial density generation: This controls the number of images used in the random density initialization.

  • GSFSC Split Resolution (A): This controls the resolution beyond which each half-map is considered independent. For smaller filaments, ideal values may lie lower (in the range of 12-15 Å) relative to the default of 20 Å.

The GSFSC split resolution is quite important for smaller filaments, and when running helical refinements without a volume input. If you notice strong oscillation in the FSC curves, especially at lower resolutions, you may want to set the split resolution to a lower number (e.g. 12-15 Å). With higher initial lowpass resolutions and lower numbers of images used for initial density generation, there is more potential for the initial density to introduce bias into the refinement, potentially causing the refinement to converge to an incorrect model. On the other hand, with lower initial lowpass resolutions and higher numbers of images used for initial density generation, the initial model may lack enough detail to reliably align particle images. While symmetric refinements are largely insensitive to these parameters (the imposition of correct symmetry overcomes most of these pitfalls), asymmetric refinements may require some experimentation to find the set of ideal parameters for initial model generation.

Using ab-initio reconstruction

Ab-initio reconstruction can sometimes be used to directly reconstruct a helical assembly from the particles. For some datasets, all projection views appear so similar that the reconstruction assigns them all to very similar poses, and hence suffers from strong preferred orientation. However, for other datasets in which views have more diversity, ab-initio reconstruction can work well for generating an initial model.

Some parameters that may increase the effectiveness of ab-initio reconstruction under these conditions are:

  • Initial/Final minibatch size: Increasing these to ≥400 and ≥1000, respectively, may increase the diversity of views used in each step of the SGD optimization (only if the dataset has enough diversity!)

  • Initial/Maximum resolution (Angstroms): Decreasing these to ≤20 and ≤10, respectively, will force the use of higher frequency information during optimization.

  • Initial structure lowpass (Fourier radius): Increasing this dramatically (i.e. from the default of 7 to between 30 and 60) can significantly impact the results, as it allows the use of higher-frequency information right from the start (which is often important for learning helical symmetry).

Reconstructions should ideally have orientation distributions that cover all or most azimuth angles. Alternately, if there are large gaps in coverage along the azimuth angle even after tweaking the above parameters, you may find more success in using an asymmetric helical refinement as detailed above.

Exploring the symmetry parameters in a reconstruction

The output of a Symmetry search utility will include a list of candidate symmetry parameters, both printed out to the streamlog, and also written to a .cs file in the job directory. They will be ranked by increasing mean squared error score, where lower scoring values are considered more plausible. The twist and rise values printed to the streamlog can be used as initial estimates for a helical refinement.

Using Helical Refinement (BETA) to refine symmetry parameters

Citations

Helical rise Δz \Delta z Δz (Å), and twist Δϕ \Delta \phi Δϕ (º)

Helical pitch, p pp (Å), number-of-subunits-per-full-turn, n nn, and hand h∈{+1,−1} h \in \{+1, -1\}h∈{+1,−1}

Conventionally, h=+1h = +1h=+1 corresponds to right-handed helices and h=−1h= -1h=−1 corresponds to left-handed helices

The helical pitch, number of subunits per turn, and hand contains the same information as the rise and twist. Note that ppp, nnn, and Δz\Delta zΔz are strictly positive, whereas hhh and Δϕ\Delta \phiΔϕ can take on either sign and thus serve to indicate the handedness of the helix. Fromp p p, nnn, and hhhvalues, the rise and twist can be calculated as:

Δϕ=360º  hn;    Δz=pn\Delta \phi = \frac{360º\;h}{n} ; \;\; \Delta z = \frac{p}{n}Δϕ=n360ºh​;Δz=np​

And similarly, from Δz\Delta zΔz and Δϕ\Delta \phiΔϕ, the pitch, number of subunits per turn, and hand can be calculated as:

n=360º∣Δϕ∣;    p=360ºΔz∣Δϕ∣;    h=sign(Δϕ)n = \frac{360º}{|\Delta \phi|} ; \; \; p = \frac{360º \Delta z}{|\Delta \phi|} ; \; \; h = sign(\Delta \phi) n=∣Δϕ∣360º​;p=∣Δϕ∣360ºΔz​;h=sign(Δϕ)

In some cases, helical assemblies will show more ambiguity in the number of subunits per turn than they do in pitch, making the (n,p,h)(n,p,h)(n,p,h) basis more useful to explore and decouple the space of symmetry parameters. On the other hand, the (Δϕ,Δz)(\Delta \phi, \Delta z)(Δϕ,Δz)basis more clearly represents the information needed to impose symmetry during a refinement, as it more directly relates to important low-level parameters such as the number of asymmetric units to perform symmetry-averaging over. For this reason, in CryoSPARC, both sets of parameters are used at different points in a helical processing workflow. In particular, helical refinements use the (Δϕ,Δz)(\Delta \phi, \Delta z)(Δϕ,Δz)basis (and this is what is stored internally), and the symmetry search utility uses either basis.

Note that even abiding by these conditions may produce multiple valid symmetry parameterizations. This is especially true when the helical assemble has an additional cyclic or dihedral point group symmetry associated with it. For example, EMPIAR-10019 with an additional C6 point group present has a helical rise (between six-fold symmetric layers) of 21.8  A˚21.8 \; Å21.8A˚. However, multiple twist values are valid due to the six-fold symmetry. Along with the reported twist value of 29.4º29.4º29.4º, a helical twist of −16(360º)+29.4º=−30.6º\frac{-1}{6}(360º) + 29.4º = -30.6º6−1​(360º)+29.4º=−30.6º will also perform the equivalent operation during reconstruction.

One way to do an asymmetric reconstruction is to launch an . This can be done by creating a Helical Refinement (BETA) job, and leaving the twist and rise parameters empty – this will mean that no symmetry will be applied during reconstruction. For such refinements, the initial density generated can impact the results. By default, helical refinement will generate an initial density by using the in-plane rotation (i.e. tilt) estimates of the filament picks, but with scrambled azimuth rotation angles. Two parameters that strongly determine how successful this method is are the following:

Once a refinement or ab-initio reconstruction has completed with a satisfactory orientation distribution, it can be passed to a job, which will search predefined range of symmetry parameters and will report the set of candidate symmetry pairs. Typically, the search ranges should be determined by inspection of the both the refined map, as well as of the 2D class averages, as these can often reveal insight as to which values of helical pitch are plausible. For more information and tips on using the search utility in your workflows, please refer to the linked job page, as well as the .

Once a pair of symmetry parameters is chosen for refinement, it can be passed to a job through the Helical twist estimate and Helical rise estimate parameters. Since these values are often not precise enough to bring the structure to very high resolution, it is recommended to enable refinement of the symmetry parameters by setting the Resolution to begin local searches of helical symmetry parameter to between 5-8 Angstroms.

[1] Doryen Bubeck. Making the most of symmetry. 2019. [ ]

[2] Shaoda He, Sjors H.W. Scheres. Helical reconstruction in RELION. Journal of Structural Biology, Volume 198, Issue 3, 2017, Pages 163-176, ISSN 1047-8477. [ ]

[3] Edward Egelman. Ambiguities in helical reconstruction. eLife 3:e04969, 2017. [ ]

asymmetric helical refinement
Symmetry search utility (BETA)
EMPIAR-10031 case study
helical refinement
http://www.cryst.bbk.ac.uk/embo2019/ppts/EMBO_2019.pdf
https://doi.org/10.1016/j.jsb.2017.02.003
https://doi.org/10.7554/eLife.04969.001
Figure A shows a continuous helix. Figure B shows a discrete helical lattice, parametrized by a twist (delta phi) and rise (delta z). Taken from [1].