CTF determination and correction for low dose tomographic tilt series

The resolution of cryo-electron tomography can be limited by the first zero of the microscope’s contrast transfer function (CTF). To achieve higher resolution, it is critical to determine the CTF and correct its phase inversions. However, the extremely low signal-to-noise ratio (SNR) and the defocus gradient in the projections of tilted specimens make this process challenging. Two programs, CTFPLOTTER and CTFPHASEFLIP, have been developed to address these issues. CTFPLOTTER obtains a 1D power spectrum by periodogram averaging and rotational averaging and it estimates the noise background with a novel approach, which uses images taken with no specimen. The background-subtracted 1D power spectra from image regions at different defocus values are then shifted to align their first zeros and averaged together. This averaging improves the SNR sufficiently that it becomes possible to determine the defocus for subsets of the tilt series rather than just the entire series. CTFPHASEFLIP corrects images line-by-line by inverting phases appropriately in thin strips of the image at nearly constant defocus. CTF correction by these methods is shown to improve the resolution of aligned, averaged particles extracted from tomograms. However, some restoration of Fourier amplitudes at high frequencies is important for seeing the benefits from CTF correction.     

Automatic CTF correction for single particles based upon multivariate statistical analysis of individual power spectra

Three-dimensional electron cryomicroscopy of randomly oriented single particles is a method that is suitable for the determination of three-dimensional structures of macromolecular complexes at molecular resolution. However, the electron-microscopical projection images are modulated by a contrast transfer function (CTF) that prevents the calculation of three-dimensional reconstructions of biological complexes at high resolution from uncorrected images. We describe here an automated method for the accurate determination and correction of the CTF parameters defocus, twofold astigmatism and amplitude-contrast proportion from single-particle images. At the same time, the method allows the frequency-dependent signal decrease (B factor) and the non-convoluted background signal to be estimated. The method involves the classification of the power spectra of single-particle images into groups with similar CTF parameters; this is done by multivariate statistical analysis (MSA) and hierarchically ascending classification (HAC). Averaging over several power spectra generates class averages with enhanced signal-to-noise ratios. The correct CTF parameters can be deduced from these class averages by applying an iterative correlation procedure with theoretical CTF functions; they are then used to correct the raw images. Furthermore, the method enables the tilt axis of the sample holder to be determined and allows the elimination of individual poor-quality images that show high drift or charging effects.     

IceBreaker: Software for high-resolution single-particle cryo-EM with non-uniform ice

1. Download : Download high-res image (263KB)
2. Download : Download full-size image

     

ASOCEM: Automatic Segmentation Of Contaminations in cryo-EM

Particle picking is currently a critical step in the cryo-electron microscopy single particle reconstruction pipeline. Contaminations in the acquired micrographs severely degrade the performance of particle pickers, resulting in many “non-particles” in the collected stack of particles. In this paper, we present ASOCEM (Automatic Segmentation Of Contaminations in cryo-EM), an automatic method to detect and segment contaminations, which requires as an input only the approximate particle size. In particular, it does not require any parameter tuning nor manual intervention. Our method is based on the observation that the statistical distribution of contaminated regions is different from that of the rest of the micrograph. This nonrestrictive assumption allows to automatically detect various types of contaminations, from the carbon edges of the supporting grid to high contrast blobs of different sizes. We demonstrate the efficiency of our algorithm using various experimental data sets containing various types of contaminations. ASOCEM is integrated as part of the KLT picker (Eldar et al., 2020) and is available at https://github.com/ShkolniskyLab/kltpicker2.     

Accurate determination of local defocus and specimen tilt in electron microscopy

