Estimating conformational landscapes from Cryo-EM particles by 3D Zernike polynomials

Single-particle cryogenic electron microscopy structure determination for membrane proteins

Membrane proteins are crucial to many cellular functions but are notoriously difficult for structural studies due to their instability outside their natural environment and their amphipathic nature with dual hydrophobic and hydrophilic regions. Single-particle cryogenic electron microscopy (cryo-EM) has emerged as a transformative approach, providing near-atomic-resolution structures without the need for crystallization. This review discusses advancements in cryo-EM, emphasizing membrane sample preparation and data processing techniques. It explores innovations in capturing membrane protein structures within native environments, analyzing their dynamics, binding partner interactions, lipid associations, and responses to electrochemical gradients. These developments continue to enhance our understanding of these vital biomolecules, advancing the contributions of structural biology for basic and translational biomedicine.

InstaMap: instant-NGP for cryo-EM density maps

Despite the parallels between problems in computer vision and cryo-electron microscopy (cryo-EM), many state-of-the-art approaches from computer vision have yet to be adapted for cryo-EM. Within the computer-vision research community, implicits such as neural radiance fields (NeRFs) have enabled the detailed reconstruction of 3D objects from few images at different camera-viewing angles. While other neural implicits, specifically density fields, have been used to map conformational heterogeneity from noisy cryo-EM projection images, most approaches represent volume with an implicit function in Fourier space, which has disadvantages compared with solving the problem in real space, complicating, for instance, masking, constraining physics or geometry, and assessing local resolution. In this work, we build on a recent development in neural implicits, a multi-resolution hash-encoding framework called instant-NGP, that we use to represent the scalar volume directly in real space and apply it to the cryo-EM density-map reconstruction problem (InstaMap). We demonstrate that for both synthetic and real data, InstaMap for homogeneous reconstruction achieves higher resolution at shorter training stages than five other real-spaced representations. We propose a solution to noise overfitting, demonstrate that InstaMap is both lightweight and fast to train, implement masking from a user-provided input mask and extend it to molecular-shape heterogeneity via bending space using a per-image vector field.

Cryo-EM heterogeneity analysis using regularized covariance estimation and kernel regression

Proteins and the complexes they form are central to nearly all cellular processes. Their flexibility, expressed through a continuum of states, provides a window into their biological functions. Cryogenic electron microscopy (cryo-EM) is an ideal tool to study these dynamic states as it captures specimens in noncrystalline conditions and enables high-resolution reconstructions. However, analyzing the heterogeneous distributions of conformations from cryo-EM data is challenging. We present RECOVAR, a method for analyzing these distributions based on principal component analysis (PCA) computed using a REgularized COVARiance estimator. RECOVAR is fast, robust, interpretable, expressive, and competitive with state-of-the-art neural network methods on heterogeneous cryo-EM datasets. The regularized covariance method efficiently computes a large number of high-resolution principal components that can encode rich heterogeneous distributions of conformations and does so robustly thanks to an automatic regularization scheme. The reconstruction method based on adaptive kernel regression resolves conformational states to a higher resolution than all other tested methods on extensive independent benchmarks while remaining highly interpretable. Additionally, we exploit favorable properties of the PCA embedding to estimate the conformational density accurately. This density allows for better interpretability of the latent space by identifying stable states and low free-energy motions. Finally, we present a scheme to navigate the high-dimensional latent space by automatically identifying these low free-energy trajectories. We make the code freely available at https://github.com/ma-gilles/recovar.

A new statistical metric for robust target detection in cryo-EM using 2D template matching

