VESPER: global and local cryo-EM map alignment using local density vectors

Xusi Han,Genki Terashi,Charles Christoffer,Siyang Chen,Daisuke Kihara
DOI: https://doi.org/10.1038/s41467-021-22401-y
IF: 16.6
2021-04-07
Nature Communications
Abstract:Abstract An increasing number of density maps of biological macromolecules have been determined by cryo-electron microscopy (cryo-EM) and stored in the public database, EMDB. To interpret the structural information contained in EM density maps, alignment of maps is an essential step for structure modeling, comparison of maps, and for database search. Here, we developed VESPER, which captures the similarity of underlying molecular structures embedded in density maps by taking local gradient directions into consideration. Compared to existing methods, VESPER achieved substantially more accurate global and local alignment of maps as well as database retrieval.
multidisciplinary sciences
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to improve the global and local alignment accuracy of cryo - EM (cryo - electron microscopy) density maps and enhance the accuracy of database retrieval. Specifically: 1. **Background problems**: - With the development of cryo - electron microscopy technology, more and more biomacromolecular structures are determined by cryo - electron microscopy and stored in public databases (such as EMDB). - In order to interpret the structural information in these density maps, the alignment of density maps is a crucial step. Alignment operations are essential for structural modeling, map comparison, and database searches. 2. **Limitations of existing methods**: - Existing methods (such as gm fit and fitmap) have limitations in global and local alignment, especially performing poorly in partial matching. - These methods also show limited accuracy when dealing with the sub - structure fitting of complex structures. 3. **VESPER's goals**: - VESPER aims to improve the global and local alignment accuracy of cryo - EM density maps by introducing a new alignment method. - VESPER uses local density vectors to capture the molecular structure information embedded in the density maps and evaluates the alignment effect by calculating the dot product of the matching vectors. - Compared with existing methods, VESPER shows higher alignment accuracy and database retrieval ability on the benchmark data set. ### Working principle of VESPER The core idea of VESPER is to represent each density map as a set of unit vectors pointing to neighboring high - density points. The directions of these vectors capture the local structural information embedded in the density map. The specific steps are as follows: 1. **Vector representation**: - Calculate the unit vectors of each density map using the mean - shift algorithm. - The vector representation captures the local molecular structure information around each voxel. 2. **Alignment scoring**: - Evaluate the alignment effect by summing the dot products of the matching vectors in the two aligned density maps. - Use the fast Fourier transform (FFT) for rotation and translation searches to find the optimal alignment position. 3. **Optimization and refinement**: - For the top 10 highest - scoring models, perform a more refined angular search (default 5° interval). - Output the highest - scoring alignment result. ### Experimental results VESPER has demonstrated its superior performance in multiple experiments, especially in global and local alignment and database retrieval. Experimental results show that VESPER exhibits high accuracy and robustness in density map alignment within different resolution ranges. In conclusion, the main goal of this paper is to develop a new method - VESPER - to improve the global and local alignment accuracy of cryo - EM density maps and improve the accuracy of database retrieval.