Simulation of transmission electron microscope images of biological specimens

X-rays, electrons, and neutrons as probes of atomic matter

X-rays, electrons, and neutrons probe different properties of matter. X-rays feel electron density (ED). Electrons sense the electrostatic potential (ESP) of electrons and nuclei. Neutrons are sensitive to nuclear coherent scattering length (NCSL). While NCSL maps are widely understood to be different, ED and ESP maps are tacitly assumed to be similar. Here, I show that the belief in ED and ESP map equivalence is mistaken, but contains a grain of truth. Using density functional theory (DFT), the Bethe-Mott (BM) relation, and the Thomas-Fermi (TF) and Cromer-Mann (CM) atomic models, I show that ED and ESP maps are indeed more similar to each other than to NCSL maps. Nevertheless, peak and integrated map values depend differently on the atomic order number and on the contributions from electrons in the inner and outer CM shells. ED and ESP maps also differ in the sign and relative magnitude of excess charge effects.

Electron videography of a lipid-protein tango

Biological phenomena, from enzymatic catalysis to synaptic transmission, originate in the structural transformations of biomolecules and biomolecular assemblies in liquid water. However, directly imaging these nanoscopic dynamics without probes or labels has been a fundamental methodological challenge. Here, we developed an approach for "electron videography"-combining liquid phase electron microscopy with molecular modeling-with which we filmed the nanoscale structural fluctuations of individual, suspended, and unlabeled membrane protein nanodiscs in liquid. Systematic comparisons with biochemical data and simulation indicate the graphene encapsulation involved can afford sufficiently reduced effects of the illuminating electron beam for these observations to yield quantitative fingerprints of nanoscale lipid-protein interactions. Our results suggest that lipid-protein interactions delineate dynamically modified membrane domains across unexpectedly long ranges. Moreover, they contribute to the molecular mechanics of the nanodisc as a whole in a manner specific to the protein within. Overall, this work illustrates an experimental approach to film, quantify, and understand biomolecular dynamics at the nanometer scale.

Cryo-forum: A framework for orientation recovery with uncertainty measure with the application in cryo-EM image analysis

In single-particle cryo-electron microscopy (cryo-EM), efficient determination of orientation parameters for particle images poses a significant challenge yet is crucial for reconstructing 3D structures. This task is complicated by the high noise levels in the datasets, which often include outliers, necessitating several time-consuming 2D clean-up processes. Recently, solutions based on deep learning have emerged, offering a more streamlined approach to the traditionally laborious task of orientation estimation. These solutions employ amortized inference, eliminating the need to estimate parameters individually for each image. However, these methods frequently overlook the presence of outliers and may not adequately concentrate on the components used within the network. This paper introduces a novel method using a 10-dimensional feature vector for orientation representation, extracting orientations as unit quaternions with an accompanying uncertainty metric. Furthermore, we propose a unique loss function that considers the pairwise distances between orientations, thereby enhancing the accuracy of our method. Finally, we also comprehensively evaluate the design choices in constructing the encoder network, a topic that has not received sufficient attention in the literature. Our numerical analysis demonstrates that our methodology effectively recovers orientations from 2D cryo-EM images in an end-to-end manner. Notably, the inclusion of uncertainty quantification allows for direct clean-up of the dataset at the 3D level. Lastly, we package our proposed methods into a user-friendly software suite named cryo-forum, designed for easy access by developers.

Small, solubilized platinum nanocrystals consist of an ordered core surrounded by mobile surface atoms

In situ structures of Platinum (Pt) nanoparticles (NPs) can be determined with graphene liquid cell transmission electron microscopy. Atomic-scale three-dimensional structural information about their physiochemical properties in solution is critical for understanding their chemical function. We here analyze eight atomic-resolution maps of small (<3 nm) colloidal Pt NPs. Their structures are composed of an ordered crystalline core surrounded by surface atoms with comparatively high mobility. 3D reconstructions calculated from cumulative doses of 8500 and 17,000 electrons/pixel, respectively, are characterized in terms of loss of atomic densities and atomic displacements. Less than 5% of the total number of atoms are lost due to dissolution or knock-on damage in five of the structures analyzed, whereas 10-16% are lost in the remaining three. Less than 5% of the atomic positions are displaced due to the increased electron irradiation in all structures. The surface dynamics will play a critical role in the diverse catalytic function of Pt NPs and must be considered in efforts to model Pt NP function computationally.