2D template matching (2DTM) can be used to detect molecules and their assemblies in cellular cryo-EM images with high positional and orientational accuracy. While 2DTM successfully detects spherical targets such as large ribosomal subunits, challenges remain in detecting smaller and more aspherical targets in various environments. In this work, a novel 2DTM metric, referred to as the 2DTM p-value, is developed to extend the 2DTM framework to more complex applications. The 2DTM p-value combines information from two previously used 2DTM metrics, namely the 2DTM signal-to-noise ratio (SNR) and z-score, which are derived from the cross-correlation coefficient between the target and the template. The 2DTM p-value demonstrates robust detection accuracies under various imaging and sample conditions and outperforms the 2DTM SNR and z-score alone. Specifically, the 2DTM p-value improves the detection of aspherical targets such as a modified artificial tubulin patch particle (500 kDa) and a much smaller clathrin monomer (193 kDa) in simulated data. It also accurately recovers mature 60S ribosomes in yeast lamellae samples, even under conditions of increased Gaussian noise. The new metric will enable the detection of a wider variety of targets in both purified and cellular samples through 2DTM.

Transient interactions between the fuzzy coat and the cross-β core of brain-derived Aβ42 filaments

Several human disorders, including Alzheimer's disease (AD), are characterized by the aberrant formation of amyloid fibrils. In many cases, the amyloid core is flanked by disordered regions, known as fuzzy coat. The structural properties of fuzzy coats, and their interactions with their environments, however, have not been fully described to date. Here, we generate conformational ensembles of two brain-derived amyloid filaments of Aβ42, corresponding respectively to the familial and sporadic forms of AD. Our approach, called metadynamic electron microscopy metainference (MEMMI), provides a characterization of the transient interactions between the fuzzy coat and the cross-β core of the filaments. These calculations indicate that the familial AD filaments are less soluble than the sporadic AD filaments, and that the fuzzy coat contributes to solubilizing both types of filament. These results illustrate how the metainference approach can help analyze cryo-EM maps for the characterization of the properties of amyloid fibrils.

Partial closure of the γ-tubulin ring complex by CDK5RAP2 activates microtubule nucleation

Microtubule nucleation is templated by the γ-tubulin ring complex (γ-TuRC), but its structure deviates from the geometry of α-/β-tubulin in the microtubule, explaining the complex's poor nucleating activity. Several proteins may activate the γ-TuRC, but the mechanisms underlying activation are not known. Here, we determined the structure of the porcine γ-TuRC purified using CDK5RAP2's centrosomin motif 1 (CM1). We identified an unexpected conformation of the γ-TuRC bound to multiple protein modules containing MZT2, GCP2, and CDK5RAP2, resulting in a long-range constriction of the γ-tubulin ring that brings it in closer agreement with the 13-protofilament microtubule. Additional CDK5RAP2 promoted γ-TuRC decoration and stimulated the microtubule-nucleating activities of the porcine γ-TuRC and a reconstituted, CM1-free human complex in single-molecule assays. Our results provide a structural mechanism for the control of microtubule nucleation by CM1 proteins and identify conformational transitions in the γ-TuRC that prime it for microtubule nucleation.

CryoSTAR: leveraging structural priors and constraints for cryo-EM heterogeneous reconstruction

Resolving conformational heterogeneity in cryogenic electron microscopy datasets remains an important challenge in structural biology. Previous methods have often been restricted to working exclusively on volumetric densities, neglecting the potential of incorporating any preexisting structural knowledge as prior or constraints. Here we present cryoSTAR, which harnesses atomic model information as structural regularization to elucidate such heterogeneity. Our method uniquely outputs both coarse-grained models and density maps, showcasing the molecular conformational changes at different levels. Validated against four diverse experimental datasets, spanning large complexes, a membrane protein and a small single-chain protein, our results consistently demonstrate an efficient and effective solution to conformational heterogeneity with minimal human bias. By integrating atomic model insights with cryogenic electron microscopy data, cryoSTAR represents a meaningful step forward, paving the way for a deeper understanding of dynamic biological processes.

Crossed laser phase plates for transmission electron microscopy

