Draft crystal structure of the vault shell at 9-A resolution

Vaults are the largest known cytoplasmic ribonucleoprotein structures and may function in innate immunity. The vault shell self-assembles from 96 copies of major vault protein and encapsulates two other proteins and a small RNA. We crystallized rat liver vaults and several recombinant vaults, all among the largest non-icosahedral particles to have been crystallized. The best crystals thus far were formed from empty vaults built from a cysteine-tag construct of major vault protein (termed cpMVP vaults), diffracting to about 9-A resolution. The asymmetric unit contains a half vault of molecular mass 4.65 MDa. X-ray phasing was initiated by molecular replacement, using density from cryo-electron microscopy (cryo-EM). Phases were improved by density modification, including concentric 24- and 48-fold rotational symmetry averaging. From this, the continuous cryo-EM electron density separated into domain-like blocks. A draft atomic model of cpMVP was fit to this improved density from 15 domain models. Three domains were adapted from a nuclear magnetic resonance substructure. Nine domain models originated in ab initio tertiary structure prediction. Three C-terminal domains were built by fitting poly-alanine to the electron density. Locations of loops in this model provide sites to test vault functions and to exploit vaults as nanocapsules.     

Draft Crystal Structure of the Vault Shell at 9-Å Resolution

Vaults are the largest known cytoplasmic ribonucleoprotein structures and may function in innate immunity. The vault shell self-assembles from 96 copies of major vault protein and encapsulates two other proteins and a small RNA. We crystallized rat liver vaults and several recombinant vaults, all among the largest non-icosahedral particles to have been crystallized. The best crystals thus far were formed from empty vaults built from a cysteine-tag construct of major vault protein (termed cpMVP vaults), diffracting to about 9-Å resolution. The asymmetric unit contains a half vault of molecular mass 4.65 MDa. X-ray phasing was initiated by molecular replacement, using density from cryo-electron microscopy (cryo-EM). Phases were improved by density modification, including concentric 24- and 48-fold rotational symmetry averaging. From this, the continuous cryo-EM electron density separated into domain-like blocks. A draft atomic model of cpMVP was fit to this improved density from 15 domain models. Three domains were adapted from a nuclear magnetic resonance substructure. Nine domain models originated in ab initio tertiary structure prediction. Three C-terminal domains were built by fitting poly-alanine to the electron density. Locations of loops in this model provide sites to test vault functions and to exploit vaults as nanocapsules.     

Draft Crystal Structure of the Vault Shell at 9-Å Resolution

Vaults are the largest known cytoplasmic ribonucleoprotein structures and may function in innate immunity. The vault shell self-assembles from 96 copies of major vault protein and encapsulates two other proteins and a small RNA. We crystallized rat liver vaults and several recombinant vaults, all among the largest non-icosahedral particles to have been crystallized. The best crystals thus far were formed from empty vaults built from a cysteine-tag construct of major vault protein (termed cpMVP vaults), diffracting to about 9-Å resolution. The asymmetric unit contains a half vault of molecular mass 4.65 MDa. X-ray phasing was initiated by molecular replacement, using density from cryo-electron microscopy (cryo-EM). Phases were improved by density modification, including concentric 24- and 48-fold rotational symmetry averaging. From this, the continuous cryo-EM electron density separated into domain-like blocks. A draft atomic model of cpMVP was fit to this improved density from 15 domain models. Three domains were adapted from a nuclear magnetic resonance substructure. Nine domain models originated in ab initio tertiary structure prediction. Three C-terminal domains were built by fitting poly-alanine to the electron density. Locations of loops in this model provide sites to test vault functions and to exploit vaults as nanocapsules.     

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.     

FitEM2EM--tools for low resolution study of macromolecular assembly and dynamics

Studies of the structure and dynamics of macromolecular assemblies often involve comparison of low resolution models obtained using different techniques such as electron microscopy or atomic force microscopy. We present new computational tools for comparing (matching) and docking of low resolution structures, based on shape complementarity. The matched or docked objects are represented by three dimensional grids where the value of each grid point depends on its position with regard to the interior, surface or exterior of the object. The grids are correlated using fast Fourier transformations producing either matches of related objects or docking models depending on the details of the grid representations. The procedures incorporate thickening and smoothing of the surfaces of the objects which effectively compensates for differences in the resolution of the matched/docked objects, circumventing the need for resolution modification. The presented matching tool FitEM2EMin successfully fitted electron microscopy structures obtained at different resolutions, different conformers of the same structure and partial structures, ranking correct matches at the top in every case. The differences between the grid representations of the matched objects can be used to study conformation differences or to characterize the size and shape of substructures. The presented low-to-low docking tool FitEM2EMout ranked the expected models at the top.     

Low-resolution structure refinement in electron microscopy