Molecular Organization of the Skin Barrier

Cryo-electron microscopy of vitreous sections allows for investigation directly in situ of the molecular architecture of skin. Recently, this technique has contributed to the elucidation of the molecular organization of the skin's permeability barrier and its stepwise formation process. The aim of this review is to provide an overview of the procedure for cryo-electron microscopy of vitreous sections, its analysis using atomic detail molecular dynamics modelling and electron microscopy simulation, and its application in the investigation of the barrier structure and formation process of the skin.

TomoTwin: generalized 3D localization of macromolecules in cryo-electron tomograms with structural data mining

Cryogenic-electron tomography enables the visualization of cellular environments in extreme detail, however, tools to analyze the full amount of information contained within these densely packed volumes are still needed. Detailed analysis of macromolecules through subtomogram averaging requires particles to first be localized within the tomogram volume, a task complicated by several factors including a low signal to noise ratio and crowding of the cellular space. Available methods for this task suffer either from being error prone or requiring manual annotation of training data. To assist in this crucial particle picking step, we present TomoTwin: an open source general picking model for cryogenic-electron tomograms based on deep metric learning. By embedding tomograms in an information-rich, high-dimensional space that separates macromolecules according to their three-dimensional structure, TomoTwin allows users to identify proteins in tomograms de novo without manually creating training data or retraining the network to locate new proteins.

A strategy combining denoising and cryo-EM single particle analysis

In cryogenic electron microscopy (cryo-EM) single particle analysis (SPA), high-resolution three-dimensional structures of biological macromolecules are determined by iteratively aligning and averaging a large number of two-dimensional projections of molecules. Since the correlation measures are sensitive to the signal-to-noise ratio, various parameter estimation steps in SPA will be disturbed by the high-intensity noise in cryo-EM. However, denoising algorithms tend to damage high frequencies and suppress mid- and high-frequency contrast of micrographs, which exactly the precise parameter estimation relies on, therefore, limiting their application in SPA. In this study, we suggest combining a cryo-EM image processing pipeline with denoising and maximizing the signal's contribution in various parameter estimation steps. To solve the inherent flaws of denoising algorithms, we design an algorithm named MScale to correct the amplitude distortion caused by denoising and propose a new orientation determination strategy to compensate for the high-frequency loss. In the experiments on several real datasets, the denoised particles are successfully applied in the class assignment estimation and orientation determination tasks, ultimately enhancing the quality of biomacromolecule reconstruction. The case study on classification indicates that our strategy not only improves the resolution of difficult classes (up to 5 Å) but also resolves an additional class. In the case study on orientation determination, our strategy improves the resolution of the final reconstructed density map by 0.34 Å compared with conventional strategy. The code is available at https://github.com/zhanghui186/Mscale.

MEPSi: A tool for simulating tomograms of membrane-embedded proteins

The throughput and fidelity of cryogenic cellular electron tomography (cryo-ET) is constantly increasing through advances in cryogenic electron microscope hardware, direct electron detection devices, and powerful image processing algorithms. However, the need for careful optimization of sample preparations and for access to expensive, high-end equipment, make cryo-ET a costly and time-consuming technique. Generally, only after the last step of the cryo-ET workflow, when reconstructed tomograms are available, it becomes clear whether the chosen imaging parameters were suitable for a specific type of sample in order to answer a specific biological question. Tools for a-priory assessment of the feasibility of samples to answer biological questions and how to optimize imaging parameters to do so would be a major advantage. Here we describe MEPSi (Membrane Embedded Protein Simulator), a simulation tool aimed at rapid and convenient evaluation and optimization of cryo-ET data acquisition parameters for studies of transmembrane proteins in their native environment. We demonstrate the utility of MEPSi by showing how to detangle the influence of different data collection parameters and different orientations in respect to tilt axis and electron beam for two examples: (1) simulated plasma membranes with embedded single-pass transmembrane αIIbβ3 integrin receptors and (2) simulated virus membranes with embedded SARS-CoV-2 spike proteins.

Skin permeability prediction with MD simulation sampling spatial and alchemical reaction coordinates