For decades since the development of phase-contrast optical microscopy, an analogous approach has been sought for maximizing the image contrast of weakly-scattering objects in transmission electron microscopy (TEM). The recent development of the laser phase plate (LPP) has demonstrated that an amplified, focused laser standing wave provides stable, tunable phase shift to the high-energy electron beam, achieving phase-contrast TEM. Building on proof-of-concept experimental demonstrations, this paper explores design improvements tailored to biological imaging. In particular, we introduce the approach of crossed laser phase plates (XLPP): two laser standing waves intersecting in the diffraction plane of the TEM, rather than a single beam as in the current LPP. We provide a theoretical model for the XLPP inside the microscope and use simulations to quantify its effect on image formation. We find that the XLPP increases information transfer at low spatial frequencies while also suppressing the ghost images formed by Kapitza-Dirac diffraction of the electron beam by the laser beam. We also demonstrate a simple acquisition scheme, enabled by the XLPP, which dramatically suppresses unwanted diffraction effects. The results of this study chart the course for future developments of LPP hardware.

Challenges, advances, and opportunities in RNA structural biology by Cryo-EM

RNAs are remarkably versatile molecules that can fold into intricate three-dimensional (3D) structures to perform diverse cellular and viral functions. Despite their biological importance, relatively few RNA 3D structures have been solved, and our understanding of RNA structure-function relationships remains in its infancy. This limitation partly arises from challenges posed by RNA's complex conformational landscape, characterized by structural flexibility, formation of multiple states, and a propensity to misfold. Recently, cryo-electron microscopy (cryo-EM) has emerged as a powerful tool for the visualization of conformationally dynamic RNA-only 3D structures. However, RNA's characteristics continue to pose challenges. We discuss experimental methods developed to overcome these hurdles, including the engineering of modular modifications that facilitate the visualization of small RNAs, improve particle alignment, and validate structural models.

DynaMight: estimating molecular motions with improved reconstruction from cryo-EM images

How to deal with continuously flexing molecules is one of the biggest outstanding challenges in single-particle analysis of proteins from cryogenic-electron microscopy (cryo-EM) images. Here, we present DynaMight, a software tool that estimates a continuous space of conformations in a cryo-EM dataset by learning three-dimensional deformations of a Gaussian pseudo-atomic model of a consensus structure for every particle image. Inversion of the learned deformations is then used to obtain an improved reconstruction of the consensus structure. We illustrate the performance of DynaMight for several experimental cryo-EM datasets. We also show how error estimates on the deformations may be obtained by independently training two variational autoencoders on half sets of the cryo-EM data, and how regularization of the three-dimensional deformations through the use of atomic models may lead to important artifacts due to model bias. DynaMight is distributed as free, open-source software, as part of RELION-5.

Cryo-EM Heterogeneity Analysis using Regularized Covariance Estimation and Kernel Regression

Proteins and the complexes they form are central to nearly all cellular processes. Their flexibility, expressed through a continuum of states, provides a window into their biological functions. Cryogenic electron microscopy (cryo-EM) is an ideal tool to study these dynamic states as it captures specimens in non-crystalline conditions and enables high-resolution reconstructions. However, analyzing the heterogeneous distributions of conformations from cryo-EM data is challenging. We present RECOVAR, a method for analyzing these distributions based on principal component analysis (PCA) computed using a REgularized COVARiance estimator. RECOVAR is fast, robust, interpretable, expressive, and competitive with the state-of-art neural network methods on heterogeneous cryo-EM datasets. The regularized covariance method efficiently computes a large number of high-resolution principal components that can encode rich heterogeneous distributions of conformations and does so robustly thanks to an automatic regularization scheme. The novel reconstruction method based on adaptive kernel regression resolves conformational states to a higher resolution than all other tested methods on extensive independent benchmarks while remaining highly interpretable. Additionally, we exploit favorable properties of the PCA embedding to estimate the conformational density accurately. This density allows for better interpretability of the latent space by identifying stable states and low free-energy motions. Finally, we present a scheme to navigate the high-dimensional latent space by automatically identifying these low free-energy trajectories. We make the code freely available at https://github.com/ma-gilles/recovar.

CryoTRANS: predicting high-resolution maps of rare conformations from self-supervised trajectories in cryo-EM

