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

Progressive Stochastic Reconstruction Technique (PSRT) for cryo electron tomography

Cryo Electron Tomography (cryoET) plays an essential role in Structural Biology, as it is the only technique that allows to study the structure of large macromolecular complexes in their close to native environment in situ. The reconstruction methods currently in use, such as Weighted Back Projection (WBP) or Simultaneous Iterative Reconstruction Technique (SIRT), deliver noisy and low-contrast reconstructions, which complicates the application of high-resolution protocols, such as Subtomogram Averaging (SA). We propose a Progressive Stochastic Reconstruction Technique (PSRT) - a novel iterative approach to tomographic reconstruction in cryoET based on Monte Carlo random walks guided by Metropolis-Hastings sampling strategy. We design a progressive reconstruction scheme to suit the conditions present in cryoET and apply it successfully to reconstructions of macromolecular complexes from both synthetic and experimental datasets. We show how to integrate PSRT into SA, where it provides an elegant solution to the region-of-interest problem and delivers high-contrast reconstructions that significantly improve template-based localization without any loss of high-resolution structural information. Furthermore, the locality of SA is exploited to design an importance sampling scheme which significantly speeds up the otherwise slow Monte Carlo approach. Finally, we design a new memory efficient solution for the specimen-level interior problem of cryoET, removing all associated artifacts.

Composition, formation, and regulation of the cytosolic c-ring, a dynamic component of the type III secretion injectisome

The injectisome is a membrane complex through which some bacteria can inject effector proteins into host cells. This study reveals that the cytosolic C-ring structure has a dynamic relationship to the rest of the injectisome, with implications for the regulation of secretion.

Many gram-negative pathogens employ a type III secretion injectisome to translocate effector proteins into eukaryotic host cells. While the structure of the distal “needle complex” is well documented, the composition and role of the functionally important cytosolic complex remain less well understood. Using functional fluorescent fusions, we found that the C-ring, an essential and conserved cytosolic component of the system, is composed of ~22 copies of SctQ (YscQ in Yersinia enterocolitica), which require the presence of YscQ_C, the product of an internal translation initiation site in yscQ, for their cooperative assembly. Photoactivated localization microscopy (PALM) reveals that in vivo, YscQ is present in both a free-moving cytosolic and a stable injectisome-bound state. Notably, fluorescence recovery after photobleaching (FRAP) shows that YscQ exchanges between the injectisome and the cytosol, with a t_½ of 68 ± 8 seconds when injectisomes are secreting. In contrast, the secretin SctC (YscC) and the major export apparatus component SctV (YscV) display minimal exchange. Under non-secreting conditions, the exchange rate of YscQ is reduced to t_½ = 134 ± 16 seconds, revealing a correlation between C-ring exchange and injectisome activity, which indicates a possible role for C-ring stability in regulation of type III secretion. The stabilization of the C-ring depends on the presence of the functional ATPase SctN (YscN). These data provide new insights into the formation and composition of the injectisome and present a novel aspect of type III secretion, the exchange of C-ring subunits, which is regulated with respect to secretion.

The type III secretion system, also known as the injectisome, is a key virulence factor in many gram-negative bacteria, and is responsible for the transmission of bacterial proteins directly into host cells. While some elements of the system are well characterized, the cytosolic components involved in substrate recognition and handling are not well understood. One of the major questions is the role of the C-ring, an essential yet enigmatic cytosolic injectisome member. We used fluorescence microscopy to analyze the architecture and behavior of the C-ring in live Y. enterocolitica bacteria, a human pathogen. We found that in vivo, the C-ring assembles cooperatively with the help of additional copies of its own C-terminal fragment and has a highly dynamic structure, with C-ring subunits exchanging between the working injectisomes and a cytosolic pool. The rate of exchange is different between secreting and non-secreting injectisomes and depends on the function of the type III secretion ATPase, indicating that the stability of the complex is altered when functioning. This dynamic behaviour raises the possibility that the C-ring is a regulator of targeted protein delivery by the type III secretion system and makes the C-ring a viable target for the development of novel anti-virulence drugs.

Structure of the immature HIV-1 capsid in intact virus particles at 8.8 Å resolution