A molecular-level understanding of skin permeation may rationalize and streamline product development, and improve quality and control, of transdermal and topical drug delivery systems. It may also facilitate toxicity and safety assessment of cosmetics and skin care products. Here, we present new molecular dynamics simulation approaches that make it possible to efficiently sample the free energy and local diffusion coefficient across the skin's barrier structure to predict skin permeability and the effects of chemical penetration enhancers. In particular, we introduce a new approach to use two-dimensional reaction coordinates in the accelerated weight histogram method, where we combine sampling along spatial coordinates with an alchemical perturbation virtual coordinate. We present predicted properties for 20 permeants, and demonstrate how our approach improves correlation with ex vivo/in vitro skin permeation data. For the compounds included in this study, the obtained log KPexp-calc mean square difference was 0.9 cm2 h-2.

Structural basis for the inhibition of IAPP fibril formation by the co-chaperonin prefoldin

Chaperones, as modulators of protein conformational states, are key cellular actors to prevent the accumulation of fibrillar aggregates. Here, we integrated kinetic investigations with structural studies to elucidate how the ubiquitous co-chaperonin prefoldin inhibits diabetes associated islet amyloid polypeptide (IAPP) fibril formation. We demonstrated that both human and archaeal prefoldin interfere similarly with the IAPP fibril elongation and secondary nucleation pathways. Using archaeal prefoldin model, we combined nuclear magnetic resonance spectroscopy with electron microscopy to establish that the inhibition of fibril formation is mediated by the binding of prefoldin’s coiled-coil helices to the flexible IAPP N-terminal segment accessible on the fibril surface and fibril ends. Atomic force microscopy demonstrates that binding of prefoldin to IAPP leads to the formation of lower amounts of aggregates, composed of shorter fibrils, clustered together. Linking structural models with observed fibrillation inhibition processes opens perspectives for understanding the interference between natural chaperones and formation of disease-associated amyloids.

Integrated kinetic and structural investigations reveal that the ubiquitous co-chaperonin prefoldin interacts with its coiled-coil helices on the islet amyloid polypeptide fibril surface and fibril ends to inhibit fibril elongation and secondary nucleation.

Atomic-number (Z)-correlated atomic sizes for deciphering electron microscopic molecular images

Atomic resolution transmission electron microscopy (TEM) has opened up a new era of molecular science by providing atomic video images of dynamic motions of single organic and inorganic molecules. However, the images often look different from the images of molecular models, because these models are designed to visualize the electronic properties of the molecule instead of nuclear electrostatic potentials that are felt by the e-beam in TEM imaging. Here, we propose a molecular model that reproduces TEM images using atomic radii correlated to atomic number (Z). The model serves to provide a priori a useful idea of how a single molecule, molecular assemblies, and thin crystals of organic or inorganic materials look in TEM.

With the advent of atomic resolution transmission electron microscopy (AR-TEM) achieving sub-Ångstrom image resolution and submillisecond time resolution, an era of cinematic molecular science where chemists can visually study the time evolution of molecular motions and reactions at atomistic precision has arrived. However, the appearance of experimental TEM images often differs greatly from that of conventional molecular models, and the images are difficult to decipher unless we know in advance the structure of the specimen molecules. The difference arises from the fundamental design of the molecular models that represent atomic connectivity and/or the electronic properties of molecules rather than the nuclear charge of atoms and electrostatic potentials that are felt by the e-beam in TEM imaging. We found a good correlation between the atomic number (Z) and the atomic size seen in TEM images when we consider shot noise in digital images. We propose Z-correlated (ZC) atomic radii for modeling AR-TEM images of single molecules and ultrathin crystals with which we can develop a good estimate of the molecular structure from the TEM image much more easily than with conventional molecular models. Two parameter sets were developed for TEM images recorded under high-noise (ZC_HN) and low-noise (ZC_LN) conditions. The molecular models will stimulate the imaginations of chemists planning to use AR-TEM for their research.

Cryo-TEM simulations of amorphous radiation-sensitive samples using multislice wave propagation

