Computational toolbox for ultrastructural quantitative analysis of filament networks in cryo-ET data

An image processing pipeline for electron cryo-tomography in RELION-5

Electron tomography of frozen, hydrated samples allows structure determination of macromolecular complexes that are embedded in complex environments. Provided that the target complexes may be localised in noisy, three-dimensional tomographic reconstructions, averaging images of multiple instances of these molecules can lead to structures with sufficient resolution for de novo atomic modelling. Although many research groups have contributed image processing tools for these tasks, a lack of standardisation and interoperability represents a barrier for newcomers to the field. Here, we present an image processing pipeline for electron tomography data in RELION-5, with functionality ranging from the import of unprocessed movies to the automated building of atomic models in the final maps. Our explicit definition of metadata items that describe the steps of our pipeline has been designed for interoperability with other software tools and provides a framework for further standardisation.

Architecture of CTPS filament networks revealed by cryo-electron tomography

The cytoophidium is a novel type of membraneless organelle, first observed in the ovaries of Drosophila using fluorescence microscopy. In vitro, purified Drosophila melanogaster CTPS (dmCTPS) can form metabolic filaments under the presence of either substrates or products, and their structures that have been analyzed using cryo-electron microscopy (cryo-EM). These dmCTPS filaments are considered the fundamental units of cytoophidia. However, due to the resolution gap between light and electron microscopy, the precise assembly pattern of cytoophidia remains unclear. In this study, we find that dmCTPS filaments can spontaneously assemble in vitro, forming network structures that reach micron-scale dimensions. Using cryo-electron tomography (cryo-ET), we reconstruct the network structures formed by dmCTPS filaments under substrate or product binding conditions and elucidate their assembly process. The dmCTPS filaments initially form structural bundles, which then further assemble into larger networks. By identifying, tracking, and statistically analyzing the filaments, we observed distinct characteristics of the structural bundles formed under different conditions. This study provides the first systematic analysis of dmCTPS filament networks, offering new insights into the relationship between cytoophidia and metabolic filaments.

Combining per-particle cryo-ET and cryo-EM single particle analysis to elucidate heterogeneous DNA-protein organization

Cryo-electron microscopy single particle analysis (cryo-EM SPA) and cryo-electron tomography (cryo-ET) have historically been employed as distinct approaches for investigating molecular structures of disparate sample types, focusing on highly purified biological macromolecules and in situ cellular contexts, respectively. However, these techniques offer inherently complementary structural insights that, when combined, provide a more comprehensive understanding of complex biological systems. For example, if both techniques are applied to the same purified biological macromolecules, cryo-ET has the ability to resolve highly flexible yet strong signal features on an individual target molecule which will not be preserved in the high-resolution cryo-EM SPA results. In this review, we highlight recent achievements utilizing such applications to unveil new insights into the chromatin assembly and activities of DNA-protein assemblies. This convergence of cryo-EM SPA and cryo-ET holds great promise for elucidating new structural aspects of these essential molecular processes.

Nonmuscle myosin 2 filaments are processive in cells

Directed transport of cellular components is often dependent on the processive movements of cytoskeletal motors. Myosin 2 motors predominantly engage actin filaments of opposing orientation to drive contractile events and are therefore not traditionally viewed as processive. However, recent in vitro experiments with purified nonmuscle myosin 2 (NM2) demonstrated myosin 2 filaments could move processively. Here, we establish processivity as a cellular property of NM2. Processive runs in central nervous system-derived CAD cells are most apparent on bundled actin in protrusions that terminate at the leading edge. We find that processive velocities in vivo are consistent with in vitro measurements. NM2 makes these processive runs in its filamentous form against lamellipodia retrograde flow, though anterograde movement can still occur in the absence of actin dynamics. Comparing the processivity of NM2 isoforms, we find that NM2A moves slightly faster than NM2B. Finally, we demonstrate that this is not a cell-specific property, as we observe processive-like movements of NM2 in the lamella and subnuclear stress fibers of fibroblasts. Collectively, these observations further broaden NM2 functionality and the biological processes in which the already ubiquitous motor can contribute.

