Normal mode-based fitting of atomic structure into electron density maps: application to sarcoplasmic reticulum Ca-ATPase

A method for the flexible docking of high-resolution atomic structures into lower resolution densities derived from electron microscopy is presented. The atomic structure is deformed by an iterative process using combinations of normal modes to obtain the best fit of the electron microscopical density. The quality of the computed structures has been evaluated by several techniques borrowed from crystallography. Two atomic structures of the SERCA1 Ca-ATPase corresponding to different conformations were used as a starting point to fit the electron density corresponding to a different conformation. The fitted models have been compared to published models obtained by rigid domain docking, and their relation to the known crystallographic structures are explored by normal mode analysis. We find that only a few number of modes contribute significantly to the transition. The associated motions involve almost exclusively rotation and translation of the cytoplasmic domains as well as displacement of cytoplasmic loops. We suggest that the movements of the cytoplasmic domains are driven by the conformational change that occurs between nonphosphorylated and phosphorylated intermediate, the latter being mimicked by the presence of vanadate at the phosphorylation site in the electron microscopy structure.     

Flexible fitting of atomic structures into electron microscopy maps using molecular dynamics

A novel method to flexibly fit atomic structures into electron microscopy (EM) maps using molecular dynamics simulations is presented. The simulations incorporate the EM data as an external potential added to the molecular dynamics force field, allowing all internal features present in the EM map to be used in the fitting process, while the model remains fully flexible and stereochemically correct. The molecular dynamics flexible fitting (MDFF) method is validated for available crystal structures of protein and RNA in different conformations; measures to assess and monitor the fitting process are introduced. The MDFF method is then used to obtain high-resolution structures of the E. coli ribosome in different functional states imaged by cryo-EM.     

YUP.SCX: coaxing atomic models into medium resolution electron density maps

The structures of large macromolecular complexes in different functional states can be determined by cryo-electron microscopy, which yields electron density maps of low to intermediate resolutions. The maps can be combined with high-resolution atomic structures of components of the complex, to produce a model for the complex that is more accurate than the formal resolution of the map. To this end, methods have been developed to dock atomic models into density maps rigidly or flexibly, and to refine a docked model so as to optimize the fit of the atomic model into the map. We have developed a new refinement method called YUP.SCX. The electron density map is converted into a component of the potential energy function to which terms for stereochemical restraints and volume exclusion are added. The potential energy function is then minimized (using simulated annealing) to yield a stereochemically-restrained atomic structure that fits into the electron density map optimally. We used this procedure to construct an atomic model of the 70S ribosome in the pre-accommodation state. Although some atoms are displaced by as much as 33 Å, they divide themselves into nearly rigid fragments along natural boundaries with smooth transitions between the fragments.     

Efficient automatic noise reduction of electron tomographic reconstructions based on iterative median filtering

A simple, fast and efficient noise-reduction protocol for three-dimensional electron tomographic reconstructions of biological material is presented. The approach is based on iterative application of median filtering and shows promise for automatic noise reduction as a pre-processor for automated data analysis tools which aim at segmentation, feature extraction and pattern recognition. The application of this algorithm produces encouraging results for a wide variety of experimental and synthetic electron tomographic reconstructions.     

Fast fitting of atomic structures to low-resolution electron density maps by surface overlap maximization

The complexities of X-ray crystallography and NMR spectroscopy for large protein complexes, and the comparative ease of approaches such as electron microscopy mean that low-resolution structures are often available long before their atomic resolution equivalents. To help bridge this gap in knowledge, we present 3SOM: an approach for finding the best fit of atomic resolution structures into lower-resolution density maps through surface overlap maximization. High-resolution templates (i.e. partial structures or models for multi-subunit complexes) and targets (lower-resolution maps) are initially represented as iso-surfaces. The latter are used first in a fast search for transformations that superimpose a significant portion of the target surface onto the template surface, which is quantified as surface overlap. The vast search space is reduced by considering key vectors that capture local surface information. The set of transformations with the highest surface overlap scores are then re-ranked by using more sophisticated scores including cross-correlation. We give a number of examples to illustrate the efficiency of the method and its restrictions. For targets for which partial complexes are available, the speed and performance of the method make it an attractive complement to existing methods, as many different hypotheses can be tested quickly on a single processor.     

Putting pharmaceuticals into the wider context of challenges to fish populations in rivers

