Quantitative fitting of atomic models into observed densities derived by electron microscopy

High-resolution yeast actin structures indicate the molecular mechanism of actin filament stiffening by cations

Actin filament assembly and the regulation of its mechanical properties are fundamental processes essential for eukaryotic cell function. Residue E167 in vertebrate actins forms an inter-subunit salt bridge with residue K61 of the adjacent subunit. Saccharomyces cerevisiae actin filaments are more flexible than vertebrate filaments and have an alanine at this position (A167). Substitution of this alanine for a glutamic acid (A167E) confers Saccharomyces cerevisiae actin filaments with salt-dependent stiffness similar to vertebrate actins. We developed an optimized cryogenic electron microscopy workflow refining sample preparation and vitrification to obtain near-atomic resolution structures of wild-type and A167E mutant Saccharomyces cerevisiae actin filaments. The difference between these structures allowed us to pinpoint the potential binding site of a filament-associated cation that controls the stiffness of the filaments in vertebrate and A167E Saccharomyces cerevisiae actins. Through an analysis of previously published high-resolution reconstructions of vertebrate actin filaments, along with a newly determined high-resolution vertebrate actin structure in the absence of potassium, we identified a unique peak near residue 167 consistent with the binding of a magnesium ion. Our findings show how magnesium can contribute to filament stiffening by directly bridging actin subunits and allosterically affecting the orientation of the DNase-I binding loop of actin, which plays a regulatory role in modulating actin filament stiffness and interactions with regulatory proteins.

CryoEM Workflow Acceleration with Feret Signatures

Common challenges in cryogenic electron microscopy, such as orientation bias, conformational diversity, and 3D misclassification, complicate single particle analysis and lead to significant resource expenditure. We previously introduced an in silico method using the maximum Feret diameter distribution, the Feret signature, to characterize sample heterogeneity of disc-shaped samples. Here, we expanded the Feret signature methodology to identify preferred orientations of samples containing arbitrary shapes with only about 1000 particles required. This method enables real-time adjustments of data acquisition parameters for optimizing data collection strategies or aiding in decisions to discontinue ineffective imaging sessions. Beyond detecting preferred orientations, the Feret signature approach can serve as an early-warning system for inconsistencies in classification during initial image processing steps, a capability that allows for strategic adjustments in data processing. These features establish the Feret signature as a valuable auxiliary tool in the context of single particle analysis, significantly accelerating the structure determination process.

The GoldX Fiducial Eraser

Gold nanoparticles with sizes in the range of 5-15 nm are a standard method of providing fiducial markers to assist with alignment during reconstruction in cryogenic electron tomography. However, due to their high electron density and resulting contrast when compared to standard cellular or biological samples, they introduce artifacts such as streaking in the reconstructed tomograms. Here, we demonstrate a tool that automatically detects these nanoparticles and suppresses them by replacing them with a local background as a post-processing step, providing a cleaner tomogram without removing any sample relevant information or introducing new artifacts or edge effects from uniform density replacements.

Novel Artificial Intelligence-Based Approaches for Ab Initio Structure Determination and Atomic Model Building for Cryo-Electron Microscopy

Single particle cryo-electron microscopy (cryo-EM) has emerged as the prevailing method for near-atomic structure determination, shedding light on the important molecular mechanisms of biological macromolecules. However, the inherent dynamics and structural variability of biological complexes coupled with the large number of experimental images generated by a cryo-EM experiment make data processing nontrivial. In particular, ab initio reconstruction and atomic model building remain major bottlenecks that demand substantial computational resources and manual intervention. Approaches utilizing recent innovations in artificial intelligence (AI) technology, particularly deep learning, have the potential to overcome the limitations that cannot be adequately addressed by traditional image processing approaches. Here, we review newly proposed AI-based methods for ab initio volume generation, heterogeneous 3D reconstruction, and atomic model building. We highlight the advancements made by the implementation of AI methods, as well as discuss remaining limitations and areas for future development.

Characterization of heterogeneity in nanodisc samples using Feret signatures

Nanodiscs have become a popular tool in structure determination of membrane proteins using cryogenic electron microscopy and single particle analysis. However, the structure determination of small membrane proteins remains challenging. When the embedded protein is in the same size range as the nanodisc, the nanodisc can significantly contribute to the alignment and classification during the structure determination process. In those cases, it is crucial to minimize the heterogeneity in the nanodisc preparations to assure maximum accuracy in the classification and alignment steps of single particle analysis. Here, we introduce a new in-silico method for the characterization of nanodisc samples that is based on analyzing the Feret diameter distribution of their particle projection as imaged in the electron microscope. We validated the method with comprehensive simulation studies and show that Feret signatures can detect subtle differences in nanodisc morphologies and composition that might otherwise go unnoticed. We used the method to identify a specific biochemical nanodisc preparation with low size variations, allowing us to obtain a structure of the 23-kDa single-span membrane protein Bcl-xL while embedded in a nanodisc. Feret signature analysis can steer experimental data collection strategies, allowing more efficient use of high-end data collection hardware, as well as image analysis investments in studies where nanodiscs significantly contribute to the total volume of the full molecular species.

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.

Novel cryo-tomography workflow reveals nanometer-scale responses of epithelial cells to matrix stiffness and dimensionality