Accurate knowledge of defocus and tilt parameters is essential for the determination of three-dimensional protein structures at high resolution using electron microscopy. We present two computer programs, CTFFIND3 and CTFTILT, which determine defocus parameters from images of untilted specimens, as well as defocus and tilt parameters from images of tilted specimens, respectively. Both programs use a simple algorithm that fits the amplitude modulations visible in a power spectrum with a calculated contrast transfer function (CTF). The background present in the power spectrum is calculated using a low-pass filter. The background is then subtracted from the original power spectrum, allowing the fitting of only the oscillatory component of the CTF. CTFTILT determines specimen tilt parameters by measuring the defocus at a series of locations on the image while constraining them to a single plane. We tested the algorithm on images of two-dimensional crystals by comparing the results with those obtained using crystallographic methods. The images also contained contrast from carbon support film that added to the visibility of the CTF oscillations. The tests suggest that the fitting procedure is able to determine the image defocus with an error of about 10nm, whereas tilt axis and tilt angle are determined with an error of about 2 degrees and 1 degrees, respectively. Further tests were performed on images of single protein particles embedded in ice that were recorded from untilted or slightly tilted specimens. The visibility of the CTF oscillations from these images was reduced due to the lack of a carbon support film. Nevertheless, the test results suggest that the fitting procedure is able to determine image defocus and tilt angle with errors of about 100 nm and 6 degrees, respectively.     

TomoAlign: A novel approach to correcting sample motion and 3D CTF in CryoET

TomoAlign is a software package that integrates tools to mitigate two important resolution limiting factors in cryoET, namely the beam-induced sample motion and the contrast transfer function (CTF) of the microscope. The package is especially focused on cryoET of thick specimens where fiducial markers are required for accurate tilt-series alignment and sample motion estimation. TomoAlign models the beam-induced sample motion undergone during the tilt-series acquisition. The motion models are used to produce motion-corrected subtilt-series centered on the particles of interest. In addition, the defocus of each particle at each tilt image is determined and can be corrected, resulting in motion-corrected and CTF-corrected subtilt-series from which the subtomograms can be computed. Alternatively, the CTF information can be passed on so that CTF correction can be carried out entirely within external packages like Relion. TomoAlign serves as a versatile tool that can streamline the cryoET workflow from initial alignment of tilt-series to final subtomogram averaging during in situ structure determination.     

A new cryo-EM single-particle ab initio reconstruction method visualizes secondary structure elements in an ATP-fueled AAA+ motor

The generation of ab initio three-dimensional (3D) models is a bottleneck in the studies of large macromolecular assemblies by single-particle cryo-electron microscopy. We describe here a novel method, in which established methods for two-dimensional image processing are combined with newly developed programs for joint rotational 3D alignment of a large number of class averages (RAD) and calculation of 3D volumes from aligned projections (VolRec). We demonstrate the power of the method by reconstructing an ∼ 660-kDa ATP-fueled AAA+ motor to 7.5 Å resolution, with secondary structure elements identified throughout the structure. We propose the method as a generally applicable automated strategy to obtain 3D reconstructions from unstained single particles imaged in vitreous ice.     

A relational database for cryoEM: experience at one year and 50 000 images

For the past year we have been using a relational database as part of an automated data collection system for cryoEM. The database is vital for keeping track of the very large number of images collected and analyzed by the automated system and essential for quantitatively evaluating the utility of methods and algorithms used in the data collection. The database can be accessed using a variety of tools including specially developed Web-based interfaces that enable a user to annotate and categorize images using a Web-based form.     

Helical inchworming: a novel translocation mechanism for a ring ATPase

Ring ATPases perform a variety of tasks in the cell. Their function involves complex communication and coordination among the often identical subunits. Translocases in this group are of particular interest as they involve both chemical and mechanical actions in their operation. We study the DNA packaging motor of bacteriophage φ29, and using single-molecule optical tweezers and single-particle cryo-electron microscopy, have discovered a novel translocation mechanism for a molecular motor.     

High-vacuum optical platform for cryo-CLEM (HOPE): A new solution for non-integrated multiscale correlative light and electron microscopy

