I am wondering if it is possible to add the functions of “align 3D volume” and rebox into particle subtraction. For particle subtraction, we sometimes only want a small part of map and would like it to be in the center and with smaller box.
Thanks,
Lei
Hi @rainfieldcn,
Thanks for the query – the recommended way to accomplish this in CryoSPARC is by performing the steps separately. For example, one possible workflow looks like:
EDIT: added two tweaks from @Hui (thanks!)
-
Volume alignment tools (calculate mask center of mass):
- Supply a volume, particles, and a mask covering the part of the map of interest; re-center to the center of this mask
- This job will report the center of mass of the mask and can be used in the next job.
-
Volume alignment tools (re-center):
- Supply the same volume and particles as in job 1. Do not input any mask. Use the calculated center of mask as the
3D coordinates of new center
.
- Supply the same volume and particles as in job 1. Do not input any mask. Use the calculated center of mask as the
-
Extract from micrographs (re-extract)
- Re-extract particles from Step 1 or 2 with “re-centering based on aligned shifts” enabled. Pick a box size large enough such that the particles are not cut off or too close to the edge. You can adjust the box size if it is too small after inspecting the output volume from step 4.
-
Homogeneous reconstruction only
- Use particles from step 3. Ensure that the output volume is not too close to the box edge
-
Create a mask around the area to subtract from the output volume from Step 4
-
Particle subtraction (subtract)
- Connect the inputs below and run the job:
- mask from Step 5,
- particles from Step 3, and
- volume from Step 4
- Connect the inputs below and run the job:
-
Downsample particles (crop)
- Set the real-space crop size to a reasonable number, based on the desired final refinement.
Best,
Michael
Hi Michael,
Thank you for reply!
Is Step #5 carried out in Step #2? Or is there another function for this cropping process?
Best,
Lei
Dear @rainfieldcn,
In step 2 you have the choice to choose a new box size. This box size must be big enough to cover both the centered region of interest, as well as the region that is being subtracted in step 4. In fact it should still be at least 2x the size of the full protein (as per our standard recommendation of the box size). This is because particle subtraction uses the CTF of the subtracted region, and if this region is near the edge of the box, signal can be cut off (see this reference for some discussion). It is only after particle subtraction occurs that we would recommend cropping to the final smallest box size, via downsample particles
Best,
Michael
Hi @mmclean,
Thank you for summarizing the necessary steps here, but I noticed a couple of potential issues of the workflow that i want to share here.
- The volume output of step1 is masked. So to get an unmasked, but translated volume for subtraction, we have to repeat step1, but this time without a mask, using the calculated number of mask center. Can cryosparc also output the unmasked volume here so we don’t have to repeat?
- If the particle box is not large enough, then the shifted volume from step 1 will be truncated. In this case, we have to do a reconstruction job using particles from step2. Therefore the box size of the particles in this step has to be big enough.
Best,
Hui
Dear @Hui,
Thanks for the report. You are right on both accounts, my apologies for the error!
Good request, you’re right that currently you have to create a new job without the mask supplied and with the manually entered re-centering coordinates. We have this recorded now to allow volume alignment tools to allow outputting the unmasked volume.
This is correct. I have updated my comment to add these two notes.
Best,
Michael