2dx_merge: data management and merging for 2D crystal images

Electron crystallography of membrane proteins determines the structure of membrane-reconstituted and two-dimensionally (2D) crystallized membrane proteins by low-dose imaging with the transmission electron microscope, and computer image processing. We have previously presented the software system 2dx, for user-friendly image processing of 2D crystal images. Its central component 2dx_image is based on the MRC program suite, and allows the optionally fully automatic processing of one 2D crystal image. We present here the program 2dx_merge, which assists the user in the management of a 2D crystal image processing project, and facilitates the merging of the data from multiple images. The merged dataset can be used as a reference to re-process all images, which usually improves the resolution of the final reconstruction. Image processing and merging can be applied iteratively, until convergence is reached. 2dx is available under the GNU General Public License at http://2dx.org.     

Single particle 3D reconstruction for 2D crystal images of membrane proteins

In cases where ultra-flat cryo-preparations of well-ordered two-dimensional (2D) crystals are available, electron crystallography is a powerful method for the determination of the high-resolution structures of membrane and soluble proteins. However, crystal unbending and Fourier-filtering methods in electron crystallography three-dimensional (3D) image processing are generally limited in their performance for 2D crystals that are badly ordered or non-flat. Here we present a single particle image processing approach, which is implemented as an extension of the 2D crystallographic pipeline realized in the 2dx software package, for the determination of high-resolution 3D structures of membrane proteins. The algorithm presented, addresses the low single-to-noise ratio (SNR) of 2D crystal images by exploiting neighborhood correlation between adjacent proteins in the 2D crystal. Compared with conventional single particle processing for randomly oriented particles, the computational costs are greatly reduced due to the crystal-induced limited search space, which allows a much finer search space compared to classical single particle processing. To reduce the considerable computational costs, our software features a hybrid parallelization scheme for multi-CPU clusters and computer with high-end graphic processing units (GPUs). We successfully apply the new refinement method to the structure of the potassium channel MloK1. The calculated 3D reconstruction shows more structural details and contains less noise than the map obtained by conventional Fourier-filtering based processing of the same 2D crystal images.     

A maximum likelihood approach to two-dimensional crystals

Maximum likelihood (ML) processing of transmission electron microscopy images of protein particles can produce reconstructions of superior resolution due to a reduced reference bias. We have investigated a ML processing approach to images centered on the unit cells of two-dimensional (2D) crystal images. The implemented software makes use of the predictive lattice node tracking in the MRC software, which is used to window particle stacks. These are then noise-whitened and subjected to ML processing. Resulting ML maps are translated into amplitudes and phases for further processing within the 2dx software package. Compared with ML processing for randomly oriented single particles, the required computational costs are greatly reduced as the 2D crystals restrict the parameter search space. The software was applied to images of negatively stained or frozen hydrated 2D crystals of different crystal order. We find that the ML algorithm is not free from reference bias, even though its sensitivity to noise correlation is lower than for pure cross-correlation alignment. Compared with crystallographic processing, the newly developed software yields better resolution for 2D crystal images of lower crystal quality, and it performs equally well for well-ordered crystal images.     

Automatic lattice determination for two-dimensional crystal images

Electron crystallography determines the structure of membrane proteins and other periodic samples by recording either images or diffraction patterns. Computer processing of recorded images requires the determination of the reciprocal lattice parameters in the Fourier transform of the image. We have developed a set of three programs 2dx_peaksearch, 2dx_findlat and 2dx_getlat, which can determine the reciprocal lattice from a Fourier transformation of a 2D crystal image automatically. 2dx_peaksearch determines a list of Fourier peak coordinates from a processed calculated diffraction pattern. These coordinates are evaluated by 2dx_findlat to determine one or more lattices, using a-priori knowledge of the real-space crystal unit cell dimensions, and the sample tilt geometry. If these are unknown, then the program 2dx_getlat can be used to obtain a guess for the unit cell dimensions. These programs are available as part of the 2dx software package for the image processing of 2D crystal images at http://2dx.org.     

Symmetry: A guide to its application in 2D electron crystallography

A defining property of a crystal is its symmetry. This mini-review sets out to summarize all aspects that define 2D crystallographic symmetry as applied to the study of macromolecular structure. It begins by defining molecular point symmetries, before covering crystallographic symmetry operations in 2D, common notation, a summary of crystallographic plane groups and theoretical methods and important considerations for the identification and application of symmetry in 2D crystal images for 3D structure determination. While many of the concepts covered here may be equally applicable to point symmetry and space group symmetry in 3D, this review has been written from the perspective of 2D electron crystallography and deals specifically with symmetry operations and crystallographic space groups in 2D crystal projection images.     