Matrix stiffness and dimensionality have been shown to be major determinants of cell behavior. However, a workflow for examining nanometer-scale responses of the associated molecular machinery is not available. Here, we describe a comprehensive, quantitative workflow that permits the analysis of cells responding to mechanical and dimensionality cues in their native state at nanometer scale by cryogenic electron tomography. Using this approach, we quantified distinct cytoskeletal nanoarchitectures and vesicle phenotypes induced in human mammary epithelial cells in response to stiffness and dimensionality of reconstituted basement membrane. Our workflow closely recapitulates the microenvironment associated with acinar morphogenesis and identified distinct differences in situ at nanometer scale. Using drug treatment, we showed that molecular events and nanometer-scale rearrangements triggered by engagement of apical cell receptors with reconstituted basement membrane correspond to changes induced by reduction of cortical tension. Our approach is fully adaptable to any kind of stiffness regime, extracellular matrix composition, and drug treatment.

Overview and applications of map and model validation tools in the CCP-EM software suite

Cryogenic electron microscopy (cryo-EM) has recently been established as a powerful technique for solving macromolecular structures. Although the best resolutions achievable are improving, a significant majority of data are still resolved at resolutions worse than 3 Å, where it is non-trivial to build or fit atomic models. The map reconstructions and atomic models derived from the maps are also prone to errors accumulated through the different stages of data processing. Here, we highlight the need to evaluate both model geometry and fit to data at different resolutions. Assessment of cryo-EM structures from SARS-CoV-2 highlights a bias towards optimising the model geometry to agree with the most common conformations, compared to the agreement with data. We present the CoVal web service which provides multiple validation metrics to reflect the quality of atomic models derived from cryo-EM data of structures from SARS-CoV-2. We demonstrate that further refinement can lead to improvement of the agreement with data without the loss of geometric quality. We also discuss the recent CCP-EM developments aimed at addressing some of the current shortcomings.

ECM dimensionality tunes actin tension to modulate endoplasmic reticulum function and spheroid phenotypes of mammary epithelial cells

Patient-derived organoids and cellular spheroids recapitulate tissue physiology with remarkable fidelity. We investigated how engagement with a reconstituted basement membrane in three dimensions (3D) supports the polarized, stress resilient tissue phenotype of mammary epithelial spheroids. Cells interacting with reconstituted basement membrane in 3D had reduced levels of total and actin-associated filamin and decreased cortical actin tension that increased plasma membrane protrusions to promote negative plasma membrane curvature and plasma membrane protein associations linked to protein secretion. By contrast, cells engaging a reconstituted basement membrane in 2D had high cortical actin tension that forced filamin unfolding and endoplasmic reticulum (ER) associations. Enhanced filamin-ER interactions increased levels of PKR-like ER kinase effectors and ER-plasma membrane contact sites that compromised calcium homeostasis and diminished cell viability. Consequently, cells with decreased cortical actin tension had reduced ER stress and survived better. Consistently, cortical actin tension in cellular spheroids regulated polarized basement membrane membrane deposition and sensitivity to exogenous stress. The findings implicate cortical actin tension-mediated filamin unfolding in ER function and underscore the importance of tissue mechanics in organoid homeostasis.

Exploring cryo-electron microscopy with molecular dynamics

Single particle analysis cryo-electron microscopy (EM) and molecular dynamics (MD) have been complimentary methods since cryo-EM was first applied to the field of structural biology. The relationship started by biasing structural models to fit low-resolution cryo-EM maps of large macromolecular complexes not amenable to crystallization. The connection between cryo-EM and MD evolved as cryo-EM maps improved in resolution, allowing advanced sampling algorithms to simultaneously refine backbone and sidechains. Moving beyond a single static snapshot, modern inferencing approaches integrate cryo-EM and MD to generate structural ensembles from cryo-EM map data or directly from the particle images themselves. We summarize the recent history of MD innovations in the area of cryo-EM modeling. The merits for the myriad of MD based cryo-EM modeling methods are discussed, as well as, the discoveries that were made possible by the integration of molecular modeling with cryo-EM. Lastly, current challenges and potential opportunities are reviewed.

Interface refinement of low- to medium-resolution Cryo-EM complexes using HADDOCK2.4

A wide range of cellular processes requires the formation of multimeric protein complexes. The rise of cryo-electron microscopy (cryo-EM) has enabled the structural characterization of these protein assemblies. The density maps produced can, however, still suffer from limited resolution, impeding the process of resolving structures at atomic resolution. In order to solve this issue, monomers can be fitted into low- to medium-resolution maps. Unfortunately, the models produced frequently contain atomic clashes at the protein-protein interfaces (PPIs), as intermolecular interactions are typically not considered during monomer fitting. Here, we present a refinement approach based on HADDOCK2.4 to remove intermolecular clashes and optimize PPIs. A dataset of 14 cryo-EM complexes was used to test eight protocols. The best-performing protocol, consisting of a semi-flexible simulated annealing refinement with centroid restraints on the monomers, was able to decrease intermolecular atomic clashes by 98% without significantly deteriorating the quality of the cryo-EM density fit.

Atomic model validation using the CCP-EM software suite

