Name: Author Spotlight: Optimizing Cryo-EM Analysis with CryoSieve for Enhanced Particle Selection Efficiency
Uploaded: 2024-05-10T14:35:00
Duration: 1 min 45 s
Description: Read this Scientific Article Protocol about Author Spotlight: Optimizing Cryo-EM Analysis with CryoSieve for Enhanced Particle Selection Efficiency at JoVE.com

This work was supported by Shenzhen Academy of Research and Translation (to M.H.), the Advanced Innovation Center for Structural Biology (to M.H.), the Beijing Frontier Research Center for Biological Structure (to M.H.), the National Key R&#38;D Program of China (No.2021YFA1001300) (to C.B.), the National Natural Science Foundation of China (No.12271291) (to C.B.), and the National Natural Science Foundation of China (No.12071244) (to Z.S.).

1. Installation
<ol>
	<li>Check and configure GPU-acceleration environment
	<ol>
		<li>Open the terminal and enter the command: nvidia-smi. Make sure the command successfully displays all information about the GPU card(s) and the CUDA version is higher than 10.2. Execute the command: conda -V to check if Conda is installed (Supplementary Figure 1).</li>
	</ol>
	</li>
	<li>Configure virtual environment
	<ol>
		<li>Enter the following command to set up the virtual environment, replacing CRYOSIEVE_ENV with your desired environment name: conda create -n CRYOSIEVE_ENV python=3.8 cudatoolkit=10.2 cupy=10.0 pytorch=1.10 -c pytorch -c conda-forge. Wait for a few minutes until the environment is successfully configured (Supplementary Figure 2). 
		NOTE: Users have the flexibility to modify the environment name as needed. The provided command is specific to CUDA 10.2. If a different CUDA version is desired, adjust the version number for cudatoolkit.</li>
	</ol>
	</li>
	<li>Install CryoSieve
	<ol>
		<li>Activate the environment by executing the command: conda activate CRYOSIEVE_ENV. Install the software by running: pip install cryosieve or conda install -c mxhulab cryosieve (Supplementary Figure 3). Enter cryosieve -h and ensure that help information is correctly displayed (Supplementary Figure 4).</li>
	</ol>
	</li>
</ol>
2. Particle sieving
<ol>
	<li>Retrieve the data
	<ol>
		<li>Download the EMPIAR-10097 final stack dataset from EMPIAR (see Table of Materials). Download the star file, the mask file (mask.mrc) and the initial model (for the re-estimate step; initial.mrc) from Github (see Table of Materials). Place all these files in a folder together (Supplementary Figure 5). 
		NOTE: The repository at https://github.com/mxhulab/cryosieve-demos utilizes Git Large File Storage (Git LFS). Installing Git LFS is essential for cloning the entire repository. Alternatively, access the file via the GitHub link and click the&#160;Download raw file button to download an individual file.</li>
	</ol>
	</li>
	<li>Process particle sieving
	<ol>
		<li>Open the terminal and use the command: cd FILEPATH to navigate to the folder where the dataset is located. Activate the Conda environment by: conda activate CRYOSIEVE_ENV.</li>
		<li>Enter the following command to start our particle sieving experiment: cryosieve --reconstruct_software relion_reconstruct --postprocess_software relion_postprocess --i T40_HA_130K-Equalized_run-data_CryoSPARC_refined.star --o output/ --mask mask.mrc --angpix 1.3099979 --num_iters 10 --frequency_start 40 --frequency_end 3 --retention_ratio 0.8 --sym C3 --num_gpus 1 --balance (Supplementary Figure 5). During the execution, the terminal will display the output logs for each iteration. 
		NOTE: Detailed instructions for each option can be found in Supplementary File 1. Processing time and the minimum requirements for execution are detailed in Supplementary File 2. T40_HA_130K-Equalized_run-data_CryoSPARC_refined.star was refined by CryoSPARC from T40_HA_130K-Equalized_run-data.star (downloaded from EMPIAR) to mitigate the effects brought about by advancements in orientation estimation techniques.</li>
	</ol>
	</li>
