Simplified Volumetric Models as an Effective Strategy for Segmenting Actin Networks in Cryo-Electron Tomograms

Summary

Here, we present a protocol to place simplified volumetric models into noisy, complex, tomographic 3D volumes. This allows the fast segmentation of actin filament densities, detection of systematic filament bending and of gaps in hair bundle filaments, as well as convenient quantification of volumetric model properties, such as distances.

Abstract

Efficient methods for the extraction of features of interest remain one of the biggest challenges for the interpretation of cryo-electron tomograms. Various automated approaches have been proposed, many of which work well for high-contrast datasets where the features of interest can be easily detected and are clearly separated from one another. Our inner ear stereocilia cryo-electron tomographic datasets are characterized by a dense array of hexagonally packed actin filaments that are frequently cross-connected. These features make automated segmentation very challenging, further aggravated by the high-noise environment of cryo-electron tomograms and the high complexity of the densely packed features. Using prior knowledge about the actin bundle organization, we have placed layers of a highly simplified ball-and-stick actin model to first obtain a global fit to the density map, followed by regional and local adjustments of the model. We show that volumetric model building not only allows us to deal with the high complexity, but also provides precise measurements and statistics about the actin bundle. Volumetric models also serve as anchoring points for local segmentation, such as in the case of the actin-actin cross connectors. Volumetric model building, particularly when further augmented by computer-based automated fitting approaches, can be a powerful alternative when conventional automated segmentation approaches are not successful.

Introduction

Cryo-electron tomography allows entire organelles or parts of cells and tissues to be visualized at nanometer resolutions in their near-native state^1,2,3 using either plunge-freezing⁴ or high-pressure freezing ultra-rapid vitrification⁵. Since only a limited electron dose can be tolerated by the cryo-preserved, unstained, frozen-hydrated sample, the tomographic 3D data is very noisy. This noise often can be significantly reduced by a variety of noise filtering algorithms^6,7, including nonlinear anisotropic diffusion⁸, bilateral filtering⁹, and recursive median filtering¹⁰.

Furthermore, tilt limitations of the microscope stage, resulting in a missing wedge of information, and the fact that the specimen thickness increases at high tilt angles, lead to 3D reconstructions with anisotropic resolution. This means that the density is smeared out in the third dimension due to lower resolution in the Z-direction. As a result, the shape of macromolecules appears distorted (i.e., less well defined and elongated in the third dimension).

Among the biggest challenges in the interpretation of tomographic data is the automated extraction of the relevant features, also known as segmentation¹¹. With sufficient unique shape features and low noise, macromolecular machines in complex 3D volumes can be identified by template matching^12,13,14; however, the success of template matching depends on tomogram resolution, a suitable search model, as well size and shape characteristics of the feature volumes. If the features of interest are sufficiently spaced apart and repeating motifs (such as large macromolecular machines) can readily be identified, tomogram subvolumes can be combined to increase the signal-to-noise ratio and to average out individual particle shape distortions. Automated segmentation of an actin filament network in electron tomograms of the thin edge of frozen-hydrated Dictyostelium discoideum cells by template matching has been reported¹⁵.

However, if features of interest are closely spaced, the data resolution anisotropy can lead to a smearing out of the map densities in the Z-direction (along the direction of the electron beam), resulting in an apparent merging of the density envelope of closely spaced macromolecular machines or supramolecular complexes. In such cases, automated approaches for segmentation, such as watershed¹⁶, boundary segmentation¹⁷, or a variety of machine learning-based classification approaches^18,19, may not be able to recognize the features of interest or establish a correct boundary around an object of interest. Often, one ends up with either a few very large pieces or with a heavily over-segmented volume, where much effort is needed to merge many small pieces until the feature of interest is perceived to be complete. Such manual curation of segmentation results can be very labor-intensive and may even fail altogether when the structure of interest is an array of closely spaced filaments that are interconnected via short linkers. In this gigantic network of filamentous structures, it can be difficult to orient oneself. This is because, due to resolution anisotropy, densities appear to blend into one another, presenting a formidable challenge for both automated as well as for interactive manual segmentation approaches. As a consequence, one can easily "jump" between filaments when only visually inspecting small regions.