Origin and arrangement of actin filaments for gliding motility in apicomplexan parasites revealed by cryo-electron tomography

The phylum Apicomplexa comprises important eukaryotic parasites that invade host tissues and cells using a unique mechanism of gliding motility. Gliding is powered by actomyosin motors that translocate host-attached surface adhesins along the parasite cell body. Actin filaments (F-actin) generated by Formin1 play a central role in this critical parasitic activity. However, their subcellular origin, path and ultrastructural arrangement are poorly understood. Here we used cryo-electron tomography to image motile Cryptosporidium parvum sporozoites and reveal the cellular architecture of F-actin at nanometer-scale resolution. We demonstrate that F-actin nucleates at the apically positioned preconoidal rings and is channeled into the pellicular space between the parasite plasma membrane and the inner membrane complex in a conoid extrusion-dependent manner. Within the pellicular space, filaments on the inner membrane complex surface appear to guide the apico-basal flux of F-actin. F-actin concordantly accumulates at the basal end of the parasite. Finally, analyzing a Formin1-depleted Toxoplasma gondii mutant pinpoints the upper preconoidal ring as the conserved nucleation hub for F-actin in Cryptosporidium and Toxoplasma. Together, we provide an ultrastructural model for the life cycle of F-actin for apicomplexan gliding motility.

Integrated Fluorescence Microscopy (iFLM) for Cryo-FIB-milling and In-situ Cryo-ET

Correlative cryo-FLM-FIB milling is a powerful sample preparation technique for in situ cryo-ET. However, correlative workflows that incorporate precise targeting remain challenging. Here, we demonstrate the development and use of an integrated Fluorescence Light Microscope (iFLM) module within a cryo-FIB-SEM to enable a coordinate-based two-point 3D correlative workflow. The iFLM guided targeting of regions of interest coupled with an automated milling process of the cryo-FIB-SEM instrument allows for the efficient preparation of 9-12 ∼200 nm thick lamellae within 24 hours. Using regular and montage-cryo-ET data collection schemes, we acquired data from FIB-milled lamellae of HeLa cells to examine cellular ultrastructure. Overall, this workflow facilitates on-the-fly targeting and automated FIB-milling of cryo-preserved cells, bacteria, and possibly high pressure frozen tissue, to produce lamellae for downstream cryo-ET data collection.

The organization of double-stranded RNA in the chikungunya virus replication organelle

Alphaviruses are mosquito-borne, positive-sense single-stranded RNA viruses. Amongst the alphaviruses, chikungunya virus is notable as a large source of human illness, especially in tropical and subtropical regions. When they invade a cell, alphaviruses generate dedicated organelles for viral genome replication, so-called spherules. Spherules form as outward-facing buds at the plasma membrane, and it has recently been shown that the thin membrane neck that connects this membrane bud with the cytoplasm is guarded by a two-megadalton protein complex that contains all the enzymatic functions necessary for RNA replication. The lumen of the spherules contains a single copy of the negative-strand template RNA, present in a duplex with newly synthesized positive-sense RNA. Less is known about the organization of this double-stranded RNA as compared to the protein components of the spherule. Here, we analyzed cryo-electron tomograms of chikungunya virus spherules in terms of the organization of the double-stranded RNA replication intermediate. We find that the double-stranded RNA has a shortened apparent persistence length as compared to unconstrained double-stranded RNA. Around half of the genome is present in either of five conformations identified by subtomogram classification, each representing a relatively straight segment of ~25-32 nm. Finally, the RNA occupies the spherule lumen at a homogeneous density, but has a preferred orientation to be perpendicular to a vector pointing from the membrane neck towards the spherule center. Taken together, this analysis lays another piece of the puzzle of the highly coordinated alphavirus genome replication.

Bringing Structure to Cell Biology with Cryo-Electron Tomography

