Dynamo: a flexible, user-friendly development tool for subtomogram averaging of cryo-EM data in high-performance computing environments

Protocols for Subtomogram Averaging of Membrane Proteins in the Dynamo Software Package

Cryo-electron tomography allows low-resolution three-dimensional (3D) viewing of cellular organelles and macromolecular complexes present as multiple copies within a tomogram. These structures are computationally extracted and averaged in order to obtain high-resolution 3D structures, and provide a map of their spatial distribution and interaction with their biological microenvironment. To do so, we apply the user-friendly Dynamo software package on a tomographic data set. Dynamo acts as a modular toolbox adaptable to different biological scenarios, allowing a custom designed pipeline for subtomogram averaging. Here, we use as a textbook example the mitochondrial docking site of the positive-strand RNA Flock house nodavirus (FHV) to describe how Dynamo coordinates several basic steps in the subtomogram averaging workflow. Our framework covers specific strategies to deal with additional issues in subtomogram averaging as tomographic data management, 3D surface visualization, automatic assignment of asymmetry and inherent loss of Fourier information in presence of preferential views.

Structure of the membrane-assembled retromer coat determined by cryo-electron tomography

Eukaryotic cells traffic proteins and lipids between different compartments using protein-coated vesicles and tubules. The retromer complex is required to generate cargo-selective tubulovesicular carriers from endosomal membranes^1-3. Conserved in eukaryotes, retromer controls the cellular localization and homeostasis of hundreds of transmembrane proteins, and its disruption is associated with major neurodegenerative disorders^4-7. How retromer is assembled and how it is recruited to form coated tubules is not known. Here we describe the structure of the retromer complex (Vps26-Vps29-Vps35) assembled on membrane tubules with the bin/amphiphysin/rvs-domain-containing sorting nexin protein Vps5, using cryo-electron tomography and subtomogram averaging. This reveals a membrane-associated Vps5 array, from which arches of retromer extend away from the membrane surface. Vps35 forms the 'legs' of these arches, and Vps29 resides at the apex where it is free to interact with regulatory factors. The bases of the arches connect to each other and to Vps5 through Vps26, and the presence of the same arches on coated tubules within cells confirms their functional importance. Vps5 binds to Vps26 at a position analogous to the previously described cargo- and Snx3-binding site, which suggests the existence of distinct retromer-sorting nexin assemblies. The structure provides insight into the architecture of the coat and its mechanism of assembly, and suggests that retromer promotes tubule formation by directing the distribution of sorting nexin proteins on the membrane surface while providing a scaffold for regulatory-protein interactions.

Fine details in complex environments: the power of cryo-electron tomography

Cryo-electron tomography (CET) is uniquely suited to obtain structural information from a wide range of biological scales, integrating and bridging knowledge from molecules to cells. In particular, CET can be used to visualise molecular structures in their native environment. Depending on the experiment, a varying degree of resolutions can be achieved, with the first near-atomic molecular structures becoming recently available. The power of CET has increased significantly in the last 5 years, in parallel with improvements in cryo-EM hardware and software that have also benefited single-particle reconstruction techniques. In this review, we cover the typical CET pipeline, starting from sample preparation, to data collection and processing, and highlight in particular the recent developments that support structural biology in situ. We provide some examples that highlight the importance of structure determination of molecules embedded within their native environment, and propose future directions to improve CET performance and accessibility.

Insights into the evolution of bacterial flagellar motors from high-throughput in situ electron cryotomography and subtomogram averaging