Cryo-correlative light and electron microscopy (cryo-CLEM) offers a unique way to analyze the high-resolution structural information of cryo-vitrified specimen by cryo-electron microscopy (cryo-EM) with the guide of the search for unique events by cryo-fluorescence microscopy (cryo-FM). To achieve cryo-FM, a trade-off must be made between the temperature and performance of objective lens. The temperature of specimen should be kept below devitrification while the distance between the objective lens and specimen should be short enough for high resolution imaging. Although special objective lens was designed in many current cryo-FM approaches, the unavoided frosting and ice contamination are still affecting the efficiency of cryo-CLEM. In addition, the correlation accuracy between cryo-FM and cryo-EM would be reduced during the current specimen transfer procedure. Here, we report an improved cryo-CLEM technique (high-vacuum optical platform for cryo-CLEM, HOPE) based on a high-vacuum optical stage and a commercial cryo-EM holder. The HOPE stage comprises of a special adapter to suit the cryo-EM holder and a high-vacuum chamber with an anti-contamination system. It provides a clean and enduring environment for cryo specimen, while the normal dry objective lens in room temperature can be used via the optical windows. The ‘touch-free’ specimen transfer via cryo-EM holder allows least specimen deformation and thus maximizes the correlation accuracy between cryo-FM and cryo-EM. Besides, we developed a software to perform semi-automatic cryo-EM acquisition of the target region localized by cryo-FM. Our work provides a new solution for cryo-CLEM and can be adapted for different commercial fluorescence microscope and electron microscope.     

Development of a new protein labeling system to map subunits and domains of macromolecular complexes for electron microscopy

Several gene fusion technologies have been successfully applied to label particular subunits or domains within macromolecular complexes to enable positional mapping of electron microscopy (EM) density maps, but exogenous fusion of a protein domain into the target polypeptide can cause unwanted structural and functional outcomes. Fab fragments from antibodies can be used as labeling reagents during EM visualization without gene manipulation of the target protein, but this method requires a panel of high-affinity antibodies that recognize a wide variety of epitopes. Linear peptide tags and their anti-tag antibodies can be used but they have a limited mapping ability as their placement is usually limited to the terminal regions of a protein. The PA dodecapeptide epitope tag (GVAMPGAEDDVV), forms a tight β-turn in the antigen binding pocket of its antibody (NZ-1). This capability allows for insertion of the PA tag into various surface-exposed loops within a multi-domain cell adhesion receptor, α_IIbβ₃ integrin. We confirmed that the purified PA-tagged integrin ectodomain fragments can form a stable complex with NZ-1 Fab. Negative stain EM of the various integrin-NZ-1 complexes revealed that a majority of the particles exhibited a clear density corresponding to the NZ-1 Fab; and the positions of the bound Fab were in good agreement with the predicted location of the inserted PA tag. The high-affinity and insertion-compatibility of the PA tag system allowed us to develop a new EM labeling methodology applicable to proteins for which good antibodies are not available.     

Two‐photon focal modulation microscopy for high‐resolution imaging in deep tissue

Two‐photon microscopy (2PM) is one of the most widely used tools for in vivo deep tissue imaging. However, the spatial resolution and penetration depth are still limited due to the strong scattering background. Here we demonstrate a two‐photon focal modulation microscopy. By utilizing the modulation and demodulation techniques, background rejection capability is enhanced, thus spatial resolution and imaging penetration depth are improved. Compared with 2PM, the transverse resolution is increased by 70%, while the axial resolution is increased to 2‐fold. Furthermore, when applied in conventional 2PM mode, it can achieve inertial‐free scanning in either transverse or axial direction with in principle unlimited scanning speed. Finally, we applied 2PFMM in thick scattering samples to further examine the imaging performance. The results show that the signal‐to‐background ratio of 2PFMM can be improved up to five times of 2PM at the depth of 500 μm. Fluorescent imaging in the mouse brain tissue. 3D Thy1‐GFP hippocampal neurons imaged by (A) 2PM compared with (B) 2PFMM; (C‐H) xy maximum‐intensity projection imaged by 2PM compared with 2PFMM. Scale bar 50 μm.      