Fortunately, in the case of the actin bundle in inner ear hair cell stereocilia, we have knowledge about the overall actin bundle organization and the directionality of the actin filaments^20,21. The actin bundle consists of hundreds of hexagonally, densely-packed actin filaments 6-8 nm in diameter, which are spaced about 12-13 nm apart from one another²².

This allowed us to take a rather different approach to segmentation that is based on simplified ball-and-stick models to represent actin filaments. The strategy involved simultaneously placing an idealized regular array of filament models into slabs of the cryo-electron tomography density maps to build up a 3D model of the actin bundle layer by layer. We ensured that the model had a general overall fit to the density map before making local adjustments to individual filament models or groups of filament models to closely match the density map. By automatic color-coding of the map density value at the filament model location, we were able to easily detect apparent gaps in the actin bundle. Volumetric models allow a quantitative analysis of volumetric properties, such as distances between actin filaments, and also lead to a simplified display of the overall 3D filamentous network organization.

In addition, models can also serve as anchoring structures for the segmentation of additional features, such as actin-actin linkers, as (portions of) individual filament models can be selected, around which appropriate radius map density zones can be generated for inspection and further segmentation.

We believe that our volumetric model-based segmentation approach is particularly useful for large filamentous structure networks that may contain gaps and inter-filament cross-connections. Segmentation algorithms tend to operate locally, whereas the human brain takes larger areas into consideration, and thus is superior to computers when it comes to recognizing filament structures, even in a complex, high-noise environment.

Protocol

The protocol follows the guidelines of Southeast University's human research ethics committee.

1. Cryo-electron tomography data source for volumetric model building

NOTE: The stereocilia cryo-electron tomographic reconstructions used for a volumetric model building have been previously published^22,23, and were obtained as described by Metlagel et al.²².

The UCSF Chimera python scripts for stereocilia modeling are provided in Supplementary File 1, Supplementary Coding File 1, Supplementary Coding File 2, Supplementary Coding File 3, Supplementary Coding File 4, and Supplementary Coding File 5.

1. Briefly, blot stereocilia from the apical surface of mouse utricle sensory epithelia onto the lacey carbon support film of an electron microscope (EM) grid. Then, perform ultra-rapid plunge-freezing vitrification and single-axis cryo-electron tomographic data collection on a cryogenic transmission electron microscopy (cryo-TEM), operated at 300 kV with a nominal defocus of 3.5-4.5 µm, using a CMOS-type camera in integration mode, with a 0.47 to 0.59 nm pixel size (see Table of Materials).
2. A typical dose for single-axis data collection is typically 80-100 electrons/Ų. Reconstruct the tomogram 3D volumes using the software package IMOD²⁴, either by weighted back-projection or the SIRT method^25,26.
3. Denoise the tomograms using recursive median or bilateral filtering in Priism²⁷ or the nonlinear anisotropic diffusion option in IMOD.
   NOTE: The key software package used for interactive visualization, volumetric model building, and quantitative analysis was the UCSF Chimera software^28,29. All software packages used are listed in the Table of Materials.

2. Cryo-electron tomography data preparation for volumetric model building

