Using Sculptor and Situs for simultaneous assembly of atomic components into low-resolution shapes

We describe an integrated software system called Sculptor that combines visualization capabilities with molecular modeling algorithms for the analysis of multi-scale data sets. Sculptor features extensive special purpose visualization techniques that are based on modern GPU programming and are capable of representing complex molecular assemblies in real-time. The integration of graphics and modeling offers several advantages. The user interface not only eases the usually steep learning curve of pure algorithmic techniques, but it also permits instant analysis and post-processing of results, as well as the integration of results from external software. Here, we implemented an interactive peak-selection strategy that enables the user to explore a preliminary score landscape generated by the colors tool of Situs. The interactive placement of components, one at a time, is advantageous for low-resolution or ambiguously shaped maps, which are sometimes difficult to interpret by the fully automatic peak selection of colors. For the subsequent refinement of the preliminary models resulting from both interactive and automatic peak selection, we have implemented a novel simultaneous multi-body docking in Sculptor and Situs that softly enforces shape complementarities between components using the normalization of the cross-correlation coefficient. The proposed techniques are freely available in Situs version 2.6 and Sculptor version 2.0.     

Evolutionary bidirectional expansion for the tracing of alpha helices in cryo-electron microscopy reconstructions

Cryo-electron microscopy (cryo-EM) enables the imaging of macromolecular complexes in near-native environments at resolutions that often permit the visualization of secondary structure elements. For example, alpha helices frequently show consistent patterns in volumetric maps, exhibiting rod-like structures of high density. Here, we introduce VolTrac (Volume Tracer) – a novel technique for the annotation of alpha-helical density in cryo-EM data sets. VolTrac combines a genetic algorithm and a bidirectional expansion with a tabu search strategy to trace helical regions. Our method takes advantage of the stochastic search by using a genetic algorithm to identify optimal placements for a short cylindrical template, avoiding exploration of already characterized tabu regions. These placements are then utilized as starting positions for the adaptive bidirectional expansion that characterizes the curvature and length of the helical region. The method reliably predicted helices with seven or more residues in experimental and simulated maps at intermediate (4–10 Å) resolution. The observed success rates, ranging from 70.6% to 100%, depended on the map resolution and validation parameters. For successful predictions, the helical axes were located within 2 Å from known helical axes of atomic structures.     

Tools for correlative cryo-fluorescence microscopy and cryo-electron tomography applied to whole mitochondria in human endothelial cells

Cryo-electron tomography (cryo-ET) allows for the visualization of biological material in a close-to-native state, in three dimensions and with nanometer scale resolution. However, due to the low signal-to-noise ratio inherent to imaging of the radiation-sensitive frozen-hydrated samples, it appears oftentimes impossible to localize structures within heterogeneous samples. Because a major potential for cryo-ET is thereby left unused, we set out to combine cryo-ET with cryo-fluorescence microscopy (cryo-FM), in order to facilitate the search for structures of interest. We describe a cryo-FM setup and workflow for correlative cryo-fluorescence and cryo-electron microscopy (cryo-CLEM) that can be easily implemented. Cells are grown on finder grids, vitally labeled with one or two fluorescent dyes, and vitrified. After a structure is located by cryo-FM (with 0.4 μm resolution), its image coordinates are translated to cryo-ET stage coordinates via a home-built software routine. We tested our workflow on whole mount primary human umbilical vein endothelial cells. The correlative routine enabled us to investigate mitochondrial ultrastructure for the first time on intact human mitochondria, and led us to find mitochondrial cristae that were connected to the intermembrane space via large slits, which challenges the current view that such connections are established exclusively via small circular pores. Taken together, this study emphasizes that cryo-CLEM can be a routinely used technique that opens up exciting new possibilities for cryo-ET.     

Parakeet: a digital twin software pipeline to assess the impact of experimental parameters on tomographic reconstructions for cryo-electron tomography