Human immunodeficiency virus type 1 (HIV-1) assembly proceeds in two stages. First, the 55 kilodalton viral Gag polyprotein assembles into a hexameric protein lattice at the plasma membrane of the infected cell, inducing budding and release of an immature particle. Second, Gag is cleaved by the viral protease, leading to internal rearrangement of the virus into the mature, infectious form. Immature and mature HIV-1 particles are heterogeneous in size and morphology, preventing high-resolution analysis of their protein arrangement in situ by conventional structural biology methods. Here we apply cryo-electron tomography and sub-tomogram averaging methods to resolve the structure of the capsid lattice within intact immature HIV-1 particles at subnanometre resolution, allowing unambiguous positioning of all α-helices. The resulting model reveals tertiary and quaternary structural interactions that mediate HIV-1 assembly. Strikingly, these interactions differ from those predicted by the current model based on in vitro-assembled arrays of Gag-derived proteins from Mason-Pfizer monkey virus. To validate this difference, we solve the structure of the capsid lattice within intact immature Mason-Pfizer monkey virus particles. Comparison with the immature HIV-1 structure reveals that retroviral capsid proteins, while having conserved tertiary structures, adopt different quaternary arrangements during virus assembly. The approach demonstrated here should be applicable to determine structures of other proteins at subnanometre resolution within heterogeneous environments.

Averaging of viral envelope glycoprotein spikes from electron cryotomography reconstructions using Jsubtomo

Enveloped viruses utilize membrane glycoproteins on their surface to mediate entry into host cells. Three-dimensional structural analysis of these glycoprotein ‘spikes’ is often technically challenging but important for understanding viral pathogenesis and in drug design. Here, a protocol is presented for viral spike structure determination through computational averaging of electron cryo-tomography data. Electron cryo-tomography is a technique in electron microscopy used to derive three-dimensional tomographic volume reconstructions, or tomograms, of pleomorphic biological specimens such as membrane viruses in a near-native, frozen-hydrated state. These tomograms reveal structures of interest in three dimensions, albeit at low resolution. Computational averaging of sub-volumes, or sub-tomograms, is necessary to obtain higher resolution detail of repeating structural motifs, such as viral glycoprotein spikes. A detailed computational approach for aligning and averaging sub-tomograms using the Jsubtomo software package is outlined. This approach enables visualization of the structure of viral glycoprotein spikes to a resolution in the range of 20-40 Å and study of the study of higher order spike-to-spike interactions on the virion membrane. Typical results are presented for Bunyamwera virus, an enveloped virus from the family Bunyaviridae. This family is a structurally diverse group of pathogens posing a threat to human and animal health.

Insights into the mechanisms of membrane curvature and vesicle scission by the small GTPase Sar1 in the early secretory pathway

The small GTPase protein Sar1 is known to be involved in both the initiation of COPII-coated vesicle formation and scission of the nascent vesicle from the endoplasmic reticulum. The molecular details for the mechanism of membrane remodeling by Sar1 remain unresolved. Here, we show that Sar1 transforms synthetic liposomes into structures of different morphologies including tubules and detached vesicles. We demonstrate that Sar1 alone is competent for vesicle scission in a manner that depends on the concentration of Sar1 molecules occupying the membrane. Sar1 molecules align on low-curvature membranes to form an extended lattice. The continuity of this lattice breaks down as the curvature locally increases. The smallest repeating unit constituting the ordered lattice is a Sar1 dimer. The three-dimensional structure of the Sar1 lattice was reconstructed by substituting spherical liposomes with galactoceramide lipid tubules of homogeneous diameter. These data suggest that Sar1 dimerization is responsible for the formation of constrictive membrane curvature. We propose a model whereby Sar1 dimers assemble into ordered arrays to promote membrane constriction and COPII-directed vesicle scission.

Structural analysis of supramolecular assemblies by cryo-electron tomography

Structural analysis of macromolecular assemblies in their physiological environment is a challenging task that is instrumental in answering fundamental questions in cellular and molecular structural biology. The continuous development of computational and analytical tools for cryo-electron tomography (cryo-ET) enables the study of these assemblies at a resolution of a few nanometers. Through the implementation of thinning procedures, cryo-ET can now be applied to the reconstruction of macromolecular structures located inside thick regions of vitrified cells and tissues, thus becoming a central tool for structural determinations in various biological disciplines. Here, we focus on the successful in situ applications of cryo-ET to reveal structures of macromolecular complexes within eukaryotic cells.

gEMpicker: a highly parallel GPU-accelerated particle picking tool for cryo-electron microscopy

Background

Picking images of particles in cryo-electron micrographs is an important step in solving the 3D structures of large macromolecular assemblies. However, in order to achieve sub-nanometre resolution it is often necessary to capture and process many thousands or even several millions of 2D particle images. Thus, a computational bottleneck in reaching high resolution is the accurate and automatic picking of particles from raw cryo-electron micrographs.

Results