1. Tomogram rotation
   NOTE: The purpose of this step is to find the "optimal angles" in which the actin core is aligned in two of the three axes (x and z), so that placing the models in 3D can be done with only one axis to worry about. The following steps are carried out within the IMOD software package.
   1. Open 3dmod (IMOD) and open the 3D image stack file in .mrc file format. Make sure the 3D stack is in 16-bit grayscale mode and displays the correct X, Y, Z dimensions in its image file header. If necessary, make corrections using the alterheader command: alterheader -d (x_pixelspacing),(y_pixelspacing),(z_pixelspacing) inputfilename.mrc in the command prompt. If the 3D stack is in TIFF file format, use the command tif2mrc inputfilename.tif outputfilename.mrc in the command prompt to create an .mrc file.
   2. Visually inspect the tomogram and open "slicer" (by pressing backslash "\") to find the best angles of rotation in the X, Y, and Z dimensions for alignment of the actin filament plane with the Z-plane. Manipulate the X rotation, Y rotation, Z rotation bar to find out the optimal angle at which the actin filaments are aligned on the X- and Z-planes. Note that changing the average thickness for the display of the tomogram (Thick : Img) can help in finding an ideal amount of averaging of the density in the Z-plane and thus increased contrast. Zoom in and out of the tomogram by utilizing the + and - keys, respectively, on the keyboard.
   3. Once the ideal angles for rotation have been identified, rotate the tomogram by typing the IMOD command rotatevol -a (Z rotation angles),(Y rotation angles),(X rotation angles) -s (x-width),(y-height),(z-depth) inputfilename.mrc outputfilename.mrc into the command prompt. Make sure to give the map plenty of space to accommodate for the rotation to avoid cutting off portions of the map upon rotation.
      NOTE: For more information on how to use the rotatevol command, use the option from the help menu, or visit the following URL (https://bio3d.colorado.edu/imod/doc/man/rotatevol.html). This command (rotatevol) does not alter the currently shown map in the IMOD window, but rather creates a new map file with a density map being rotated.
   4. Using the slicer window, identify the X, Y, and Z coordinates of two opposite corners of a map to identify the cropping region containing the stereocilium using the trimvol command.
   5. Use the IMOD command Trimvol to crop the map at the previously obtained coordinates, by typing trimvol -x (x-coordiante 1),( x-coordiante 2) -y (y-coordiante 1),( y-coordiante 2) -z (z-coordiante 1),( z-coordiante 2) inputfilename.mrc outputfilename.mrc into the command prompt. Now, the map is much smaller, and thus easier and faster to work with from this point onward.
      NOTE: For more information on how to use the trimvol command, use the option from the help menu, or visit the following URL (https://bio3d.colorado.edu/imod/doc/man/trimvol.html). These steps create a .mrc file that has been rotated and trimmed. 
2. Tomogram filtering
   NOTE: This step uses a nonlinear anisotropic diffusion (NAD) filter provided by IMOD for noise reduction. We used as initial value the ones recommended in the IMOD help page for non-linear anisotropic diffusion filtering, as listed in the following URL (https://bio3d.colorado.edu/imod/doc/NADexample.html), and was used as the initial values.
   1. Type etomo in the command prompt to initiate IMOD's etomo graphic user interface.
   2. Select the Nonlinear Anisotropic Diffusion option from the etomo menu and select the map file (Pick a volume).
   3. In order to find a K value and the number of iterations that provide the best filtering results, as judged by visual inspection, apply the NAD filter to a smaller test volume to extract the specified NAD-filter test volume. Click the Rubberband Tool on top of the 3dmod interface, left-click drag to select test volume, and click the Hi and Lo buttons to indicate the Z-slice to be the boundary of the test volume (Get Test Volume Range from 3dmod). Extract the test volume by clicking Extract Test Volume.
   4. Filter the test volume using various K values (Find a K Value for Test Volume > List of K values), such as 0.1, 1,5, 10, 15, 25, 30, 50, and 75. Perform NAD filtering with each of the specified K values (Run with Different K Values > View Different K Values Test Results). Use the slicer window to evaluate the performance of the NAD filter and select a K value that will be used for filtering the entire volume.
   5. Use the K value identified in step 2.2.4 and test different numbers of iterations, such as 2, 5, 8, 11, 15, and 21 (Find an Iteration Number for Test Volume > List of iterations). Evaluate the NAD filter performance for different iterations by either clicking View Different Iteration Test Results or using a slicer tool. Choose a value for a different number of iterations.
   6. Use the identified K value and the number of iterations (Filter Full Volume > K value and Iterations) and filter the entire volume (Filter Full Volume). A newly filtered volume will appear in the same directory where the loaded full volume is located with a .nad file extension attached to the .mrc file.
      NOTE: These steps create a .mrc file that has been filtered with a NAD filter and displays a lower noise and thus an improved signal-to-noise ratio.

3. Volumetric model building

1. 3D model placement preparation
   NOTE: The purpose of this step is to prepare and analyze the map for placing 3D volumetric models into the rotated and denoised density map using the UCSF Chimera program.
   1. Using the slicer window, obtain a cross-sectional view (slab of 30 slices/28.4 nm) to identify the Z-coordinates of the center of the density slab, into which the actin filament model is to be placed. Click on the center of the filament density and note the Z-coordinates in the IMOD main window.
   2. Keep note of the coordinates identified for future model placement in UCSF Chimera.
      NOTE: The Z-coordinates noted from this step is used during step 4 to easily find and place pre-made actin filament models at the correct Z-heights, making the model placement faster.
2. 3D model placement
   NOTE: The purpose of this step is to place a 3D volumetric model into the prepared density map using UCSF Chimera.
   1. In UCSF Chimera, open the filtered, rotated map (File > Open)
   2. Check whether the parameters of the map are set correctly by selecting Features > Planes and Feature > Coordinates in the Volume Viewer to open the feature and coordinates interface. Check Voxel Size to see if the correct voxel spacing is set by using the header command in the command prompt (in IMOD). If not, correct the Voxel Size interface to the correct voxel spacing. Center the map to set the camera at the center of the viewing window (Origin Index > Center).
   3. Open the Camera Control window (Tools > Viewing Control > Camera). Then, set the camera viewing to the orthographic view (projection > orthographic)
   4. Load the pre-made python script model (File > Open > ActinFilamentPlane.py) and open as many actin filament plane scripts as there are actin filament planes in the map.
   5. Relocate the loaded models by turning on the Movement Mouse Mode (Tools > Movement > Movement Mouse Mode). Enable the movement mouse option (Set Movement Mouse Mode > Move Selection). Use Ctrl + Drag to multi-select any models present in the session and right-click (or any other key that is specified in Preference > Mouse) to move the models in the desired direction.
      NOTE: The models have been created with balls and sticks with certain thickness (radius of ActinFilamentPlane.py script) and with certain spacings (Lattice_Spacing of ActinFilamentPlane.py script).
   6. To visualize just the model that one is interested in, open Model Panel (Tools > General Controls > Model Panel). Click checkboxes under S (Show) to either enable or disable the visibility of a given model.
   7. Activate the command line panel (Tools > General Controls > Command Line).
   8. After placing each actin filament plane model to its proper Z-height, delete any excess actin filaments in each actin filament plane model by selecting the excess portion of the models, and by typing del sel (Delete Selected) in the command line at the bottom of the UCSF Chimera main view window.
   9. Only visualize one actin filament plane at a time, by going to the Model Panel window and selecting and hiding all other actin filament models (Ctrl + Click > Hide).
   10. To visualize portions of the map density that are in the immediate proximity of the single actin filament plane model of interest, select the model (Model Panel > select) and open the zone control panel in Volume Viewer (Volume Viewer > Features > Zone). Set the zone radius to 100 Å, then click Zone to zone the map in a radius of 100 Å (= 10 nm) from the selected model. Note that the dimension in Chimera is set in Angstrøm (1 Å = 0.1 nm).
   11. If needed, adjust the model to fit to the density map by turning on the Movement Mouse Mode. Multi-select any models that are not in the proper position to move the models to be properly placed inside the zoned density map.
   12. To overcome a programming bug in ActinFilamentPlane.py, run the FixingMarkerID.py script, which fixes the missing MarkerID of the "atoms" (balls in the ball-and-stick model). Select the entire model and type runscript FixingMarkerID.py in the command line (the script should be specified by correct directory, such as C: \ directory \ Script.py).
   13. If the density map indicates that the actin filament is curved, make sure to divide the straight actin filament into multiple parts to accommodate for the curvature by selecting the desired models and typing runscript dividelinks.py number. Additional markers (specified by the number in the command) are placed inside the actin model equidistant from one another, making "joints" where the model can be bent.
   14. Adjust the joints to reflect the curvature by moving additional marker "joints" along the zoned map by using Movement Mouse Mode. Add additional markers to the model as needed by using the dividelinks.py script.
   15. To verify that all the all the actin filaments are properly placed inside the actin densities, use clipping bars by clicking Tools > Viewing Control > Camera in the UCSF Chimera drop-down menu to open the Camera Control window, then select Side View to see the two clipping bars (near- and far-clipping planes) of the camera. Drag two yellow bars into a narrow slit by left-clicking and drag the clipping portion by using the Middle Mouse Button to inspect the models one small bit at a time.
       NOTE: The stereocilium actin core model is now complete.
3. Membrane segmentation
   NOTE: The purpose of this step is to create a surface model of the stereocilia membrane.
   1. The segmentation of surfaces utilizes the Volume Tracer tool (Tools > Volume Data > Volume Tracer).
   2. In the Volume Tracer window, click Mouse, and only check the options Place Markers on Data Planes and Link new marker to selected marker. This allows precisely placing the markers on the same Z-plane (when a single Z-plane is shown). Also check Place Markers Using the Middle Mouse Button in the Volume Tracer window.
   3. Create a new marker set (File > New Marker set). Marker set 1 should be created and selected on the Volume Tracer window.
   4. Use Volume Viewer to display single Z-planes in the main viewing window (Axis to Z > One.)
   5. Using the Middle Mouse Button, trace out the membrane from one end of the visible membrane to another. All markers should be connected by links from one end to another.
   6. De-select the last-placed marker by pressing Ctrl + Click on the background where not a single selectable item is located. This ensures that the next marker to be placed is connected with links in the new plane, and not with a previous membrane trace of another Z-plane.
   7. Repeat the membrane tracing for every 10 nm in the same direction of tracing from one end of the visible membrane to another side.
   8. Once many parallel bands of membranes are segmented, go to the Volume Tracer window and click Features > Surfaces. This enables the surface interface in the Volume Tracer window. Click Create next to the Surfaces interface to create surfaces in between bands of membrane segmentation. This creates a thin flap of the surface model that fills the spaced-out membrane traces.
      NOTE: The stereocilium membrane surface model creation is now complete.
4. Crosslinker modelling
   NOTE: The purpose of this step is to model the crosslinkers between actin filaments.
   1. Crosslinkers can also be placed by utilizing the Volume Tracer tool.
   2. In the Volume Tracer window, left-click to check only Place Markers on Surfaces and Link new marker to selected marker. Then, check Place Markers Using the Middle Mouse Button in the Volume Tracer window.
   3. Only visualize densities around the single actin filament plane model of interest by using the zone function. Set the zone radius to 100 Å to visualize the map within a radius of 100 Å from the selected model.
   4. Place markers across visible densities of crosslinkers by clicking the Middle Mouse Button. Once two markers are placed (and thus a single crosslinker is modeled) make sure to de-select the second marker to prevent links continuously being connected to the next crosslinker model.
   5. Repeat the steps above for every actin filament plane and repeat the entire above process to all three major directions in which crosslinkers may be located. When zoning actin filament planes in the diagonal direction, manually select actin filament models individually.
      NOTE: For each of the actin crosslinker models, we can record the position on the actin filaments to which the putative crosslinkers connect. This information could, in principle, be extracted for a follow-up analysis using sub-tomogram averaging, which is beyond the scope of this study.

4. Quantitative analysis of the 3D model

1. Gap analysis
   NOTE: The purpose of this step is to detect gaps in the actin filaments.
   1. Selected items can be analyzed. Specify their parameters, like their color, using the Inspect Selection tool (Action > Inspect).
   2. The color of the selected item can be specified by selecting any selectable item in UCSF Chimera, and by changing the color next to the Color section of the Inspect Selection window.
   3. Specify the size of the markers and links by selecting any markers and links in UCSF Chimera, and by changing the values next to the radius section of the Inspect Selection window, where the units are in Angstrøm.
   4. Perform gap analysis by auto-detection of density gaps, utilizing the Values at Atom Positions function (Tools > Volume Data > Values at Atom Positions).
   5. Select the map in which the values will be read by selecting the desired map next to the Volume Data section of the Values at Atom Positions window. Select the model in which all the markers will be assigned a parameter by selecting the desired model next to the Molecule section.
   6. With the Histogram button at the bottom of the Values at Atom Positions window, a histogram is created in a new window termed Render/Select by Attribute. Ctrl + Click on the histogram will add or remove cutoff points, for which different colors can be selected. Uncheck Keep Opaque to apply a transparent color to the scene. Markers with a density below a given threshold, which represents gaps, can not only be colored but can also be selected, and thus allow us to count the number markers in the gap region. To select, click on Select > Attribute > value_mapname.
2. Determining inter-actin distances
   1. For measuring inter-actin distances, place markers on a single XZ-plane (referred to in UCSF Chimera as a Y-plane) by using the Volume Tracer window, clicking Mouse, and checking only Place Markers on Data Planes and unchecking Link new marker to selected marker. This allows precisely placing the markers on the same Y-plane (when a single Z-plane is shown). Also check Place Markers Using Middle Mouse Button in the Volume Tracer window and place the markers on a target Y-plane of interest for inter-actin distance measurement.
   2. Once single Y-plane's actin positions are successfully marked, select all the markers and type the command findclash #model-spec test self overlapCutoff -200 hbondAllowance 0 log true Linewidth 10 pbColor black in the command line at the bottom of the actin filament plane. This finds all inter-marker clashes within a distance of 200 Å (overlapCutoff -200), and creates black lines (pbColor black) with a thickness of 10 Å (Linewidth 10) between the markers.
   3. Since there are overlapping markers within the given range of 200 Å, remove redundant measurements that are not links between only the closest, adjacent markers. Do this by using the RemoveCross.py script provided on UCSF Chimera website, which removes all the overlapping links that are not the shortest link with its adjacent markers.
   4. To measure the distance of the bonds that are remaining, select all pseudobonds (links that are created by findclash) and use the python script pblengths.py by typing runscript pblengths.py. See the measured values from the Reply Log (Tools > Utilities > Reply Log). Then, copy the obtained values from the Reply Log and transfer to a spreadsheet or any other suitable tool to create a histogram.

Discussion

We have shown that automated approaches for segmentation, such as watershed segmentation, can fail in the high-noise and high-complexity environment of hair cell stereocilia cryo-electron tomograms. Distinguishing which part of this filamentous network represents actin filaments and what constitutes crosslinks on a local environment level seems challenging at best when just inspecting small tomographic subvolumes. The model building approach used in this study benefits from prior knowledge of the large-scale order of the actin bundle, which helps in developing an expectation about the orientation of actin filaments and the crosslinker densities. Perhaps even more significant is that a human brain can easily find patterns by considering the larger context beyond the local density distribution, whereas a computer algorithm only works for a relatively small region that is considered by the algorithm; hence, larger-scale trends cannot easily be taken into consideration. By fitting a model globally to a layer of density, we avoided the confusion that can occur when attempting to create a model for small portions of single actin filaments at a time. Of course, such global fitting assumes an order that extends over large distances. However, since we had an unexpected small but significant gradual bending of the actin filaments, the global fit was only an initial approximation, and required local adjustments of the model to fit the density map. Since the initial model was a good starting point, the adjustments could be made with high confidence. One big advantage of our approach was that we could choose to display only a defined zone of density, which helped to reduce the complexity of the scenery. Furthermore, viewing of the map density slab along the filament model axis helped to identify the unexpected curvature, which we would have most likely missed when simply displaying smaller subvolumes. Placement of the initial model also facilitated quick zooming in and out, to alternate between an overall view of the respective layer of actin filaments and the detailed views to make model adjustments.

The critical steps within the protocol included the rotation of the map after visual inspection, the creation and placement of the model into the density map, as well as the division of the filament model into smaller segments. The atom position of the segments could then be spatially adjusted to fit the density map, and/or color-coded to detect gaps.

This approach of actin model building can also be modified by placing a set of "atoms" (i.e., the balls of the ball-and-stick model) into filament densities by using a cross-sectional view of a 10-30 slice/9.47-28.4 nm averaged slab of density, which then can be connected by bonds (i.e., the sticks of the ball-and-stick model). We have used this approach, which is a modification from the protocol described here in detail, for the volumetric model building in the taper region of hair cell stereocilia²³. Furthermore, as we have described here, our volumetric model building approach is also well suited for the segmentation and model building of membranes.

While volumetric model building can be applied to any density map that shows filamentous features, the technique we have described here is most efficient when we have an array of regularly spaced filaments, for which a global fit of a volumetric model can be obtained. It also depends on the filamentous features to alter their directionality in a gradual manner. If there are sudden kinks and sharp turns in the filamentous structures, our approach may not be particularly helpful for segmentation.

In the meantime, our collaborators have developed an automated approach for automated filament tracing that follows a similar concept used here for manual segmentation^30,31. Going forward, the best approach may well be a hybrid of manual identification and placement of an initial sparse model (even only a few balls) into the density as a starting point, and then letting a search and fitting algorithm finish the tracing of the filaments.

Simplified volumetric models reduce the complexity of a system and allows certain patterns, such as actin filament bending near the tip, to be better appreciated. Also, the volumetric model can be used as an "anchor" to display a zone of density around the chosen anchoring ball-and-stick model, which allows the detection and visualization of crosslinker densities between adjacent actin filaments. The ability to select individual filaments and to set appropriate radii as a zone in which density is displayed again allows the overwhelming complexity of the scenery to be reduced to a manageable level.

One advantage of this volumetric model building approach of global fitting, followed by local adjustments, was that we were able to identify regions where actin filaments appeared to be interrupted, and significant gaps in actin filaments were indicated by the absence of map density. Since we had placed a volumetric ball-and-stick model, we could make use of a routine in the UCSF Chimera software package that color-codes each model ball position according to the map density value at that location. This approach allowed a fast detection and visualization of actin filament gaps in the actin bundle, which is a biologically significant feature we found in our cryo-electron tomogram, and which would have been very difficult to detect and visualize with traditional segmentation approaches. Yet another advantage of our volumetric model is that volumetric properties including lengths and distances can be easily obtained, which allows for actual numbers and thus a statistical analysis to be performed.

In summary, interactive manual model point placement, possibly further augmented by subsequent automated local fitting and filament tracing capabilities, is a rather promising approach for the visualization and quantitative analysis of electron tomographic subcellular volumes. This is because it uses the power of the human brain for pattern recognition and the power of computer science for model optimization.

Materials

Name	Company	Catalog Number	Comments
Chimera	RBVI		Version 1.16 https://www.cgl.ucsf.edu/chimera/download.html
Chimera	RBVI		Version 1.16 https://www.cgl.ucsf.edu/chimera/cgi-bin/secure/chimera-get.py?file=win64/chimera-1.16-win64.exe
Chimera	RBVI		Version 1.16 https://www.cgl.ucsf.edu/chimera/cgi-bin/secure/chimera-get.py?file=mac64/chimera-1.16-mac64.dmg
Chimera	RBVI		Version 1.16 https://www.cgl.ucsf.edu/chimera/cgi-bin/secure/chimera-get.py?file=linux_x86_64/chimera-1.16-linux_x86_64.bin
Excel	Microsoft		Version 2211 https://www.office.com/?auth=1
Falcon II	Thermofisher		https://www.thermofisher.com/de/de/home/electron-microscopy/products/accessories-em/falcon-detector.html
IMOD	University of Colorado		Version 4.11.1 https://bio3d.colorado.edu/imod/download.html
PC Desktop	Intel		Windows 10, ver. 22H2
PC Laptop	Gigabyte		Windows 10, ver. 22H2
Powerpoint	Microsoft		Version 2211 https://www.office.com/?auth=1
Titan Krios Electron Microscope	Thermofisher		https://www.thermofisher.com/de/de/home/electron-microscopy/products/transmission-electron-microscopes/krios-g4-cryo-tem.html
Word	Microsoft		Version 2211 https://www.office.com/?auth=1