The natural range of fish species in our rivers is related to flow, elevation, temperature, local habitat and connectivity. For over 2000 years, humans have altered to varying degrees the river habitat. In the past 200 years, we added to the environmental disruption by discharging poorly treated sewage, nutrients and industrial waste into our rivers. For many rivers, the low point arrived during the period of 1950s–1970s, when rapid economic development overrode environmental concerns and dissolved oxygen concentrations dropped to zero. In these more enlightened times, gross river pollution is a thing of the past in the Developed World. However, persistent legacy chemical contaminants can be found in fish long after their discharge ceased. Changes in habitat quality and morphology caused and continue to cause the disappearance of fish species. The range of fish stressors has now increased as temperatures rise, and non-native fish introductions bring new diseases. The threat from pharmaceuticals to fish populations remains hypothetical, and no studies have yet linked change in fish populations to exposure.     

The structure of fibrin. An electron microscopic investigation

An electron microscopic investigation on the structure of fibrin is reported. Fibrin morphology was investigated in a wide variety of experimental conditions, and by carefully controlled staining procedures. Two main band patterns A (230-Å-spaced main dark bands) and B (230-Å-spaced main light bands) are observed for stained fibrin; it is shown that the former results from the superimposition of both positive and negative staining, and the latter is given by positive staining. By suitable denaturation experiments, it was found that the fiber is composed of a regular alternation of lose and dense regions along the axis. We have assumed that the monomer of fibrin is described by the three-nodular model of Hall and Slayter, as supported by recent investigations. The monomers are arrayed according to a head-to-tail sequence along the fiber, and to a staggered lateral association. This model accounts for all the experimental observations, and predicts well the high-resolution band pattern of fibrin. It further agrees with the results of a recent work on the early stages of the fibrinogen-to-fibrin conversion.     

Polysome structure in sea urchin eggs and embryos: an electron microscopic analysis

Primary cultures of cardiac myocytes from newborn normal and genetically cardiomyopathic (strain UM-X7.1) hamsters were analyzed by electron microscopy and immunofluorescent staining for myosin, actin, tropomyosin, and alpha-actinin. Antibody staining of these contractile proteins demonstrates that both normal and cardiomyopathic (CM) myocytes contain prominent myofibrils after 3 days in culture, although the CM myofibrils are disarrayed and not aligned as those in normal cells. The disarray becomes even more pronounced in CM cells after 5 days in culture. The immunofluorescent staining patterns of individual myofibrils in normal and CM cells were similar for myosin, actin, and tropomyosin. However, alpha-actinin staining reveals that the CM myofibrils have abnormally wide and irregularly shaped Z bands. Electron microscopy confirms the irregular Z-band appearance as well as the myofibril disarray. Thus, CM cardiomyocytes clearly show an aberrant pattern of myofibril structure and organization in culture.     

A core-weighted fitting method for docking atomic structures into low-resolution maps: application to cryo-electron microscopy

Cryo-electron microscopy of "single particles" is a powerful method to analyze structures of large macromolecular assemblies that are not amenable to investigation by traditional X-ray crystallographic methods. A key step in these studies is to obtain atomic interpretations of multiprotein complexes by fitting atomic structures of individual components into maps obtained from electron microscopic data. Here, we report the use of a "core-weighting" method, combined with a grid-threading Monte Carlo (GTMC) approach for this purpose. The "core" of an individual structure is defined to represent the part where the density distribution is least likely to be altered by other components that comprise the macromolecular assembly of interest. The performance of the method has been evaluated by its ability to determine the correct fit of (i) the alpha-chain of the T-cell receptor variable domain into a simulated map of the alphabeta complex at resolutions between 5 and 40 A, and (ii) the E2 catalytic domain of the pyruvate dehydrogenase into an experimentally determined map, at 14 A resolution, of the icosahedral complex formed by 60 copies of this enzyme. Using the X-ray structures of the two test cases as references, we demonstrate that, in contrast to more traditional methods, the combination of the core-weighting method and the grid-threading Monte Carlo approach can identify the correct fit reliably and rapidly from the low-resolution maps that are typical of structures determined with the use of single-particle electron microscopy.     

UCSF tomography: an integrated software suite for real-time electron microscopic tomographic data collection, alignment, and reconstruction