Surface topography of the p3 and p6 annexin V crystal forms determined by atomic force microscopy

Annexin V is a member of a family of structurally homologous proteins sharing the ability to bind to negatively charged phospholipid membranes in a Ca(2+)-dependent manner. The structure of the soluble form of annexin V has been solved by X-ray crystallography, while electron crystallography of two-dimensional (2D) crystals has been used to reveal the structure of its membrane-bound form. Two 2D crystal forms of annexin V have been reported to date, with either p6 or p3 symmetry. Atomic force microscopy has previously been used to investigate the growth and the topography of the p6 crystal form on supported phospholipid bilayers (Reviakine et al., 1998). The surface structure of the second crystal form, p3, is presented in this study, along with an improved topographic map of the p6 crystal form. The observed topography is correlated with the structure determined by X-ray crystallography.     

Electron crystallography of macromolecular periodic arrays on phospholipid monolayers

Electron crystallography has the potential of yielding structural information equivalent to x-ray diffraction. The major difficulty has been preparing specimens with the required structural order and size for diffraction and imaging in the electron microscope. 2D crystallization on phospholipid monolayers is capable of fulfilling both of these requirements. Crystals can form as a result of specific interactions with a protein's ligand or an analog, suitably linked to a lipid tail; or on a surface of complementary head-group charge. With such choices, the availability of a suitable lipid is limited only by synthetic chemistry. Ultimately, it is the quality and regularity of the protein-protein interactions that determine the crystalline order, as it is with any protein crystal. In the case of streptavidin, the monolayer crystal diffracts beyond 2.5 Å. A 3 Å projection map reconstructed from electron diffraction amplitudes and phases from images shows density which can be interpreted as β-sheets and clusters of side chains. It remains to be shown that the monolayer crystals are flat and diffract as well at high tilt angle as untilted. Technological issues such as charging must be resolved. With parallel advances in data collection and processing, electron crystallography of monolayer macromolecular crystals will eventually take its place beside x-ray crystallography and NMR as a routine and efficient structural technique.     

Two-dimensional crystallization of the ryanodine receptor Ca2+ release channel on lipid membranes

The ryanodine receptor (RyR) is the largest known membrane protein with a total molecular mass of 2.3 x 10(3) kDa. Well ordered, two-dimensional (2D) crystals are an essential prerequisite to enable RyR structure determination by electron crystallography. Conventionally, the 2D crystallization of membrane proteins is based on a 'trial-and-error' strategy, which is both time-consuming and chance-directed. By adopting a new strategy that utilizes protein sequence information and predicted transmembrane topology, we successfully crystallized the RyR on positively charged lipid membranes. Image processing of negatively stained crystals reveals that they are well ordered, with diffraction spots of IQ < or = 4 extending to approximately 20 angstroms, the resolution attainable in negative stain. The RyR crystals obtained on the charged lipid membrane have characteristics consistent with 2D arrays that have been observed in native sarcoplasmic reticulum of muscle tissues. These crystals provide ideal materials to enable structural analysis of RyR by high-resolution electron crystallography. Moreover, the reconstituted native-like 2D array provides an ideal model system to gain structural insights into the mechanism of RyR-mediated Ca2+ signaling processes, in which the intrinsic ability of RyR oligomers to organize into a 2D array plays a crucial role.     

Automated macromolecular model building for X-ray crystallography using ARP/wARP version 7

ARP/wARP is a software suite to build macromolecular models in X-ray crystallography electron density maps. Structural genomics initiatives and the study of complex macromolecular assemblies and membrane proteins all rely on advanced methods for 3D structure determination. ARP/wARP meets these needs by providing the tools to obtain a macromolecular model automatically, with a reproducible computational procedure. ARP/wARP 7.0 tackles several tasks: iterative protein model building including a high-level decision-making control module; fast construction of the secondary structure of a protein; building flexible loops in alternate conformations; fully automated placement of ligands, including a choice of the best-fitting ligand from a 'cocktail'; and finding ordered water molecules. All protocols are easy to handle by a nonexpert user through a graphical user interface or a command line. The time required is typically a few minutes although iterative model building may take a few hours.     

An Open-Source Mesh Generation Platform for Biophysical Modeling Using Realistic Cellular Geometries