Cryogenic electron microscopy (cryo-EM) has revolutionized structural biology, enabling efficient determination of structures at near-atomic resolutions. However, a common challenge arises from the severe imbalance among various conformations of vitrified particles, leading to low-resolution reconstructions in rare conformations due to a lack of particle images in these quasi-stable states. We introduce CryoTRANS, a method that predicts high-resolution maps of rare conformations by constructing a self-supervised pseudo-trajectory between density maps of varying resolutions. This trajectory is represented by an ordinary differential equation parameterized by a deep neural network, ensuring retention of detailed structures from high-resolution density maps. By leveraging a single high-resolution density map, CryoTRANS significantly improves the reconstruction of rare conformations and has been validated on four real-world datasets: alpha-2-macroglobulin, actin-binding protein complexes, SARS-CoV-2 spike glycoprotein, and the 70S ribosome. CryoTRANS can also predict high-resolution structures in cryogenic electron tomography maps using a high-resolution cryo-EM map.

Accurate Prediction of Protein Structural Flexibility by Deep Learning Integrating Intricate Atomic Structures and Cryo-EM Density Information

The dynamics of proteins are crucial for understanding their mechanisms. However, computationally predicting protein dynamic information has proven challenging. Here, we propose a neural network model, RMSF-net, which outperforms previous methods and produces the best results in a large-scale protein dynamics dataset; this model can accurately infer the dynamic information of a protein in only a few seconds. By learning effectively from experimental protein structure data and cryo-electron microscopy (cryo-EM) data integration, our approach is able to accurately identify the interactive bidirectional constraints and supervision between cryo-EM maps and PDB models in maximizing the dynamic prediction efficacy. Rigorous 5-fold cross-validation on the dataset demonstrates that RMSF-net achieves test correlation coefficients of 0.746 ± 0.127 at the voxel level and 0.765 ± 0.109 at the residue level, showcasing its ability to deliver dynamic predictions closely approximating molecular dynamics simulations. Additionally, it offers real-time dynamic inference with minimal storage overhead on the order of megabytes. RMSF-net is a freely accessible tool and is anticipated to play an essential role in the study of protein dynamics.

Use of Zernike moments to characterize dose conformity for radiotherapy treatment plans

Dose conformity is an essential parameter used in radiotherapy and radiosurgery that measures the correspondence of the dose distribution derived from a Treatment Planning System (TPS) with the actual volume to be treated, the Planning Treatment Volume (PTV). The present work uses a method based on the expansion of dose distributions and PTVs by three-dimensional Zernike polynomials and further comparison of their moments to define a general criterion of dose conformity. To carry on this study, data coming from 20 patients comprising 80 datasets exported from the TPS, which included imaging data (PTVs) and dose distributions corresponding to different treatment modalities: three-dimensional conformal radiotherapy, intensity-modulated radiotherapy (IMRT) and volumetric modulated arc therapy (VMAT), were used. The expansions in Zernike polynomials were obtained up to order 6 and reconstructed dose distributions and PTVs were obtained and compared, and several definitions for a general dose conformity index were proposed. Results indicate agreement between the proposed dose conformity index and the Conformation Number CN. The proposed method allows for a systematic approach to the analysis of dose distributions with further extensions in AI applications.

Accurate model and ensemble refinement using cryo-electron microscopy maps and Bayesian inference

Converting cryo-electron microscopy (cryo-EM) data into high-quality structural models is a challenging problem of outstanding importance. Current refinement methods often generate unbalanced models in which physico-chemical quality is sacrificed for excellent fit to the data. Furthermore, these techniques struggle to represent the conformational heterogeneity averaged out in low-resolution regions of density maps. Here we introduce EMMIVox, a Bayesian inference approach to determine single-structure models as well as structural ensembles from cryo-EM maps. EMMIVox automatically balances experimental information with accurate physico-chemical models of the system and the surrounding environment, including waters, lipids, and ions. Explicit treatment of data correlation and noise as well as inference of accurate B-factors enable determination of structural models and ensembles with both excellent fit to the data and high stereochemical quality, thus outperforming state-of-the-art refinement techniques. EMMIVox represents a flexible approach to determine high-quality structural models that will contribute to advancing our understanding of the molecular mechanisms underlying biological functions.