</ol>
3. Finding the optimal iteration
<ol>
	<li>Check resolutions
	<ol>
		<li>Use the command: grep &#34;+ FINAL RESOLUTION:&#34; output/_postprocess*.txt to print resolution results for the 10 iterations of sieving (Figure 1). Since the particle stack filtered in the 7th iteration has the highest resolution with the fewest particles, it is likely to provide the optimal result. 
		NOTE: To avoid unintentional information transfer from discarded to retained particles15&#160;and to ensure that the particle stack posts the 7th iteration is indeed optimal, users are required to execute a re-estimation step for nearby iterations. In this protocol, iterations 4, 5, 6, 7, and 8&#160;are subjected to verification.</li>
	</ol>
	</li>
	<li>Import sieved particles
	<ol>
		<li>Open the CryoSPARC web interface and follow these steps: Enter a workspace and click on the Builder button at the top right of the panel. In the panel, select and click on the Import Particle Stack option. In the Parameters section of the Particle Stack Import panel, specify the Particle meta path as the _iter{n}.star file located in the output folder of the completed results and the particle data path to the folder where the mrcs file is stored. Click the Queue Job button, then click the Queue button to initiate the process. Use the same way to import the remaining iterations that need re-estimation (Supplementary Figure 6A).</li>
	</ol>
	</li>
	<li>Import initial model
	<ol>
		<li>Click on the Builder button at the top right of the panel. In the panel, select and click on the Import 3D Volumes option.</li>
		<li>Specify the Volume data path as the initial.mrc file. Click the Queue Job button, then click the Queue button to initiate the process (Supplementary Figure 6B). 
		NOTE: The initial model can also be generated through ab initio reconstruction (Supplementary File 3).</li>
	</ol>
	</li>
	<li>Homogeneous refinement (Build job)
	<ol>
		<li>Click on the Builder button at the top right of the panel. In the panel, select and click on the Homogeneous Refinement option. 
		NOTE: Non-uniform refinement is also applicable.</li>
	</ol>
	</li>
	<li>Homogeneous refinement (Import particles)
	<ol>
		<li>In the main panel on the left, open the job for importing the particle stack of the 5th iteration (or the desired iteration). Drag the imported particles module from the right side of the main panel and drop it into the Particle stacks section of the Builder on the right. Close the Import Particle Stack job by clicking on the red X in the top right corner of the main panel.</li>
		<li>Open the job for importing 3D volumes. Drag the imported volumes module from the right side of the main panel and drop it into the Initial volume section of the Builder on the right.</li>
	</ol>
	</li>
	<li>Homogeneous refinement (Modify the parameters)
	<ol>
		<li>Under the Parameters fold, locate the Symmetry option and set it to C3. Find the Force re-do GS split option and disable it. Click the Queue Job button, then click the Queue button to initiate the Homogeneous Refinement. Perform Homogeneous Refinement for the remaining iterations using the same method (Supplementary Figure 6C-D). 
		NOTE: The Force re-do GS split option is critical. Disabling this option ensures that CryoSPARC retains the gold-standard split given by the star file, avoiding overfitting. A detailed rationale for disabling Force Re-do GS Split can be found in Supplementary File 4.</li>
	</ol>
	</li>
	<li>Wait until all the jobs finish running to obtain the results. Based on the results, it is confirmed that the particle stack filtered in the 6th iteration is&#160;the actual optimal result. 
	NOTE: It is normal for the results obtained to have minor random deviations from the results provided in this protocol. These deviations do not affect the overall conclusion.</li>
</ol>

Enhancing Cryo-EM Map Quality with CryoSieve Particle Selection

