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

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

Radiation dose reduction and image enhancement in biological imaging through equally-sloped tomography

Electron tomography is currently the highest resolution imaging modality available to study the 3D structures of pleomorphic macromolecular assemblies, viruses, organelles and cells. Unfortunately, the resolution is currently limited to 3–5 nm by several factors including the dose tolerance of biological specimens and the inaccessibility of certain tilt angles. Here we report the first experimental demonstration of equally-sloped tomography (EST) to alleviate these problems. As a proof of principle, we applied EST to reconstructing frozen-hydrated keyhole limpet hemocyanin molecules from a tilt-series taken with constant slope increments. In comparison with weighted back-projection (WBP), the algebraic reconstruction technique (ART) and the simultaneous algebraic reconstruction technique (SART), EST reconstructions exhibited higher contrast, less peripheral noise, more easily detectable molecular boundaries and reduced missing wedge effects. More importantly, EST reconstructions including only two-thirds the original images appeared to have the same resolution as full WBP reconstructions, suggesting that EST can either reduce the dose required to reach a given resolution or allow higher resolutions to be achieved with a given dose. EST was also applied to reconstructing a frozen-hydrated bacterial cell from a tilt-series taken with constant angular increments. The results confirmed similar benefits when standard tilts are utilized

A Bayesian approach to single-particle electron cryo-tomography in RELION-4.0

We present a new approach for macromolecular structure determination from multiple particles in electron cryo-tomography (cryo-ET) data sets. Whereas existing subtomogram averaging approaches are based on 3D data models, we propose to optimise a regularised likelihood target that approximates a function of the 2D experimental images. In addition, analogous to Bayesian polishing and contrast transfer function (CTF) refinement in single-particle analysis, we describe the approaches that exploit the increased signal-to-noise ratio in the averaged structure to optimise tilt-series alignments, beam-induced motions of the particles throughout the tilt-series acquisition, defoci of the individual particles, as well as higher-order optical aberrations of the microscope. Implementation of our approaches in the open-source software package RELION aims to facilitate their general use, particularly for those researchers who are already familiar with its single-particle analysis tools. We illustrate for three applications that our approaches allow structure determination from cryo-ET data to resolutions sufficient for de novo atomic modelling.This work was funded by the UK Research and Innovation (UKRI) Medical Research Council (MC_UP_A025_1013 to SHWS; and MC_UP_1201/16 to JAGB), the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (ERC-CoG-2014, grant 648432, MEMBRANEFUSION to JAGB and ERC StG-2019, grant 852915 CRYTOCOP to GZ); the Swiss National Science Foundation (grant 205321_179041/1 to DC-D), the Max Planck Society (to JAGB) and the UKRI Biotechnology and Biological Sciences Research Council (grant BB/T002670/1 to GZ). TAMB is a recipient of a Sir Henry Dale Fellowship, jointly funded by the Wellcome Trust and the Royal Society (202231/Z/16/Z). JZ was partially funded by the European Union’s Horizon 2020 research and innovation program (ERC-ADG-2015, grant 692726, GlobalBioIm to Michael Unser)

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

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

Tools for Subtomogram Averaging for Cryo Electron Tomography

Trabajo presentado en el 44º Congreso Nacional de la Sociedad Española de Bioquímica y Biología Molecular (SEBBM), celebrado en Málaga (España) del 06 al 09 de septiembre de 2022

An adaptive wavelet viscosity method for hyperbolic conservation laws

We extend the multiscale finite element viscosity method for hyperbolic conservation laws developed in terms of hierarchical finite element bases to a (pre-orthogonal spline-)wavelet basis. Depending on an appropriate error criterion, the multiscale framework allows for a controlled adaptive resolution of discontinuities of the solution. The nonlinearity in the weak form is treated by solving a least-squares data fitting problem