Dynamics in cryo EM reconstructions visualized with maximum-likelihood derived variance maps

CryoEM data capture the dynamic character associated with biological macromolecular assemblies by preserving the various conformations of the individual specimens at the moment of flash freezing. Regions of high variation in the data set are apparent in the image reconstruction due to the poor density that results from the lack of superposition of these regions. These observations are qualitative and, to date, only preliminary efforts have been made to quantitate the heterogeneity in the ensemble of particles that are individually imaged. We developed and tested a quantitative method for simultaneously computing a reconstruction of the particle and a map of the space-varying heterogeneity of the particle based on an entire data set. The method uses a maximum likelihood algorithm that explicitly takes into account the continuous variability from one instance to another instance of the particle. The result describes the heterogeneity of the particle as a variance to be plotted at every voxel of the reconstructed density. The test, employing time resolved data sets of virus maturation, not only recapitulated local variations obtained with difference map analysis, but revealed a remarkable time dependent reduction in the overall particle dynamics that was unobservable with classical methods of analysis.