Data-driven regularization lowers the size barrier of cryo-EM structure determination

Macromolecular structure determination by electron cryo-microscopy (cryo-EM) is limited by the alignment of noisy images of individual particles. Because smaller particles have weaker signals, alignment errors impose size limitations on its applicability. Here, we explore how image alignment is improved by the application of deep learning to exploit prior knowledge about biological macromolecular structures that would otherwise be difficult to express mathematically. We train a denoising convolutional neural network on pairs of half-set reconstructions from the electron microscopy data bank (EMDB) and use this denoiser as an alternative to a commonly used smoothness prior. We demonstrate that this approach, which we call Blush regularization, yields better reconstructions than do existing algorithms, in particular for data with low signal-to-noise ratios. The reconstruction of a protein-nucleic acid complex with a molecular weight of 40 kDa, which was previously intractable, illustrates that denoising neural networks will expand the applicability of cryo-EM structure determination for a wide range of biological macromolecules.

Confronting heterogeneity in cryogenic electron microscopy data: Innovative strategies and future perspectives with data-driven methods

The surge in the influx of data from cryogenic electron microscopy (cryo-EM) experiments has intensified the demand for robust algorithms capable of autonomously managing structurally heterogeneous datasets. This presents a wealth of exciting opportunities from a data science viewpoint, inspiring the development of numerous innovative, application-specific methods, many of which leverage contemporary data-driven techniques. However, addressing the challenges posed by heterogeneous datasets remains a paramount yet unresolved issue in the field. Here, we explore the subtleties of this challenge and the array of strategies devised to confront it. We pinpoint the shortcomings of existing methodologies and deliberate on prospective avenues for improvement. Specifically, our discussion focuses on strategies to mitigate model overfitting and manage data noise, as well as the effects of constraints, priors, and invariances on the optimization process.

Accurate global and local 3D alignment of cryo-EM density maps using local spatial structural features

Advances in cryo-electron microscopy (cryo-EM) imaging technologies have led to a rapidly increasing number of cryo-EM density maps. Alignment and comparison of density maps play a crucial role in interpreting structural information, such as conformational heterogeneity analysis using global alignment and atomic model assembly through local alignment. Here, we present a fast and accurate global and local cryo-EM density map alignment method called CryoAlign, that leverages local density feature descriptors to capture spatial structure similarities. CryoAlign is a feature-based cryo-EM map alignment tool, in which the employment of feature-based architecture enables the rapid establishment of point pair correspondences and robust estimation of alignment parameters. Extensive experimental evaluations demonstrate the superiority of CryoAlign over the existing methods in terms of both alignment accuracy and speed.

Improving resolution and resolvability of single-particle cryoEM structures using Gaussian mixture models

Cryogenic electron microscopy is widely used in structural biology, but its resolution is often limited by the dynamics of the macromolecule. Here we developed a refinement protocol based on Gaussian mixture models that integrates particle orientation and conformation estimation and improves the alignment for flexible domains of protein structures. We demonstrated this protocol on multiple datasets, resulting in improved resolution and resolvability, locally and globally, by visual and quantitative measures.

MDSPACE and MDTOMO Software for Extracting Continuous Conformational Landscapes from Datasets of Single Particle Images and Subtomograms Based on Molecular Dynamics Simulations: Latest Developments in ContinuousFlex Software Package