Image simulation plays a central role in the development and practice of high-resolution electron microscopy, including transmission electron microscopy of frozen-hydrated specimens (cryo-EM). Simulating images with contrast that matches the contrast observed in experimental images remains challenging, especially for amorphous samples. Current state-of-the-art simulators apply post hoc scaling to approximate empirical solvent contrast, attenuated image intensity due to specimen thickness and amplitude contrast. This practice fails for images that require spatially variable scaling, e.g. simulations of a crowded or cellular environment. Modeling both the signal and the noise accurately is necessary to simulate images of biological specimens with contrast that is correct on an absolute scale. The 'frozen plasmon' method is introduced to explicitly model spatially variable inelastic scattering processes in cryo-EM specimens. This approach produces amplitude contrast that depends on the atomic composition of the specimen, reproduces the total inelastic mean free path as observed experimentally and allows for the incorporation of radiation damage in the simulation. These improvements are quantified using the matched filter concept to compare simulation and experiment. The frozen plasmon method, in combination with a new mathematical formulation for accurately sampling the tabulated atomic scattering potentials onto a Cartesian grid, is implemented in the open-source software package cisTEM.

SINGLE: Atomic-resolution structure identification of nanocrystals by graphene liquid cell EM

Analysis of the three-dimensional (3D) structures of nanocrystals with solution-phase transmission electron microscopy is beginning to reveal their unique physiochemical properties. We developed a "one-particle Brownian 3D reconstruction method" based on imaging of ensembles of colloidal nanocrystals using graphene liquid cell electron microscopy. Projection images of differently rotated nanocrystals are acquired using a direct electron detector with high temporal (<2.5 ms) resolution and analyzed to obtain an ensemble of 3D reconstructions. Here, we introduce computational methods required for successful atomic-resolution 3D reconstruction: (i) tracking of the individual particles throughout the time series, (ii) subtraction of the interfering background of the graphene liquid cell, (iii) identification and rejection of low-quality images, and (iv) tailored strategies for 2D/3D alignment and averaging that differ from those used in biological cryo-electron microscopy. Our developments are made available through the open-source software package SINGLE.

Machine Learning to Reveal Nanoparticle Dynamics from Liquid-Phase TEM Videos

Liquid-phase transmission electron microscopy (TEM) has been recently applied to materials chemistry to gain fundamental understanding of various reaction and phase transition dynamics at nanometer resolution. However, quantitative extraction of physical and chemical parameters from the liquid-phase TEM videos remains bottlenecked by the lack of automated analysis methods compatible with the videos' high noisiness and spatial heterogeneity. Here, we integrate, for the first time, liquid-phase TEM imaging with our customized analysis framework based on a machine learning model called U-Net neural network. This combination is made possible by our workflow to generate simulated TEM images as the training data with well-defined ground truth. We apply this framework to three typical systems of colloidal nanoparticles, concerning their diffusion and interaction, reaction kinetics, and assembly dynamics, all resolved in real-time and real-space by liquid-phase TEM. A diversity of properties for differently shaped anisotropic nanoparticles are mapped, including the anisotropic interaction landscape of nanoprisms, curvature-dependent and staged etching profiles of nanorods, and an unexpected kinetic law of first-order chaining assembly of concave nanocubes. These systems representing properties at the nanoscale are otherwise experimentally inaccessible. Compared to the prevalent image segmentation methods, U-Net shows a superior capability to predict the position and shape boundary of nanoparticles from highly noisy and fluctuating background-a challenge common and sometimes inevitable in liquid-phase TEM videos. We expect our framework to push the potency of liquid-phase TEM to its full quantitative level and to shed insights, in high-throughput and statistically significant fashion, on the nanoscale dynamics of synthetic and biological nanomaterials.

Could Egg White Lysozyme be Solved by Single Particle Cryo-EM?

The combination of high-end cryogenic transmission electron microscopes (cryo-EM), direct electron detectors, and advanced image algorithms allows researchers to obtain the 3D structures of much smaller macromolecules than years ago. However, there are still major challenges for the single-particle cryo-EM method to achieve routine structure determinations for macromolecules much smaller than 100 kDa, which are the majority of all plant and animal proteins. These challenges include sample characteristics such as sample heterogeneity, beam damage, ice layer thickness, stability, and quality, as well as hardware limitations such as detector performance, beam, and phase plate quality. Here, single particle data sets were simulated for samples that were ideal in terms of homogeneity, distribution, and stability, but with realistic parameters for ice layer, dose, detector performance, and beam characteristics. Reference data were calculated for human apo-ferritin using identical parameters reported for an experimental data set downloaded from EMPIAR. Processing of the simulated data set resulted in a value of 1.86 Å from 20 214 particles, similar to a 2 Å density map obtained from 29 224 particles selected from real micrographs. Simulated data sets were then generated for a 14 kDa protein, hen egg white lysozyme (HEWL), with and without an ideal phase plate (PP). Whereas we could not obtain a high-resolution 3D reconstruction of HEWL for the data set without PP, the one with PP resulted in a 2.78 Å resolution density map from 225 751 particles. Our simulator and simulations could help in pushing the size limits of cryo-EM.