In cryo-electron tomography (cryo-ET) of biological samples, the quality of tomographic reconstructions can vary depending on the transmission electron microscope (TEM) instrument and data acquisition parameters. In this paper, we present Parakeet, a ‘digital twin’ software pipeline for the assessment of the impact of various TEM experiment parameters on the quality of three-dimensional tomographic reconstructions. The Parakeet digital twin is a digital model that can be used to optimize the performance and utilization of a physical instrument to enable in silico optimization of sample geometries, data acquisition schemes and instrument parameters. The digital twin performs virtual sample generation, TEM image simulation, and tilt series reconstruction and analysis within a convenient software framework. As well as being able to produce physically realistic simulated cryo-ET datasets to aid the development of tomographic reconstruction and subtomogram averaging programs, Parakeet aims to enable convenient assessment of the effects of different microscope parameters and data acquisition parameters on reconstruction quality. To illustrate the use of the software, we present the example of a quantitative analysis of missing wedge artefacts on simulated planar and cylindrical biological samples and discuss how data collection parameters can be modified for cylindrical samples where a full 180° tilt range might be measured.     

ContinuousFlex: Software package for analyzing continuous conformational variability of macromolecules in cryo electron microscopy and tomography data

ContinuousFlex is a user-friendly open-source software package for analyzing continuous conformational variability of macromolecules in cryo electron microscopy (cryo-EM) and cryo electron tomography (cryo-ET) data. In 2019, ContinuousFlex became available as a plugin for Scipion, an image processing software package extensively used in the cryo-EM field. Currently, ContinuousFlex contains software for running (1) recently published methods HEMNMA-3D, TomoFlow, and NMMD; (2) earlier published methods HEMNMA and StructMap; and (3) methods for simulating cryo-EM and cryo-ET data with conformational variability and methods for data preprocessing. It also includes external software for molecular dynamics simulation (GENESIS) and normal mode analysis (ElNemo), used in some of the mentioned methods. The HEMNMA software has been presented in the past, but not the software of other methods. Besides, ContinuousFlex currently also offers a deep learning extension of HEMNMA, named DeepHEMNMA. In this article, we review these methods in the context of the ContinuousFlex package, developed to facilitate their use by the community.     

Cryo-electron tomography of cells: connecting structure and function

Cryo-electron tomography (cryo-ET) allows the visualization of cellular structures under close-to-life conditions and at molecular resolution. While it is inherently a static approach, yielding structural information about supramolecular organization at a certain time point, it can nevertheless provide insights into function of the structures imaged, in particular, when supplemented by other approaches. Here, we review the use of experimental methods that supplement cryo-ET imaging of whole cells. These include genetic and pharmacological manipulations, as well as correlative light microscopy and cryo-ET. While these methods have mostly been used to detect and identify structures visualized in cryo-ET or to assist the search for a feature of interest, we expect that in the future they will play a more important role in the functional interpretation of cryo-tomograms.     

Locating macromolecules and determining structures inside bacterial cells using electron cryotomography

Electron cryotomography (cryo-ET) is an imaging technique uniquely suited to the study of bacterial ultrastructure and cell biology. Recent years have seen a surge in structural and cell biology research on bacteria using cryo-ET. This research has driven major technical developments in the field, with applications emerging to address a wide range of biological questions. In this review, we explore the diversity of cryo-ET approaches used for structural and cellular microbiology, with a focus on in situ localization and structure determination of macromolecules. The first section describes strategies employed to locate target macromolecules within large cellular volumes. Next, we explore methods to study thick specimens by sample thinning. Finally, we review examples of macromolecular structure determination in a cellular context using cryo-ET. The examples outlined serve as powerful demonstrations of how the cellular location, structure, and function of any bacterial macromolecule of interest can be investigated using cryo-ET.     

An integrative simulation model linking major biochemical reactions of actin-polymerization to structural properties of actin filaments