A real-time alignment and reconstruction scheme for electron microscopic tomography (EMT) has been developed and integrated within our UCSF tomography data collection software. This newly integrated software suite provides full automation from data collection to real-time reconstruction by which the three-dimensional (3D) reconstructed volume is immediately made available at the end of each data collection. Real-time reconstruction is achieved by calculating a weighted back-projection on a small Linux cluster (five dual-processor compute nodes) concurrently with the UCSF tomography data collection running on the microscope's computer, and using the fiducial-marker free alignment data generated during the data collection process. The real-time reconstructed 3D volume provides users with immediate feedback to fully asses all aspects of the experiment ranging from sample choice, ice thickness, experimental parameters to the quality of specimen preparation. This information can be used to guide subsequent data collections. Access to the reconstruction is especially useful in low-dose cryo EMT where such information is very difficult to obtain due to extraordinary low signal to noise ratio in each 2D image. In our environment, we generally collect 2048 x 2048 pixel images which are subsequently computationally binned four-fold for the on-line reconstruction. Based upon experiments performed with thick and cryo specimens at various CCD magnifications (50000x-80000x), alignment accuracy is sufficient to support this reduced resolution but should be refined before calculating a full resolution reconstruction. The reduced resolution has proven to be quite adequate to assess sample quality, or to screen for the best data set for full-resolution reconstruction, significantly improving both productivity and efficiency of system resources. The total time from start of data collection to a final reconstructed volume (512 x 512 x 256 pixels) is about 50 min for a +/-70 degrees 2k x 2k pixel tilt series acquired at every 1 degrees.     

An emerging view of presynaptic structure from electron microscopic studies

In response to calcium influx, some of the synaptic vesicles in presynaptic terminals fuse rapidly with the presynaptic membrane, allowing fast synaptic transmission. The regulated recycling of synaptic vesicles at the terminals is required for a sustained release of neurotransmitters. Localization of 'ready to be released' vesicles in close vicinities to voltage-gated calcium channels enables the rapid release of neurotransmitters. Thus, recycling vesicles must translocate from the sites of endocytosis to these release sites. However, the sub-cellular organization that supports this local vesicular traffic remains poorly understood. We will review the results of various electron microscopy studies, which have begun to unveil the structure of presynaptic terminals.     

BCL::EM-Fit: rigid body fitting of atomic structures into density maps using geometric hashing and real space refinement

Cryo-electron microscopy (cryoEM) can visualize large macromolecular assemblies at resolutions often below 10? and recently as good as 3.8-4.5 ?. These density maps provide important insights into the biological functioning of molecular machineries such as viruses or the ribosome, in particular if atomic-resolution crystal structures or models of individual components of the assembly can be placed into the density map. The present work introduces a novel algorithm termed BCL::EM-Fit that accurately fits atomic-detail structural models into medium resolution density maps. In an initial step, a "geometric hashing" algorithm provides a short list of likely placements. In a follow up Monte Carlo/Metropolis refinement step, the initial placements are optimized by their cross correlation coefficient. The resolution of density maps for a reliable fit was determined to be 10 ? or better using tests with simulated density maps. The algorithm was applied to fitting of capsid proteins into an experimental cryoEM density map of human adenovirus at a resolution of 6.8 and 9.0 ?, and fitting of the GroEL protein at 5.4 ?. In the process, the handedness of the cryoEM density map was unambiguously identified. The BCL::EM-Fit algorithm offers an alternative to the established Fourier/Real space fitting programs. BCL::EM-Fit is free for academic use and available from a web server or as downloadable binary file at http://www.meilerlab.org.     

Probing the structure of F-actin: cross-links constrain atomic models and modify actin dynamics.

Cross-links between protomers in F-actin can be used as a very sensitive probe of both the dynamics and structure of F-actin. We have characterized filaments formed from a previously described yeast actin Q41C mutant, where disulfide bonds can be formed between the Cys41 that is introduced into subdomain-2 and Cys374 on an adjacent protomer. We find that the distribution of cross-linked n-mers shows no cooperativity and corresponds to a random probability cross-linking reaction. The random distribution suggests that disulfide formation does not cause a significant perturbation of the F-actin structure. Consistent with this lack of perturbation, three-dimensional reconstructions of extensively cross-linked filaments, using a new approach to helical image analysis, show very small structural changes with respect to uncross-linked filaments. This finding is in conflict with refined models but in agreement with the original Holmes et al. model for F-actin. Under conditions where 94 % of the protomers are linked by disulfide bonds, the distribution of filament twist becomes more heterogeneous with respect to control filaments. A molecular model suggests that strain, introduced by the disulfide, is relieved by increasing the twist of the long-pitch actin helices. Disulfide formation makes yeast actin filaments approximately three times less flexible in terms of bending and similar, in this respect, to vertebrate skeletal muscle F-actin. These observations support previous reports that the rigidity of F-actin can be controlled by the position of subdomain-2, and that this region is more flexible in yeast F-actin than in skeletal muscle F-actin.