Structural Transitions in Ceramide Cubic Phases during Formation of the Human Skin Barrier

The stratum corneum is the outermost layer of human skin and the primary barrier toward the environment. The barrier function is maintained by stacked layers of saturated long-chain ceramides, free fatty acids, and cholesterol. This structure is formed through a reorganization of glycosylceramide-based bilayers with cubic-like symmetry into ceramide-based bilayers with stacked lamellar symmetry. The process is accompanied by deglycosylation of glycosylceramides and dehydration of the skin barrier lipid structure. Using coarse-grained molecular dynamics simulation, we show the effects of deglycosylation and dehydration on bilayers of human skin glycosylceramides and ceramides, folded in three dimensions with cubic (gyroid) symmetry. Deglycosylation of glycosylceramides destabilizes the cubic lipid bilayer phase and triggers a cubic-to-lamellar phase transition. Furthermore, subsequent dehydration of the deglycosylated lamellar ceramide system closes the remaining pores between adjacent lipid layers and locally induces a ceramide chain transformation from a hairpin-like to a splayed conformation.

Human skin barrier structure and function analyzed by cryo-EM and molecular dynamics simulation

In the present study we have analyzed the molecular structure and function of the human skin's permeability barrier using molecular dynamics simulation validated against cryo-electron microscopy data from near native skin. The skin's barrier capacity is located to an intercellular lipid structure embedding the cells of the superficial most layer of skin - the stratum corneum. According to the splayed bilayer model (Iwai et al., 2012) the lipid structure is organized as stacked bilayers of ceramides in a splayed chain conformation with cholesterol associated with the ceramide sphingoid moiety and free fatty acids associated with the ceramide fatty acid moiety. However, knowledge about the lipid structure's detailed molecular organization, and the roles of its different lipid constituents, remains circumstantial. Starting from a molecular dynamics model based on the splayed bilayer model, we have, by stepwise structural and compositional modifications, arrived at a thermodynamically stable molecular dynamics model expressing simulated electron microscopy patterns matching original cryo-electron microscopy patterns from skin extremely closely. Strikingly, the closer the individual molecular dynamics models' lipid composition was to that reported in human stratum corneum, the better was the match between the models' simulated electron microscopy patterns and the original cryo-electron microscopy patterns. Moreover, the closest-matching model's calculated water permeability and thermotropic behaviour were found compatible with that of human skin. The new model may facilitate more advanced physics-based skin permeability predictions of drugs and toxicants. The proposed procedure for molecular dynamics based analysis of cellular cryo-electron microscopy data might be applied to other biomolecular systems.

Tunable Ampere phase plate for low dose imaging of biomolecular complexes

A novel device that can be used as a tunable support-free phase plate for transmission electron microscopy of weakly scattering specimens is described. The device relies on the generation of a controlled phase shift by the magnetic field of a segment of current-carrying wire that is oriented parallel or antiparallel to the electron beam. The validity of the concept is established using both experimental electron holographic measurements and a theoretical model based on Ampere's law. Computer simulations are used to illustrate the resulting contrast enhancement for studies of biological cells and macromolecules.

Polymersomes with asymmetric membranes and self-assembled superstructures using pentablock quintopolymers resolved by electron tomography

Polystyrene-block-poly(1,4-isoprene)-block-poly(dimethyl siloxane)-block-poly(tert-butyl methacrylate)-block-poly(2-vinyl pyridine), PS-b-PI-b-PDMS-b-PtBMA-b-P2VP, self-assembles in acetone into polymersomes with asymmetric (directional) PI-b-PDMS membranes.

A classical view on nonclassical nucleation