Recently, there has been a dramatic improvement in the quality and quantity of data derived using cryogenic electron microscopy (cryo-EM). This is also associated with a large increase in the number of atomic models built. Although the best resolutions that are achievable are improving, often the local resolution is variable, and a significant majority of data are still resolved at resolutions worse than 3 Å. Model building and refinement is often challenging at these resolutions, and hence atomic model validation becomes even more crucial to identify less reliable regions of the model. Here, a graphical user interface for atomic model validation, implemented in the CCP-EM software suite, is presented. It is aimed to develop this into a platform where users can access multiple complementary validation metrics that work across a range of resolutions and obtain a summary of evaluations. Based on the validation estimates from atomic models associated with cryo-EM structures from SARS-CoV-2, it was observed that models typically favor adopting the most common conformations over fitting the observations when compared with the model agreement with data. At low resolutions, the stereochemical quality may be favored over data fit, but care should be taken to ensure that the model agrees with the data in terms of resolvable features. It is demonstrated that further re-refinement can lead to improvement of the agreement with data without the loss of geometric quality. This also highlights the need for improved resolution-dependent weight optimization in model refinement and an effective test for overfitting that would help to guide the refinement process.

Rapid tool for cell nanoarchitecture integrity assessment

With the rapid increase and accessibility of high-resolution imaging technologies of cells, the interpretation of results relies more and more on the assumption that the three-dimensional integrity of the surrounding cellular landscape is not compromised by the experimental setup. However, the only available technology for directly probing the structural integrity of whole-cell preparations at the nanoscale is electron cryo-tomography, which is time-consuming, costly, and complex. We devised an accessible, inexpensive and reliable screening assay to quickly report on the compatibility of experimental protocols with preserving the structural integrity of whole-cell preparations at the nanoscale. Our Rapid Cell Integrity Assessment (RCIA) assay is executed at room temperature and relies solely on light microscopy imaging. Using cellular electron cryo-tomography as a benchmark, we verify that RCIA accurately unveils the adverse impact of reagents and/or protocols such as those used for virus inactivation or to arrest dynamic processes on the cellular nanoarchitecture.

The actomyosin interface contains an evolutionary conserved core and an ancillary interface involved in specificity

Plasmodium falciparum, the causative agent of malaria, moves by an atypical process called gliding motility. Actomyosin interactions are central to gliding motility. However, the details of these interactions remained elusive until now. Here, we report an atomic structure of the divergent Plasmodium falciparum actomyosin system determined by electron cryomicroscopy at the end of the powerstroke (Rigor state). The structure provides insights into the detailed interactions that are required for the parasite to produce the force and motion required for infectivity. Remarkably, the footprint of the myosin motor on filamentous actin is conserved with respect to higher eukaryotes, despite important variability in the Plasmodium falciparum myosin and actin elements that make up the interface. Comparison with other actomyosin complexes reveals a conserved core interface common to all actomyosin complexes, with an ancillary interface involved in defining the spatial positioning of the motor on actin filaments.

Plasmodium falciparum moves by an atypical process called gliding motility which comprises of atypical myosin A (PfMyoA) and filaments of the dynamic and divergent PfActin-1 (PfAct1). Here authors present the cryo-EM structure of PfMyoA bound to filamentous PfAct1 stabilized with jasplakinolide and provide insights into the interactions that are required for the parasite to produce the force and motion required for infectivity.

Structural basis of αE-catenin-F-actin catch bond behavior

Cell-cell and cell-matrix junctions transmit mechanical forces during tissue morphogenesis and homeostasis. α-Catenin links cell-cell adhesion complexes to the actin cytoskeleton, and mechanical load strengthens its binding to F-actin in a direction-sensitive manner. Specifically, optical trap experiments revealed that force promotes a transition between weak and strong actin-bound states. Here, we describe the cryo-electron microscopy structure of the F-actin-bound αE-catenin actin-binding domain, which in solution forms a five-helix bundle. In the actin-bound structure, the first helix of the bundle dissociates and the remaining four helices and connecting loops rearrange to form the interface with actin. Deletion of the first helix produces strong actin binding in the absence of force, suggesting that the actin-bound structure corresponds to the strong state. Our analysis explains how mechanical force applied to αE-catenin or its homolog vinculin favors the strongly bound state, and the dependence of catch bond strength on the direction of applied force.

Simulation-Based Methods for Model Building and Refinement in Cryoelectron Microscopy

Advances in cryoelectron microscopy (cryo-EM) have revolutionized the structural investigation of large macromolecular assemblies. In this review, we first provide a broad overview of modeling methods used for flexible fitting of molecular models into cryo-EM density maps. We give special attention to approaches rooted in molecular simulations-atomistic molecular dynamics and Monte Carlo. Concise descriptions of the methods are given along with discussion of their advantages, limitations, and most popular alternatives. We also describe recent extensions of the widely used molecular dynamics flexible fitting (MDFF) method and discuss how different model-building techniques could be incorporated into new hybrid modeling schemes and simulation workflows. Finally, we provide two illustrative examples of model-building and refinement strategies employing MDFF, cascade MDFF, and RosettaCM. These examples come from recent cryo-EM studies that elucidated transcription preinitiation complexes and shed light on the functional roles of these assemblies in gene expression and gene regulation.

Advances in Structure Modeling Methods for Cryo-Electron Microscopy Maps