Recent advances in cryo-electron microscopy have marked only the beginning of the potential of this technique. To bring structure into cell biology, the modality of cryo-electron tomography has fast developed into a bona fide in situ structural biology technique where structures are determined in their native environment, the cell. Nearly every step of the cryo-focused ion beam-assisted electron tomography (cryo-FIB-ET) workflow has been improved upon in the past decade, since the first windows were carved into cells, unveiling macromolecular networks in near-native conditions. By bridging structural and cell biology, cryo-FIB-ET is advancing our understanding of structure-function relationships in their native environment and becoming a tool for discovering new biology.

Unleashed Actin Assembly in Capping Protein-Deficient B16-F1 Cells Enables Identification of Multiple Factors Contributing to Filopodium Formation

Background: Filopodia are dynamic, finger-like actin-filament bundles that overcome membrane tension by forces generated through actin polymerization at their tips to allow extension of these structures a few microns beyond the cell periphery. Actin assembly of these protrusions is regulated by accessory proteins including heterodimeric capping protein (CP) or Ena/VASP actin polymerases to either terminate or promote filament growth. Accordingly, the depletion of CP in B16-F1 melanoma cells was previously shown to cause an explosive formation of filopodia. In Ena/VASP-deficient cells, CP depletion appeared to result in ruffling instead of inducing filopodia, implying that Ena/VASP proteins are absolutely essential for filopodia formation. However, this hypothesis was not yet experimentally confirmed. Methods: Here, we used B16-F1 cells and CRISPR/Cas9 technology to eliminate CP either alone or in combination with Ena/VASP or other factors residing at filopodia tips, followed by quantifications of filopodia length and number. Results: Unexpectedly, we find massive formations of filopodia even in the absence of CP and Ena/VASP proteins. Notably, combined inactivation of Ena/VASP, unconventional myosin-X and the formin FMNL3 was required to markedly impair filopodia formation in CP-deficient cells. Conclusions: Taken together, our results reveal that, besides Ena/VASP proteins, numerous other factors contribute to filopodia formation.

Deciphering the molecular mechanisms of actin cytoskeleton regulation in cell migration using cryo-EM

The actin cytoskeleton plays a key role in cell migration and cellular morphodynamics in most eukaryotes. The ability of the actin cytoskeleton to assemble and disassemble in a spatiotemporally controlled manner allows it to form higher-order structures, which can generate forces required for a cell to explore and navigate through its environment. It is regulated not only via a complex synergistic and competitive interplay between actin-binding proteins (ABP), but also by filament biochemistry and filament geometry. The lack of structural insights into how geometry and ABPs regulate the actin cytoskeleton limits our understanding of the molecular mechanisms that define actin cytoskeleton remodeling and, in turn, impact emerging cell migration characteristics. With the advent of cryo-electron microscopy (cryo-EM) and advanced computational methods, it is now possible to define these molecular mechanisms involving actin and its interactors at both atomic and ultra-structural levels in vitro and in cellulo. In this review, we will provide an overview of the available cryo-EM methods, applicable to further our understanding of the actin cytoskeleton, specifically in the context of cell migration. We will discuss how these methods have been employed to elucidate ABP- and geometry-defined regulatory mechanisms in initiating, maintaining, and disassembling cellular actin networks in migratory protrusions.

ArpC5 isoforms regulate Arp2/3 complex-dependent protrusion through differential Ena/VASP positioning

Regulation of the Arp2/3 complex is required for productive nucleation of branched actin networks. An emerging aspect of regulation is the incorporation of subunit isoforms into the Arp2/3 complex. Specifically, both ArpC5 subunit isoforms, ArpC5 and ArpC5L, have been reported to fine-tune nucleation activity and branch junction stability. We have combined reverse genetics and cellular structural biology to describe how ArpC5 and ArpC5L differentially affect cell migration. Both define the structural stability of ArpC1 in branch junctions and, in turn, by determining protrusion characteristics, affect protein dynamics and actin network ultrastructure. ArpC5 isoforms also affect the positioning of members of the Ena/Vasodilator-stimulated phosphoprotein (VASP) family of actin filament elongators, which mediate ArpC5 isoform-specific effects on the actin assembly level. Our results suggest that ArpC5 and Ena/VASP proteins are part of a signaling pathway enhancing cell migration.