Electron microscopy holdings of the Protein Data Bank: the impact of the resolution revolution,new validation tools,and implications for the future

As a discipline, structural biology has been transformed by the three-dimensional electron microscopy (3DEM) “Resolution Revolution” made possible by convergence of robust cryo-preservation of vitrified biological materials, sample handling systems, and measurement stages operating a liquid nitrogen temperature, improvements in electron optics that preserve phase information at the atomic level, direct electron detectors (DEDs), high-speed computing with graphics processing units, and rapid advances in data acquisition and processing software. 3DEM structure information (atomic coordinates and related metadata) are archived in the open-access Protein Data Bank (PDB), which currently holds more than 11,000 3DEM structures of proteins and nucleic acids, and their complexes with one another and small-molecule ligands (~ 6% of the archive). Underlying experimental data (3DEM density maps and related metadata) are stored in the Electron Microscopy Data Bank (EMDB), which currently holds more than 21,000 3DEM density maps. After describing the history of the PDB and the Worldwide Protein Data Bank (wwPDB) partnership, which jointly manages both the PDB and EMDB archives, this review examines the origins of the resolution revolution and analyzes its impact on structural biology viewed through the lens of PDB holdings. Six areas of focus exemplifying the impact of 3DEM across the biosciences are discussed in detail (icosahedral viruses, ribosomes, integral membrane proteins, SARS-CoV-2 spike proteins, cryogenic electron tomography, and integrative structure determination combining 3DEM with complementary biophysical measurement techniques), followed by a review of 3DEM structure validation by the wwPDB that underscores the importance of community engagement.     

Weighted averaging partial least squares regression (WA-PLS): an improved method for reconstructing environmental variables from species assemblages

Weighted averaging regression and calibration form a simple, yet powerful method for reconstructing environmental variables from species assemblages. Based on the concepts of niche-space partitioning and ecological optima of species (indicator values), it performs well with noisy, species-rich data that cover a long ecological gradient (>3 SD units). Partial least squares regression is a linear method for multivariate calibration that is popular in chemometrics as a robust alternative to principal component regression. It successively selects linear components so as to maximize predictive power. In this paper the ideas of the two methods are combined. It is shown that the weighted averaging method is a form of partial least squares regression applied to transformed data that uses the first PLS-component only. The new combined method, ast squares, consists of using further components, namely as many as are useful in terms of predictive power. The further components utilize the residual structure in the species data to improve the species parameters (optima) in the final weighted averaging predictor. Simulations show that the new method can give 70% reduction in prediction error in data sets with low noise, but only a small reduction in noisy data sets. In three real data sets of diatom assemblages collected for the reconstruction of acidity and salinity, the reduction in prediction error was zero, 19% and 32%.     

绿色荧光蛋白——照亮生命科学的一盏明灯

绿色荧光蛋白的发现及应用具有划时代的重要意义,它不仅为当代生物学研究提供了极为实用的基本研究手段,并且在此基础上改造发展和发现了一系列荧光蛋白,拓展了应用范围。这使得对微观生物学的研究也可以进入一个时空结合,研究鲜活动态过程的新时代。本文主要回顾总结了绿色荧光蛋白的发现、优化改造及其应用。     

Nautilus pompilius hemocyanin: 9 A cryo-EM structure and molecular model reveal the subunit pathway and the interfaces between the 70 functional units