Understanding and controlling nucleation is important for many crystallization applications. Calcium carbonate (CaCO₃) is often used as a model system to investigate nucleation mechanisms. Despite its great importance in geology, biology, and many industrial applications, CaCO₃ nucleation is still a topic of intense discussion, with new pathways for its growth from ions in solution proposed in recent years. These new pathways include the so-called nonclassical nucleation mechanism via the assembly of thermodynamically stable prenucleation clusters, as well as the formation of a dense liquid precursor phase via liquid-liquid phase separation. Here, we present results from a combined experimental and computational investigation on the precipitation of CaCO₃ in dilute aqueous solutions. We propose that a dense liquid phase (containing 4-7 H₂O per CaCO₃ unit) forms in supersaturated solutions through the association of ions and ion pairs without significant participation of larger ion clusters. This liquid acts as the precursor for the formation of solid CaCO₃ in the form of vaterite, which grows via a net transfer of ions from solution according to z Ca²⁺ + z CO₃^2- → z CaCO₃ The results show that all steps in this process can be explained according to classical concepts of crystal nucleation and growth, and that long-standing physical concepts of nucleation can describe multistep, multiphase growth mechanisms.

Single-protein detection in crowded molecular environments in cryo-EM images

We present an approach to study macromolecular assemblies by detecting component proteins' characteristic high-resolution projection patterns, calculated from their known 3D structures, in single electron cryo-micrographs. Our method detects single apoferritin molecules in vitreous ice with high specificity and determines their orientation and location precisely. Simulations show that high spatial-frequency information and-in the presence of protein background-a whitening filter are essential for optimal detection, in particular for images taken far from focus. Experimentally, we could detect small viral RNA polymerase molecules, distributed randomly among binding locations, inside rotavirus particles. Based on the currently attainable image quality, we estimate a threshold for detection that is 150 kDa in ice and 300 kDa in 100 nm thick samples of dense biological material.

Computation in electron microscopy

Some uses of the computer and computation in high-resolution transmission electron microscopy are reviewed. The theory of image calculation using Bloch wave and multislice methods with and without aberration correction is reviewed and some applications are discussed. The inverse problem of reconstructing the specimen structure from an experimentally measured electron microscope image is discussed. Some future directions of software development are given.

Phase contrast STEM for thin samples: Integrated differential phase contrast

It has been known since the 1970s that the movement of the center of mass (COM) of a convergent beam electron diffraction (CBED) pattern is linearly related to the (projected) electrical field in the sample. We re-derive a contrast transfer function (CTF) for a scanning transmission electron microscopy (STEM) imaging technique based on this movement from the point of view of image formation and continue by performing a two-dimensional integration on the two images based on the two components of the COM movement. The resulting integrated COM (iCOM) STEM technique yields a scalar image that is linear in the phase shift caused by the sample and therefore also in the local (projected) electrostatic potential field of a thin sample. We confirm that the differential phase contrast (DPC) STEM technique using a segmented detector with 4 quadrants (4Q) yields a good approximation for the COM movement. Performing a two-dimensional integration, just as for the COM, we obtain an integrated DPC (iDPC) image which is approximately linear in the phase of the sample. Beside deriving the CTFs of iCOM and iDPC, we clearly point out the objects of the two corresponding imaging techniques, and highlight the differences to objects corresponding to COM-, DPC-, and (HA) ADF-STEM. The theory is validated with simulations and we present first experimental results of the iDPC-STEM technique showing its capability for imaging both light and heavy elements with atomic resolution and a good signal to noise ratio (SNR).

PF2fit: Polar Fast Fourier Matched Alignment of Atomistic Structures with 3D Electron Microscopy Maps

There continue to be increasing occurrences of both atomistic structure models in the PDB (possibly reconstructed from X-ray diffraction or NMR data), and 3D reconstructed cryo-electron microscopy (3D EM) maps (albeit at coarser resolution) of the same or homologous molecule or molecular assembly, deposited in the EMDB. To obtain the best possible structural model of the molecule at the best achievable resolution, and without any missing gaps, one typically aligns (match and fits) the atomistic structure model with the 3D EM map. We discuss a new algorithm and generalized framework, named PF(2) fit (Polar Fast Fourier Fitting) for the best possible structural alignment of atomistic structures with 3D EM. While PF(2) fit enables only a rigid, six dimensional (6D) alignment method, it augments prior work on 6D X-ray structure and 3D EM alignment in multiple ways: Scoring. PF(2) fit includes a new scoring scheme that, in addition to rewarding overlaps between the volumes occupied by the atomistic structure and 3D EM map, rewards overlaps between the volumes complementary to them. We quantitatively demonstrate how this new complementary scoring scheme improves upon existing approaches. PF(2) fit also includes two scoring functions, the non-uniform exterior penalty and the skeleton-secondary structure score, and implements the scattering potential score as an alternative to traditional Gaussian blurring. Search. PF(2) fit utilizes a fast polar Fourier search scheme, whose main advantage is the ability to search over uniformly and adaptively sampled subsets of the space of rigid-body motions. PF(2) fit also implements a new reranking search and scoring methodology that considerably improves alignment metrics in results obtained from the initial search.