Ena/VASP clustering at microspike tips involves lamellipodin but not I-BAR proteins, and absolutely requires unconventional myosin-X

Sheet-like membrane protrusions at the leading edge, termed lamellipodia, drive 2D-cell migration using active actin polymerization. Microspikes comprise actin-filament bundles embedded within lamellipodia, but the molecular mechanisms driving their formation and their potential functional relevance have remained elusive. Microspike formation requires the specific activity of clustered Ena/VASP proteins at their tips to enable processive actin assembly in the presence of capping protein, but the factors and mechanisms mediating Ena/VASP clustering are poorly understood. Systematic analyses of B16-F1 melanoma mutants lacking potential candidate proteins revealed that neither inverse BAR-domain proteins, nor lamellipodin or Abi is essential for clustering, although they differentially contribute to lamellipodial VASP accumulation. In contrast, unconventional myosin-X (MyoX) identified here as proximal to VASP was obligatory for Ena/VASP clustering and microspike formation. Interestingly, and despite the invariable distribution of other relevant marker proteins, the width of lamellipodia in MyoX-KO mutants was significantly reduced as compared with B16-F1 control, suggesting that microspikes contribute to lamellipodium stability. Consistently, MyoX removal caused marked defects in protrusion and random 2D-cell migration. Strikingly, Ena/VASP-deficiency also uncoupled MyoX cluster dynamics from actin assembly in lamellipodia, establishing their tight functional association in microspike formation.

Spaghetti Tracer: A Framework for Tracing Semiregular Filamentous Densities in 3D Tomograms

Within cells, cytoskeletal filaments are often arranged into loosely aligned bundles. These fibrous bundles are dense enough to exhibit a certain regularity and mean direction, however, their packing is not sufficient to impose a symmetry between-or specific shape on-individual filaments. This intermediate regularity is computationally difficult to handle because individual filaments have a certain directional freedom, however, the filament densities are not well segmented from each other (especially in the presence of noise, such as in cryo-electron tomography). In this paper, we develop a dynamic programming-based framework, Spaghetti Tracer, to characterizing the structural arrangement of filaments in the challenging 3D maps of subcellular components. Assuming that the tomogram can be rotated such that the filaments are oriented in a mean direction, the proposed framework first identifies local seed points for candidate filament segments, which are then grown from the seeds using a dynamic programming algorithm. We validate various algorithmic variations of our framework on simulated tomograms that closely mimic the noise and appearance of experimental maps. As we know the ground truth in the simulated tomograms, the statistical analysis consisting of precision, recall, and F1 scores allows us to optimize the performance of this new approach. We find that a bipyramidal accumulation scheme for path density is superior to straight-line accumulation. In addition, the multiplication of forward and backward path densities provides for an efficient filter that lifts the filament density above the noise level. Resulting from our tests is a robust method that can be expected to perform well (F1 scores 0.86-0.95) under experimental noise conditions.

Quantitative Cryo-Electron Tomography

The three-dimensional organization of biomolecules important for the functioning of all living systems can be determined by cryo-electron tomography imaging under native biological contexts. Cryo-electron tomography is continually expanding and evolving, and the development of new methods that use the latest technology for sample thinning is enabling the visualization of ever larger and more complex biological systems, allowing imaging across scales. Quantitative cryo-electron tomography possesses the capability of visualizing the impact of molecular and environmental perturbations in subcellular structure and function to understand fundamental biological processes. This review provides an overview of current hardware and software developments that allow quantitative cryo-electron tomography studies and their limitations and how overcoming them may allow us to unleash the full power of cryo-electron tomography.

Capturing actin assemblies in cells using in situ cryo-electron tomography

Actin contributes to an exceptionally wide range of cellular processes through the assembly and disassembly of highly dynamic and ordered structures. Visualizing these structures in cells can help us understand how the molecular players of the actin machinery work together to produce force-generating systems. In recent years, cryo-electron tomography (cryo-ET) has become the method of choice for structural analysis of the cell interior at the molecular scale. Here we review advances in cryo-ET workflows that have enabled this transformation, especially the automation of sample preparation procedures, data collection, and processing. We discuss new structural analyses of dynamic actin assemblies in cryo-preserved cells, which have provided mechanistic insights into actin assembly and function at the nanoscale. Finally, we highlight the latest visual proteomics studies of actin filaments and their interactors reaching sub-nanometer resolutions in cells.