Cryo electron microscopy (cryo-EM) instrumentation allows obtaining 3D reconstruction of the structure of biomolecular complexes in vitro (purified complexes studied by single particle analysis) and in situ (complexes studied in cells by cryo electron tomography). Standard cryo-EM approaches allow high-resolution reconstruction of only a few conformational states of a molecular complex, as they rely on data classification into a given number of classes to increase the resolution of the reconstruction from the most populated classes while discarding all other classes. Such discrete classification approaches result in a partial picture of the full conformational variability of the complex, due to continuous conformational transitions with many, uncountable intermediate states. In this article, we present the software with a user-friendly graphical interface for running two recently introduced methods, namely, MDSPACE and MDTOMO, to obtain continuous conformational landscapes of biomolecules by analyzing in vitro and in situ cryo-EM data (single particle images and subtomograms) based on molecular dynamics simulations of an available atomic model of one of the conformations. The MDSPACE and MDTOMO software is part of the open-source ContinuousFlex software package (starting from version 3.4.2 of ContinuousFlex), which can be run as a plugin of the Scipion software package (version 3.1 and later), broadly used in the cryo-EM field.

Virus structures revealed by advanced cryoelectron microscopy methods

Before the resolution revolution, cryoelectron microscopy (cryo-EM) single-particle analysis (SPA) already achieved resolutions beyond 4 Å for certain icosahedral viruses, enabling ab initio atomic model building of these viruses. As the only samples that achieved such high resolution at that time, cryo-EM method development was closely intertwined with the improvement of reconstructions of symmetrical viruses. Viral morphology exhibits significant diversity, ranging from small to large, uniform to non-uniform, and from containing single symmetry to multiple symmetries. Furthermore, viruses undergo conformational changes during their life cycle. Several methods, such as asymmetric reconstruction, Ewald sphere correction, cryoelectron tomography (cryo-ET), and sub-tomogram averaging (STA), have been developed and applied to determine virus structures in vivo and in vitro. This review outlines current advanced cryo-EM methods for high-resolution structure determination of viruses and summarizes accomplishments obtained with these approaches. Moreover, persisting challenges in comprehending virus structures are discussed and we propose potential solutions.

Scipion-EM-ProDy: A Graphical Interface for the ProDy Python Package within the Scipion Workflow Engine Enabling Integration of Databases, Simulations and Cryo-Electron Microscopy Image Processing

Macromolecular assemblies, such as protein complexes, undergo continuous structural dynamics, including global reconfigurations critical for their function. Two fast analytical methods are widely used to study these global dynamics, namely elastic network model normal mode analysis and principal component analysis of ensembles of structures. These approaches have found wide use in various computational studies, driving the development of complex pipelines in several software packages. One common theme has been conformational sampling through hybrid simulations incorporating all-atom molecular dynamics and global modes of motion. However, wide functionality is only available for experienced programmers with limited capabilities for other users. We have, therefore, integrated one popular and extensively developed software for such analyses, the ProDy Python application programming interface, into the Scipion workflow engine. This enables a wider range of users to access a complete range of macromolecular dynamics pipelines beyond the core functionalities available in its command-line applications and the normal mode wizard in VMD. The new protocols and pipelines can be further expanded and integrated into larger workflows, together with other software packages for cryo-electron microscopy image analysis and molecular simulations. We present the resulting plugin, Scipion-EM-ProDy, in detail, highlighting the rich functionality made available by its development.

Novel Artificial Intelligence-Based Approaches for Ab Initio Structure Determination and Atomic Model Building for Cryo-Electron Microscopy

Single particle cryo-electron microscopy (cryo-EM) has emerged as the prevailing method for near-atomic structure determination, shedding light on the important molecular mechanisms of biological macromolecules. However, the inherent dynamics and structural variability of biological complexes coupled with the large number of experimental images generated by a cryo-EM experiment make data processing nontrivial. In particular, ab initio reconstruction and atomic model building remain major bottlenecks that demand substantial computational resources and manual intervention. Approaches utilizing recent innovations in artificial intelligence (AI) technology, particularly deep learning, have the potential to overcome the limitations that cannot be adequately addressed by traditional image processing approaches. Here, we review newly proposed AI-based methods for ab initio volume generation, heterogeneous 3D reconstruction, and atomic model building. We highlight the advancements made by the implementation of AI methods, as well as discuss remaining limitations and areas for future development.