A real-space structure refinement method, originally developed for macromolecular X-ray crystallography, has been applied to protein structure analysis by electron microscopy (EM). This method simultaneously optimizes the fit of an atomic model to a density map and the stereo-chemical properties of the model by minimizing an energy function. The performance of this method is characterized at different resolution and signal-to-noise ratio conditions typical for EM electron density maps. A multi-resolution scheme is devised to improve the convergence of the refinement on the global energy minimum. Applications of the method to various model systems are demonstrated here. The first case is the arrangement of FlgE molecules in the helical filament of flagellar hook, in which refinement with segmented rigid bodies improves the density correlation and reduces severe van der Waals contacts among the symmetry-related subunits. The second case is a conformational analysis of the NSF AAA ATPase in which a multi-conformer model is used in the refinement to investigate the arrangement of the two ATPase domains in the molecule. The third case is a docking simulation in which the crystal structure of actin and the NOE data from NMR experiments on the dematin headpiece are combined with a low-resolution EM density map to generate an atomic model of the F-actin-dematin headpiece structure.     

Exploring global distortions of biological macromolecules and assemblies from low-resolution structural information and elastic network theory

A theory of elastic normal modes is described for the exploration of global distortions of biological structures and their assemblies based upon low-resolution image data. Structural information at low resolution, e.g. from density maps measured by cryogenic electron microscopy (cryo-EM), is used to construct discrete multi-resolution models for the electron density using the techniques of vector quantization. The elastic normal modes computed based on these discretized low-resolution models are found to compare well with the normal modes obtained at atomic resolution. The quality of the normal modes describing global displacements of the molecular system is found to depend on the resolution of the synthetic EM data and the extent of reductionism in the discretized representation. However, models that reproduce the functional rearrangements of our test set of molecules are achieved for realistic values of experimental resolution. Thus large conformational changes as occur during the functioning of biological macromolecules and assemblies can be elucidated directly from low-resolution structural data through the application of elastic normal mode theory and vector quantization.     

NMMD: Efficient Cryo-EM Flexible Fitting Based on Simultaneous Normal Mode and Molecular Dynamics atomic displacements

Atomic models of cryo electron microscopy (cryo-EM) maps of biomolecular conformations are often obtained by flexible fitting of the maps with available atomic structures of other conformations (e.g., obtained by X-ray crystallography). This article presents a new flexible fitting method, NMMD, which combines normal mode analysis (NMA) and molecular dynamics simulation (MD). Given an atomic structure and a cryo-EM map to fit, NMMD simultaneously estimates global atomic displacements based on NMA and local displacements based on MD. NMMD was implemented by modifying EMfit, a flexible fitting method using MD only, in GENESIS 1.4. As EMfit, NMMD can be run with replica exchange umbrella sampling procedure. The new method was tested using a variety of EM maps (synthetic and experimental, with different noise levels and resolutions). The results of the tests show that adding normal modes to MD-based fitting makes the fitting faster (40% in average) and, in the majority of cases, more accurate.     

Using situs for flexible and rigid-body fitting of multiresolution single-molecule data

We describe here a set of multiresolution visualization and docking procedures that we refer to as the Situs package. The package was developed to provide an efficient and robust method for the fitting of atomic structures into low-resolution data. The current release was optimized specifically for the visualization and docking of single molecules. A novel 3D graphics viewer, volslice3d, was developed for the package to provide an immersive virtual reality environment for measuring and rendering volumetric data sets. The precision of single-molecule, rigid-body docking was tested on simulated (noise-free) low-resolution density maps. For spatial resolutions near 20 A typically arising in electron microscopy image reconstructions, a docking precision on the order of 1 A can be achieved. The shape-matching score captured the correct solutions in all 10 trial cases and was sufficiently stringent to yield unique matches in 8 systems. Novel routines were developed for the flexible docking of atomic structures whose shape deviates from the corresponding low-resolution shape. Test calculations on isoforms of actin and lactoferrin demonstrate that the flexible docking faithfully reproduces conformational differences with a precision < 2 A if atomic structures are locally conserved.     

The architecture of CopA from Archeaoglobus fulgidus studied by cryo-electron microscopy and computational docking

CopA uses ATP to pump Cu(+) across cell membranes. X-ray crystallography has defined atomic structures of several related P-type ATPases. We have determined a structure of CopA at 10?? resolution by cryo-electron microscopy of a new crystal form and used computational molecular docking to study the interactions between the N-terminal metal-binding domain (NMBD) and other elements of the molecule. We found that the shorter-chain lipids used to produce these crystals are associated with movements of the cytoplasmic domains, with a novel dimer interface and with disordering of the NMBD, thus offering evidence for the transience of its interaction with the other cytoplasmic domains. Docking identified a binding site that matched the location of the NMBD in our previous structure by cryo-electron microscopy, allowing a more detailed view of its binding configuration and further support for its role in autoinhibition.     