We have developed “gEMpicker”, a highly parallel correlation-based particle picking tool. To our knowledge, gEMpicker is the first particle picking program to use multiple graphics processor units (GPUs) to accelerate the calculation. When tested on the publicly available keyhole limpet hemocyanin dataset, we find that gEMpicker gives similar results to the FindEM program. However, compared to calculating correlations on one core of a contemporary central processor unit (CPU), running gEMpicker on a modern GPU gives a speed-up of about 27 ×. To achieve even higher processing speeds, the basic correlation calculations are accelerated considerably by using a hierarchy of parallel programming techniques to distribute the calculation over multiple GPUs and CPU cores attached to multiple nodes of a computer cluster. By using a theoretically optimal reduction algorithm to collect and combine the cluster calculation results, the speed of the overall calculation scales almost linearly with the number of cluster nodes available.

Conclusions

The very high picking throughput that is now possible using GPU-powered workstations or computer clusters will help experimentalists to achieve higher resolution 3D reconstructions more rapidly than before.

Mechanism of membranous tunnelling nanotube formation in viral genome delivery

In internal membrane-containing viruses, a lipid vesicle enclosed by the icosahedral capsid protects the genome. It has been postulated that this internal membrane is the genome delivery device of the virus. Viruses built with this architectural principle infect hosts in all three domains of cellular life. Here, using a combination of electron microscopy techniques, we investigate bacteriophage PRD1, the best understood model for such viruses, to unveil the mechanism behind the genome translocation across the cell envelope. To deliver its double-stranded DNA, the icosahedral protein-rich virus membrane transforms into a tubular structure protruding from one of the 12 vertices of the capsid. We suggest that this viral nanotube exits from the same vertex used for DNA packaging, which is biochemically distinct from the other 11. The tube crosses the capsid through an aperture corresponding to the loss of the peripentonal P3 major capsid protein trimers, penton protein P31 and membrane protein P16. The remodeling of the internal viral membrane is nucleated by changes in osmolarity and loss of capsid-membrane interactions as consequence of the de-capping of the vertices. This engages the polymerization of the tail tube, which is structured by membrane-associated proteins. We have observed that the proteo-lipidic tube in vivo can pierce the gram-negative bacterial cell envelope allowing the viral genome to be shuttled to the host cell. The internal diameter of the tube allows one double-stranded DNA chain to be translocated. We conclude that the assembly principles of the viral tunneling nanotube take advantage of proteo-lipid interactions that confer to the tail tube elastic, mechanical and functional properties employed also in other protein-membrane systems.

The structure of the COPII transport-vesicle coat assembled on membranes

Coat protein complex II (COPII) mediates formation of the membrane vesicles that export newly synthesised proteins from the endoplasmic reticulum. The inner COPII proteins bind to cargo and membrane, linking them to the outer COPII components that form a cage around the vesicle. Regulated flexibility in coat architecture is essential for transport of a variety of differently sized cargoes, but structural data on the assembled coat has not been available. We have used cryo-electron tomography and subtomogram averaging to determine the structure of the complete, membrane-assembled COPII coat. We describe a novel arrangement of the outer coat and find that the inner coat can assemble into regular lattices. The data reveal how coat subunits interact with one another and with the membrane, suggesting how coordinated assembly of inner and outer coats can mediate and regulate packaging of vesicles ranging from small spheres to large tubular carriers. DOI:http://dx.doi.org/10.7554/eLife.00951.001.

In situ structural analysis of the Yersinia enterocolitica injectisome

Injectisomes are multi-protein transmembrane machines allowing pathogenic bacteria to inject effector proteins into eukaryotic host cells, a process called type III secretion. Here we present the first three-dimensional structure of Yersinia enterocolitica and Shigella flexneri injectisomes in situ and the first structural analysis of the Yersinia injectisome. Unexpectedly, basal bodies of injectisomes inside the bacterial cells showed length variations of 20%. The in situ structures of the Y. enterocolitica and S. flexneri injectisomes had similar dimensions and were significantly longer than the isolated structures of related injectisomes. The crystal structure of the inner membrane injectisome component YscD appeared elongated compared to a homologous protein, and molecular dynamics simulations documented its elongation elasticity. The ring-shaped secretin YscC at the outer membrane was stretched by 30–40% in situ, compared to its isolated liposome-embedded conformation. We suggest that elasticity is critical for some two-membrane spanning protein complexes to cope with variations in the intermembrane distance.

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

Humans and other animals can use the five senses—touch, taste, sight, smell, and hearing—to interpret the world around them. Single-celled organisms, however, must rely on molecular cues to understand their immediate surroundings. In particular, bacteria gather information about external conditions, including potential hosts nearby, by secreting protein sensors that can relay messages back to the cell.