We report on an advanced universal Monte Carlo simulation model of actin polymerization processes offering a broad application panel. The model integrates major actin-related reactions, such as assembly of actin nuclei, association/dissociation of monomers to filament ends, ATP-hydrolysis via ADP-Pi formation and ADP-ATP exchange, filament branching, fragmentation and annealing or the effects of regulatory proteins. Importantly, these reactions are linked to information on the nucleotide state of actin subunits in filaments (ATP hydrolysis) and the distribution of actin filament lengths. The developed stochastic simulation modelling schemes were validated on: i) synthetic theoretical data generated by a deterministic model and ii) sets of our and published experimental data obtained from fluorescence pyrene-actin experiments. Build on an open-architecture principle, the designed model can be extended for predictive evaluation of the activities of other actin-interacting proteins and can be applied for the analysis of experimental pyrene actin-based or fluorescence microscopy data. We provide a user-friendly, free software package ActinSimChem that integrates the implemented simulation algorithms and that is made available to the scientific community for modelling in silico any specific actin-polymerization system.     

Revealing bacterial cell biology using cryo-electron tomography

Visualizing macromolecules inside bacteria at a high spatial resolution has remained a challenge owing to their small size and limited resolution of optical microscopy techniques. Recent advances in cryo-electron tomography (cryo-ET) imaging methods have revealed the spatial and temporal assemblies of many macromolecules involved in different cellular processes in bacteria at a resolution of a few nanometers in their native milieu. Specifically, the application of cryo-focused ion beam (cryo-FIB) milling to thin bacterial specimens makes them amenable for high-resolution cryo-ET data collection. In this review, we highlight recent research in three emerging areas of bacterial cell biology that have benefited from the cryo-FIB-ET technology - cytoskeletal filament assembly, intracellular organelles, and multicellularity.     

Multiscale imaging of neurons grown in culture: from light microscopy to cryo-electron tomography

Cryo-electron tomography (cryo-ET) allows the visualization of supramolecular architecture in cells preserved in a close-to-physiological state. In order to supplement the structural information obtained by cryo-ET with the functional state of the molecules involved based on fluorescent labeling we developed a method of correlating light microscopy and cryo-ET. This method is suitable for investigating complicated cellular landscapes such as mature neurons grown in culture. It has the advantage that a correlation is obtained without exposing a feature of interest to additional electron irradiation, and that it does not rely on visual recognition of features. Different modes of correlation are presented here: a feature identified on a light microscopy image is used to guide the cryo-ET investigation, and cryo-tomograms are correlated to light microscopy images. Cryo-tomograms of a neuronal synapse and of an isolated presynaptic terminal are shown as examples of the correlative method. The correlation method presented here can be expected to provide new insights into the structure-function relationship of supramolecular organization in neurons.     

Regulation of the actin-activated ATPase and in vitro motility activities of monomeric and filamentous Acanthamoeba myosin II.

The actin-activated Mg(2+)-ATPase activity of filamentous Acanthamoeba myosin II is inhibited by phosphorylation of 3 serine residues at the tip of the tail of each heavy chain. From previous studies, it had been concluded that the activity of each molecule in the filament was regulated by the global state of phosphorylation of the filament and was independent of its own phosphorylation state. The actin-activated Mg(2+)-ATPase activity of monomeric phosphorylated myosin II was not known because it polymerizes under the ionic conditions necessary for the expression of this activity. We have now found conditions to maintain myosin II monomeric and active during the enzyme assay. The actin-activated Mg(2+)-ATPase activities of monomeric dephosphorylated and phosphorylated myosin II were found to be the same as the activity of filamentous dephosphorylated myosin II. These results support the conclusion that phosphorylation regulates filamentous myosin II by affecting filament conformation. Consistent with their equivalent enzymatic activities, monomeric and filamentous dephosphorylated myosin II were equally active in an in vitro motility assay in which myosin adsorbed to a surface drives the movement of F-actin. In contrast to their very different enzymatic activities, however, filamentous and monomeric phosphorylated myosin II had similar activities in the in vitro motility assay; both were much less active than monomeric and filamentous dephosphorylated myosin II. One interpretation of these results is that the rate-limiting steps in the two assays are different and that, while the rate-limiting step for actin-activated Mg(2+)-ATPase activity is regulated only at the level of the filament, the rate-limiting step for motility can also be regulated at the level of the monomer.     