Cryo-electron microscopy (cryo-EM) has now become a widely used technique for structure determination of macromolecular complexes. For modeling molecular structures from density maps of different resolutions, many algorithms have been developed. These algorithms can be categorized into rigid fitting, flexible fitting, and de novo modeling methods. It is also observed that machine learning (ML) techniques have been increasingly applied following the rapid progress of the ML field. Here, we review these different categories of macromolecule structure modeling methods and discuss their advances over time.

Simultaneous Determination of Protein Structure and Dynamics Using Cryo-Electron Microscopy

Cryo-electron microscopy is rapidly emerging as a powerful technique to determine the structures of complex macromolecular systems elusive to other techniques. Because many of these systems are highly dynamical, characterizing their movements is also a crucial step to unravel their biological functions. To achieve this goal, we report an integrative modeling approach to simultaneously determine structure and dynamics of macromolecular systems from cryo-electron microscopy density maps. By quantifying the level of noise in the data and dealing with their ensemble-averaged nature, this approach enables the integration of multiple sources of information to model ensembles of structures and infer their populations. We illustrate the method by characterizing structure and dynamics of the integral membrane receptor STRA6, thus providing insights into the mechanisms by which it interacts with retinol binding protein and translocates retinol across the membrane.

Bayesian Weighing of Electron Cryo-Microscopy Data for Integrative Structural Modeling