In situ structural information on molecular machines can be invaluable in understanding their assembly, mechanism and evolution. Here, the use of electron cryotomography (ECT) to obtain significant insights into how an archetypal molecular machine, the bacterial flagellar motor, functions and how it has evolved is described. Over the last decade, studies using a high-throughput, medium-resolution ECT approach combined with genetics, phylogenetic reconstruction and phenotypic analysis have revealed surprising structural diversity in flagellar motors. Variations in the size and the number of torque-generating proteins in the motor visualized for the first time using ECT has shown that these variations have enabled bacteria to adapt their swimming torque to the environment. Much of the structural diversity can be explained in terms of scaffold structures that facilitate the incorporation of additional motor proteins, and more recent studies have begun to infer evolutionary pathways to higher torque-producing motors. This review seeks to highlight how the emerging power of ECT has enabled the inference of ancestral states from various bacterial species towards understanding how, and `why', flagellar motors have evolved from an ancestral motor to a diversity of variants with adapted or modified functions.

Routine single particle CryoEM sample and grid characterization by tomography

Single particle cryo-electron microscopy (cryoEM) is often performed under the assumption that particles are not adsorbed to the air-water interfaces and in thin, vitreous ice. In this study, we performed fiducial-less tomography on over 50 different cryoEM grid/sample preparations to determine the particle distribution within the ice and the overall geometry of the ice in grid holes. Surprisingly, by studying particles in holes in 3D from over 1000 tomograms, we have determined that the vast majority of particles (approximately 90%) are adsorbed to an air-water interface. The implications of this observation are wide-ranging, with potential ramifications regarding protein denaturation, conformational change, and preferred orientation. We also show that fiducial-less cryo-electron tomography on single particle grids may be used to determine ice thickness, optimal single particle collection areas and strategies, particle heterogeneity, and de novo models for template picking and single particle alignment.

In Situ Imaging of Bacterial Secretion Systems by Electron Cryotomography

The unique property of electron cryotomography (ECT) is its capability to resolve the structure of macromolecular machines in their cellular context. The integration of ECT data with high-resolution structures of purified subcomplexes and live-cell fluorescence light microscopy can generate pseudo-atomic models that lead to a mechanistic understanding across size and time scales. Recent advances in electron detection, sample thinning, data acquisition, and data processing have significantly enhanced the applicability and performance of ECT. Here we describe a detailed workflow for an ECT experiment, including cell culture, vitrification, data acquisition, data reconstruction, tomogram analysis, and subtomogram averaging. This protocol provides an entry point to the technique for students and researchers and indicates the many possible variations arising from specific target properties and the available instrumentation.

An Open-Source Storage Solution for Cryo-Electron Microscopy Samples

Cryo-electron microscopy (cryo-EM) enables the study of biological structures in situ in great detail and to solve protein structures at Ångstrom level resolution. Due to recent advances in instrumentation and data processing, the field of cryo-EM is a rapidly growing. Access to facilities and national centers that house the state-of-the-art microscopes is limited due to the ever-rising demand, resulting in long wait times between sample preparation and data acquisition. To improve sample storage, we have developed a cryo-storage system with an efficient, high storage capacity that enables sample storage in a highly organized manner. This system is simple to use, cost-effective and easily adaptable for any type of grid storage box and dewar and any size cryo-EM laboratory.

Shielding and activation of a viral membrane fusion protein

Entry of enveloped viruses relies on insertion of hydrophobic residues of the viral fusion protein into the host cell membrane. However, the intermediate conformations during fusion remain unknown. Here, we address the fusion mechanism of Rift Valley fever virus. We determine the crystal structure of the Gn glycoprotein and fit it with the Gc fusion protein into cryo-electron microscopy reconstructions of the virion. Our analysis reveals how the Gn shields the hydrophobic fusion loops of the Gc, preventing premature fusion. Electron cryotomography of virions interacting with membranes under acidic conditions reveals how the fusogenic Gc is activated upon removal of the Gn shield. Repositioning of the Gn allows extension of Gc and insertion of fusion loops in the outer leaflet of the target membrane. These data show early structural transitions that enveloped viruses undergo during host cell entry and indicate that analogous shielding mechanisms are utilized across diverse virus families.

Viral fusion proteins undergo extensive conformational changes during entry but intermediate conformations often remain unknown. Here, the authors show how Gn of Rift Valley fever virus fusion protein shields hydrophobic fusion loops of Gc and how these loops embed in the target membrane at acidic conditions.

ICON-MIC: Implementing a CPU/MIC Collaboration Parallel Framework for ICON on Tianhe-2 Supercomputer

Electron tomography (ET) is an important technique for studying the three-dimensional structures of the biological ultrastructure. Recently, ET has reached sub-nanometer resolution for investigating the native and conformational dynamics of macromolecular complexes by combining with the sub-tomogram averaging approach. Due to the limited sampling angles, ET reconstruction typically suffers from the "missing wedge" problem. Using a validation procedure, iterative compressed-sensing optimized nonuniform fast Fourier transform (NUFFT) reconstruction (ICON) demonstrates its power in restoring validated missing information for a low-signal-to-noise ratio biological ET dataset. However, the huge computational demand has become a bottleneck for the application of ICON. In this work, we implemented a parallel acceleration technology ICON-many integrated core (MIC) on Xeon Phi cards to address the huge computational demand of ICON. During this step, we parallelize the element-wise matrix operations and use the efficient summation of a matrix to reduce the cost of matrix computation. We also developed parallel versions of NUFFT on MIC to achieve a high acceleration of ICON by using more efficient fast Fourier transform (FFT) calculation. We then proposed a hybrid task allocation strategy (two-level load balancing) to improve the overall performance of ICON-MIC by making full use of the idle resources on Tianhe-2 supercomputer. Experimental results using two different datasets show that ICON-MIC has high accuracy in biological specimens under different noise levels and a significant acceleration, up to 13.3 × , compared with the CPU version. Further, ICON-MIC has good scalability efficiency and overall performance on Tianhe-2 supercomputer.

Structure and assembly of the Ebola virus nucleocapsid

Ebola and Marburg viruses are filoviruses: filamentous, enveloped viruses that cause haemorrhagic fever. Filoviruses are within the order Mononegavirales, which also includes rabies virus, measles virus, and respiratory syncytial virus. Mononegaviruses have non-segmented, single-stranded negative-sense RNA genomes that are encapsidated by nucleoprotein and other viral proteins to form a helical nucleocapsid. The nucleocapsid acts as a scaffold for virus assembly and as a template for genome transcription and replication. Insights into nucleoprotein-nucleoprotein interactions have been derived from structural studies of oligomerized, RNA-encapsidating nucleoprotein, and cryo-electron microscopy of nucleocapsid or nucleocapsid-like structures. There have been no high-resolution reconstructions of complete mononegavirus nucleocapsids. Here we apply cryo-electron tomography and subtomogram averaging to determine the structure of Ebola virus nucleocapsid within intact viruses and recombinant nucleocapsid-like assemblies. These structures reveal the identity and arrangement of the nucleocapsid components, and suggest that the formation of an extended α-helix from the disordered carboxy-terminal region of nucleoprotein-core links nucleoprotein oligomerization, nucleocapsid condensation, RNA encapsidation, and accessory protein recruitment.

RAZA: A Rapid 3D z-crossings algorithm to segment electron tomograms and extract organelles and macromolecules

Resolving the 3D architecture of cells to atomic resolution is one of the most ambitious challenges of cellular and structural biology. Central to this process is the ability to automate tomogram segmentation to identify sub-cellular components, facilitate molecular docking and annotate detected objects with associated metadata. Here we demonstrate that RAZA (Rapid 3D z-crossings algorithm) provides a robust, accurate, intuitive, fast, and generally applicable segmentation algorithm capable of detecting organelles, membranes, macromolecular assemblies and extrinsic membrane protein domains. RAZA defines each continuous contour within a tomogram as a discrete object and extracts a set of 3D structural fingerprints (major, middle and minor axes, surface area and volume), enabling selective, semi-automated segmentation and object extraction. RAZA takes advantage of the fact that the underlying algorithm is a true 3D edge detector, allowing the axes of a detected object to be defined, independent of its random orientation within a cellular tomogram. The selectivity of object segmentation and extraction can be controlled by specifying a user-defined detection tolerance threshold for each fingerprint parameter, within which segmented objects must fall and/or by altering the number of search parameters, to define morphologically similar structures. We demonstrate the capability of RAZA to selectively extract subgroups of organelles (mitochondria) and macromolecular assemblies (ribosomes) from cellular tomograms. Furthermore, the ability of RAZA to define objects and their contours, provides a basis for molecular docking and rapid tomogram annotation.

The structure and flexibility of conical HIV-1 capsids determined within intact virions

HIV-1 contains a cone-shaped capsid encasing the viral genome. This capsid is thought to follow fullerene geometry-a curved hexameric lattice of the capsid protein, CA, closed by incorporating 12 CA pentamers. Current models for core structure are based on crystallography of hexameric and cross-linked pentameric CA, electron microscopy of tubular CA arrays, and simulations. Here, we report subnanometer-resolution cryo-electron tomography structures of hexameric and pentameric CA within intact HIV-1 particles. Whereas the hexamer structure is compatible with crystallography studies, the pentamer forms using different interfaces. Determining multiple structures revealed how CA flexes to form the variably curved core shell. We show that HIV-1 CA assembles both aberrant and perfect fullerene cones, supporting models in which conical cores assemble de novo after maturation.

Accelerating electron tomography reconstruction algorithm ICON with GPU

Electron tomography (ET) plays an important role in studying in situ cell ultrastructure in three-dimensional space. Due to limited tilt angles, ET reconstruction always suffers from the "missing wedge" problem. With a validation procedure, iterative compressed-sensing optimized NUFFT reconstruction (ICON) demonstrates its power in the restoration of validated missing information for low SNR biological ET dataset. However, the huge computational demand has become a major problem for the application of ICON. In this work, we analyzed the framework of ICON and classified the operations of major steps of ICON reconstruction into three types. Accordingly, we designed parallel strategies and implemented them on graphics processing units (GPU) to generate a parallel program ICON-GPU. With high accuracy, ICON-GPU has a great acceleration compared to its CPU version, up to 83.7×, greatly relieving ICON's dependence on computing resource.

Cryo-electron tomography reveals novel features of a viral RNA replication compartment

Positive-strand RNA viruses, the largest genetic class of viruses, include numerous important pathogens such as Zika virus. These viruses replicate their RNA genomes in novel, membrane-bounded mini-organelles, but the organization of viral proteins and RNAs in these compartments has been largely unknown. We used cryo-electron tomography to reveal many previously unrecognized features of Flock house nodavirus (FHV) RNA replication compartments. These spherular invaginations of outer mitochondrial membranes are packed with electron-dense RNA fibrils and their volumes are closely correlated with RNA replication template length. Each spherule’s necked aperture is crowned by a striking cupped ring structure containing multifunctional FHV RNA replication protein A. Subtomogram averaging of these crowns revealed twelve-fold symmetry, concentric flanking protrusions, and a central electron density. Many crowns were associated with long cytoplasmic fibrils, likely to be exported progeny RNA. These results provide new mechanistic insights into positive-strand RNA virus replication compartment structure, assembly, function and control.

DOI:http://dx.doi.org/10.7554/eLife.25940.001

9Å structure of the COPI coat reveals that the Arf1 GTPase occupies two contrasting molecular environments

COPI coated vesicles mediate trafficking within the Golgi apparatus and between the Golgi and the endoplasmic reticulum. Assembly of a COPI coated vesicle is initiated by the small GTPase Arf1 that recruits the coatomer complex to the membrane, triggering polymerization and budding. The vesicle uncoats before fusion with a target membrane. Coat components are structurally conserved between COPI and clathrin/adaptor proteins. Using cryo-electron tomography and subtomogram averaging, we determined the structure of the COPI coat assembled on membranes in vitro at 9 Å resolution. We also obtained a 2.57 Å resolution crystal structure of βδ-COP. By combining these structures we built a molecular model of the coat. We additionally determined the coat structure in the presence of ArfGAP proteins that regulate coat dissociation. We found that Arf1 occupies contrasting molecular environments within the coat, leading us to hypothesize that some Arf1 molecules may regulate vesicle assembly while others regulate coat disassembly.

DOI:http://dx.doi.org/10.7554/eLife.26691.001

The Dynamo package for tomography and subtomogram averaging: components for MATLAB, GPU computing and EC2 Amazon Web Services

Dynamo is a package for the processing of tomographic data. As a tool for subtomogram averaging, it includes different alignment and classification strategies. Furthermore, its data-management module allows experiments to be organized in groups of tomograms, while offering specialized three-dimensional tomographic browsers that facilitate visualization, location of regions of interest, modelling and particle extraction in complex geometries. Here, a technical description of the package is presented, focusing on its diverse strategies for optimizing computing performance. Dynamo is built upon mbtools (middle layer toolbox), a general-purpose MATLAB library for object-oriented scientific programming specifically developed to underpin Dynamo but usable as an independent tool. Its structure intertwines a flexible MATLAB codebase with precompiled C++ functions that carry the burden of numerically intensive operations. The package can be delivered as a precompiled standalone ready for execution without a MATLAB license. Multicore parallelization on a single node is directly inherited from the high-level parallelization engine provided for MATLAB, automatically imparting a balanced workload among the threads in computationally intense tasks such as alignment and classification, but also in logistic-oriented tasks such as tomogram binning and particle extraction. Dynamo supports the use of graphical processing units (GPUs), yielding considerable speedup factors both for native Dynamo procedures (such as the numerically intensive subtomogram alignment) and procedures defined by the user through its MATLAB-based GPU library for three-dimensional operations. Cloud-based virtual computing environments supplied with a pre-installed version of Dynamo can be publicly accessed through the Amazon Elastic Compute Cloud (EC2), enabling users to rent GPU computing time on a pay-as-you-go basis, thus avoiding upfront investments in hardware and longterm software maintenance.

The advent of structural biology in situ by single particle cryo-electron tomography

Single particle tomography (SPT), also known as subtomogram averaging, is a powerful technique uniquely poised to address questions in structural biology that are not amenable to more traditional approaches like X-ray crystallography, nuclear magnetic resonance, and conventional cryoEM single particle analysis. Owing to its potential for in situ structural biology at subnanometer resolution, SPT has been gaining enormous momentum in the last five years and is becoming a prominent, widely used technique. This method can be applied to unambiguously determine the structures of macromolecular complexes that exhibit compositional and conformational heterogeneity, both in vitro and in situ. Here we review the development of SPT, highlighting its applications and identifying areas of ongoing development.

TomoMiner and TomoMinerCloud: A Software Platform for Large-Scale Subtomogram Structural Analysis

Cryo-electron tomography (cryo-ET) captures the 3D electron density distribution of macromolecular complexes in close to native state. With the rapid advance of cryo-ET acquisition technologies, it is possible to generate large numbers (>100,000) of subtomograms, each containing a macromolecular complex. Often, these subtomograms represent a heterogeneous sample due to variations in the structure and composition of a complex in situ form or because particles are a mixture of different complexes. In this case subtomograms must be classified. However, classification of large numbers of subtomograms is a time-intensive task and often a limiting bottleneck. This paper introduces an open source software platform, TomoMiner, for large-scale subtomogram classification, template matching, subtomogram averaging, and alignment. Its scalable and robust parallel processing allows efficient classification of tens to hundreds of thousands of subtomograms. In addition, TomoMiner provides a pre-configured TomoMinerCloud computing service permitting users without sufficient computing resources instant access to TomoMiners high-performance features.

Challenges of Integrating Stochastic Dynamics and Cryo-Electron Tomograms in Whole-Cell Simulations

Cryo-electron tomography (cryo-ET) has rapidly emerged as a powerful tool to investigate the internal, three-dimensional spatial organization of the cell. In parallel, the GPU-based technology to perform spatially resolved stochastic simulations of whole cells has arisen, allowing the simulation of complex biochemical networks over cell cycle time scales using data taken from -omics, single molecule experiments, and in vitro kinetics. By using real cell geometry derived from cryo-ET data, we have the opportunity to imbue these highly detailed structural data-frozen in time-with realistic biochemical dynamics and investigate how cell structure affects the behavior of the embedded chemical reaction network. Here we present two examples to illustrate the challenges and techniques involved in integrating structural data into stochastic simulations. First, a tomographic reconstruction of Saccharomyces cerevisiae is used to construct the geometry of an entire cell through which a simple stochastic model of an inducible genetic switch is studied. Second, a tomogram of the nuclear periphery in a HeLa cell is converted directly to the simulation geometry through which we study the effects of cellular substructure on the stochastic dynamics of gene repression. These simple chemical models allow us to illustrate how to build whole-cell simulations using cryo-ET derived geometry and the challenges involved in such a process.

Focus: The interface between data collection and data processing in cryo-EM

We present a new software package called Focus that interfaces cryo-transmission electron microscopy (cryo-EM) data collection with computer image processing. Focus creates a user-friendly environment to import and manage data recorded by direct electron detectors and perform elemental image processing tasks in a high-throughput manner while new data is being acquired at the microscope. It provides the functionality required to remotely monitor the progress of data collection and data processing, which is essential now that automation in cryo-EM allows a steady flow of images of single particles, two-dimensional crystals, or electron tomography data to be recorded in overnight sessions. The rapid detection of any errors that may occur greatly increases the productivity of recording sessions at the electron microscope.

Cell-free reconstitution reveals centriole cartwheel assembly mechanisms

How cellular organelles assemble is a fundamental question in biology. The centriole organelle organizes around a nine-fold symmetrical cartwheel structure typically ∼100 nm high comprising a stack of rings that each accommodates nine homodimers of SAS-6 proteins. Whether nine-fold symmetrical ring-like assemblies of SAS-6 proteins harbour more peripheral cartwheel elements is unclear. Furthermore, the mechanisms governing ring stacking are not known. Here we develop a cell-free reconstitution system for core cartwheel assembly. Using cryo-electron tomography, we uncover that the Chlamydomonas reinhardtii proteins CrSAS-6 and Bld10p together drive assembly of the core cartwheel. Moreover, we discover that CrSAS-6 possesses autonomous properties that ensure self-organized ring stacking. Mathematical fitting of reconstituted cartwheel height distribution suggests a mechanism whereby preferential addition of pairs of SAS-6 rings governs cartwheel growth. In conclusion, we have developed a cell-free reconstitution system that reveals fundamental assembly principles at the root of centriole biogenesis.

The centriole is an organelle composed of rings of SAS-6 proteins that form a cartwheel structure. Here the authors develop a cell-free system to examine core cartwheel assembly of C. reinhardtii proteins and discover that CrSAS-6 has autonomous properties that facilitates self-organized stacking of pairs of rings.

Assembly of a nucleus-like structure during viral replication in bacteria

We observed the assembly of a nucleus-like structure in bacteria during viral infection. Using fluorescence microscopy and cryo-electron tomography, we showed that Pseudomonas chlororaphis phage 201φ2-1 assembled a compartment that separated viral DNA from the cytoplasm. The phage compartment was centered by a bipolar tubulin-based spindle, and it segregated phage and bacterial proteins according to function. Proteins involved in DNA replication and transcription localized inside the compartment, whereas proteins involved in translation and nucleotide synthesis localized outside. Later during infection, viral capsids assembled on the cytoplasmic membrane and moved to the surface of the compartment for DNA packaging. Ultimately, viral particles were released from the compartment and the cell lysed. These results demonstrate that phages have evolved a specialized structure to compartmentalize viral replication.

Correlated fluorescence microscopy and cryo-electron tomography of virus-infected or transfected mammalian cells

Correlative light and electron microscopy (CLEM) combines spatiotemporal information from fluorescence light microscopy (fLM) with high-resolution structural data from cryo-electron tomography (cryo-ET). These technologies provide opportunities to bridge knowledge gaps between cell and structural biology. Here we describe our protocol for correlated cryo-fLM, cryo-electron microscopy (cryo-EM), and cryo-ET (i.e., cryo-CLEM) of virus-infected or transfected mammalian cells. Mammalian-derived cells are cultured on EM substrates, using optimized conditions that ensure that the cells are spread thinly across the substrate and are not physically disrupted. The cells are then screened by fLM and vitrified before acquisition of cryo-fLM and cryo-ET images, which is followed by data processing. A complete session from grid preparation through data collection and processing takes 5-15 d for an individual experienced in cryo-EM.

The Architecture of the Cytoplasmic Region of Type III Secretion Systems

Type III secretion systems (T3SSs) are essential devices in the virulence of many Gram-negative bacterial pathogens. They mediate injection of protein effectors of virulence from bacteria into eukaryotic host cells to manipulate them during infection. T3SSs involved in virulence (vT3SSs) are evolutionarily related to bacterial flagellar protein export apparatuses (fT3SSs), which are essential for flagellar assembly and cell motility. The structure of the external and transmembrane parts of both fT3SS and vT3SS is increasingly well-defined. However, the arrangement of their cytoplasmic and inner membrane export apparatuses is much less clear. Here we compare the architecture of the cytoplasmic regions of the vT3SSs of Shigella flexneri and the vT3SS and fT3SS of Salmonella enterica serovar Typhimurium at ~5 and ~4 nm resolution using electron cryotomography and subtomogram averaging. We show that the cytoplasmic regions of vT3SSs display conserved six-fold symmetric features including pods, linkers and an ATPase complex, while fT3SSs probably only display six-fold symmetry in their ATPase region. We also identify other morphological differences between vT3SSs and fT3SSs, such as relative disposition of their inner membrane-attached export platform, C-ring/pods and ATPase complex. Finally, using classification, we find that both types of apparatuses can loose elements of their cytoplasmic region, which may therefore be dynamic.

Resolving macromolecular structures from electron cryo-tomography data using subtomogram averaging in RELION

Electron cryo-tomography (cryo-ET) is a technique that is used to produce 3D pictures (tomograms) of complex objects such as asymmetric viruses, cellular organelles or whole cells from a series of tilted electron cryo-microscopy (cryo-EM) images. Averaging of macromolecular complexes found within tomograms is known as subtomogram averaging, and this technique allows structure determination of macromolecular complexes in situ. Subtomogram averaging is also gaining in popularity for the calculation of initial models for single-particle analysis. We describe herein a protocol for subtomogram averaging from cryo-ET data using the RELION software (http://www2.mrc-lmb.cam.ac.uk/relion). RELION was originally developed for cryo-EM single-particle analysis, and the subtomogram averaging approach presented in this protocol has been implemented in the existing workflow for single-particle analysis so that users may conveniently tap into existing capabilities of the RELION software. We describe how to calculate 3D models for the contrast transfer function (CTF) that describe the transfer of information in the imaging process, and we illustrate the results of classification and subtomogram averaging refinement for cryo-ET data of purified hepatitis B capsid particles and Saccharomyces cerevisiae 80S ribosomes. Using the steps described in this protocol, along with the troubleshooting and optimization guidelines, high-resolution maps can be obtained in which secondary structure elements are resolved subtomogram.

An atomic model of HIV-1 capsid-SP1 reveals structures regulating assembly and maturation

Immature HIV-1 assembles at and buds from the plasma membrane before proteolytic cleavage of the viral Gag polyprotein induces structural maturation. Maturation can be blocked by maturation inhibitors (MIs), thereby abolishing infectivity. The CA (capsid) and SP1 (spacer peptide 1) region of Gag is the key regulator of assembly and maturation and is the target of MIs. We applied optimized cryo-electron tomography and subtomogram averaging to resolve this region within assembled immature HIV-1 particles at 3.9 angstrom resolution and built an atomic model. The structure reveals a network of intra- and intermolecular interactions mediating immature HIV-1 assembly. The proteolytic cleavage site between CA and SP1 is inaccessible to protease. We suggest that MIs prevent CA-SP1 cleavage by stabilizing the structure, and MI resistance develops by destabilizing CA-SP1.

Cryo-Electron Tomography and Subtomogram Averaging

Cryo-electron tomography (cryo-ET) allows 3D volumes to be reconstructed from a set of 2D projection images of a tilted biological sample. It allows densities to be resolved in 3D that would otherwise overlap in 2D projection images. Cryo-ET can be applied to resolve structural features in complex native environments, such as within the cell. Analogous to single-particle reconstruction in cryo-electron microscopy, structures present in multiple copies within tomograms can be extracted, aligned, and averaged, thus increasing the signal-to-noise ratio and resolution. This reconstruction approach, termed subtomogram averaging, can be used to determine protein structures in situ. It can also be applied to facilitate more conventional 2D image analysis approaches. In this chapter, we provide an introduction to cryo-ET and subtomogram averaging. We describe the overall workflow, including tomographic data collection, preprocessing, tomogram reconstruction, subtomogram alignment and averaging, classification, and postprocessing. We consider theoretical issues and practical considerations for each step in the workflow, along with descriptions of recent methodological advances and remaining limitations.

A Molecular-Level Account of the Antigenic Hantaviral Surface

Hantaviruses, a geographically diverse group of zoonotic pathogens, initiate cell infection through the concerted action of Gn and Gc viral surface glycoproteins. Here, we describe the high-resolution crystal structure of the antigenic ectodomain of Gn from Puumala hantavirus (PUUV), a causative agent of hemorrhagic fever with renal syndrome. Fitting of PUUV Gn into an electron cryomicroscopy reconstruction of intact Gn-Gc spike complexes from the closely related but non-pathogenic Tula hantavirus localized Gn tetramers to the membrane-distal surface of the virion. The accuracy of the fitting was corroborated by epitope mapping and genetic analysis of available PUUV sequences. Interestingly, Gn exhibits greater non-synonymous sequence diversity than the less accessible Gc, supporting a role of the host humoral immune response in exerting selective pressure on the virus surface. The fold of PUUV Gn is likely to be widely conserved across hantaviruses.

•

We describe the high-resolution crystal structure of a hantaviral Gn ectodomain

•

Electron cryotomography analysis reveals the ultrastructure of Gn-Gc assembly

•

X-ray fitting and mapping analysis reveals the antigenic hantavirus surface

•

The Gn fold is likely to be widely conserved across this group of viruses

Nucleic Acid Binding by Mason-Pfizer Monkey Virus CA Promotes Virus Assembly and Genome Packaging

The Gag polyprotein of retroviruses drives immature virus assembly by forming hexameric protein lattices. The assembly is primarily mediated by protein-protein interactions between capsid (CA) domains and by interactions between nucleocapsid (NC) domains and RNA. Specific interactions between NC and the viral RNA are required for genome packaging. Previously reported cryoelectron microscopy analysis of immature Mason-Pfizer monkey virus (M-PMV) particles suggested that a basic region (residues RKK) in CA may serve as an additional binding site for nucleic acids. Here, we have introduced mutations into the RKK region in both bacterial and proviral M-PMV vectors and have assessed their impact on M-PMV assembly, structure, RNA binding, budding/release, nuclear trafficking, and infectivity using in vitro and in vivo systems. Our data indicate that the RKK region binds and structures nucleic acid that serves to promote virus particle assembly in the cytoplasm. Moreover, the RKK region appears to be important for recruitment of viral genomic RNA into Gag particles, and this function could be linked to changes in nuclear trafficking. Together these observations suggest that in M-PMV, direct interactions between CA and nucleic acid play important functions in the late stages of the viral life cycle.

IMPORTANCE

Assembly of retrovirus particles is driven by the Gag polyprotein, which can self-assemble to form virus particles and interact with RNA to recruit the viral genome into the particles. Generally, the capsid domains of Gag contribute to essential protein-protein interactions during assembly, while the nucleocapsid domain interacts with RNA. The interactions between the nucleocapsid domain and RNA are important both for identifying the genome and for self-assembly of Gag molecules. Here, we show that a region of basic residues in the capsid protein of the betaretrovirus Mason-Pfizer monkey virus (M-PMV) contributes to interaction of Gag with nucleic acid. This interaction appears to provide a critical scaffolding function that promotes assembly of virus particles in the cytoplasm. It is also crucial for packaging the viral genome and thus for infectivity. These data indicate that, surprisingly, interactions between the capsid domain and RNA play an important role in the assembly of M-PMV.

Alignment algorithms and per-particle CTF correction for single particle cryo-electron tomography

Single particle cryo-electron tomography (cryoSPT) extracts features from cryo-electron tomograms, followed by 3D classification, alignment and averaging to generate improved 3D density maps of such features. Robust methods to correct for the contrast transfer function (CTF) of the electron microscope are necessary for cryoSPT to reach its resolution potential. Many factors can make CTF correction for cryoSPT challenging, such as lack of eucentricity of the specimen stage, inherent low dose per image, specimen charging, beam-induced specimen motions, and defocus gradients resulting both from specimen tilting and from unpredictable ice thickness variations. Current CTF correction methods for cryoET make at least one of the following assumptions: that the defocus at the center of the image is the same across the images of a tiltseries, that the particles all lie at the same Z-height in the embedding ice, and/or that the specimen, the cryo-electron microscopy (cryoEM) grid and/or the carbon support are flat. These experimental conditions are not always met. We have developed a CTF correction algorithm for cryoSPT without making any of the aforementioned assumptions. We also introduce speed and accuracy improvements and a higher degree of automation to the subtomogram averaging algorithms available in EMAN2. Using motion-corrected images of isolated virus particles as a benchmark specimen, recorded with a DE20 direct detection camera, we show that our CTF correction and subtomogram alignment routines can yield subtomogram averages close to 4/5 Nyquist frequency of the detector under our experimental conditions.

Recent progress in structure and dynamics of dual-membrane-spanning bacterial nanomachines

Advances in hard-ware and soft-ware for electron cryo-microscopy and tomography have provided unprecedented structural insights into large protein complexes in bacterial membranes. Tomographic volumes of native complexes in situ, combined with other structural and functional data, reveal functionally important conformational changes. Here, we review recent progress in elucidating the structure and mechanism of dual-membrane-spanning nanomachines involved in bacterial motility, adhesion, pathogenesis and biofilm formation, including the type IV pilus assembly machinery and the type III and VI secretions systems. We highlight how these new structural data shed light on the assembly and action of such machines and discuss future directions for more detailed mechanistic understanding of these massive, fascinating complexes.

Acidic pH-Induced Conformations and LAMP1 Binding of the Lassa Virus Glycoprotein Spike

Lassa virus is an enveloped, bi-segmented RNA virus and the most prevalent and fatal of all Old World arenaviruses. Virus entry into the host cell is mediated by a tripartite surface spike complex, which is composed of two viral glycoprotein subunits, GP1 and GP2, and the stable signal peptide. Of these, GP1 binds to cellular receptors and GP2 catalyzes fusion between the viral envelope and the host cell membrane during endocytosis. The molecular structure of the spike and conformational rearrangements induced by low pH, prior to fusion, remain poorly understood. Here, we analyzed the three-dimensional ultrastructure of Lassa virus using electron cryotomography. Sub-tomogram averaging yielded a structure of the glycoprotein spike at 14-Å resolution. The spikes are trimeric, cover the virion envelope, and connect to the underlying matrix. Structural changes to the spike, following acidification, support a viral entry mechanism dependent on binding to the lysosome-resident receptor LAMP1 and further dissociation of the membrane-distal GP1 subunits.

Using Tomoauto: A Protocol for High-throughput Automated Cryo-electron Tomography

Cryo-electron tomography (Cryo-ET) is a powerful three-dimensional (3-D) imaging technique for visualizing macromolecular complexes in their native context at a molecular level. The technique involves initially preserving the sample in its native state by rapidly freezing the specimen in vitreous ice, then collecting a series of micrographs from different angles at high magnification, and finally computationally reconstructing a 3-D density map. The frozen-hydrated specimen is extremely sensitive to the electron beam and so micrographs are collected at very low electron doses to limit the radiation damage. As a result, the raw cryo-tomogram has a very low signal to noise ratio characterized by an intrinsically noisy image. To better visualize subjects of interest, conventional imaging analysis and sub-tomogram averaging in which sub-tomograms of the subject are extracted from the initial tomogram and aligned and averaged are utilized to improve both contrast and resolution. Large datasets of tilt-series are essential to understanding and resolving the complexes at different states, conditions, or mutations as well as obtaining a large enough collection of sub-tomograms for averaging and classification. Collecting and processing this data can be a major obstacle preventing further analysis. Here we describe a high-throughput cryo-ET protocol based on a computer-controlled 300kV cryo-electron microscope, a direct detection device (DDD) camera and a highly effective, semi-automated image-processing pipeline software wrapper library tomoauto developed in-house. This protocol has been effectively utilized to visualize the intact type III secretion system (T3SS) in Shigella flexneri minicells. It can be applicable to any project suitable for cryo-ET.

The Structure of the Mouse Serotonin 5-HT3 Receptor in Lipid Vesicles

The function of membrane proteins is best understood if their structure in the lipid membrane is known. Here, we determined the structure of the mouse serotonin 5-HT3 receptor inserted in lipid bilayers to a resolution of 12 Å without stabilizing antibodies by cryo electron tomography and subtomogram averaging. The reconstruction reveals protein secondary structure elements in the transmembrane region, the extracellular pore, and the transmembrane channel pathway, showing an overall similarity to the available X-ray model of the truncated 5-HT3 receptor determined in the presence of a stabilizing nanobody. Structural analysis of the 5-HT3 receptor embedded in a lipid bilayer allowed the position of the membrane to be determined. Interactions between the densely packed receptors in lipids were visualized, revealing that the interactions were maintained by the short horizontal helices. In combination with methodological improvements, our approach enables the structural analysis of membrane proteins in response to voltage and ligand gating.

Structure of a bacterial type III secretion system in contact with a host membrane in situ

Many bacterial pathogens of animals and plants use a conserved type III secretion system (T3SS) to inject virulence effector proteins directly into eukaryotic cells to subvert host functions. Contact with host membranes is critical for T3SS activation, yet little is known about T3SS architecture in this state or the conformational changes that drive effector translocation. Here we use cryo-electron tomography and sub-tomogram averaging to derive the intact structure of the primordial Chlamydia trachomatis T3SS in the presence and absence of host membrane contact. Comparison of the averaged structures demonstrates a marked compaction of the basal body (4 nm) occurs when the needle tip contacts the host cell membrane. This compaction is coupled to a stabilization of the cytosolic sorting platform–ATPase. Our findings reveal the first structure of a bacterial T3SS from a major human pathogen engaged with a eukaryotic host, and reveal striking ‘pump-action' conformational changes that underpin effector injection.

Bacterial type III secretion systems (T3SSs) inject virulence effector proteins into eukaryotic cells and are activated by host membrane contact. Here the authors report the in situ structure of the Chlamydia trachomatis T3SS in the presence or absence of host membrane, and observe compaction of the basal body embedded in the bacterial envelope.

Automated batch fiducial-less tilt-series alignment in Appion using Protomo

The field of electron tomography has benefited greatly from manual and semi-automated approaches to marker-based tilt-series alignment that have allowed for the structural determination of multitudes of in situ cellular structures as well as macromolecular structures of individual protein complexes. The emergence of complementary metal-oxide semiconductor detectors capable of detecting individual electrons has enabled the collection of low dose, high contrast images, opening the door for reliable correlation-based tilt-series alignment. Here we present a set of automated, correlation-based tilt-series alignment, contrast transfer function (CTF) correction, and reconstruction workflows for use in conjunction with the Appion/Leginon package that are primarily targeted at automating structure determination with cryogenic electron microscopy.

The Structure of Immature Virus-Like Rous Sarcoma Virus Gag Particles Reveals a Structural Role for the p10 Domain in Assembly

The polyprotein Gag is the primary structural component of retroviruses. Gag consists of independently folded domains connected by flexible linkers. Interactions between the conserved capsid (CA) domains of Gag mediate formation of hexameric protein lattices that drive assembly of immature virus particles. Proteolytic cleavage of Gag by the viral protease (PR) is required for maturation of retroviruses from an immature form into an infectious form. Within the assembled Gag lattices of HIV-1 and Mason-Pfizer monkey virus (M-PMV), the C-terminal domain of CA adopts similar quaternary arrangements, while the N-terminal domain of CA is packed in very different manners. Here, we have used cryo-electron tomography and subtomogram averaging to study in vitro-assembled, immature virus-like Rous sarcoma virus (RSV) Gag particles and have determined the structure of CA and the surrounding regions to a resolution of ∼8 Å. We found that the C-terminal domain of RSV CA is arranged similarly to HIV-1 and M-PMV, whereas the N-terminal domain of CA adopts a novel arrangement in which the upstream p10 domain folds back into the CA lattice. In this position the cleavage site between CA and p10 appears to be inaccessible to PR. Below CA, an extended density is consistent with the presence of a six-helix bundle formed by the spacer-peptide region. We have also assessed the affect of lattice assembly on proteolytic processing by exogenous PR. The cleavage between p10 and CA is indeed inhibited in the assembled lattice, a finding consistent with structural regulation of proteolytic maturation.

IMPORTANCE

Retroviruses first assemble into immature virus particles, requiring interactions between Gag proteins that form a protein layer under the viral membrane. Subsequently, Gag is cleaved by the viral protease enzyme into separate domains, leading to rearrangement of the virus into its infectious form. It is important to understand how Gag is arranged within immature retroviruses, in order to understand how virus assembly occurs, and how maturation takes place. We used the techniques cryo-electron tomography and subtomogram averaging to obtain a detailed structural picture of the CA domains in immature assembled Rous sarcoma virus Gag particles. We found that part of Gag next to CA, called p10, folds back and interacts with CA when Gag assembles. This arrangement is different from that seen in HIV-1 and Mason-Pfizer monkey virus, illustrating further structural diversity of retroviral structures. The structure provides new information on how the virus assembles and undergoes maturation.

Insight into the Assembly of Viruses with Vertical Single β-barrel Major Capsid Proteins

Archaeal viruses constitute the least explored niche within the virosphere. Structure-based approaches have revealed close relationships between viruses infecting organisms from different domains of life. Here, using biochemical and cryo-electron microscopy techniques, we solved the structure of euryarchaeal, halophilic, internal membrane-containing Haloarcula hispanica icosahedral virus 2 (HHIV-2). We show that the density of the two major capsid proteins (MCPs) recapitulates vertical single β-barrel proteins and that disulfide bridges stabilize the capsid. Below, ordered density is visible close to the membrane and at the five-fold vertices underneath the host-interacting vertex complex underpinning membrane-protein interactions. The HHIV-2 structure exemplifies the division of conserved architectural elements of a virion, such as the capsid, from those that evolve rapidly due to selective environmental pressure such as host-recognizing structures. We propose that in viruses with two vertical single β-barrel MCPs the vesicle is indispensable, and membrane-protein interactions serve as protein-railings for guiding the assembly.

Three-dimensional structure of the human breast cancer resistance protein (BCRP/ABCG2) in an inward-facing conformation

ABCG2 is an efflux drug transporter that plays an important role in drug resistance and drug disposition. In this study, the first three-dimensional structure of human full-length ABCG2 analysed by electron crystallography from two-dimensional crystals in the absence of nucleotides and transported substrates is reported at 2 nm resolution. In this state, ABCG2 forms a symmetric homodimer with a noncrystallographic twofold axis perpendicular to the two-dimensional crystal plane, as confirmed by subtomogram averaging. This configuration suggests an inward-facing configuration similar to murine ABCB1, with the nucleotide-binding domains (NBDs) widely separated from each other. In the three-dimensional map, densities representing the long cytoplasmic extensions from the transmembrane domains that connect the NBDs are clearly visible. The structural data have allowed the atomic model of ABCG2 to be refined, in which the two arms of the V-shaped ABCG2 homodimeric complex are in a more closed and narrower conformation. The structural data and the refined model of ABCG2 are compatible with the biochemical analysis of the previously published mutagenesis studies, providing novel insight into the structure and function of the transporter.

Isolation, cryotomography, and three-dimensional reconstruction of centrioles

Centrioles and basal bodies (referred to hereafter as centrioles for simplicity) are microtubule-based cylindrical organelles that are typically ∼450-nm long and ∼250nm in diameter. The centriole is composed of three distinct regions: the distal part characterized by microtubule doublets, the central core that harbors microtubule triplets, which are also present in the proximal part that also contains the cartwheel, a structure crucial for centriole assembly. The cartwheel was initially revealed by conventional electron microscopy of resin-embedded samples and is thought to impart the near universal ninefold symmetry of centrioles. Deciphering the native architecture of the cartwheel has proven challenging owing to its small dimensions and the difficulties in isolating it. Here, we present a method to purify and analyze the structure of the exceptionally long Trichonympha centriole by cryotomography and subtomogram averaging. Using this method, we revealed the native architecture of the proximal cartwheel-containing region at ∼40Å-resolution. This method can be applied as a general strategy for uncovering the structure of centrioles in other species.

Single particle tomography in EMAN2

Single particle tomography (SPT or subtomogram averaging) offers a powerful alternative to traditional 2-D single particle reconstruction for studying conformationally or compositionally heterogeneous macromolecules. It can also provide direct observation (without labeling or staining) of complexes inside cells at nanometer resolution. The development of computational methods and tools for SPT remains an area of active research. Here we present the EMAN2.1 SPT toolbox, which offers a full SPT processing pipeline, from particle picking to post-alignment analysis of subtomogram averages, automating most steps. Different algorithm combinations can be applied at each step, providing versatility and allowing for procedural cross-testing and specimen-specific strategies. Alignment methods include all-vs-all, binary tree, iterative single-model refinement, multiple-model refinement, and self-symmetry alignment. An efficient angular search, Graphic Processing Unit (GPU) acceleration and both threaded and distributed parallelism are provided to speed up processing. Finally, automated simulations, per particle reconstruction of subtiltseries, and per-particle Contrast Transfer Function (CTF) correction have been implemented. Processing examples using both real and simulated data are shown for several structures.

Localize.pytom: a modern webserver for cryo-electron tomography

Localize.pytom, available through http://localize.pytom.org is a webserver for the localize module in the PyTom package. It is a free website and open to all users and there is no login requirement. The server accepts tomograms as they are imaged and reconstructed by Cryo-Electron Tomography (CET) and returns densities and coordinates of candidate-macromolecules in the tomogram. Localization of macromolecules in cryo-electron tomograms is one of the key procedures to unravel structural features of imaged macromolecules. Positions of localized molecules are further used for structural analysis by single particle procedures such as fine alignment, averaging and classification. Accurate localization can be furthermore used to generate molecular atlases of whole cells. Localization uses a cross-correlation-based score and requires a reference volume as input. A reference can either be a previously detected macromolecular structure or extrapolated on the server from a specific PDB chain. Users have the option to use either coarse or fine angular sampling strategies based on uniformly distributed rotations and to accurately compensate for the CET common ‘Missing Wedge’ artefact during sampling. After completion, all candidate macromolecules cut out from the tomogram are available for download. Their coordinates are stored and available in XML format, which can be easily integrated into successive analysis steps in other software. A pre-computed average of the first one hundred macromolecules is also available for immediate download, and the user has the option to further analyse the average, based on the detected score distribution in a novel web-density viewer.

Structure of the type VI secretion system contractile sheath

Bacteria use rapid contraction of a long sheath of the type VI secretion system (T6SS) to deliver effectors into a target cell. Here, we present an atomic-resolution structure of a native contracted Vibrio cholerae sheath determined by cryo-electron microscopy. The sheath subunits, composed of tightly interacting proteins VipA and VipB, assemble into a six-start helix. The helix is stabilized by a core domain assembled from four β strands donated by one VipA and two VipB molecules. The fold of inner and middle layers is conserved between T6SS and phage sheaths. However, the structure of the outer layer is distinct and suggests a mechanism of interaction of the bacterial sheath with an accessory ATPase, ClpV, that facilitates multiple rounds of effector delivery. Our results provide a mechanistic insight into assembly of contractile nanomachines that bacteria and phages use to translocate macromolecules across membranes.