Intermediate filament systems are collapsed onto the nuclear surface after isolation of nuclei from tissue culture cells

Standard procedures for the biochemical isolation and purification of tissue culture cell nuclei are very similar to recently described methods for the preparation of detergent-resistant cytoskeletons. But the latter in addition to extracted nuclei contain cytoplasmic filament systems. We, therefore, investigated the distribution of these filamentous systems during nuclear isolation both by immunofluorescence microscopy and by SDS-gel electrophoresis. The intermediate filaments copurify with each single isolated nucleus. The majority of the filaments is collapsed onto isolated nuclei and still constitutes a filamentous system. This can be shown by partially unfolding the collapsed filament systems.     

i-ADHoRe 3.0--fast and sensitive detection of genomic homology in extremely large data sets

Comparative genomics is a powerful means to gain insight into the evolutionary processes that shape the genomes of related species. As the number of sequenced genomes increases, the development of software to perform accurate cross-species analyses becomes indispensable. However, many implementations that have the ability to compare multiple genomes exhibit unfavorable computational and memory requirements, limiting the number of genomes that can be analyzed in one run. Here, we present a software package to unveil genomic homology based on the identification of conservation of gene content and gene order (collinearity), i-ADHoRe 3.0, and its application to eukaryotic genomes. The use of efficient algorithms and support for parallel computing enable the analysis of large-scale data sets. Unlike other tools, i-ADHoRe can process the Ensembl data set, containing 49 species, in 1?h. Furthermore, the profile search is more sensitive to detect degenerate genomic homology than chaining pairwise collinearity information based on transitive homology. From ultra-conserved collinear regions between mammals and birds, by integrating coexpression information and protein-protein interactions, we identified more than 400 regions in the human genome showing significant functional coherence. The different algorithmical improvements ensure that i-ADHoRe 3.0 will remain a powerful tool to study genome evolution.     

Discussion on common data analysis strategies used in MS-based proteomics

Current proteomics technology is limited in resolving the proteome complexity of biological systems. The main issue at stake is to increase throughput and spectra quality so that spatiotemporal dimensions, population parameters and the complexity of protein modifications on a quantitative scale can be considered. MS-based proteomics and protein arrays are the main players in large-scale proteome analysis and an integration of these two methodologies is powerful but presently not sufficient for detailed quantitative and spatiotemporal proteome characterization. Improvements of instrumentation for MS-based proteomics have been achieved recently resulting in data sets of approximately one million spectra which is a large step in the right direction. The corresponding raw data range from 50 to 100?Gb and are frequently made available. Multidimensional LC-MS data sets have been demonstrated to identify and quantitate 2000-8000 proteins from whole cell extracts. The analysis of the resulting data sets requires several steps from raw data processing, to database-dependent search, statistical evaluation of the search result, quantitative algorithms and statistical analysis of quantitative data. A large number of software tools have been proposed for the above-mentioned tasks. However, it is not the aim of this review to cover all software tools, but rather discuss common data analysis strategies used by various algorithms for each of the above-mentioned steps in a non-redundant approach and to argue that there are still some areas which need improvements.     

Self-assembly incompetence of synemin is related to the property of its head and rod domains

The mechanisms regulating the intermediate filament (IF) protein assembly are complex and not yet fully understood. All vertebrate cytoplasmic IF proteins have a central alpha-helical rod domain flanked by variable head and tail domains. The IF protein synemin cannot homopolymerize to form filament networks; it needs an appropriate copolymerization partner. To elucidate the roles of the vimentin head domain, the TAAL motif in the 2A region, and the TYRKLLEGEE motif in the 2B region of the rod domain in synemin filament formation, we have prepared a series of synemin constructs by site-directed mutagenesis and chimeric synemins having the vimentin head domain. The assembly properties of synemin constructs were assessed by the immunofluorescence of transient transfection into cultured SW13 cells without endogenous IFs. Our data showed that the formation of a filamentous network required at least the vimentin-like head domain and both the 2A and 2B regions of the rod domain.     