Conformational heterogeneity and probability distributions from single-particle cryo-electron microscopy

Single-particle cryo-electron microscopy (cryo-EM) is a technique that takes projection images of biomolecules frozen at cryogenic temperatures. A major advantage of this technique is its ability to image single biomolecules in heterogeneous conformations. While this poses a challenge for data analysis, recent algorithmic advances have enabled the recovery of heterogeneous conformations from the noisy imaging data. Here, we review methods for the reconstruction and heterogeneity analysis of cryo-EM images, ranging from linear-transformation-based methods to nonlinear deep generative models. We overview the dimensionality-reduction techniques used in heterogeneous 3D reconstruction methods and specify what information each method can infer from the data. Then, we review the methods that use cryo-EM images to estimate probability distributions over conformations in reduced subspaces or predefined by atomistic simulations. We conclude with the ongoing challenges for the cryo-EM community.

MDTOMO method for continuous conformational variability analysis in cryo electron subtomograms based on molecular dynamics simulations

Cryo electron tomography (cryo-ET) allows observing macromolecular complexes in their native environment. The common routine of subtomogram averaging (STA) allows obtaining the three-dimensional (3D) structure of abundant macromolecular complexes, and can be coupled with discrete classification to reveal conformational heterogeneity of the sample. However, the number of complexes extracted from cryo-ET data is usually small, which restricts the discrete-classification results to a small number of enough populated states and, thus, results in a largely incomplete conformational landscape. Alternative approaches are currently being investigated to explore the continuity of the conformational landscapes that in situ cryo-ET studies could provide. In this article, we present MDTOMO, a method for analyzing continuous conformational variability in cryo-ET subtomograms based on Molecular Dynamics (MD) simulations. MDTOMO allows obtaining an atomic-scale model of conformational variability and the corresponding free-energy landscape, from a given set of cryo-ET subtomograms. The article presents the performance of MDTOMO on a synthetic ABC exporter dataset and an in situ SARS-CoV-2 spike dataset. MDTOMO allows analyzing dynamic properties of molecular complexes to understand their biological functions, which could also be useful for structure-based drug discovery.

Methods for Cryo-EM Single Particle Reconstruction of Macromolecules Having Continuous Heterogeneity

Macromolecules change their shape (conformation) in the process of carrying out their functions. The imaging by cryo-electron microscopy of rapidly-frozen, individual copies of macromolecules (single particles) is a powerful and general approach to understanding the motions and energy landscapes of macromolecules. Widely-used computational methods already allow the recovery of a few distinct conformations from heterogeneous single-particle samples, but the treatment of complex forms of heterogeneity such as the continuum of possible transitory states and flexible regions remains largely an open problem. In recent years there has been a surge of new approaches for treating the more general problem of continuous heterogeneity. This paper surveys the current state of the art in this area.

ZART: A novel multiresolution reconstruction algorithm with motion-blur correction for single particle analysis

One of the main purposes of CryoEM Single Particle Analysis is to reconstruct the three-dimensional structure of a macromolecule thanks to the acquisition of many particle images representing different poses of the sample. By estimating the orientation of each projected particle, it is possible to recover the underlying 3D volume by multiple 3D reconstruction methods, usually working either in Fourier or in real space. However, the reconstruction from the projected images works under the assumption that all particles in the dataset correspond to the same conformation of the macromolecule. Although this requisite holds for some macromolecules, it is not true for flexible specimens, leading to motion-induced artefacts in the reconstructed CryoEM maps. In this work, we introduce a new Algebraic Reconstruction Technique called ZART, which is able to include continuous flexibility information during the reconstruction process to improve local resolution and reduce motion blurring. The conformational changes are modelled through Zernike3D polynomials. Our implementation allows for a multiresolution description of the macromolecule adapting itself to the local resolution of the reconstructed map. In addition, ZART has also proven to be a useful algorithm in cases where flexibility is not so dominant, as it improves the overall aspect of the reconstructed maps by improving their local and global resolution.