<ol>
	<li>Bai, X. C., Fernandez, I. S., Mcmullan, G., Scheres, S. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Ribosome+structures+to+near-atomic+resolution+from+thirty+thousand+cryo-em+particles.">Ribosome structures to near-atomic resolution from thirty thousand cryo-em particles.</a> elife. 2, 00461 (2013).</li><li>Campbell, M. G., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Movies+of+ice-embedded+particles+enhance+resolution+in+electron+cryo-microscopy.">Movies of ice-embedded particles enhance resolution in electron cryo-microscopy.</a> Structure. 20 (11), 1823-1828 (2012).</li><li>Li, X., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Electron+counting+and+beam-induced+motion+correction+enable+near-atomic-resolution+single-particle+cryo-em.">Electron counting and beam-induced motion correction enable near-atomic-resolution single-particle cryo-em.</a> Nat Meth. 10 (6), 584-590 (2013).</li><li>Grant, T., Rohou, A., Grigorieff, N. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Cis+tem,+user-friendly+software+for+single-particle+image+processing.">Cis tem, user-friendly software for single-particle image processing.</a> eLife. 7, e35383 (2018).</li><li>Scheres, S. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Relion:+Implementation+of+a+bayesian+approach+to+cryo-em+structure+determination.">Relion: Implementation of a bayesian approach to cryo-em structure determination.</a> J Str Biol. 180 (3), 519-530 (2012).</li><li>Punjani, A., Rubinstein, J. L., Fleet, D. J., Brubaker, M. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Cryosparc:+Algorithms+for+rapid+unsupervised+cryo-em+structure+determination.">Cryosparc: Algorithms for rapid unsupervised cryo-em structure determination.</a> Nat Meth. 14 (3), 290-296 (2017).</li><li>K&#252;hlbrandt, W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+resolution+revolution.">The resolution revolution.</a> Science. 343 (6178), 1443-1444 (2014).</li><li>Dubochet, J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Cryo-electron+microscopy+of+vitrified+specimens.">Cryo-electron microscopy of vitrified specimens.</a> Quart Rev Biophys. 21 (2), 129-228 (1988).</li><li>Wagner, T., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Sphire-cryolo+is+a+fast+and+accurate+fully+automated+particle+picker+for+cryo-EM.">Sphire-cryolo is a fast and accurate fully automated particle picker for cryo-EM.</a> Comm Biol. 2 (1), 218 (2019).</li><li>Bepler, T., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Positive-unlabeled+convolutional+neural+networks+for+particle+picking+in+cryo-electron+micrographs.">Positive-unlabeled convolutional neural networks for particle picking in cryo-electron micrographs.</a> Nat Meth. 16 (11), 1153-1160 (2019).</li><li>Wang, F., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Deeppicker:+A+deep+learning+approach+for+fully+automated+particle+picking+in+cryo-em.">Deeppicker: A deep learning approach for fully automated particle picking in cryo-em.</a> J Str Biol. 195 (3), 325-336 (2016).</li><li>Heimowitz, A., And&#233;n, J., Singer, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Apple+picker:+Automatic+particle+picking,+a+low-effort+cryo-em+framework.">Apple picker: Automatic particle picking, a low-effort cryo-em framework.</a> J Str Biol. 204 (2), 215-227 (2018).</li><li>Glaeser, R. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=How+good+can+single-particle+cryo-em+become?+What+remains+before+it+approaches+its+physical+limits.">How good can single-particle cryo-em become? What remains before it approaches its physical limits.</a> Ann Rev Biophys. 48, 45-61 (2019).</li><li>Diiorio, M. C., Kulczyk, A. W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+robust+single-particle+cryo-electron+microscopy+(cryo-em)+processing+workflow+with+cryosparc,+relion,+and+scipion.">A robust single-particle cryo-electron microscopy (cryo-em) processing workflow with cryosparc, relion, and scipion.</a> J Vis Exp. (179), e63387 (2022).</li><li>Zhu, J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+minority+of+final+stacks+yields+superior+amplitude+in+single-particle+cryo-em.">A minority of final stacks yields superior amplitude in single-particle cryo-em.</a> Nat Comm. 14 (1), 7822 (2023).</li><li>Zhou, Y., Moscovich, A., Bendory, T., Bartesaghi, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Unsupervised+particle+sorting+for+high-resolution+single-particle+cryo-em.">Unsupervised particle sorting for high-resolution single-particle cryo-em.</a> Inv Probl. 36 (4), 044002 (2020).</li><li>M&#233;ndez, J., Garduno, E., Carazo, J. M., Sorzano, C. O. S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Identification+of+incorrectly+oriented+particles+in+cryo-em+single+particle+analysis.">Identification of incorrectly oriented particles in cryo-em single particle analysis.</a> J Str Biol. 213 (3), 107771 (2021).</li><li>Tan, Y. Z., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Addressing+preferred+specimen+orientation+in+single-particle+cryo-em+through+tilting.">Addressing preferred specimen orientation in single-particle cryo-em through tilting.</a> Nat Meth. 14 (8), 793-796 (2017).</li><li>Rosenthal, P. B., Henderson, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Optimal+determination+of+particle+orientation,+absolute+hand,+and+contrast+loss+in+single-particle+electron+cryomicroscopy.">Optimal determination of particle orientation, absolute hand, and contrast loss in single-particle electron cryomicroscopy.</a> J Mol Biol. 333 (4), 721-745 (2003).</li><li>Kaur, S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Local+computational+methods+to+improve+the+interpretability+and+analysis+of+cryo-em+maps.">Local computational methods to improve the interpretability and analysis of cryo-em maps.</a> Nat Comm. 12 (1), 1240 (2021).</li><li>Stagg, S. M., Noble, A. J., Spilman, M., Chapman, M. S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Reslog+plots+as+an+empirical+metric+of+the+quality+of+cryo-em+reconstructions.">Reslog plots as an empirical metric of the quality of cryo-em reconstructions.</a> J Str Biol. 185 (3), 418-426 (2014).</li><li>Vargas, J., Ot&#243;n, J., Marabini, R., Carazo, J. M., Sorzano, C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Particle+alignment+reliability+in+single+particle+electron+cryomicroscopy:+A+general+approach.">Particle alignment reliability in single particle electron cryomicroscopy: A general approach.</a> Sci Rep. 6 (1), 21626 (2016).</li><li>Vargas, J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Particle+quality+assessment+and+sorting+for+automatic+and+semiautomatic+particle-picking+techniques.">Particle quality assessment and sorting for automatic and semiautomatic particle-picking techniques.</a> J Str Biol. 183 (3), 342-353 (2013).</li><li>Vargas, J., Melero, R., Gomez-Blanco, J., Carazo, J. -. M., Sorzano, C. O. S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Quantitative+analysis+of+3d+alignment+quality:+Its+impact+on+soft-validation,+particle+pruning+and+homogeneity+analysis.">Quantitative analysis of 3d alignment quality: Its impact on soft-validation, particle pruning and homogeneity analysis.</a> Sci Rep. 7 (1), 6307 (2017).</li><li>Carragher, B., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Current+outcomes+when+optimizing+'standard'sample+preparation+for+single-particle+cryo-em.">Current outcomes when optimizing 'standard'sample preparation for single-particle cryo-em.</a> J Microsc. 276 (1), 39-45 (2019).</li><li>Drulyte, I., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Approaches+to+altering+particle+distributions+in+cryo-electron+microscopy+sample+preparation.">Approaches to altering particle distributions in cryo-electron microscopy sample preparation.</a> Acta Crystallographica Sec D: Str Biol. 74 (6), 560-571 (2018).</li><li>Glaeser, R. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=How+good+can+cryo-em+become.">How good can cryo-em become.</a> Nat Meth. 13 (1), 28-32 (2016).</li><li>Kim, L. Y., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Benchmarking+cryo-em+single+particle+analysis+workflow.">Benchmarking cryo-em single particle analysis workflow.</a> Front Mol Biosci. 5, 50 (2018).</li><li>Weissenberger, G., Henderikx, R. J., Peters, P. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Understanding+the+invisible+hands+of+sample+preparation+for+cryo-em.">Understanding the invisible hands of sample preparation for cryo-em.</a> Nat Meth. 18 (5), 463-471 (2021).</li></ol>