Advances in imaging methods such as electron microscopy, tomography, and other modalities are enabling high-resolution reconstructions of cellular and organelle geometries. Such advances pave the way for using these geometries for biophysical and mathematical modeling once these data can be represented as a geometric mesh, which, when carefully conditioned, enables the discretization and solution of partial differential equations. In this work, we outline the steps for a naïve user to approach the Geometry-preserving Adaptive MeshER software version 2, a mesh generation code written in C++ designed to convert structural data sets to realistic geometric meshes while preserving the underlying shapes. We present two example cases: 1) mesh generation at the subcellular scale as informed by electron tomography and 2) meshing a protein with a structure from x-ray crystallography. We further demonstrate that the meshes generated by the Geometry-preserving Adaptive MeshER software are suitable for use with numerical methods. Together, this collection of libraries and tools simplifies the process of constructing realistic geometric meshes from structural biology data.     

A helical processing pipeline for EM structure determination of membrane proteins

Electron crystallography plays a key role in the structural biology of integral membrane proteins (IMPs) by offering one of the most direct means of providing insight into the functional state of these molecular machines in their lipid-associated forms, and also has the potential to facilitate examination of physiologically relevant transitional states and complexes. Helical or tubular crystals, which are the natural product of proteins crystallizing on the surface of a cylindrical vesicle, offer some unique advantages, such as three-dimensional (3D) information from a single view, compared to other crystalline forms. While a number of software packages are available for processing images of helical crystals to produce 3D electron density maps, widespread exploitation of helical image reconstruction is limited by a lack of standardized approaches and the initial effort and specialized expertise required. Our goal is to develop an integrated pipeline to enable structure determination by transmission electron microscopy (TEM) of IMPs in the form of tubular crystals. We describe here the integration of standard Fourier-Bessel helical analysis techniques into Appion, an integrated, database-driven pipeline.     

Organization of transmembrane helices in photosystem II: comparison of plants and cyanobacteria

Electron microscopy and X-ray crystallography are revealing the structure of photosystem II. Electron crystallography has yielded a 3D structure at sufficient resolution to identify subunit positioning and transmembrane organization of the reaction-centre core complex of spinach. Single-particle analyses are providing 3D structures of photosystem II-light-harvesting complex II supercomplexes that can be used to incorporate high-resolution structural data emerging from electron and X-ray crystallography. The positions of the chlorins and metal centres within photosystem II are now available. It can be concluded that photosystem II is a dimeric complex with the transmembrane helices of CP47/D2 proteins related to those of the CP43/D1 proteins by a twofold axis within each monomer. Further, both electron microscopy and X-ray analyses show that P(680) is not a 'special pair' and that cytochrome b559 is located on the D2 side of the reaction centres some distance from P(680). However, although comparison of the electron microscopy and X-ray models for spinach and Synechococcus elongatus show considerable similarities, there seem to be differences in the number and positioning of some small subunits.     

Projection maps of three members of the KdgM outer membrane protein family

We present the projection structures of the three outer membrane porins KdgM and KdgN from Erwinia chrysanthemi and NanC from Escherichia coli, based on 2D electron crystallography. A wide screening of 2D crystallization conditions yielded tubular crystals of a suitable size and quality to perform high-resolution electron microscopy. Data processing of untilted samples allowed us to separate the information of the two crystalline layers and resulted in projection maps to a resolution of up to 7 Å. All three proteins exhibit a similar putative β-barrel structure and the three crystal forms have the same symmetry. However, there are differences in the packing arrangements of the monomers as well as the densities of the projections. To interpret these projections, secondary structure prediction was performed using β-barrel specific prediction algorithms. The predicted transmembrane β-barrels have a high similarity in the arrangement of the putative β-strands and the loops, but do not match those of OmpG, a related protein porin whose structure was solved.     

SwarmPS: rapid, semi-automated single particle selection software

Single particle analysis (SPA) coupled with high-resolution electron cryo-microscopy is emerging as a powerful technique for the structure determination of membrane protein complexes and soluble macromolecular assemblies. Current estimates suggest that approximately 10(4)-10(5) particle projections are required to attain a 3A resolution 3D reconstruction (symmetry dependent). Selecting this number of molecular projections differing in size, shape and symmetry is a rate-limiting step for the automation of 3D image reconstruction. Here, we present Swarm(PS), a feature rich GUI based software package to manage large scale, semi-automated particle picking projects. The software provides cross-correlation and edge-detection algorithms. Algorithm-specific parameters are transparently and automatically determined through user interaction with the image, rather than by trial and error. Other features include multiple image handling (approximately 10(2)), local and global particle selection options, interactive image freezing, automatic particle centering, and full manual override to correct false positives and negatives. Swarm(PS) is user friendly, flexible, extensible, fast, and capable of exporting boxed out projection images, or particle coordinates, compatible with downstream image processing suites.     