Cryo-electron microscopy (cryo-EM) has become a mainstream technique for determining the structures of complex biological systems. However, accurate integrative structural modeling has been hampered by the challenges in objectively weighing cryo-EM data against other sources of information due to the presence of random and systematic errors, as well as correlations, in the data. To address these challenges, we introduce a Bayesian scoring function that efficiently and accurately ranks alternative structural models of a macromolecular system based on their consistency with a cryo-EM density map as well as other experimental and prior information. The accuracy of this approach is benchmarked using complexes of known structure and illustrated in three applications: the structural determination of the GroEL/GroES, RNA polymerase II, and exosome complexes. The approach is implemented in the open-source Integrative Modeling Platform (http://integrativemodeling.org), thus enabling integrative structure determination by combining cryo-EM data with other sources of information.

Local Tension on Talin in Focal Adhesions Correlates with F-Actin Alignment at the Nanometer Scale

Cellular force transmission and mechanotransduction are critical in embryogenesis, normal physiology, and many diseases. Talin plays a key role in these processes by linking integrins to force-generating actomyosin. Using the previously characterized FRET-based talin tension sensor, we observed variations of tension both between and within individual focal adhesions in the same cell. Assembling and sliding adhesions showed gradients with higher talin tension toward the cell center, whereas mature, stable adhesions had uniform talin tension. Total talin accumulation was maximal in high-tension regions; by contrast, vinculin intensity was flat or maximal at the adhesion center, and actin intensity was maximal toward the cell center. To investigate mechanism, we combined talin tension imaging with cellular cryotomography to visualize the correlated actin organization at nanometer resolution. Regions of high talin tension had highly aligned linear actin filaments, whereas regions of low tension had less-well-aligned F-actin. These results reveal an orchestrated spatiotemporal relationship between talin tension, actin/vinculin localization, local actin organization, and focal adhesion dynamics.

Analytical Approaches to Improve Accuracy in Solving the Protein Topology Problem

To take advantage of recent advances in genomics and proteomics it is critical that the three-dimensional physical structure of biological macromolecules be determined. Cryo-Electron Microscopy (cryo-EM) is a promising and improving method for obtaining this data, however resolution is often not sufficient to directly determine the atomic scale structure. Despite this, information for secondary structure locations is detectable. De novo modeling is a computational approach to modeling these macromolecular structures based on cryo-EM derived data. During de novo modeling a mapping between detected secondary structures and the underlying amino acid sequence must be identified. DP-TOSS (Dynamic Programming for determining the Topology Of Secondary Structures) is one tool that attempts to automate the creation of this mapping. By treating the correspondence between the detected structures and the structures predicted from sequence data as a constraint graph problem DP-TOSS achieved good accuracy in its original iteration. In this paper, we propose modifications to the scoring methodology of DP-TOSS to improve its accuracy. Three scoring schemes were applied to DP-TOSS and tested: (i) a skeleton-based scoring function; (ii) a geometry-based analytical function; and (iii) a multi-well potential energy-based function. A test of 25 proteins shows that a combination of these schemes can improve the performance of DP-TOSS to solve the topology determination problem for macromolecule proteins.

Marker-free method for accurate alignment between correlated light, cryo-light, and electron cryo-microscopy data using sample support features

Combining fluorescence microscopy with electron cryo-tomography allows, in principle, spatial localization of tagged macromolecular assemblies and structural features within the cellular environment. To allow precise localization and scale integration between the two disparate imaging modalities, accurate alignment procedures are needed. Here, we describe a marker-free method for aligning images from light or cryo-light fluorescence microscopy and from electron cryo-microscopy that takes advantage of sample support features, namely the holes in the carbon film. We find that the accuracy of this method, as judged by prediction errors of the hole center coordinates, is better than 100 nm.

Structure of anthrax lethal toxin prepore complex suggests a pathway for efficient cell entry

Anthrax toxin comprises three soluble proteins: protective antigen (PA), lethal factor (LF), and edema factor (EF). PA must be cleaved by host proteases before it oligomerizes and forms a prepore, to which LF and EF bind. After endocytosis of this tripartite complex, the prepore transforms into a narrow transmembrane pore that delivers unfolded LF and EF into the host cytosol. Here, we find that translocation of multiple 90-kD LF molecules is rapid and efficient. To probe the molecular basis of this translocation, we calculated a three-dimensional map of the fully loaded (PA63)7-(LF)3 prepore complex by cryo-electron microscopy (cryo-EM). The map shows three LFs bound in a similar way to one another, via their N-terminal domains, to the surface of the PA heptamer. The model also reveals contacts between the N- and C-terminal domains of adjacent LF molecules. We propose that this molecular arrangement plays an important role in the maintenance of translocation efficiency through the narrow PA pore.

PEM-fitter: A Coarse-Grained Method to Validate Protein Candidate Models

The volumetric images produced by Cryo-Electron Microscopy (cryo-EM) technique are used to model macromolecular assemblies and machines. De novo protein modeling uses these images to computationally model the structure of the molecules. Many candidate conformations are usually generated during the intermediate step. Conventionally, each of these candidates is evaluated by time-consuming approaches such as potential energy. We introduce an initial version of a geometrical screening method that uses the skeleton of the cryo-EM images to evaluate candidate structures. The aim of this method is to reduce the number of native-like candidate conformations and, therefore, reduce the time required for structural evaluation by energy calculations. A test of two datasets was performed. The first dataset contains 10 proteins and shows that our method can successfully detect the correct native structure for the given skeleton among a set of different protein structures. The second dataset contains 12 proteins and shows that our method can filter slightly modified decoy conformations of the same protein. The efficiency of the method is also reported.

A novel fast Fourier transform accelerated off-grid exhaustive search method for cryo-electron microscopy fitting

This paper presents a novel fast Fourier transform (FFT)-based exhaustive search method extended to off-grid translational and rotational degrees of freedom. The method combines the advantages of the FFT-based exhaustive search, which samples all the conformations of a system under study on a grid, with a local optimization technique that guarantees to find the nearest optimal off-grid conformation. The method is demonstrated on a fitting problem and can be readily applied to a docking problem. The algorithm first samples a scoring function on a six-dimensional grid of sizeN6using the FFT. This operation has an asymptotic complexity ofO(N6logN). Then, the method performs the off-grid search using a local quadratic approximation of the cost function and the trust-region optimization algorithm. The computation of the quadratic approximation is also accelerated by FFT at the same additional asymptotic cost ofO(N6logN). The method is demonstrated by fitting atomic protein models into several simulated and experimental maps from cryo-electron microscopy. The method is available at https://team.inria.fr/nano-d/software/offgridfit.

Improved metrics for comparing structures of macromolecular assemblies determined by 3D electron-microscopy

Recent developments in 3-dimensional electron microcopy (3D-EM) techniques and a concomitant drive to look at complex molecular structures, have led to a rapid increase in the amount of volume data available for biomolecules. This creates a demand for better methods to analyse the data, including improved scores for comparison, classification and integration of data at different resolutions. To this end, we developed and evaluated a set of scoring functions that compare 3D-EM volumes. To test our scores we used a benchmark set of volume alignments derived from the Electron Microscopy Data Bank. We find that the performance of different scores vary with the map-type, resolution and the extent of overlap between volumes. Importantly, adding the overlap information to the local scoring functions can significantly improve their precision and accuracy in a range of resolutions. A combined score involving the local mutual information and overlap (LMI_OV) performs best overall, irrespective of the map category, resolution or the extent of overlap, and we recommend this score for general use. The local mutual information score itself is found to be more discriminatory than cross-correlation coefficient for intermediate-to-low resolution maps or when the map size and density distribution differ significantly. For comparing map surfaces, we implemented two filters to detect the surface points, including one based on the 'extent of surface exposure'. We show that scores that compare surfaces are useful at low resolutions and for maps with evident surface features. All the scores discussed are implemented in TEMPy (http://tempy.ismb.lon.ac.uk/).

Nano-scale actin-network characterization of fibroblast cells lacking functional Arp2/3 complex

Arp2/3 complex is thought to be the primary protrusive force generator in cell migration by controlling the assembly and turnover of the branched filament network that pushes the leading edge of moving cells forward. However, mouse fibroblasts without functional Arp2/3 complex migrate at rates similar to wild-type cells, contradicting this paradigm. We show by correlative fluorescence and large-scale cryo-tomography studies combined with automated actin-network analysis that the absence of functional Arp2/3 complex has profound effects on the nano-scale architecture of actin networks. Our quantitative analysis at the single-filament level revealed that cells lacking functional Arp2/3 complex fail to regulate location-dependent fine-tuning of actin filament growth and organization that is distinct from its role in the formation and regulation of dendritic actin networks.

Three-Dimensional Structures of Full-Length, Membrane-Embedded Human α(IIb)β(3) Integrin Complexes

Integrins are bidirectional, allosteric transmembrane receptors that play a central role in hemostasis and arterial thrombosis. Using cryo-electron microscopy, multireference single-particle reconstruction methods, and statistics-based computational fitting approaches, we determined three-dimensional structures of human integrin αIIbβ3 embedded in a lipid bilayer (nanodiscs) while bound to domains of the cytosolic regulator talin and to extracellular ligands. We also determined the conformations of integrin in solution by itself to localize the membrane and the talin-binding site. To our knowledge, our data provide unprecedented three-dimensional information about the conformational states of intact, full-length integrin within membrane bilayers under near-physiological conditions and in the presence of cytosolic activators and extracellular ligands. We show that αIIbβ3 integrins exist in a conformational equilibrium clustered around four main states. These conformations range from a compact bent nodule to two partially extended intermediate conformers and finally to a fully upright state. In the presence of nanodiscs and the two ligands, the equilibrium is significantly shifted toward the upright conformation. In this conformation, the receptor extends ∼20 nm upward from the membrane. There are no observable contacts between the two subunits other than those in the headpiece near the ligand-binding pocket, and the α- and β-subunits are well separated with their cytoplasmic tails ∼8 nm apart. Our results indicate that extension of the ectodomain is possible without separating the legs or extending the hybrid domain, and that the ligand-binding pocket is not occluded by the membrane in any conformations of the equilibrium. Further, they suggest that integrin activation may be influenced by equilibrium shifts.

Use of evolutionary information in the fitting of atomic level protein models in low resolution cryo-EM map of a protein assembly improves the accuracy of the fitting

Protein-protein interface residues, especially those at the core of the interface, exhibit higher conservation than residues in solvent exposed regions. Here, we explore the ability of this differential conservation to evaluate fittings of atomic models in low-resolution cryo-EM maps and select models from the ensemble of solutions that are often proposed by different model fitting techniques. As a prelude, using a non-redundant and high-resolution structural dataset involving 125 permanent and 95 transient complexes, we confirm that core interface residues are conserved significantly better than nearby non-interface residues and this result is used in the cryo-EM map analysis. From the analysis of inter-component interfaces in a set of fitted models associated with low-resolution cryo-EM maps of ribosomes, chaperones and proteasomes we note that a few poorly conserved residues occur at interfaces. Interestingly a few conserved residues are not in the interface, though they are close to the interface. These observations raise the potential requirement of refitting the models in the cryo-EM maps. We show that sampling an ensemble of models and selection of models with high residue conservation at the interface and in good agreement with the density helps in improving the accuracy of the fit. This study indicates that evolutionary information can serve as an additional input to improve and validate fitting of atomic models in cryo-EM density maps.

Refinement of atomic models in high resolution EM reconstructions using Flex-EM and local assessment

As the resolutions of Three Dimensional Electron Microscopic reconstructions of biological macromolecules are being improved, there is a need for better fitting and refinement methods at high resolutions and robust approaches for model assessment. Flex-EM/MODELLER has been used for flexible fitting of atomic models in intermediate-to-low resolution density maps of different biological systems. Here, we demonstrate the suitability of the method to successfully refine structures at higher resolutions (2.5-4.5Å) using both simulated and experimental data, including a newly processed map of Apo-GroEL. A hierarchical refinement protocol was adopted where the rigid body definitions are relaxed and atom displacement steps are reduced progressively at successive stages of refinement. For the assessment of local fit, we used the SMOC (segment-based Manders' overlap coefficient) score, while the model quality was checked using the Qmean score. Comparison of SMOC profiles at different stages of refinement helped in detecting regions that are poorly fitted. We also show how initial model errors can have significant impact on the goodness-of-fit. Finally, we discuss the implementation of Flex-EM in the CCP-EM software suite.

The neuronal porosome complex in health and disease

Cup-shaped secretory portals at the cell plasma membrane called porosomes mediate the precision release of intravesicular material from cells. Membrane-bound secretory vesicles transiently dock and fuse at the base of porosomes facing the cytosol to expel pressurized intravesicular contents from the cell during secretion. The structure, isolation, composition, and functional reconstitution of the neuronal porosome complex have greatly progressed, providing a molecular understanding of its function in health and disease. Neuronal porosomes are 15 nm cup-shaped lipoprotein structures composed of nearly 40 proteins, compared to the 120 nm nuclear pore complex composed of >500 protein molecules. Membrane proteins compose the porosome complex, making it practically impossible to solve its atomic structure. However, atomic force microscopy and small-angle X-ray solution scattering studies have provided three-dimensional structural details of the native neuronal porosome at sub-nanometer resolution, providing insights into the molecular mechanism of its function. The participation of several porosome proteins previously implicated in neurotransmission and neurological disorders, further attest to the crosstalk between porosome proteins and their coordinated involvement in release of neurotransmitter at the synapse.

Electron Microscopy and Image Processing: Essential Tools for Structural Analysis of Macromolecules

Macromolecular electron microscopy typically depicts the structures of macromolecular complexes ranging from ∼200 kDa to hundreds of MDa. The amount of specimen required, a few micrograms, is typically 100 to 1000 times less than needed for X-ray crystallography or nuclear magnetic resonance spectroscopy. Micrographs of frozen-hydrated (cryogenic) specimens portray native structures, but the original images are noisy. Computational averaging reduces noise, and three-dimensional reconstructions are calculated by combining different views of free-standing particles ("single-particle analysis"). Electron crystallography is used to characterize two-dimensional arrays of membrane proteins and very small three-dimensional crystals. Under favorable circumstances, near-atomic resolutions are achieved. For structures at somewhat lower resolution, pseudo-atomic models are obtained by fitting high-resolution components into the density. Time-resolved experiments describe dynamic processes. Electron tomography allows reconstruction of pleiomorphic complexes and subcellular structures and modeling of macromolecules in their cellular context. Significant information is also obtained from metal-coated and dehydrated specimens.

Integrative Modeling of Macromolecular Assemblies from Low to Near-Atomic Resolution

While conventional high-resolution techniques in structural biology are challenged by the size and flexibility of many biological assemblies, recent advances in low-resolution techniques such as cryo-electron microscopy (cryo-EM) and small angle X-ray scattering (SAXS) have opened up new avenues to define the structures of such assemblies. By systematically combining various sources of structural, biochemical and biophysical information, integrative modeling approaches aim to provide a unified structural description of such assemblies, starting from high-resolution structures of the individual components and integrating all available information from low-resolution experimental methods. In this review, we describe integrative modeling approaches, which use complementary data from either cryo-EM or SAXS. Specifically, we focus on the popular molecular dynamics flexible fitting (MDFF) method, which has been widely used for flexible fitting into cryo-EM maps. Second, we describe hybrid molecular dynamics, Rosetta Monte-Carlo and minimum ensemble search (MES) methods that can be used to incorporate SAXS into pseudoatomic structural models. We present concise descriptions of the two methods and their most popular alternatives, along with select illustrative applications to protein/nucleic acid assemblies involved in DNA replication and repair.

Cell adhesion. The minimal cadherin-catenin complex binds to actin filaments under force

Linkage between the adherens junction (AJ) and the actin cytoskeleton is required for tissue development and homeostasis. In vivo findings indicated that the AJ proteins E-cadherin, β-catenin, and the filamentous (F)-actin binding protein αE-catenin form a minimal cadherin-catenin complex that binds directly to F-actin. Biochemical studies challenged this model because the purified cadherin-catenin complex does not bind F-actin in solution. Here, we reconciled this difference. Using an optical trap-based assay, we showed that the minimal cadherin-catenin complex formed stable bonds with an actin filament under force. Bond dissociation kinetics can be explained by a catch-bond model in which force shifts the bond from a weakly to a strongly bound state. These results may explain how the cadherin-catenin complex transduces mechanical forces at cell-cell junctions.

Neuronal Porosome-The Secretory Portal at the Nerve Terminal: It's Structure-Function, Composition, and Reconstitution

Cup-shaped secretory portals at the cell plasma membrane called porosomes mediate secretion from cells. Membrane bound secretory vesicles transiently dock and fuse at the cytosolic compartment of the porosome base to expel intravesicular contents to the outside during cell secretion. In the past decade, the structure, isolation, composition, and functional reconstitution of the neuronal porosome complex has been accomplished providing a molecular understanding of its structure-function. Neuronal porosomes are 15 nm cup-shaped lipoprotein structures composed of nearly 40 proteins. Being a membrane-associated supramolecular complex has precluded determination of the atomic structure of the porosome. However recent studies using small-angle X-ray solution scattering (SAXS), provide at sub-nanometer resolution, the native 3D structure of the neuronal porosome complex associated with docked synaptic vesicle at the nerve terminal. Additionally, results from the SAXS study and earlier studies using atomic force microscopy, provide the possible molecular mechanism involved in porosome-mediated neurotransmitter release at the nerve terminal.

Holoenzyme structures of endothelial nitric oxide synthase - an allosteric role for calmodulin in pivoting the FMN domain for electron transfer

While the three-dimensional structures of heme- and flavin-binding domains of the NOS isoforms have been determined, the structures of the holoenzymes remained elusive. Application of electron cryo-microscopy and structural modeling of the bovine endothelial nitric oxide synthase (eNOS) holoenzyme produced detailed models of the intact holoenzyme in the presence and absence of Ca(2+)/calmodulin (CaM). These models accommodate the cross-electron transfer from the reductase in one monomer to the heme in the opposite monomer. The heme domain acts as the anchoring dimeric structure for the entire enzyme molecule, while the FMN domain is activated by CaM to move flexibly to bridge the distance between the reductase and oxygenase domains. Our results indicate that the key regulatory role of CaM involves the stabilization of structural intermediates and precise positioning of the pivot for the FMN domain tethered shuttling motion to accommodate efficient and rapid electron transfer in the homodimer of eNOS.

HermiteFit: fast-fitting atomic structures into a low-resolution density map using three-dimensional orthogonal Hermite functions

HermiteFit, a novel algorithm for fitting a protein structure into a low-resolution electron-density map, is presented. The algorithm accelerates the rotation of the Fourier image of the electron density by using three-dimensional orthogonal Hermite functions. As part of the new method, an algorithm for the rotation of the density in the Hermite basis and an algorithm for the conversion of the expansion coefficients into the Fourier basis are presented. HermiteFit was implemented using the cross-correlation or the Laplacian-filtered cross-correlation as the fitting criterion. It is demonstrated that in the Hermite basis the Laplacian filter has a particularly simple form. To assess the quality of density encoding in the Hermite basis, an analytical way of computing the crystallographic R factor is presented. Finally, the algorithm is validated using two examples and its efficiency is compared with two widely used fitting methods, ADP_EM and colores from the Situs package. HermiteFit will be made available at http://nano-d.inrialpes.fr/software/HermiteFit or upon request from the authors.

Porosome in Cystic Fibrosis

Macromolecular structures embedded in the cell plasma membrane called ‘porosomes’, are involved in the regulated fractional release of intravesicular contents from cells during secretion. Porosomes range in size from 15 nm in neurons and astrocytes to 100-180 nm in the exocrine pancreas and neuroendocrine cells. Porosomes have been isolated from a number of cells, and their morphology, composition, and functional reconstitution well documented. The 3D contour map of the assembly of proteins within the porosome complex, and its native X-ray solution structure at sub-nm resolution has also advanced. This understanding now provides a platform to address diseases that may result from secretory defects. Water and ion binding to mucin impart hydration, critical for regulating viscosity of the mucus in the airways epithelia. Appropriate viscosity is required for the movement of mucus by the underlying cilia. Hence secretion of more viscous mucus prevents its proper transport, resulting in chronic and fatal airways disease such as cystic fibrosis (CF). CF is caused by the malfunction of CF transmembrane conductance regulator (CFTR), a chloride channel transporter, resulting in viscous mucus in the airways. Studies in mice lacking functional CFTR secrete highly viscous mucous that adhered to the epithelium. Since CFTR is known to interact with the t-SNARE protein syntaxin-1A, and with the chloride channel CLC-3, which are also components of the porosome complex, the interactions between CFTR and the porosome complex in the mucin-secreting human airway epithelial cell line Calu-3 was hypothesized and tested. Results from the study demonstrate the presence of approximately 100 nm in size porosome complex composed of 34 proteins at the cell plasma membrane in Calu-3 cells, and the association of CFTR with the complex. In comparison, the nuclear pore complex measures 120 nm and is comprised of over 500 protein molecules. The involvement of CFTR in porosome-mediated mucin secretion is hypothesized, and is currently being tested.

Conformational States of macromolecular assemblies explored by integrative structure calculation

A detailed description of macromolecular assemblies in multiple conformational states can be very valuable for understanding cellular processes. At present, structural determination of most assemblies in different biologically relevant conformations cannot be achieved by a single technique and thus requires an integrative approach that combines information from multiple sources. Different techniques require different computational methods to allow efficient and accurate data processing and analysis. Here, we summarize the latest advances and future challenges in computational methods that help the interpretation of data from two techniques—mass spectrometry and three-dimensional cryo-electron microscopy (with focus on alignment and classification of heterogeneous subtomograms from cryo-electron tomography). We evaluate how new developments in these two broad fields will lead to further integration with atomic structures to broaden our picture of the dynamic behavior of assemblies in their native environment.

Three-dimensional reconstructions of actin filaments capped by Arp2/3 complex

The primary function of Arp2/3 complex is the generation of free barbed ends by nucleating new filaments from the sides of pre-existing filaments. The pathway of branch formation is complex and involves nucleation promoting factors, actin monomers and nucleotides. A less prominent function of Arp2/3 complex is capping of actin filament pointed ends. Here we show, using electron microscopy, electron tomography, and image reconstruction of negatively-stained samples at ∼2-3nm resolution, that Arp2/3 complex bound to the pointed ends of actin filaments has a conformation similar to that in the branch junction with the Arps arranged in an actin-filament like configuration. This is direct evidence for the existence of two distinct activation pathways for Arp2/3 complex, one in the context of branch formation, one in the context of pointed-end capping, with essentially the same conformational end point.

The joys and perils of flexible fitting

While performing their functions, biological macromolecules often form large, dynamically changing macromolecular assemblies. Only a relatively small number of such assemblies have been accessible to the atomic-resolution techniques X-ray crystallography and NMR. Electron microscopy in conjunction with image reconstruction has become the preferred alternative for revealing the structures of such macromolecular complexes. However, for most assemblies the achievable resolution is too low to allow accurate atomic modeling directly from the data. Yet, useful models often can be obtained by fitting atomic models of individual components into a low-resolution reconstruction of the entire assembly. Several algorithms for achieving optimal fits in this context were developed recently, many allowing considerable degrees of flexibility to account for binding-induced conformational changes of the assembly components. This chapter describes the advantages and potential pitfalls of these methods and puts them into perspective with alternative approaches such as iterative modular fitting of rigid-body domains.

Structural studies on full-length talin1 reveal a compact auto-inhibited dimer: implications for talin activation

Talin is a large adaptor protein that activates integrins and couples them to cytoskeletal actin. Talin contains an N-terminal FERM (band 4.1, ezrin, radixin, moesin) domain (the head) linked to a flexible rod comprised of 13 amphipathic helical bundles (R1-R13) that terminate in a C-terminal helix (DD) that forms an anti-parallel dimer. We derived a three-dimensional structural model of full-length talin at a resolution of approximately 2.5nm using EM reconstruction of full-length talin and the known shapes of the individual domains and inter-domain angles as derived from small angle X-ray scattering. Talin adopts a compact conformation consistent with a dimer in which the two talin rods form a donut-shaped structure, with the two talin heads packed side by side occupying the hole at the center of this donut. In this configuration, the integrin binding site in the head domain and the actin-binding site at the carboxy-terminus of the rod are masked, implying that talin must unravel before it can support integrin activation and engage the actin cytoskeleton.