A Workflow for Visualizing Protein Conformational Dynamics from Cryo-EM Maps Using Human Asparagine Synthetase

Summary

This protocol outlines how to perform 3D variability analysis (3DVA) in CryoSPARC and variability refinement in Phenix to analyze conformational heterogeneity in cryo-EM data. Using human asparagine synthetase as an example, it outlines how to generate variability maps, visualize principal components, and refine corresponding atomic models.

Abstract

This article presents a comprehensive workflow for analyzing conformational heterogeneity in single-particle cryo-electron microscopy (cryo-EM) data using 3D variability analysis (3DVA) in CryoSPARC and variability refinement in Phenix. The protocol describes how to compute variability maps along principal components derived from particle images, visualize principal modes of structural variation, and refine corresponding ensembles of atomic models against individual variability frames. The resulting outputs include a series of 3D variability maps, refined multi-model atomic structures representing discrete conformational states, and visualization files suitable for comparative structural analysis and movie generation. This workflow is broadly applicable to cryo-EM reconstructions in the resolution range (typically ~2.0-6.0 Å) where conformational heterogeneity can be meaningfully interpreted at the backbone and side-chain levels. It is particularly suited for capturing continuous domain motions, hinge-bending movements, and coordinated local rearrangements that are not easily resolved by discrete classification alone. The combined use of variability maps and refined model ensembles enables direct structural comparison of conformational extremes, supports quantitative analyses such as root-mean-square fluctuation (RMSF)  calculations, and provides a framework for linking dynamic structural features to biochemical or functional hypotheses. Human asparagine synthetase (ASNS) is used as an example to demonstrate the practical implementation of the workflow.

Introduction

Characterizing conformational heterogeneity in macromolecular structures is an increasingly important aspect of single-particle cryo-electron microscopy (cryo-EM) analysis. The goal of this protocol is to provide a practical workflow for extracting, visualizing, and refining dominant modes of structural variability directly from cryo-EM datasets using three-dimensional variability analysis (3DVA) in CryoSPARC¹ combined with variability refinement in Phenix². This article is intended to guide readers in determining whether this computational approach is appropriate for their data and in understanding the types of outputs ....

Protocol

1. Data preparation

1. Generation of consensus EM map
   1. Use single-particle cryo-electron microscopy (cryo-EM) to determine the consensus 3D structure of the macromolecular complex.
   2. Vitrify samples on cryo-grids and image using a high-voltage transmission electron microscope equipped with a direct electron detector.
   3. Process the collected movies to correct beam-induced motion, estimate imaging parameters, and select particle images.
   4. Perform iterative classification and refinement of the particles to generate a consensus EM map.
   5. Perform data processing using CryoSPARC (v4 or higher)

Results

This section presents representative outputs generated by the 3D variability analysis (3DVA) and variability refinement workflow and provides guidance on their interpretation. Human asparagine synthetase (ASNS) is used as an illustrative dataset to demonstrate successful outcomes, common artifacts, and comparative analyses enabled by the protocol. The complete set of 3DVA-derived variability maps and the corresponding refined models generated from the ASNS datasets have been deposited in the Zenodo public repository. Det.......

Discussion

Critical protocol steps and decision points
Effective use of this protocol depends on several key steps in the workflow. The quality of the particle stack and mask used for final cryo-EM refinement is critical, as 3D variability analysis (3DVA) operates directly on these inputs rather than on a post-processed map. Selection of the number of principal components (PCs) should reflect dataset size and expected complexity; too few components may underrepresent variability, whereas too many may introduc.......

Materials

List of materials used in this article

Name	Company	Catalog Number	Comments
CryoSPARC v4.4+	Structura Biotechnology	https://cryosparc.com	Software for cryo-EM processing
Phenix v1.24+	Phenix Consortium	https://phenix-online.org	Software for Variability refinement
UCSF Chimera	Resource for Biocomputing, Visualization, and Informatics, UCSF	https://www.cgl.ucsf.edu/chimera/	Software for Molecular Visualization
UCSF ChimeraX	Resource for Biocomputing, Visualization, and Informatics, UCSF	https://www.cgl.ucsf.edu/chimerax/	Software for Molecular Visualization
CPPTRAJ	AMBER molecular dynamics software suite	https://amberhub.chpc.utah.edu/cpptraj/	Python tools for RMSF calculation
Linux Workstation (in-house assembled GPU-based system)
Workstation Case	Supermicro	https://www.supermicro.com	Workstation case: Supermicro CSE-743TQ-1200B-SQ 1200W 4U Server Super Chassis
CPU	Intel	https://www.intel.com	CPU: XEON E5-2690V3 x2
Mother Board	Supermicro	https://www.supermicro.com	Mother Board: Supermicro Extended ATX DDR4 LGA 2011 Motherboard X10DAI-O
Memory	Crucial	https://www.crucial.com	Memory: Crucial 64GB Kit (16GBx4) DDR4 2133 (PC4-2133) DR x4 ECC Registered. 288-Pin Server Memory CT4K16G4RFD4213 / CT4C16G4RFD4213 x2
Hard drive	Western Digital	https://www.westerndigital.com	HDD: WD RE 4TB Enterprise Hard Drive
RAID CONTROLLER	LSI Logic	https://www.broadcom.com	RAID CONTROLLER: LSI LOGIC MegaRAID SAS 9271-8i Kit
GPU	NVIDIA	https://www.nvidia.com	Quadro M6000 graphics card with 24 GB GDDR5 memory