Image formation modeling in cryo-electron microscopy

Accurate modeling of image formation in cryo-electron microscopy is an important requirement for quantitative image interpretation and optimization of the data acquisition strategy. Here we present a forward model that accounts for the specimen's scattering properties, microscope optics, and detector response. The specimen interaction potential is calculated with the isolated atom superposition approximation (IASA) and extended with the influences of solvent's dielectric and ionic properties as well as the molecular electrostatic distribution. We account for an effective charge redistribution via the Poisson-Boltzmann approach and find that the IASA-based potential forms the dominant part of the interaction potential, as the contribution of the redistribution is less than 10%. The electron wave is propagated through the specimen by a multislice approach and the influence of the optics is included via the contrast transfer function. We incorporate the detective quantum efficiency of the camera due to the difference between signal and noise transfer characteristics, instead of using only the modulation transfer function. The full model was validated against experimental images of 20S proteasome, hemoglobin, and GroEL. The simulations adequately predict the effects of phase contrast, changes due to the integrated electron flux, thickness, inelastic scattering, detective quantum efficiency and acceleration voltage. We suggest that beam-induced specimen movements are relevant in the experiments whereas the influence of the solvent amorphousness can be neglected. All simulation parameters are based on physical principles and, when necessary, experimentally determined.

Shape-based regularization of electron tomographic reconstruction

We introduce a tomographic reconstruction method implemented using a shape-based regularization technique. Spatial models of known features in the structure being reconstructed are integrated into the reconstruction process as regularizers. Our regularization scheme is driven locally through shape information obtained from segmentation and compared with a known spatial model. We demonstrated our method on tomography data from digital phantoms, simulated data, and experimental electron tomography (ET) data of virus complexes. Our reconstruction showed reduced blurring and an improvement in the resolution of the reconstructed volume was also measured. This method also produced improved demarcation of spike boundaries in viral membranes when compared with popular techniques like weighted back projection and the algebraic reconstruction technique. Improved ET reconstructions will provide better structure elucidation and improved feature visualization, which can aid in solving key biological issues. Our method can also be generalized to other tomographic modalities.

Hydration-layer models for cryo-EM image simulation

To compare cryo-EM images and 3D reconstructions with atomic structures in a quantitative way it is essential to model the electron scattering by solvent (water or ice) that surrounds protein assemblies. The most rigorous method for determining the density of solvating water atoms for this purpose has been to perform molecular-dynamics (MD) simulations of the protein-water system. In this paper we adapt the ideas of bulk-water modeling that are used in the refinement of X-ray crystal structures to the cryo-EM solvent-modeling problem. We present a continuum model for solvent density which matches MD-based results to within sampling errors. However, we also find that the simple binary-mask model of Jiang and Brünger (1994) performs nearly as well as the new model. We conclude that several methods are now available for rapid and accurate modeling of cryo-EM images and maps of solvated proteins.

The human skin barrier is organized as stacked bilayers of fully extended ceramides with cholesterol molecules associated with the ceramide sphingoid moiety

The skin barrier is fundamental to terrestrial life and its evolution; it upholds homeostasis and protects against the environment. Skin barrier capacity is controlled by lipids that fill the extracellular space of the skin's surface layer--the stratum corneum. Here we report on the determination of the molecular organization of the skin's lipid matrix in situ, in its near-native state, using a methodological approach combining very high magnification cryo-electron microscopy (EM) of vitreous skin section defocus series, molecular modeling, and EM simulation. The lipids are organized in an arrangement not previously described in a biological system-stacked bilayers of fully extended ceramides (CERs) with cholesterol molecules associated with the CER sphingoid moiety. This arrangement rationalizes the skin's low permeability toward water and toward hydrophilic and lipophilic substances, as well as the skin barrier's robustness toward hydration and dehydration, environmental temperature and pressure changes, stretching, compression, bending, and shearing.