Single particle analysis of filamentous and highly elongated macromolecular assemblies

The application of single particle techniques to the three-dimensional analysis of electron microscope images of elongated or filamentous macromolecular assemblies is evaluated, taking as an example the muscle thin filament. Although the thin filament contains local helical symmetry, because of the inherent variable twist along it, the helical coherence does not extend for large enough distances to allow the symmetry to be used for full reconstruction of the tropomyosin/troponin repeat along the filament. The muscle thin filament therefore represents a general case of a filamentous object in that it is not possible to exploit symmetry in a full analysis. Due to the nature of the imaging process in the electron microscope, only projections of the thin filament around its long axis are available without tilting the grid. Crucially, projection images around a single axis do not provide enough information to assign Euler angles ab initio using current methods. Tests with a model thin filament structure indicated that an out-of-plane tilt of approximately 20 degrees was needed for ab initio angular assignment of sufficient accuracy to calculate a 3D structure to a resolution of approximately 25 A. If no out-of-plane views are available, an alternative approach is to use a prior 3D model as a reference for the initial angle assignment. Tests with the thin filament model indicated that reasonably accurate angular assignment can be made using a reference containing actin, but lacking the regulatory proteins tropomyosin and troponin. We also found that an adaptation of the exact filtered back projection method is required to allow the correct weighting of projection images in which the particle has a very large axial ratio. This adaptation resulted in significant improvements in the reconstruction.     

Ca2+-calmodulin-dependent polymerization of actin by myelin basic protein

The interaction between myelin basic protein (MBP) and G-actin was studied under nonpolymerizing conditions, i.e.,2mM HEPES, pH 7.5, 0.1 mM CaCl2 and 0.2 mM ATP. Fluorescence studies using pyrenyl-actin and the measurements of ATP hydrolysis rate show that MBP induces changes in the structure of the actin monomer similar to those occurring during polymerization by salt. Electron microscope observations of the MBP-G-actin complex reveal the presence of filamentous structures which appear as separate filaments or as bundles of filaments in lateral association. These filaments are polar as visualized by attachment of heavy meromyosin. The biochemical data together with electron microscope observations suggest that the binding of MBP to G-actin under non-polymerizing conditions induces an interaction between actin monomers leading to the formation of filamentous structures which may be similar to F-actin filaments. The effects of MBP on G-actin can be reversed by calmodulin in the presence of Ca2+.     

Comparative Map and Trait Viewer (CMTV): an integrated bioinformatic tool to construct consensus maps and compare QTL and functional genomics data across genomes and experiments

In the past few decades, a wealth of genomic data has been produced in a wide variety of species using a diverse array of functional and molecular marker approaches. In order to unlock the full potential of the information contained in these independent experiments, researchers need efficient and intuitive means to identify common genomic regions and genes involved in the expression of target phenotypic traits across diverse conditions. To address this need, we have developed a Comparative Map and Trait Viewer (CMTV) tool that can be used to construct dynamic aggregations of a variety of types of genomic datasets. By algorithmically determining correspondences between sets of objects on multiple genomic maps, the CMTV can display syntenic regions across taxa, combine maps from separate experiments into a consensus map, or project data from different maps into a common coordinate framework using dynamic coordinate translations between source and target maps. We present a case study that illustrates the utility of the tool for managing large and varied datasets by integrating data collected by CIMMYT in maize drought tolerance research with data from public sources. This example will focus on one of the visualization features for Quantitative Trait Locus (QTL) data, using likelihood ratio (LR) files produced by generic QTL analysis software and displaying the data in a unique visual manner across different combinations of traits, environments and crosses. Once a genomic region of interest has been identified, the CMTV can search and display additional QTLs meeting a particular threshold for that region, or other functional data such as sets of differentially expressed genes located in the region; it thus provides an easily used means for organizing and manipulating data sets that have been dynamically integrated under the focus of the researchers specific hypothesis.