Bacteria export these sensors via secretion systems that enable the organism both to receive information about the environment and to invade a host cell. A total of seven separate secretion systems, known as types I–VII, have been identified. These different secretion systems handle distinct cargoes, allowing the bacterial cell to respond to a range of feedback from the external milieu.

The type III secretion system, also known as the ‘injectisome’, is found in bacterial species that are enclosed by two membranes separated by a periplasmic space. The injectisome comprises different components that combine to form the basal body, which spans the inner and outer membranes, and a projection from the basal body, called the hollow needle, that mediates the export of cargo from a bacterium to its host or the local environment.

The distance between the inner and outer membranes may vary across species or according to environmental conditions, so the basal body must be able to accommodate these changes. However, no mechanism has yet been established that might introduce such elasticity into the injectisome. Now, Kudryashev et al. have generated three-dimensional structures for the injectisomes of two species of bacteria, Shigella flexneri and Yersinia enterocolitica, and shown that the size of the basal body can fluctuate by up to 20%.

Kudryashev et al. imaged whole injectisomes in these two species and found that the height of the basal body was proportional to the distance between the inner and outer membranes. To probe how this could occur, the properties of two proteins that are important components of the basal body were studied in greater detail. YscD, a protein that extends across the periplasmic space, was crystallized and its structure was then determined and used to develop a computer model to assess its compressibility: this model indicated that YscD could stretch or contract by up to 50% of its total length. The outer membrane component YscC also appeared elastic: when the protein was isolated and introduced into synthetic membranes, its length was reduced 30–40% relative to that observed in intact bacterial membranes.

A further experiment confirmed the adaptability of the basal body: when the separation of the membranes was deliberately increased by placing bacteria in a high-salt medium, the basal body extended approximately 10% in length. Cumulatively, therefore, these experiments suggest that the in-built flexibility of the basal body of the injectisome allows bacteria to adjust to environmental changes while maintaining their sensory abilities and host-invasion potential.

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

Orthobunyavirus ultrastructure and the curious tripodal glycoprotein spike

The genus Orthobunyavirus within the family Bunyaviridae constitutes an expanding group of emerging viruses, which threaten human and animal health. Despite the medical importance, little is known about orthobunyavirus structure, a prerequisite for understanding virus assembly and entry. Here, using electron cryo-tomography, we report the ultrastructure of Bunyamwera virus, the prototypic member of this genus. Whilst Bunyamwera virions are pleomorphic in shape, they display a locally ordered lattice of glycoprotein spikes. Each spike protrudes 18 nm from the viral membrane and becomes disordered upon introduction to an acidic environment. Using sub-tomogram averaging, we derived a three-dimensional model of the trimeric pre-fusion glycoprotein spike to 3-nm resolution. The glycoprotein spike consists mainly of the putative class-II fusion glycoprotein and exhibits a unique tripod-like arrangement. Protein–protein contacts between neighbouring spikes occur at membrane-proximal regions and intra-spike contacts at membrane-distal regions. This trimeric assembly deviates from previously observed fusion glycoprotein arrangements, suggesting a greater than anticipated repertoire of viral fusion glycoprotein oligomerization. Our study provides evidence of a pH-dependent conformational change that occurs during orthobunyaviral entry into host cells and a blueprint for the structure of this group of emerging pathogens.

Orthobunyaviruses comprise a group of emerging arboviruses within the Bunyaviridae, the largest family of membrane-containing viruses. In spite of the continued medical impact upon human and animal health, little is known about orthobunyavirus structure or the process of host cell entry. Here, we address this paucity of information through electron cryo-microscopy analysis of Bunyamwera virus, the prototypic representative of this genus. We reveal that Bunyamwera virions are pleomorphic and display locally-ordered lattices of viral glycoprotein spikes on the envelope surface. The three-dimensional structure of the glycoprotein spike was resolved to 3.0-nm resolution. The spike is composed of the attachment and fusion glycoproteins and comprises a unique tripodal organization. This glycoprotein arrangement contrasts those observed in other virus families. Consistent with the established pH-dependent mechanism of membrane fusion during host cell entry, we provide evidence for the disruption of this tripodal assembly upon exposure to acidic environments. These data constitute a blueprint for orthobunyavirus architecture and support a case for broadened structural diversity within the Bunyaviridae family.

The respiratory syncytial virus nucleoprotein-RNA complex forms a left-handed helical nucleocapsid