A computational comparison of the atomic models of the actomyosin interface

Several atomic models of the actomyosin interface have been proposed based on the docking together of their component structures using electron microscopy and resonance energy-transfer measurements. Although these models are in approximate agreement in the location of the binding interfaces when myosin is tightly bound to actin, their relationships to molecular docking simulations based on computational free-energy calculations are investigated here. Both rigid-docking and flexible-docking conformational search strategies were used to identify free-energy minima at the interfaces between atomic models of myosin and actin. These results suggest that the docking model produced by resonance energy-transfer data is closer to a free-energy minimum at the interface than are the available atomic models based on electron microscopy. The conformational searches were performed using both scallop and chicken skeletal muscle myosins and identified similarly oriented actin-binding interfaces that serve to validate that these models are at the global minimum. These results indicate that the existing docking models are close to but not precisely at the lowest-energy initial contact site for strong binding between myosin and actin that should represent an initial contact between the two proteins; therefore, conformational changes are likely to be important during the transition to a strongly bound complex.     

Three-dimensional structure of a voltage-gated potassium channel at 2.5 nm resolution

BACKGROUND: The voltage-gated potassium channel Shaker from Drosophila consists of a tetramer of identical subunits, each containing six transmembrane segments. The atomic structure of a bacterial homolog, the potassium channel KcsA, is much smaller than Shaker. It does not have a voltage sensor and other important domains like the N-terminal tetramerization (T1) domain. The structure of these additional elements has to be studied in the more complex voltage-gated channels. RESULTS: We determined the three-dimensional structure of the entire Shaker channel at 2.5 nm resolution using electron microscopy. The four-fold symmetric structure shows a large and a small domain linked by thin 2 nm long connectors. To interpret the structure, we used the crystal structures of the isolated T1 domain and the KcsA channel. A unique density assignment was made based on the symmetry and dimensions of the crystal structures and domains, identifying the smaller domain as the cytoplasmic mass of Shaker containing T1 and the larger domain as embedded in the membrane. CONCLUSIONS: The two-domain architecture of the Shaker channel is consistent with the recently proposed "hanging gondola" model for the T1 domain, putting the T1 domain at a distance from the membrane domain but attached to it by thin connectors. The space between the two domains is sufficient to permit cytoplasmic access of ions and the N-terminal inactivation domain to the pore region. A hanging gondola architecture has also been observed in the nicotinic acetylcholine receptor and the KcsA structure, suggesting that it is a common element of ion channels.     

Normal mode based flexible fitting of high-resolution structure into low-resolution experimental data from cryo-EM

A new method for the flexible fitting of high-resolution structures into low-resolution maps of macromolecular complexes from electron microscopy has been recently described in applications to simulated electron density maps. This method uses a linear combination of low-frequency normal modes in an iterative manner to deform the structure optimally to conform to the low-resolution electron density map. Gradient-following techniques in the coordinate space of collective normal modes are used to optimize the overall correlation coefficient between computed and measured electron densities. With this approach, multi-scale flexible fitting can be performed using all-atoms or Calpha atoms. In this paper, illustrative studies of normal mode based flexible fitting to experimental cryo-EM maps are presented for three different systems. Large, functionally relevant conformational changes for elongation factor G bound to the ribosome, Escherichia coli RNA polymerase and cowpea chlorotic mottle virus are elucidated as the result of the application of NMFF from high-resolution structures to cryo-electron microscopy maps.     

Single particle image reconstruction of the human recombinant Kv2.1 channel

Kv2.1 channels are widely expressed in neuronal and endocrine cells and generate slowly activating K⁺ currents, which contribute to repolarization in these cells. Kv2.1 is expressed at high levels in the mammalian brain and is a major component of the delayed rectifier current in the hippocampus. In addition, Kv2.1 channels have been implicated in the regulation of membrane repolarization, cytoplasmic calcium levels, and insulin secretion in pancreatic β-cells. They are therefore an important drug target for the treatment of Type II diabetes mellitus. We used electron microscopy and single particle image analysis to derive a three-dimensional density map of recombinant human Kv2.1. The tetrameric channel is egg-shaped with a diameter of ∼80 Å and a long axis of ∼120 Å. Comparison to known crystal structures of homologous domains allowed us to infer the location of the cytoplasmic and transmembrane assemblies. There is a very good fit of the Kv1.2 crystal structure to the assigned transmembrane assembly of Kv2.1. In other low-resolution maps of K⁺ channels, the cytoplasmic N-terminal and transmembrane domains form separate rings of density. In contrast, Kv2.1 displays contiguous density that connects the rings, such that there are no large windows between the channel interior and the cytoplasmic space. The crystal structure of KcsA is thought to be in a closed conformation, and the good fit of the KcsA crystal structure to the Kv2.1 map suggests that our preparations of Kv2.1 may also represent a closed conformation. Substantial cytoplasmic density is closely associated with the T1 tetramerization domain and is ascribed to the ∼184 kDa C-terminal regulatory domains within each tetramer.     