Hemocyanins are giant extracellular oxygen carriers in the hemolymph of many molluscs. Nautilus pompilius (Cephalopoda) hemocyanin is a cylindrical decamer of a 350 kDa polypeptide subunit that in turn is a “pearl-chain” of seven different functional units (FU-a to FU-g). Each globular FU has a binuclear copper centre that reversibly binds one O₂ molecule, and the 70-FU decamer is a highly allosteric protein. Its primary structure and an 11 Å cryo-electron microscopy (cryo-EM) structure have recently been determined, and the crystal structures of two related FU types are available in the databanks. However, in molluscan hemocyanin, the precise subunit pathway within the decamer, the inter-FU interfaces, and the allosteric unit are still obscure, but this knowledge is crucial to understand assembly and allosterism of these proteins. Here we present the cryo-EM structure of Nautilus hemocyanin at 9.1 Å resolution (FSC_1/2-bit criterion), and its molecular model obtained by rigid-body fitting of the individual FUs. In this model we identified the subunit dimer, the subunit pathway, and 15 types of inter-FU interface. Four interface types correspond to the association mode of the two protomers in the published Octopus FU-g crystal. Other interfaces explain previously described morphological structures such as the fenestrated wall (which shows D5 symmetry), the three horizontal wall tiers, the major and minor grooves, the anchor structure and the internal collar (which unexpectedly has C5 symmetry). Moreover, the potential calcium/magnesium and N-glycan binding sites have emerged. Many interfaces have amino acid constellations that might transfer allosteric interaction between FUs. From their topologies we propose that the prime allosteric unit is the oblique segment between major and minor groove, consisting of seven FUs from two different subunits. Thus, the 9 Å structure of Nautilus hemocyanin provides fundamentally new insight into the architecture and function of molluscan hemocyanins.     

A Structural Model for Apolipoprotein C-II Amyloid Fibrils: Experimental Characterization and Molecular Dynamics Simulations

The self-assembly of specific proteins to form insoluble amyloid fibrils is a characteristic feature of a number of age-related and debilitating diseases. Lipid-free human apolipoprotein C-II (apoC-II) forms characteristic amyloid fibrils and is one of several apolipoproteins that accumulate in amyloid deposits located within atherosclerotic plaques. X-ray diffraction analysis of aligned apoC-II fibrils indicated a simple cross-β-structure composed of two parallel β-sheets. Examination of apoC-II fibrils using transmission electron microscopy, scanning transmission electron microscopy, and atomic force microscopy indicated that the fibrils are flat ribbons composed of one apoC-II molecule per 4.7-Å rise of the cross-β-structure. Cross-linking results using single-cysteine substitution mutants are consistent with a parallel in-register structural model for apoC-II fibrils. Fluorescence resonance energy transfer analysis of apoC-II fibrils labeled with specific fluorophores provided distance constraints for selected donor-acceptor pairs located within the fibrils. These findings were used to develop a simple ‘letter-G-like’ β-strand-loop-β-strand model for apoC-II fibrils. Fully solvated all-atom molecular dynamics (MD) simulations showed that the model contained a stable cross-β-core with a flexible connecting loop devoid of persistent secondary structure. The time course of the MD simulations revealed that charge clusters in the fibril rearrange to minimize the effects of same-charge interactions inherent in parallel in-register models. Our structural model for apoC-II fibrils suggests that apoC-II monomers fold and self-assemble to form a stable cross-β-scaffold containing relatively unstructured connecting loops.     

The quaternary structure of a calcium-permeable AMPA receptor: conservation of shape and symmetry across functionally distinct subunit assemblies