SCWRL and MolIDE: computer programs for side-chain conformation prediction and homology modeling

SCWRL and MolIDE are software applications for prediction of protein structures. SCWRL is designed specifically for the task of prediction of side-chain conformations given a fixed backbone usually obtained from an experimental structure determined by X-ray crystallography or NMR. SCWRL is a command-line program that typically runs in a few seconds. MolIDE provides a graphical interface for basic comparative (homology) modeling using SCWRL and other programs. MolIDE takes an input target sequence and uses PSI-BLAST to identify and align templates for comparative modeling of the target. The sequence alignment to any template can be manually modified within a graphical window of the target-template alignment and visualization of the alignment on the template structure. MolIDE builds the model of the target structure on the basis of the template backbone, predicted side-chain conformations with SCWRL and a loop-modeling program for insertion-deletion regions with user-selected sequence segments. SCWRL and MolIDE can be obtained at (http://dunbrack.fccc.edu/Software.php).     

Use of negative stain and single-particle image processing to explore dynamic properties of flexible macromolecules

Flexible macromolecules pose special difficulties for structure determination by crystallography or NMR. Progress can be made by electron microscopy, but electron cryo-microscopy of unstained, hydrated specimens is limited to larger macromolecules because of the inherently low signal-to-noise ratio. For three-dimensional structure determination, the single particles must be invariant in structure. Here, we describe how we have used negative staining and single-particle image processing techniques to explore the structure and flexibility of single molecules of two motor proteins: myosin and dynein. Critical for the success of negative staining is a hydrophilic, thin carbon film, because it produces a low noise background around each molecule, and stabilises the molecule against damage by the stain. The strategy adopted for single-particle image processing exploits the flexibility available within the SPIDER software suite. We illustrate the benefits of successive rounds of image alignment and classification, and the use of whole molecule averages and movies to analyse and display both structure and flexibility within the dynein motor.     

XMIPP: a new generation of an open-source image processing package for electron microscopy

X-windows based microscopy image processing package (Xmipp) is a specialized suit of image processing programs, primarily aimed at obtaining the 3D reconstruction of biological specimens from large sets of projection images acquired by transmission electron microscopy. This public-domain software package was introduced to the electron microscopy field eight years ago, and since then it has changed drastically. New methodologies for the analysis of single-particle projection images have been added to classification, contrast transfer function correction, angular assignment, 3D reconstruction, reconstruction of crystals, etc. In addition, the package has been extended with functionalities for 2D crystal and electron tomography data. Furthermore, its current implementation in C++, with a highly modular design of well-documented data structures and functions, offers a convenient environment for the development of novel algorithms. In this paper, we present a general overview of a new generation of Xmipp that has been re-engineered to maximize flexibility and modularity, potentially facilitating its integration in future standardization efforts in the field. Moreover, by focusing on those developments that distinguish Xmipp from other packages available, we illustrate its added value to the electron microscopy community.     

Rapid prototyping of arbitrary 2D and 3D wireframe DNA origami

Wireframe DNA origami assemblies can now be programmed automatically from the top-down using simple wireframe target geometries, or meshes, in 2D and 3D, using either rigid, six-helix bundle (6HB) or more compliant, two-helix bundle (DX) edges. While these assemblies have numerous applications in nanoscale materials fabrication due to their nanoscale spatial addressability and high degree of customization, no easy-to-use graphical user interface software yet exists to deploy these algorithmic approaches within a single, standalone interface. Further, top-down sequence design of 3D DX-based objects previously enabled by DAEDALUS was limited to discrete edge lengths and uniform vertex angles, limiting the scope of objects that can be designed. Here, we introduce the open-source software package ATHENA with a graphical user interface that automatically renders single-stranded DNA scaffold routing and staple strand sequences for any target wireframe DNA origami using DX or 6HB edges, including irregular, asymmetric DX-based polyhedra with variable edge lengths and vertices demonstrated experimentally, which significantly expands the set of possible 3D DNA-based assemblies that can be designed. ATHENA also enables external editing of sequences using caDNAno, demonstrated using asymmetric nanoscale positioning of gold nanoparticles, as well as providing atomic-level models for molecular dynamics, coarse-grained dynamics with oxDNA, and other computational chemistry simulation approaches.