Respiratory syncytial virus (RSV) is an important human pathogen. Its nucleocapsid (NC), which comprises the negative sense RNA viral genome coated by the viral nucleoprotein N, is a critical assembly that serves as template for both mRNA synthesis and genome replication. We have previously described the X-ray structure of an NC-like structure: a decameric ring formed of N-RNA that mimics one turn of the helical NC. In the absence of experimental data we had hypothesized that the NC helix would be right-handed, as the N–N contacts in the ring appeared to more easily adapt to that conformation. We now unambiguously show that the RSV NC is a left-handed helix. We further show that the contacts in the ring can be distorted to maintain key N–N-protein interactions in a left-handed helix, and discuss the implications of the resulting atomic model of the helical NC for viral replication and transcription.

Structural biology in situ--the potential of subtomogram averaging

Cryo-electron tomography provides low-resolution 3D views of cells, organelles, or viruses. Macromolecular complexes present in multiple copies can be subsequently identified within the 3D reconstruction (the tomogram), computationally extracted, and averaged to obtain higher resolution 3D structures, as well as a map of their spatial distribution. This method, called subtomogram averaging or subvolume averaging, allows structures of macromolecular complexes to be resolved in situ. Recent applications have provided in situ structural data at resolutions of 2-4 nm on samples including polysomes, nuclear pores, vesicle coats, and viral surface proteins. Here I describe the method and discuss limitations, advances and recent applications. I speculate how the method will solve more structures at higher resolution, allowing in situ structural biology.

Combined approaches to study virus structures

A virus particle must work as a safe box for protecting its genome, but at the same time it has to undergo dramatic conformational changes in order to preserve itself by propagating in a cell infection. Thus, viruses are miniaturized wonders whose structural complexity requires them to be investigated by a combination of different techniques that can tackle both static and dynamic processes. In this chapter we will illustrate how major structural techniques such as X-ray crystallography and electron microscopy have been and can be combined with other techniques to determine the structure of complex viruses. The power of these hybrid method approaches are revealed through the various examples provided.

Cryo-electron microscopy--a primer for the non-microscopist

Cryo-electron microscopy (cryo-EM) is increasingly becoming a mainstream technology for studying the architecture of cells, viruses and protein assemblies at molecular resolution. Recent developments in microscope design and imaging hardware, paired with enhanced image processing and automation capabilities, are poised to further advance the effectiveness of cryo-EM methods. These developments promise to increase the speed and extent of automation, and to improve the resolutions that may be achieved, making this technology useful to determine a wide variety of biological structures. Additionally, established modalities for structure determination, such as X-ray crystallography and nuclear magnetic resonance spectroscopy, are being routinely integrated with cryo-EM density maps to achieve atomic-resolution models of complex, dynamic molecular assemblies. In this review, which is directed towards readers who are not experts in cryo-EM methodology, we provide an overview of emerging themes in the application of this technology to investigate diverse questions in biology and medicine. We discuss the ways in which these methods are being used to study structures of macromolecular assemblies that range in size from whole cells to small proteins. Finally, we include a description of how the structural information obtained by cryo-EM is deposited and archived in a publicly accessible database.

High precision alignment of cryo-electron subtomograms through gradient-based parallel optimization

BACKGROUND

Cryo-electron tomography emerges as an important component for structural system biology. It not only allows the structural characterization of macromolecular complexes, but also the detection of their cellular localizations in near living conditions. However, the method is hampered by low resolution, missing data and low signal-to-noise ratio (SNR). To overcome some of these difficulties and enhance the nominal resolution one can align and average a large set of subtomograms. Existing methods for obtaining the optimal alignments are mostly based on an exhaustive scanning of all but discrete relative rigid transformations (i.e. rotations and translations) of one subtomogram with respect to the other.

RESULTS

In this paper, we propose gradient-guided alignment methods based on two popular subtomogram similarity measures, a real space as well as a Fourier-space constrained score. We also propose a stochastic parallel refinement method that increases significantly the efficiency for the simultaneous refinement of a set of alignment candidates. We estimate that our stochastic parallel refinement is on average about 20 to 40 fold faster in comparison to the standard independent refinement approach. Results on simulated data of model complexes and experimental structures of protein complexes show that even for highly distorted subtomograms and with only a small number of very sparsely distributed initial alignment seeds, our combined methods can accurately recover true transformations with a substantially higher precision than the scanning based alignment methods.

CONCLUSIONS

Our methods increase significantly the efficiency and accuracy for subtomogram alignments, which is a key factor for the systematic classification of macromolecular complexes in cryo-electron tomograms of whole cells.