AMPA receptors (AMPA-Rs) are formed as heterotetrameric combinations of subunits known as GluR1-GluR4. The calcium permeability of AMPA-Rs is controlled by the identity of the amino-acid side chain contributed by each subunit at a key position in the conductance pathway, which can be either a glutamine (Q) or an arginine (R). Tetramers assembled only from Q-containing subunits are calcium permeable. In contrast, tetramers that incorporate R-containing subunits are calcium impermeable. Both forms play key roles in physiological and pathophysiological processes in the central nervous system. Here, using electron microscopy, we present the first quaternary structure of a calcium-permeable Q-homomeric AMPA-R. The receptor is elongated, with overall 2-fold symmetry and a large central vestibule. It is thus similar to the structure previously reported for an AMPA-R assembled exclusively from R-subunits. Both structures differ from those reported for brain-derived but urea-washed “native” AMPA-Rs, which exhibited multiple asymmetrical conformations. However, even transient exposure of our Q-homomeric AMPA-Rs to urea significantly attenuates the binding of a conformationally specific antibody. As a result, we propose a model in which all AMPA-Rs share a 2-fold symmetrical structure and in which subunit-dependent differences in assembly, trafficking, and electrophysiology are mediated within the framework of fundamentally similar quaternary conformations.     

A modern pollen–climate calibration set from northern Europe: developing and testing a tool for palaeoclimatological reconstructions

Aims and location The potential of pollen records in quantitative climate reconstructions has been widely debated but seldom tested. Our aim is to develop a pollen–climate transfer function for northern Europe and test its performance and inference power by numerical cross‐validation with modern climate data. Annual mean temperature (T_ann) was assessed as the critical climatic variable because T_ann has a distinct south–north gradient (5.5 to ?4.7 °C) in the study region with a corresponding zonal vegetation gradient from the hemiboreal zone in the south to the northern boreal zone in the north. Methods We collected 137 pollen surface samples from small‐ to medium size lakes from southern Estonia to northern Finland. The transfer function for T_ann was developed with weighted averaging partial least squares (WA‐PLS) regression. All 102 terrestrial pollen and spore types were included in the calculation sum and all 137 surface samples and all 102 taxa were included in the transfer function. The performance of the WA‐PLS transfer function was evaluated by leave‐one‐out cross‐validation. Results A cross‐validated root mean square error of prediction (RMSEP) of our model is 0.89 °C and the coefficient of determination (r²) between the observed meteorological T_ann values and those predicted by the model in leave‐one‐out cross‐validation is 0.88. The RMSEP as a percentage of the gradient length of T_ann is 8.8%. These figures indicate high performance statistics for our transfer function compared with other inference models. This is probably because of standardization of our surface‐sampling and pollen‐analytical procedures, careful selection of the surface sample sites with consideration of the relevant pollen source area, the simple patterns of vegetation zones and climate in the study area, and the mostly natural floristic composition of the forests in northern Europe. However, we also demonstrate the limitations of our model in reliably detecting fine‐scale climatic variability. Main conclusions The study shows the strong influence of T_ann on modern pollen composition and demonstrates the potential of pollen data for long‐term climate reconstructions in northern Europe. It also provides evidence against simple interpretations of fine‐scale variations in a single climate reconstruction. In particular, our results highlight the importance of careful study design and implementation in the construction of pollen–climate transfer functions.     

Sex determination of African Penguins Spheniscus demersus using bill measurements: method comparisons and implications for use

African Penguins Spheniscus demersus are sexually dimorphic; on average, males are larger than females but measurements overlap making sex determination difficult through observations alone. We developed a discriminant function, using bill length and depth from a sample of birds sexed from gonad visualisation during post-mortem, which correctly classified 93% of the individuals. Cross-validation correctly assigned 90% of DNA-sexed birds and 91% of birds sexed by partner measurement comparisons. The use of discriminant function score cutpoints, while leaving 16% and 29% of birds unclassified, improved accuracy of birds sexed by DNA to 97% and of those sexed by partner comparison to 99%. Bill depth was found to be a discriminating variable. However, two techniques for measuring bill depth are currently in use for African Penguins. While these measurements are correlated (r = 0.85), they differ on average by 1.4?mm hindering accuracy of sex determination when using a discriminant function developed from the other bill depth measurement. Exploration of adult bill morphology of birds sexed from DNA at different colonies suggests the discriminant functions can be applied throughout the African Penguins’ South African range.