Building the atomic model for the bacterial flagellar filament by electron cryomicroscopy and image analysis

The bacterial flagellar filament is a helical propeller for bacterial locomotion. It is a well-ordered helical assembly of a single protein, flagellin, and its tubular structure is formed by 11 protofilaments, each in either of the two distinct conformations, L- and R-type, for supercoiling. We have been studying the three-dimensional structures of the flagellar filaments by electron cryomicroscopy and recently obtained a density map of the R-type filament up to 4 angstroms resolution from an image data set containing only about 41,000 molecular images. The density map showed the features of the alpha-helical backbone and some large side chains, which allowed us to build the complete atomic model as one of the first atomic models of macromolecules obtained solely by electron microscopy image analysis (Yonekura et al., 2003a). We briefly review the structure and the structure analysis, and point out essential techniques that have made this analysis possible.     

Fingerprint-based structure retrieval using electron density

We present a computational approach that can quickly search a large protein structural database to identify structures that fit a given electron density, such as determined by cryo-electron microscopy. We use geometric invariants (fingerprints) constructed using 3D Zernike moments to describe the electron density, and reduce the problem of fitting of the structure to the electron density to simple fingerprint comparison. Using this approach, we are able to screen the entire Protein Data Bank and identify structures that fit two experimental electron densities determined by cryo-electron microscopy.     

Structural characterization of components of protein assemblies by comparative modeling and electron cryo-microscopy

We explore structural characterization of protein assemblies by a combination of electron cryo-microscopy (cryoEM) and comparative protein structure modeling. Specifically, our method finds an optimal atomic model of a given assembly subunit and its position within an assembly by fitting alternative comparative models into a cryoEM map. The alternative models are calculated by MODELLER [J. Mol. Biol. 234 (1993) 313] from different sequence alignments between the modeled protein and its template structures. The fitting of these models into a cryoEM density map is performed either by FOLDHUNTER [J. Mol. Biol. 308 (2001) 1033] or by a new density fitting module of MODELLER (Mod-EM). Identification of the most accurate model is based on the correlation between the model accuracy and the quality of fit into the cryoEM density map. To quantify this correlation, we created a benchmark consisting of eight proteins of different structural folds with corresponding density maps simulated at five resolutions from 5 to 15 angstroms, with three noise levels each. Each of the proteins in the set was modeled based on 300 different alignments to their remotely related templates (12-32% sequence identity), spanning the range from entirely inaccurate to essentially accurate alignments. The benchmark revealed that one of the most accurate models can usually be identified by the quality of its fit into the cryoEM density map, even for noisy maps at 15 angstroms resolution. Therefore, a cryoEM density map can be helpful in improving the accuracy of a comparative model. Moreover, a pseudo-atomic model of a component in an assembly may be built better with comparative models of the native subunit sequences than with experimentally determined structures of their homologs.     

Real space refinement of acto-myosin structures from sectioned muscle

We have adapted a real space refinement protocol originally developed for high-resolution crystallographic analysis for use in fitting atomic models of actin filaments and myosin subfragment 1 (S1) to 3-D images of thin-sectioned, plastic-embedded whole muscle. The rationale for this effort is to obtain a refinement protocol that will optimize the fit of the model to the density obtained by electron microscopy and correct for poor geometry introduced during the manual fitting of a high-resolution atomic model into a lower resolution 3-D image. The starting atomic model consisted of a rigor acto-S1 model obtained by X-ray crystallography and helical reconstruction of electron micrographs. This model was rebuilt to fit 3-D images of rigor insect flight muscle at a resolution of 7 nm obtained by electron tomography and image averaging. Our highly constrained real space refinement resulted in modest improvements in the agreement of model and reconstruction but reduced the number of conflicting atomic contacts by 70% without loss of fit to the 3-D density. The methodology seems to be well suited to the derivation of stereochemically reasonable atomic models that are consistent with experimentally determined 3-D reconstructions computed from electron micrographs.     

The three-dimensional structure of P2 myelin protein.

The three-dimensional structure of P2 protein from peripheral nervous system myelin has been determined at 2.7 A resolution by X-ray crystallography. The single isomorphous replacement/anomalous map was interpreted using skeletonized electron density on a computer graphics system. An atomic model was built using fragment fitting. The structure forms a compact 10-stranded up-and-down beta-barrel which encapsulates residual electron density that we interpret as a fatty acid molecule. This beta-barrel shows some similarity to, but is different from, the retinol binding protein family of structures. The relationship of the P2 structure to a family of cytoplasmic, lipid binding proteins is described.