1
|
Hamitouche I, Jonic S. DeepHEMNMA: ResNet-based hybrid analysis of continuous conformational heterogeneity in cryo-EM single particle images. Front Mol Biosci 2022; 9:965645. [PMID: 36158571 PMCID: PMC9493108 DOI: 10.3389/fmolb.2022.965645] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2022] [Accepted: 08/08/2022] [Indexed: 11/13/2022] Open
Abstract
Single-particle cryo-electron microscopy (cryo-EM) is a technique for biomolecular structure reconstruction from vitrified samples containing many copies of a biomolecular complex (known as single particles) at random unknown 3D orientations and positions. Cryo-EM allows reconstructing multiple conformations of the complexes from images of the same sample, which usually requires many rounds of 2D and 3D classifications to disentangle and interpret the combined conformational, orientational, and translational heterogeneity. The elucidation of different conformations is the key to understand molecular mechanisms behind the biological functions of the complexes and the key to novel drug discovery. Continuous conformational heterogeneity, due to gradual conformational transitions giving raise to many intermediate conformational states of the complexes, is both an obstacle for high-resolution 3D reconstruction of the conformational states and an opportunity to obtain information about multiple coexisting conformational states at once. HEMNMA method, specifically developed for analyzing continuous conformational heterogeneity in cryo-EM, determines the conformation, orientation, and position of the complex in each single particle image by image analysis using normal modes (the motion directions simulated for a given atomic structure or EM map), which in turn allows determining the full conformational space of the complex but at the price of high computational cost. In this article, we present a new method, referred to as DeepHEMNMA, which speeds up HEMNMA by combining it with a residual neural network (ResNet) based deep learning approach. The performance of DeepHEMNMA is shown using synthetic and experimental single particle images.
Collapse
Affiliation(s)
| | - Slavica Jonic
- IMPMC - UMR 7590 CNRS, Sorbonne Université, MNHN, Paris, France
| |
Collapse
|
2
|
Krieger JM, Sorzano COS, Carazo JM, Bahar I. Protein dynamics developments for the large scale and cryoEM: case study of ProDy 2.0. Acta Crystallogr D Struct Biol 2022; 78:399-409. [PMID: 35362464 PMCID: PMC8972803 DOI: 10.1107/s2059798322001966] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2021] [Accepted: 02/18/2022] [Indexed: 11/24/2022] Open
Abstract
Cryo-electron microscopy (cryoEM) has become a well established technique with the potential to produce structures of large and dynamic supramolecular complexes that are not amenable to traditional approaches for studying structure and dynamics. The size and low resolution of such molecular systems often make structural modelling and molecular dynamics simulations challenging and computationally expensive. This, together with the growing wealth of structural data arising from cryoEM and other structural biology methods, has driven a trend in the computational biophysics community towards the development of new pipelines for analysing global dynamics using coarse-grained models and methods. At the centre of this trend has been a return to elastic network models, normal mode analysis (NMA) and ensemble analyses such as principal component analysis, and the growth of hybrid simulation methodologies that make use of them. Here, this field is reviewed with a focus on ProDy, the Python application programming interface for protein dynamics, which has been developed over the last decade. Two key developments in this area are highlighted: (i) ensemble NMA towards extracting and comparing the signature dynamics of homologous structures, aided by the recent SignDy pipeline, and (ii) pseudoatom fitting for more efficient global dynamics analyses of large and low-resolution supramolecular assemblies from cryoEM, revisited in the CryoDy pipeline. It is believed that such a renewal and extension of old models and methods in new pipelines will be critical for driving the field forward into the next cryoEM revolution.
Collapse
Affiliation(s)
- James Michael Krieger
- Biocomputing Unit, Centro Nacional de Biotecnología (CSIC), Calle Darwin 3, 28049 Madrid, Spain
| | - Carlos Oscar S. Sorzano
- Biocomputing Unit, Centro Nacional de Biotecnología (CSIC), Calle Darwin 3, 28049 Madrid, Spain
| | - Jose Maria Carazo
- Biocomputing Unit, Centro Nacional de Biotecnología (CSIC), Calle Darwin 3, 28049 Madrid, Spain
| | - Ivet Bahar
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, 800 Murdoch Building, 3420 Forbes Avenue, Pittsburgh, PA 15213, USA
| |
Collapse
|
3
|
Advances in Xmipp for Cryo-Electron Microscopy: From Xmipp to Scipion. Molecules 2021; 26:molecules26206224. [PMID: 34684805 PMCID: PMC8537808 DOI: 10.3390/molecules26206224] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2021] [Revised: 09/28/2021] [Accepted: 09/29/2021] [Indexed: 11/21/2022] Open
Abstract
Xmipp is an open-source software package consisting of multiple programs for processing data originating from electron microscopy and electron tomography, designed and managed by the Biocomputing Unit of the Spanish National Center for Biotechnology, although with contributions from many other developers over the world. During its 25 years of existence, Xmipp underwent multiple changes and updates. While there were many publications related to new programs and functionality added to Xmipp, there is no single publication on the Xmipp as a package since 2013. In this article, we give an overview of the changes and new work since 2013, describe technologies and techniques used during the development, and take a peek at the future of the package.
Collapse
|
4
|
Harastani M, Eltsov M, Leforestier A, Jonic S. HEMNMA-3D: Cryo Electron Tomography Method Based on Normal Mode Analysis to Study Continuous Conformational Variability of Macromolecular Complexes. Front Mol Biosci 2021; 8:663121. [PMID: 34095222 PMCID: PMC8170028 DOI: 10.3389/fmolb.2021.663121] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2021] [Accepted: 04/09/2021] [Indexed: 12/28/2022] Open
Abstract
Cryogenic electron tomography (cryo-ET) allows structural determination of biomolecules in their native environment (in situ). Its potential of providing information on the dynamics of macromolecular complexes in cells is still largely unexploited, due to the challenges of the data analysis. The crowded cell environment and continuous conformational changes of complexes make difficult disentangling the data heterogeneity. We present HEMNMA-3D, which is, to the best of our knowledge, the first method for analyzing cryo electron subtomograms in terms of continuous conformational changes of complexes. HEMNMA-3D uses a combination of elastic and rigid-body 3D-to-3D iterative alignments of a flexible 3D reference (atomic structure or electron microscopy density map) to match the conformation, orientation, and position of the complex in each subtomogram. The elastic matching combines molecular mechanics simulation (Normal Mode Analysis of the 3D reference) and experimental, subtomogram data analysis. The rigid-body alignment includes compensation for the missing wedge, due to the limited tilt angle of cryo-ET. The conformational parameters (amplitudes of normal modes) of the complexes in subtomograms obtained through the alignment are processed to visualize the distribution of conformations in a space of lower dimension (typically, 2D or 3D) referred to as space of conformations. This allows a visually interpretable insight into the dynamics of the complexes, by calculating 3D averages of subtomograms with similar conformations from selected (densest) regions and by recording movies of the 3D reference's displacement along selected trajectories through the densest regions. We describe HEMNMA-3D and show its validation using synthetic datasets. We apply HEMNMA-3D to an experimental dataset describing in situ nucleosome conformational variability. HEMNMA-3D software is available freely (open-source) as part of ContinuousFlex plugin of Scipion V3.0 (http://scipion.i2pc.es).
Collapse
Affiliation(s)
- Mohamad Harastani
- IMPMC-UMR 7590 CNRS, Sorbonne Université, Muséum National d'Histoire Naturelle, Paris, France
| | - Mikhail Eltsov
- Department of Integrated Structural Biology, Institute of Genetics and Molecular and Cellular Biology, Illkirch, France
| | - Amélie Leforestier
- Laboratoire de Physique des Solides, UMR 8502 CNRS, Université Paris-Saclay, Paris, France
| | - Slavica Jonic
- IMPMC-UMR 7590 CNRS, Sorbonne Université, Muséum National d'Histoire Naturelle, Paris, France
| |
Collapse
|
5
|
Sanchez Sorzano CO, Alvarez-Cabrera AL, Kazemi M, Carazo JM, Jonić S. StructMap: Elastic Distance Analysis of Electron Microscopy Maps for Studying Conformational Changes. Biophys J 2017; 110:1753-1765. [PMID: 27119636 DOI: 10.1016/j.bpj.2016.03.019] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2015] [Revised: 03/14/2016] [Accepted: 03/16/2016] [Indexed: 02/06/2023] Open
Abstract
Single-particle electron microscopy (EM) has been shown to be very powerful for studying structures and associated conformational changes of macromolecular complexes. In the context of analyzing conformational changes of complexes, distinct EM density maps obtained by image analysis and three-dimensional (3D) reconstruction are usually analyzed in 3D for interpretation of structural differences. However, graphic visualization of these differences based on a quantitative analysis of elastic transformations (deformations) among density maps has not been done yet due to a lack of appropriate methods. Here, we present an approach that allows such visualization. This approach is based on statistical analysis of distances among elastically aligned pairs of EM maps (one map is deformed to fit the other map), and results in visualizing EM maps as points in a lower-dimensional distance space. The distances among points in the new space can be analyzed in terms of clusters or trajectories of points related to potential conformational changes. The results of the method are shown with synthetic and experimental EM maps at different resolutions.
Collapse
Affiliation(s)
- Carlos Oscar Sanchez Sorzano
- Biocomputing Unit, Centro Nacional de Biotecnología, Consejo Superior de Investigaciones Científicas, Campus de Cantoblanco, Madrid, Spain
| | - Ana Lucia Alvarez-Cabrera
- Biocomputing Unit, Centro Nacional de Biotecnología, Consejo Superior de Investigaciones Científicas, Campus de Cantoblanco, Madrid, Spain
| | - Mohsen Kazemi
- Biocomputing Unit, Centro Nacional de Biotecnología, Consejo Superior de Investigaciones Científicas, Campus de Cantoblanco, Madrid, Spain
| | - Jose María Carazo
- Biocomputing Unit, Centro Nacional de Biotecnología, Consejo Superior de Investigaciones Científicas, Campus de Cantoblanco, Madrid, Spain
| | - Slavica Jonić
- IMPMC, Sorbonne Universités, CNRS UMR 7590, Université Pierre et Marie Curie, Muséum National d'Histoire Naturelle, IRD UMR 206, Paris, France.
| |
Collapse
|
6
|
Versatility of Approximating Single-Particle Electron Microscopy Density Maps Using Pseudoatoms and Approximation-Accuracy Control. BIOMED RESEARCH INTERNATIONAL 2016; 2016:7060348. [PMID: 28097146 PMCID: PMC5209604 DOI: 10.1155/2016/7060348] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/30/2016] [Accepted: 11/03/2016] [Indexed: 01/09/2023]
Abstract
Three-dimensional Gaussian functions have been shown useful in representing electron microscopy (EM) density maps for studying macromolecular structure and dynamics. Methods that require setting a desired number of Gaussian functions or a maximum number of iterations may result in suboptimal representations of the structure. An alternative is to set a desired error of approximation of the given EM map and then optimize the number of Gaussian functions to achieve this approximation error. In this article, we review different applications of such an approach that uses spherical Gaussian functions of fixed standard deviation, referred to as pseudoatoms. Some of these applications use EM-map normal mode analysis (NMA) with elastic network model (ENM) (applications such as predicting conformational changes of macromolecular complexes or exploring actual conformational changes by normal-mode-based analysis of experimental data) while some other do not use NMA (denoising of EM density maps). In applications based on NMA and ENM, the advantage of using pseudoatoms in EM-map coarse-grain models is that the ENM springs are easily assigned among neighboring grains thanks to their spherical shape and uniformed size. EM-map denoising based on the map coarse-graining was so far only shown using pseudoatoms as grains.
Collapse
|
7
|
Jonić S, Vargas J, Melero R, Gómez-Blanco J, Carazo JM, Sorzano COS. Denoising of high-resolution single-particle electron-microscopy density maps by their approximation using three-dimensional Gaussian functions. J Struct Biol 2016; 194:423-33. [PMID: 27085420 DOI: 10.1016/j.jsb.2016.04.007] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2016] [Revised: 04/12/2016] [Accepted: 04/12/2016] [Indexed: 12/22/2022]
Abstract
Cryo-electron microscopy (cryo-EM) of frozen-hydrated preparations of isolated macromolecular complexes is the method of choice to obtain the structure of complexes that cannot be easily studied by other experimental methods due to their flexibility or large size. An increasing number of macromolecular structures are currently being obtained at subnanometer resolution but the interpretation of structural details in such EM-derived maps is often difficult because of noise at these high-frequency signal components that reduces their contrast. In this paper, we show that the method for EM density-map approximation using Gaussian functions can be used for denoising of single-particle EM maps of high (typically subnanometer) resolution. We show its denoising performance using simulated and experimental EM density maps of several complexes.
Collapse
Affiliation(s)
- S Jonić
- IMPMC, Sorbonne Universités - CNRS UMR 7590, UPMC Univ Paris 6, MNHN, IRD UMR 206, 75005 Paris, France.
| | - J Vargas
- Biocomputing Unit, Centro Nacional de Biotecnología - CSIC, Campus de Cantoblanco, Darwin 3, 28049 Madrid, Spain
| | - R Melero
- Biocomputing Unit, Centro Nacional de Biotecnología - CSIC, Campus de Cantoblanco, Darwin 3, 28049 Madrid, Spain
| | - J Gómez-Blanco
- Biocomputing Unit, Centro Nacional de Biotecnología - CSIC, Campus de Cantoblanco, Darwin 3, 28049 Madrid, Spain
| | - J M Carazo
- Biocomputing Unit, Centro Nacional de Biotecnología - CSIC, Campus de Cantoblanco, Darwin 3, 28049 Madrid, Spain
| | - C O S Sorzano
- Biocomputing Unit, Centro Nacional de Biotecnología - CSIC, Campus de Cantoblanco, Darwin 3, 28049 Madrid, Spain
| |
Collapse
|
8
|
Joubert P, Habeck M. Bayesian inference of initial models in cryo-electron microscopy using pseudo-atoms. Biophys J 2016; 108:1165-75. [PMID: 25762328 DOI: 10.1016/j.bpj.2014.12.054] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2014] [Revised: 12/11/2014] [Accepted: 12/23/2014] [Indexed: 11/28/2022] Open
Abstract
Single-particle cryo-electron microscopy is widely used to study the structure of macromolecular assemblies. Tens of thousands of noisy two-dimensional images of the macromolecular assembly viewed from different directions are used to infer its three-dimensional structure. The first step is to estimate a low-resolution initial model and initial image orientations. This is a challenging global optimization problem with many unknowns, including an unknown orientation for each two-dimensional image. Obtaining a good initial model is crucial for the success of the subsequent refinement step. We introduce a probabilistic algorithm for estimating an initial model. The algorithm is fast, has very few algorithmic parameters, and yields information about the precision of estimated model parameters in addition to the parameters themselves. Our algorithm uses a pseudo-atomic model to represent the low-resolution three-dimensional structure, with isotropic Gaussian components as moveable pseudo-atoms. This leads to a significant reduction in the number of parameters needed to represent the three-dimensional structure, and a simplified way of computing two-dimensional projections. It also contributes to the speed of the algorithm. We combine the estimation of the unknown three-dimensional structure and image orientations in a Bayesian framework. This ensures that there are very few parameters to set, and specifies how to combine different types of prior information about the structure with the given data in a systematic way. To estimate the model parameters we use Markov chain Monte Carlo sampling. The advantage is that instead of just obtaining point estimates of model parameters, we obtain an ensemble of models revealing the precision of the estimated parameters. We demonstrate the algorithm on both simulated and real data.
Collapse
Affiliation(s)
- Paul Joubert
- Felix-Bernstein Institute for Mathematical Statistics, Georg-August-Universität Göttingen, Göttingen, Germany.
| | - Michael Habeck
- Felix-Bernstein Institute for Mathematical Statistics, Georg-August-Universität Göttingen, Göttingen, Germany; Max Planck Institute for Biophysical Chemistry, Göttingen, Germany.
| |
Collapse
|
9
|
Kim J, Kim JG, Yun G, Lee PS, Kim DN. Toward Modular Analysis of Supramolecular Protein Assemblies. J Chem Theory Comput 2015; 11:4260-72. [PMID: 26575921 DOI: 10.1021/acs.jctc.5b00329] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Despite recent advances in molecular simulation technologies, analysis of high-molecular-weight structures is still challenging. Here, we propose an automated model reduction procedure aiming to enable modular analysis of these structures. It employs a component mode synthesis for the reduction of finite element protein models. Reduced models may consist of real biological subunits or artificial partitions whose dynamics is described using the degrees of freedom at the substructural interfaces and a small set of dominant vibrational modes only. Notably, the proper number of dominant modes is automatically determined using a novel estimator for eigenvalue errors without calculating the reference eigensolutions of the full model. The performance of the proposed approach is thoroughly investigated by analyzing 50 representative structures including a crystal structure of GroEL and an electron density map of a ribosome.
Collapse
Affiliation(s)
| | - Jin-Gyun Kim
- Mechanical Systems Safety Research Division, Korea Institute of Machinery & Materials , Gajeongbuk-ro 156, Yuseong-gu, Daejeon 34103, Republic of Korea
| | | | - Phill-Seung Lee
- Department of Mechanical Engineering, Korea Advanced Institute of Science and Technology , Daehak-ro 291, Yuseong-gu, Daejeon 34141, Republic of Korea
| | | |
Collapse
|
10
|
Sorzano COS, de la Rosa-Trevín JM, Tama F, Jonić S. Hybrid Electron Microscopy Normal Mode Analysis graphical interface and protocol. J Struct Biol 2014; 188:134-41. [PMID: 25268657 DOI: 10.1016/j.jsb.2014.09.005] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2014] [Revised: 09/18/2014] [Accepted: 09/20/2014] [Indexed: 10/24/2022]
Abstract
This article presents an integral graphical interface to the Hybrid Electron Microscopy Normal Mode Analysis (HEMNMA) approach that was developed for capturing continuous motions of large macromolecular complexes from single-particle EM images. HEMNMA was shown to be a good approach to analyze multiple conformations of a macromolecular complex but it could not be widely used in the EM field due to a lack of an integral interface. In particular, its use required switching among different software sources as well as selecting modes for image analysis was difficult without the graphical interface. The graphical interface was thus developed to simplify the practical use of HEMNMA. It is implemented in the open-source software package Xmipp 3.1 (http://xmipp.cnb.csic.es) and only a small part of it relies on MATLAB that is accessible through the main interface. Such integration provides the user with an easy way to perform the analysis of macromolecular dynamics and forms a direct connection to the single-particle reconstruction process. A step-by-step HEMNMA protocol with the graphical interface is given in full details in Supplementary material. The graphical interface will be useful to experimentalists who are interested in studies of continuous conformational changes of macromolecular complexes beyond the modeling of continuous heterogeneity in single particle reconstruction.
Collapse
Affiliation(s)
| | | | - Florence Tama
- RIKEN, Advanced Institute for Computational Sciences, Kobe, Hyogo 650-0047, Japan
| | - Slavica Jonić
- IMPMC, Sorbonne Universités - CNRS UMR 7590, UPMC Univ Paris 6, MNHN, IRD UMR 206, 75005 Paris, France.
| |
Collapse
|
11
|
Jin Q, Sorzano C, de la Rosa-Trevín J, Bilbao-Castro J, Núñez-Ramírez R, Llorca O, Tama F, Jonić S. Iterative Elastic 3D-to-2D Alignment Method Using Normal Modes for Studying Structural Dynamics of Large Macromolecular Complexes. Structure 2014; 22:496-506. [DOI: 10.1016/j.str.2014.01.004] [Citation(s) in RCA: 50] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2013] [Revised: 12/28/2013] [Accepted: 01/02/2014] [Indexed: 10/25/2022]
|