Mechanistic insights into actin force generation during vesicle formation from cryo-electron tomography

Actin assembly provides force for a multitude of cellular processes. Compared to actin-assembly-based force production during cell migration, relatively little is understood about how actin assembly generates pulling forces for vesicle formation. Here, cryo-electron tomography identified actin filament number, organization, and orientation during clathrin-mediated endocytosis in human SK-MEL-2 cells, showing that force generation is robust despite variance in network organization. Actin dynamics simulations incorporating a measured branch angle indicate that sufficient force to drive membrane internalization is generated through polymerization and that assembly is triggered from ∼4 founding "mother" filaments, consistent with tomography data. Hip1R actin filament anchoring points are present along the entire endocytic invagination, where simulations show that it is key to pulling force generation, and along the neck, where it targets filament growth and makes internalization more robust. Actin organization described here allowed direct translation of structure to mechanism with broad implications for other actin-driven processes.

Cryo-electron tomography of the onion cell wall shows bimodally oriented cellulose fibers and reticulated homogalacturonan networks

One hallmark of plant cells is their cell wall. They protect cells against the environment and high turgor and mediate morphogenesis through the dynamics of their mechanical and chemical properties. The walls are a complex polysaccharidic structure. Although their biochemical composition is well known, how the different components organize in the volume of the cell wall and interact with each other is not well understood and yet is key to the wall's mechanical properties. To investigate the ultrastructure of the plant cell wall, we imaged the walls of onion (Allium cepa) bulbs in a near-native state via cryo-focused ion beam milling (cryo-FIB milling) and cryo-electron tomography (cryo-ET). This allowed the high-resolution visualization of cellulose fibers in situ. We reveal the coexistence of dense fiber fields bathed in a reticulated matrix we termed "meshing," which is more abundant at the inner surface of the cell wall. The fibers adopted a regular bimodal angular distribution at all depths in the cell wall and bundled according to their orientation, creating layers within the cell wall. Concomitantly, employing homogalacturonan (HG)-specific enzymatic digestion, we observed changes in the meshing, suggesting that it is-at least in part-composed of HG pectins. We propose the following model for the construction of the abaxial epidermal primary cell wall: the cell deposits successive layers of cellulose fibers at -45° and +45° relative to the cell's long axis and secretes the surrounding HG-rich meshing proximal to the plasma membrane, which then migrates to more distal regions of the cell wall.

Quantitative mapping of keratin networks in 3D

Mechanobiology requires precise quantitative information on processes taking place in specific 3D microenvironments. Connecting the abundance of microscopical, molecular, biochemical, and cell mechanical data with defined topologies has turned out to be extremely difficult. Establishing such structural and functional 3D maps needed for biophysical modeling is a particular challenge for the cytoskeleton, which consists of long and interwoven filamentous polymers coordinating subcellular processes and interactions of cells with their environment. To date, useful tools are available for the segmentation and modeling of actin filaments and microtubules but comprehensive tools for the mapping of intermediate filament organization are still lacking. In this work, we describe a workflow to model and examine the complete 3D arrangement of the keratin intermediate filament cytoskeleton in canine, murine, and human epithelial cells both, in vitro and in vivo. Numerical models are derived from confocal airyscan high-resolution 3D imaging of fluorescence-tagged keratin filaments. They are interrogated and annotated at different length scales using different modes of visualization including immersive virtual reality. In this way, information is provided on network organization at the subcellular level including mesh arrangement, density and isotropic configuration as well as details on filament morphology such as bundling, curvature, and orientation. We show that the comparison of these parameters helps to identify, in quantitative terms, similarities and differences of keratin network organization in epithelial cell types defining subcellular domains, notably basal, apical, lateral, and perinuclear systems. The described approach and the presented data are pivotal for generating mechanobiological models that can be experimentally tested.