All other authors declare no competing interests.

Cryo-EM stands as a pivotal technique for elucidating the structures of biological molecules. In this process, after data collection via microscopy, particle extraction from micrographs is essential, followed by their classification in multiple stages to compile the final stack. A common challenge is the predominance of damaged or undesirably conformed particles, underscoring the need for repeated particle selection to attain high-resolution density maps. This makes particle selection a critical step in cryo-EM SPA for achieving high-quality density maps. Existing particle selection techniques include the statistical non-tilt validation algorithm22, the z-score-based approach23, and the angular accuracy estimation method24.
CryoSieve emerges as a valuable tool in this context, adept at eliminating a significant number of extraneous particles from the final stack. This reduction not only enhances the reconstruction&#39;s computational efficiency but also streamlines the process. It offers a comprehensive suite for particle selection, where the extent of particle discard and the consequent improvement in resolution largely hinge on the initial data quality and the methodologies employed in data processing.
In this manuscript, we have presented a complete workflow of particle sieving using the real case dataset of influenza hemagglutinin trimer (EMPIAR entry: 10097). The steps covered and discussed here can be summarized as particle sieving and pose re-estimation. The final 3D reconstructed volume achieved a resolution of 3.62 &#197;, and side chains in alpha-helices were clearer in the post-processed volume compared to the published density map.
CryoSieve is an open-source method which is available on GitHub (https://github.com/mxhulab/cryosieve). A detailed tutorial can also be found on its homepage. Users can install and use it by following the tutorial. Additionally, two modules, cryosieve-csrefine and cryosieve-csrhbfactor, are provided. The cryosieve-csrefine module is specifically crafted to automate the sequential execution of various operations within CryoSPARC (Supplementary File 5). These operations include importing particle stacks and conducting ab initio, homogeneous refinement, or non-uniform refinement jobs. On the other hand, the cryosieve-csrhbfactor module is designed to automate the determination of the Rosenthal-Henderson B-factor by leveraging the capabilities of cryosieve-csrefine (Supplementary File 6).
Presently, this method&#39;s application is confined to single conformation scenarios. Consequently, in instances where particles represent multiple conformations, their capabilities are limited. Users are advised to initially engage in 3D classification to segregate particles of disparate conformations before employing it for refined particle selection. Moreover, although the method demonstrates proficiency in filtering out over 50% of particles from the final stack, the origins of these discarded particles and the underlying reasons for their negligible contribution to reconstruction quality remain unclear. This gap in understanding necessitates additional research to comprehensively address and potentially rectify this limitation.
There are three possible existing methods of particle sorting or particle sieving. First of all, cisTEM4&#160;can report a score for each single particle image after 3D refinement. Users could sort particles using the cisTEM score to discard particles. The angular graph consistency (AGC) approach17&#160;is also a method to discard misaligned particles. Furthermore, the non-alignment classification5&#160;is a traditional way to discard particles using 3D classification. We compared the quality of particles retained by these methods with CryoSieve and found that the retained particles of CryoSieve are of higher quality15. The method presented here significantly outperforms alternative methods and achieves the smallest number of particles at the same resolution.
As demonstrated in the result, the majority of particles in a cryo-EM final stack do not contribute to density map reconstruction. In other words, among all particles gathered during image acquisition, only a select few, namely the finest subset, actually contribute to the final reconstruction. Consequently, the ratio of this final subset to the total number of collected particles could serve as a quantitative metric for assessing sample quality. The higher this ratio, the better the sample quality. Despite technical advancements that have made cryo-EM more accessible to structural biologists, sample preparation remains a major bottleneck in the workflow. Scientists and engineers are thus focusing their efforts on this challenge25. In single-particle analysis (SPA), sample preparation consists of two crucial steps: sample optimization and grid preparation. The former involves purifying the specimen while maintaining its optimal biochemical state. The latter entails preparing the sample for analysis in the microscope, including chemical or plasma treatment of the grid, sample deposition, and vitrification. Numerous techniques have been proposed to address macromolecular instability, but the efficacy of one approach over another depends on the sample&#39;s characteristics25,26. Currently, grid preparation results are heavily influenced by the user&#39;s expertise and experience, which can make the process time-consuming and challenging27,28. The numerous variables encountered in sample and grid preparation pose challenges in establishing cause-and-effect relationships, as researchers can only assess the sample at the molecular level using the microscope. As a result, quantitative statistics from comparisons of different sample and grid preparation protocols are still lacking, and a systematic approach is necessary to investigate trends and comprehend the fundamental mechanisms of sample behavior29.

Cryogenic electron microscopy (cryo-EM) single-particle analysis (SPA) has become a dominant method to determine high-resolution three-dimensional density maps of biological macromolecules. Due to a series of technological innovations1,2,3,4,5,6, named resolution revolution7, cryo-EM has the capability to determine the structures of biological macromolecules with up to atomic resolution at an unprecedented rate. This breakthrough marks the beginning of a new era in molecular insights, overtaking X-ray crystallography as the predominant technique and answering longstanding biological questions.
Cryo-EM SPA diverges from X-ray crystallography by not requiring the crystallization of biological macromolecules. Instead, a solution containing the target biological macromolecules is rapidly frozen in vitreous ice. It is then imaged with an electron beam to produce a series of micrographs, bypassing the need for crystallization8. Subsequently, particle-picking algorithms are utilized to extract individual raw particles from these micrographs4,9,10,11,12. As cryo-EM does not depend on crystallization, it is natural that extracted particles are predominantly damaged or in undesired conformational states, necessitating multiple rounds of particle selection to achieve a high-resolution density map. In cryo-EM SPA image processing, particle selection is therefore crucial for obtaining high-resolution density maps13.
In cryo-EM SPA, standard particle selection methods include two-dimensional (2D) and three-dimensional (3D) classification14. 2D classification categorizes particles into a predefined number of groups, yielding an average image and an estimated 2D resolution for each class. Researchers can then visually inspect these classes, removing particles from lower resolution groups to use the remaining ones in reconstructions aimed at achieving higher resolution. Once particle poses are established using refinement algorithms, researchers will proceed with 3D classification, clustering particles into multiple classes. This enables visual inspection of the reconstructed density map for each class, allowing for the exclusion of undesirable particles, such as those from undesired conformations. Following multiple rounds of classification, a final stack comprising relatively high-quality particles is obtained. These final stacks are instrumental in producing atomic or near-atomic resolution density maps.
Zhu and her colleagues have demonstrated that further particle selection can be conducted on these final stacks15. CryoSieve15, an innovative iterative method for particle selection, can be applied to enhance the quality of the final density map by significantly reducing the number of particles. While other particle sorting criteria and software, such as the normalized cross-correlation (NCC) method16, the angular graph consistency (AGC) approach17, and non-alignment classification5, are currently in use within the field, this method has been shown to outperform these algorithms in terms of effectiveness.
In this study, we present a detailed guide to the entire process. As a case study, we applied this new method to the dataset of the influenza hemagglutinin trimer (EMPIAR entry: 10097)18, which includes 130,000 particles in its final stack. Our procedure successfully discarded about 73.8% of the particles from the final stack of this dataset, improving the resolution of the reconstructed density map from 4.11 &#197; to 3.62 &#197;. In addition to the influenza hemagglutinin trimer, results from multiple datasets are presented in earlier publication15, showcasing a variety of resolutions and molecular weights of biomolecules.

<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<tbody>
		<tr>
			<td>CryoSPARC</td>
			<td>Structura Biotechnology Inc. Toronto, Canada</td>
			<td></td>
			<td>CryoSPARC (Cryo-EM Single Particle Ab-Initio Reconstruction and Classification) is a state of the art HPC software solution for complete processing of single-particle cryo-electron microscopy (cryo-EM) data. CryoSPARC is useful for solving cryo-EM structures of membrane proteins, viruses, complexes, flexible molecules, small particles, phase plate data and negative stain data.</td>
		</tr>
		<tr>
			<td>EMPIAR-10097 Dataset</td>
			<td></td>
			<td>https://ftp.ebi.ac.uk/empiar/world_availability/10097/data/Particle-Stack/T40_HA_130K-Equalized-Particle-Stack.mrcs</td>
			<td>This dataset comprises single-particle cryo-EM data of the Influenza Hemagglutinin trimer, characterized by its highly preferred orientation, collected using a 40-degree tilted collection strategy.</td>
		</tr>
		<tr>
			<td>initial.mrc</td>
			<td></td>
			<td>https://github.com/mxhulab/cryosieve-demos/tree/master/EMPIAR-10097</td>
			<td></td>
		</tr>
		<tr>
			<td>mask.mrc</td>
			<td></td>
			<td>https://github.com/mxhulab/cryosieve-demos/tree/master/EMPIAR-10097</td>
			<td></td>
		</tr>
		<tr>
			<td>RELION</td>
			<td></td>
			<td>4.0-beta-2</td>
			<td>RELION (REgularised LIkelihood OptimisatioN) is an open-source software for cryo-electron microscopy (cryo-EM) data processing, particularly for refining macromolecular structures. Utilizing a Bayesian approach, it excels in separating signal from noise, enabling high-resolution structure determination. RELION supports single-particle analysis, tomography, and sub-tomogram averaging, and has become widely used in structural biology due to its effectiveness and user-friendly interface.</td>
		</tr>
		<tr>
			<td>T40_HA_130K-Equalized_run-data_CryoSPARC_refined.star</td>
			<td></td>
			<td>https://github.com/mxhulab/cryosieve-demos/tree/master/EMPIAR-10097</td>
			<td>Metadata file for the final stack of particles from EMPIAR-10097</td>
		</tr>
	</tbody>
</table>

enhancing density maps by removing the majority of particles in single particle cryogenic electron microscopy final stacks

Over the past decade, advancements in technology and methodology within the field of cryogenic electron microscopy (cryo-EM) single-particle analysis (SPA) have substantially improved our capacity for high-resolution structural examination of biological macromolecules. This advancement has ushered in a new era of molecular insights, replacing X-ray crystallography as the dominant method and providing answers to longstanding questions in biology. Since cryo-EM does not depend on crystallization, which is a significant limitation of X-ray crystallography, it captures particles of varying quality. Consequently, the selection of particles is crucial, as the quality of the selected particles directly influences the resolution of the reconstructed density map. An innovative iterative approach for particle selection, termed CryoSieve, significantly improves the quality of reconstructed density maps by effectively reducing the number of particles in the final stack. Experimental evidence shows that this method can eliminate the majority of particles in final stacks, resulting in a notable enhancement in the quality of density maps. This article outlines the detailed workflow of this approach and showcases its application on a real-world dataset.

Author Spotlight: Optimizing Cryo-EM Analysis with CryoSieve for Enhanced Particle Selection Efficiency

Enhancing Density Maps by Removing the Majority of Particles in Single Particle Cryogenic Electron Microscopy Final Stacks

In this protocol, we utilized the influenza hemagglutinin trimer dataset (EMPIAR entry: 10097) as a demonstration of the efficacy of this process. Due to the preferred orientation of the sample, data acquisition required tilting at 40&#176;. The protein exhibits C3 symmetry and has a molecular weight of 150 kDa.
We have implemented the protocol described earlier to process the final particle stack. It progressively removed 20% of the particles in each iteration, resulting in a retention ratio of 80.0%, 64.0%, 51.2%, and so on. As depicted in Figure 1&#160;and Figure 2, the resolution of the retained particles initially improved but eventually decreased. Among the iterations, the 6th iteration was identified as the most optimal subset, containing the fewest particles yet achieving the highest resolution. Our algorithm successfully identified a subset of particles comprising only 26.2% of the original stack, resulting in an improved resolution from 4.19 &#197; to 3.62 &#197; (re-estimated by CryoSPARC), shown in Figure 2. Furthermore, density maps before and after using CryoSieve were compared in Figure 3. Model-to-map Fourier Shell Correlation (FSC) curve and half-maps FSC curve of the reconstructed density maps before and after the method are also shown (Figure 3A-B). Raw density maps and sharp density maps obtained were also compared, with the equivalent contour level applied (Figure 3C). The side chains of sharp density maps were compared, showing the enhancement of reconstructed density maps. The estimated Rosenthal-Henderson B-factor was also adopted for the criteria of particle quality19. After removing the majority of particles in the final stack, the Rosenthal-Henderson B-factor raised from 226.9 &#197;2 to 146.2 &#197;2&#160;(Figure 3D). Local resolution, local B-factor20, and ResLog21&#160;were also utilized for comparison, indicating that CryoSieve indeed enhances both the quality of the density maps and the particles (Figure 4).
<img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/66617/66617fig01.jpg" /> 
Figure 1: Resolutions of each iteration. Resolutions that have been reported are highlighted in red boxes. <a href="https://www.jove.com/files/ftp_upload/66617/66617fig01large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 2" class="xfigimg" src="/files/ftp_upload/66617/66617fig02.jpg" /> 
Figure 2: Resolutions of each iteration. Resolutions identified by homogeneous refinement jobs are highlighted in red boxes. <a href="https://www.jove.com/files/ftp_upload/66617/66617fig02large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 3" class="xfigimg" src="/files/ftp_upload/66617/66617fig03.jpg" /> 
Figure 3: Density maps. (A) Comparison of model-to-map FSC curve of reconstructed density maps before and after using CryoSieve. The y-axis represents FSC, while the x-axis represents resolution. The red dashed line marks the threshold of 0.5 for the FSC. The vertical dashed line illustrates the resolution of the density maps obtained under a threshold of 0.5. (B) Half-maps FSC curve were obtained from reconstructed density maps before and after using CryoSieve via CryoSPARC. The y-axis represents FSC, while the x-axis represents resolution. (C) Raw density maps and sharp density maps were shown for both the CryoSieve-retained particles and the complete set of particles in the final stacks. The equivalent contour level of 0.65 was applied for raw density maps. The equivalent contour level of 0.84 was applied for sharp density maps. Sharp density maps were directly obtained by CryoSPARC. The sharp density maps were auto-postprocessed, first FSC-weighted (based on FSCs given by CryoSPARC). Then, the B-factor was sharpened using the auto-determined B-factors (232.0 &#197;2 for all particles in the final stack and 160.8 &#197;2 for CryoSieve). The side chains in the sharp density maps were compared, incorporating atomic models for reference. Red arrows highlight the improved regions. (D) The estimated Rosenthal-Henderson B-factor was shown for both the CryoSieve-retained particles and the complete set of particles in the final stacks. The y-axis represents the number of particles used, and the x-axis represents the reciprocal of the square of the resolution. Moving from top to bottom, each point represents half the particles of the previous one. The resolutions were determined by refinement. B-factors were determined using a least-squares approximation of the measured points, as shown by the fitting curves. The estimated Rosenthal and Henderson&#39;s B-factors are indicated in the legends: orange represents particles retained by CryoSieve, while blue denotes all particles in the final stack. <a href="https://www.jove.com/files/ftp_upload/66617/66617fig03large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 4" class="xfigimg" src="/files/ftp_upload/66617/66617fig04.jpg" /> 
Figure 4: Comparison of various metrics of density maps. (A) Comparison of local resolution maps before and after using CryoSieve obtained by CryoSPARC. The local resolution ranges between 7 &#197; (red) and 3.5 &#197; (blue). (B) Comparison of density maps before and after using CryoSieve, colored with the local B-factor map obtained by LocBFactor using a resolution range of [20-3.5] &#197;. (C), Comparison of ResLog plots before and after using CryoSieve obtained by CryoSPARC. <a href="https://www.jove.com/files/ftp_upload/66617/66617fig04large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
Supplementary Figure 1: Using commands nvidia-smi and conda -V to verify the prerequisites. If the prerequisites are met, typing the command nvidia-smi will display the GPU driver version, the CUDA version, and the status of the GPU cards. Similarly, entering the command conda -V should correctly display the installed version of Conda. <a href="https://www.jove.com/files/ftp_upload/66617/Supplementary_Figure_1.zip" target="_blank">Please click here to download this File.</a>
Supplementary Figure 2: The process of creating new GPU-acceleration environments. The screen displays the output generated by the command used to create the Conda environment. <a href="https://www.jove.com/files/ftp_upload/66617/Supplementary_Figure_2.zip" target="_blank">Please click here to download this File.</a>
Supplementary Figure 3: Installation of CryoSieve in the GPU-acceleration environment. After activating the newly created Conda environment, the screen displays the output resulting from executing the command to install CryoSieve using Pip. <a href="https://www.jove.com/files/ftp_upload/66617/Supplementary_Figure_3.zip" target="_blank">Please click here to download this File.</a>
Supplementary Figure 4: Help information. <a href="https://www.jove.com/files/ftp_upload/66617/Supplementary_Figure_4.zip" target="_blank">Please click here to download this File.</a>
Supplementary Figure 5: Running process. Upon executing CryoSieve through the command line, the screen then displays information regarding the running process. <a href="https://www.jove.com/files/ftp_upload/66617/Supplementary_Figure_5.zip" target="_blank">Please click here to download this File.</a>
Supplementary Figure 6: The configuration of CryoSPARC&#39;s jobs. (A) Import particle stack. (B) Import 3D volumes. (C-D) Homogeneous refinement. <a href="https://www.jove.com/files/ftp_upload/66617/Supplementary_Figure_6.zip" target="_blank">Please click here to download this File.</a>
Supplementary File 1: Options of CryoSieve. <a href="https://www.jove.com/files/ftp_upload/66617/supplementary_1.pdf" target="_blank">Please click here to download this File.</a>
Supplementary File 2: Processing time and minimal requirement for running Cryosieve. <a href="https://www.jove.com/files/ftp_upload/66617/supplementary_2.pdf" target="_blank">Please click here to download this File.</a>
Supplementary File 3: Generation of initial model by CryoSPARC. <a href="https://www.jove.com/files/ftp_upload/66617/supplementary_3.pdf" target="_blank">Please click here to download this File.</a>
Supplementary File 4: Rationale for disabling force re-do GS split. <a href="https://www.jove.com/files/ftp_upload/66617/supplementary_4.pdf" target="_blank">Please click here to download this File.</a>
Supplementary File 5: Options of cryosieve-csrefine. <a href="https://www.jove.com/files/ftp_upload/66617/supplementary_5.pdf" target="_blank">Please click here to download this File.</a>
Supplementary File 6: Options of cryosieve-csrhbfactor. <a href="https://www.jove.com/files/ftp_upload/66617/supplementary_6.pdf" target="_blank">Please click here to download this File.</a>

Read this Scientific Article Protocol about Author Spotlight: Optimizing Cryo-EM Analysis with CryoSieve for Enhanced Particle Selection Efficiency at JoVE.com

An advanced particle selection method for cryo-EM, namely CryoSieve, improves density map resolution by removing a majority of particles in final stacks, as demonstrated through its application on a real-world dataset.

Over the past decade, advancements in technology and methodology within the field of cryogenic electron microscopy (cryo-EM) single-particle analysis ...