research papers
Conventions and workflows for using Situs
^{a}Department of Physiology and Biophysics and Institute for Computational Biomedicine, Weill Medical College of Cornell University, 1300 York Avenue, New York, NY 10065, USA
^{*}Correspondence email: wriggers@biomachina.org
Situs is a modular program package for the multiscale modeling of atomic resolution structures and lowresolution biophysical data from tomography or smallangle Xray scattering. This article provides an overview of recent developments in the Situs package, with an emphasis on workflows and conventions that are important for practical applications. The modular design of the programs facilitates scripting in the bash shell that allows specific programs to be combined in creative ways that go beyond the original intent of the developers. Several scriptingenabled functionalities, such as flexible transformations of data type, the use of symmetry constraints or the creation of twodimensional projection images, are described. The processing of lowresolution biophysical maps in such workflows follows not only first principles but often relies on implicit conventions. Situs conventions related to map formats, resolution, correlation functions and feature detection are reviewed and summarized. The compatibility of the Situs workflow with CCP4 conventions and programs is discussed.
Keywords: Situs; multiresolution data.
1. Introduction
Situs is a modular commandlinebased opensource package written in C/C++ and is available under the GNU GPL License. Originally designed in 1998/9 to assist in the visualization and interpretation of cryoelectron microscopy (cryoEM) density maps (Wriggers et al., 1999), its scope has been extended over the years to model multiresolution data from a variety of biophysical sources, including tomography and smallangle Xray scattering (SAXS; Wriggers & Chacón, 2001b). The usefulness of hybrid multiscale methodologies that combine atomic structures with lower resolution density maps or coarsegrained models has been well established (Mendelson & Morris, 1997; Lindert et al., 2009) and the historical evolution of Situs with application examples was reviewed in Wriggers (2010).
The focus of this paper is on practical applications, specifically on workflows and conventions that are implicit to Situs programs. Over the years, owing to the modular design, the number of possible combinations of programs has increased to a point where it has become difficult to document possible workflows in our online tutorials. Based on specific modeling tasks, some new usage examples are provided here to inspire users to experiment on their own. Also, in an effort to bridge Situs to other software, many of the implicit conventions are described for the first time.
Our online tutorials (http://situs.biomachina.org) now include Unix bashshell scripts for the automatic generation of tutorial solutions. Much of the workflow complexity originates from the commandlinebased scripting that allows programs to be combined in creative ways. Fig. 1 shows a typical example. The hypothetical problem is that one would like to bring two volumetric density maps into register. The maps can be format converted (§2) with the map2map tool or processed with `volume algebra' tools (§8) such as voledit. For technical reasons, the rigidbody matching tools collage, colores and matchp(oin)t require an atomic PDB file for docking to a target map (Fig. 1). Therefore, the second map must be intermittently transformed to the atomic (PDB) format using vol2pdb so it is free from the cubic lattice (for the rotation and translation in the docking). After the matching of the pseudoPDB map, it can be interpolated back into the volumetric format through projection onto the original lattice with pdb2vol.
Many of the Situs tools rely on implicit conventions for setting parameters. It is perhaps surprising to readers from the crystallographic community that important parameters of volumetric density maps such as resolution, density levels and even map formats are not strictly defined in the hybrid modeling community. For example, the surface isolevel used for visualizing a volume map is an intuitive concept that is surprisingly difficult to solve computationally. Although one can attempt to set the isolevel based on the enclosed volume (Harpaz et al., 1994), the resolution lowering leads to a shift in density, eroding convex features and filling up concave features of the atomic structure. To prevent convex features from protruding from the lowresolution surface, the isolevel of cryoEM maps is empirically set to enclose 120–150% of the molecular volume depending on the overall shape of the system. This is just one example where empirical `fudge factors' trump first principles. Over time, software developers have implemented conventions for a multitude of such quantities as they were breaking new ground. Although an effort is under way to standardize such conventions (Heymann et al., 2005), it is still often necessary to investigate the source code when sharing data between different software packages. In an effort to create more transparency, the most important Situs conventions are documented here.
The remainder of this paper is organized as follows. §2 describes the evolution of the Situs and CCP4based map formats. §3 exemplifies the conversion between lowresolution structure types using smallangle Xrayrelated bead models. §4 contains a comparison of resolution conventions used for multiscale biophysical data. §5 describes conventions for coarsegrained models used in structure matching. §6 presents correlationbased fitting approaches and the conventions used for computing the crosscorrelation. §7 presents workflows enabled by shell scripting, such as the implementation of symmetry constraints and twodimensional projection. The paper concludes with a discussion of extended and supplemental Situs functions in §8.
2. Map format conventions
In this section, we describe threedimensional densitymap formats directly supported by Situs, such as the original Situs map format and the conventions regarding various CCP4derived formats developed at MRC Cambridge, England.
Volumetric maps (Fig. 1) come in a variety of map formats, many of which were inspired by crystallographic formats. For example, the map2map formatconversion utility of Situs accepts cryoEM and crystallographic density files in ASCII (text), CCP4, MRC, Situs, SPIDER and XPLOR formats, automatically adjusts to the machine architecture (endianism) and supports permutation of axes as well as nonorthogonal unit cells (this is accomplished by trilinear interpolation to a cubic lattice, if necessary). In the late 1990s there were no universally accepted standards in the cryoEM community for setting the origin of the map coordinate system, which is critical for the docking of atomic structures. The minimalist Situs format was conceived to keep track of the coordinate system, to enforce a cubic lattice and to be independent of the everchanging mapformat standards in the community. Although this ASCIIbased (and thus easily readable/editable) format was initially meant to be of limited use only within the Situs package, it is now supported by the moleculargraphics programs VMD (Humphrey et al., 1996), Chimera (Pettersen et al., 2004) and Sculptor (Birmanns et al., 2011), by the EMAN2 reconstruction package (Ludtke et al., 1999) and also by the em2em formatconversion tool (http://www.imagescience.de/em2em). In the map format, a short header holds the voxel spacing WIDTH, the map origin as defined by the threedimensional coordinates of the first voxel ORIGX, ORIGY, ORIGZ and the map dimensions (number of increments) NX, NY, NZ. This minimalist header is followed by the data fields such that x increments change fastest and z increments change slowest.
To take advantage of a more compact binary data storage, CCP4 and MRC file formats were recently adopted for direct use within Situs programs, eliminating the map2map formatconversion step from the workflow for many users. The rendering of the map formats was coordinated with the developers of em2em (Michael Schatz), Chimera (Tom Goddard), VMD (John Stone) and Sculptor (Stefan Birmanns). We have made sure that the Situs format matches all other format conventions, such that different map formats and PDB structures are rendered consistently. This effort required, in particular, a revisiting of the CCP4derived map formats developed at MRC Cambridge, which have changed and diversified over the years. The MRC file format was once identical to the CCP4 format; however, incremental changes and the lack of utility of some CCP4 features in the EM community have caused them to become incompatible. Details of the various CCP4 and MRC map formats, and the Situs conventions for reading them, are provided as Supplementary Material^{1}. Our recent efforts were timely because the CCP4 format, as generated by em2em, has become the official format of the EMDB map data bank (Tagari et al., 2002). The detailed conventions we arrived at are described in the Supplementary Material^{1} (pp. 7–8). We hope that our conventions will be more widely adopted and that they help to ensure that maps display correctly.
3. Structuretype conversion
Situs supports three multiresolution structure types: atomic structures (PDB format), volumetric density maps (§2) and SAXS bead models (PDB format). Fig. 1 shows how to interconvert between atomic and volume data.

Situs also supports threedimensional bead models from SAXS (Chacón et al., 2000). Fig. 2 shows an application of SAXSrelated tools in the visualization and atomic interpretation of SAXSderived shapes (Wriggers & Chacón, 2001b). To test the docking accuracy, the pdb2saxs tool was created, which projects atomic structures onto a hexagonal closepacked lattice, with userdefined bead radii written to the PDB occupancy field. The resulting models (Fig. 2b) served as `simulated' lowresolution data in Wriggers & Chacón (2001b).
The bead models in PDB format can be transformed into density maps for subsequent docking using a hardsphere kernel in pdb2vol. The SAXS modeler then has access to all docking strategies supported by Situs, including correlationbased docking (collage or colores), featurepoint based matching (matchpt) and even flexible realspace fitting (§8).
One problem in the interpretation of SAXS data is the visualization of the beads. It is useful to render not the densely packed beads themselves (Fig. 2b), but a transparent wire mesh or envelope (Fig. 2c) that shows the fitted atomic structure. This envelope was created by isocontouring an intermediate volumetric map, which was generated from the beads by convolution with a soft kernel such as a Gaussian (using pdb2vol).
Our approach to the rendering and interpretation of SAXS data has been adopted by other groups in the structural biology community (Lipfert et al., 2007; Fagan et al., 2009; Alvarado et al., 2009). More details, including a workflow and a listing of programs that are used to dock an atomic structure into lowresolution SAXS models, are provided in an online tutorial (http://situs.biomachina.org/tutorial_saxs.html).
4. Resolution convention
This section relates the resolution convention used in Situs with the crystallographic resolution and the well known Rayleigh point resolution in optics.
Resolution R is a quantity in Fourier space and has dimension Å^{−1}, but the realspace quantity r = R^{−1} is also often termed `resolution' in biophysical parlance. The symbols r and R are used here to differentiate between the spatial and frequency domain.
In crystallography, let R_{c} be the radius of a circular region within which Fourier terms contribute to the crystallographic synthesis; one can then say that r_{c} = R_{c}^{−1} is the crystallographic resolution (Frank, 2006). In contrast, Rayleigh considered the resolvability of two points in real space (Stenkamp & Jensen, 1984). The realspace image of one point, whose Fourier transform is limited by a disk of radius R_{c}, is the Airy pattern [J_{1}(2πrR_{c})/(2πrR_{c})]^{2}, where J_{1} is the firstorder Bessel function. The Rayleigh criterion is satisfied when the first minimum of one point's Airy pattern coincides with the central maximum of the other. This critical pointtopoint distance r_{p} turns out to be r_{p} = 0.610r_{c} (see Appendix A3 in Radermacher, 1988). The Airy pattern is dominated by a bellshaped central Airy disk, which measures a `fullwidth at halfmaximum' FWHM = 0.514r_{c}. Ignoring the weak outer rings (resulting from the hard limit in Fourier space), the Airy pattern can be approximated by a Gaussian function that matches the Airy disk (http://en.wikipedia.org/wiki/Airy_disk).
In the development of the pdb2vol tool of Situs in the 1990s, the width of the Gaussian convolution kernel was set empirically to mimic cryoEM maps at published resolution values (for resolution estimation in cryoEM, see Frank, 2006). The empirical Situs resolution r_{s} was set to r_{s} = 2σ (or 1.471 FWHM), where σ is the threedimensional standard deviation of the Gaussian kernel, exp(−3r^{2}/2σ^{2}). One can relate the Situs resolution r_{s} to the crystallographic resolution r_{c} and to the point resolution r_{p} by matching the Gaussian to the Airy disk of the same FWHM. It follows that r_{s} = 1.239r_{p} = 0.756r_{c}. Note that the Sculptor moleculargraphics software (Birmanns et al., 2011) from our laboratory also follows this convention.
The pdb2mrc tool of the EMAN package (Ludtke et al., 1999; renamed e2pdb2mrc.py in EMAN2) is another popular resolutionlowering tool which established its own resolution measure, r_{e}. In EMAN, the functional form of the Gaussian realspace kernel is exp(−π^{2}r^{2}/r_{e}^{2}). The Fourier transform exp(−r_{e}^{2}/R^{2}) is not strictly limited to a hard radius, as in the crystallographic case, but a reasonable radius of the soft Gaussian in Fourier space is set by EMAN at the 1/e cutoff where R = R_{e}. After matching the Gaussian realspace kernel to the Airy disk, it follows that r_{e} = 1.886 FWHM = 1.589r_{p} = 1.282r_{s} = 0.969r_{c}.
In summary, both the EMAN and Situs resolution values are larger than the point resolution and smaller than the crystallographic resolution: r_{p} < r_{s} < r_{e} < r_{c}. However, the Situs resolution r_{s} is closer to the Rayleigh criterion r_{p}, whereas the EMAN resolution r_{e} is almost identical to the crystallographic resolution r_{c} as defined by Frank (2006). When comparing similar maps between Situs (or Sculptor) and EMAN, the resolution values r_{s} (in units of Å) are smaller by a factor of 1.282 compared with r_{e}.
5. Nextgeneration featurepointbased matching
This section describes a recent paradigm shift in setting the level of detail of coarsegrained models to better represent resolved features of the structural data. Prior to Situs v.2.5, the level of detail of the coarsegrained model was limited by the algorithm, whereas now the level of detail is matched to the spatial resolution of the structures.
The idea of using coarsegrained models of feature points for the docking of multiscale structures is a classic idea that predates Situs (Wriggers et al., 1998). We have shown in earlier work that this approach is advantageous at resolutions below 10 Å because the points provide `interior features' (and an encoding of the molecular shape) even in the absence of interior density variations from the secondary structure. However, the earlier Situs tools had a limited number of feature points and required both point clouds to be of equal size. We have recently published a new anchorpoint registration technique that overcomes these combinatorial limitations and is able to dock smaller point clouds to larger ones (Birmanns & Wriggers, 2007).
Fig. 3 shows a typical application of featurepointbased rigidbody matching. The matchpt utility is a commandline program for matching arbitrarysized threedimensional point sets (coarsegrained models), which can be generated on the fly or by using the output of the Situs programs qpdb and qvol. The utility can dock a subunit into a larger target map, i.e. find N feature points within another set with M points, N < M, and match them. To solve this problem, matchpt uses a heuristic and investigates only a subset of all possible permutations of feature points (Birmanns & Wriggers, 2007).
The idea of matching point sets was based on the observation that for many lowresolution maps numeric values of the crosscorrelation (CC) are often in a narrow range and less discriminatory compared with the r.m.s.d. values of the feature points (Wriggers et al., 1999). This is a consequence of the fact that feature points can reliably and reproducibly encode the molecular shape, even in the absence of interior (secondarystructure) density features. Therefore, it makes sense for difficult lowresolution maps to use matchpt as an alternative to the CCbased tools colores and collage. In the default mode, a user would explore the quality of the match of the point clouds by minimizing their r.m.s.d. Alternatively, the minimum of the statistical variability (here the sum of average variabilities of both point sets) can be used to select an optimum N and M, since this variability was found to be a good estimator for the docking accuracy (Wriggers & Birmanns, 2001). Finally, a user may wish to explore the standard crosscorrelation (CC), which is discretely sampled by the solutions of the pointcloud matching.
In Situs v.2.6, the matchpt tool was improved to replace all the functionality of the classic Situs tools qdock and qrange, and the online tutorials were updated accordingly. This introduced a paradigm shift for the way the level of detail of the coarsegrained models (N and M) is estimated. In situations where a smaller structure is docked into a larger density (e.g. an oligomeric map), the units parameter defines the fraction of occupied volume (which may be noninteger), i.e. it estimates how many atomic input structures fit into the target volume. M is then defined as units × N. To estimate the number N of feature points, one can divide the volume of the atomic structure by the volume of a resolution element, r_{s}^{3} (where r_{s} is the resolution value of the target map in Å). This calculation gives an upper bound for the number of features N contained in the structure (and, via M, the number of features contained in the volume) at the given map resolution. To avoid overfitting and to find an optimal number, it is useful in practical applications to bracket N between the 30 to 50% level of this upper bound.
6. Nextgeneration correlationbased matching
This section highlights recent advances in correlationbased fitting technology and the underlying convention in the calculation of the CC.
The Situs crosscorrelation coefficient (CC) for volumetric maps ρ_{vol}(r) and ρ_{calc}(r) is defined as
where ρ_{vol}(r) is a lowresolution map and ρ_{calc}(r) is an atomic structure subject to rigidbody movements, projection to the volume V and convolution with a Gaussian kernel (Wriggers & Chacón, 2001a). The CC values are normalized, but unlike in the Pearson the averages are not subtracted from the densities. Consequently, the calculation is slightly more efficient than that of the Pearson CC. The convention takes into account that ρ_{calc}(r) often has the physical meaning of a density with well defined positive amplitude that corresponds to the represented lowresolution structure. This way, CC ∈ [0, 1] for ρ_{vol}(r), ρ_{calc}(r) > 0. One can show that maximizing the Situs CC value also maximizes the Pearson so dropping the subtraction of averages in (1) yields no performance penalty. The Situs CC convention has also been adopted by the Sculptor visualization program (Birmanns et al., 2011).
Chacón & Wriggers (2002) introduced colores, a widely used CCbased fitting tool that takes advantage of Fourier correlation theory to rapidly scan the translational degrees of freedom of a probe molecule relative to a (fixed) targetdensity map (whereas the rotations are sampled exhaustively by the enumeration of a list of homogeneously distributed Euler angles). The performance of the standard CC (1) is limited to resolutions higher than 10 Å, where densities exhibit the internal (i.e. secondary) structure. The major advantage of colores is that it extends the viable resolution range to ∼30 Å by means of an (optional) Laplacian operator [applied to both ρ_{vol}(r) and ρ_{calc}(r) in equation 1] that emphasizes shapecontour information in addition to the traditional volume correlation (also, a masking filter was implemented that suppresses singularities of the Laplacian of ρ_{vol} at density edges resulting from thresholding or segmentation).
Recently, the new collage was introduced which performs a conjugategradient optimization of the same scoring functions known from colores. The main innovation is the simultaneous optimization of multiple rigid fragments that `see' each other and avoid steric clashes by means of the normalization in (1). Birmanns et al. (2011) showed that this approach yields more accurate fits. Also, if all density is accounted for by the fragments it is no longer necessary to use the Laplacian filter option, even at low resolution.
toolFig. 4 shows an application example of collage using the simultaneous optimization of six monomers. A single run of offlattice Powell optimization is applied that refines a preliminary multifragment model (consisting here of six input PDB files) to the nearest maximum of the CC. The start model of fragments could, for example, be derived manually by eye in a graphics program (as was performed here) or it could be based on colores or matchpt solutions. We have created a new multifragment online tutorial to explain this powerful new approach.
7. Scriptingbased workflows
This section describes the benefits of bashshell scripting to tie together multiple tasks. Example scripts are provided here for the implementation of symmetry constraints in multifragment fitting and for the creation of twodimensional projections from threedimensional maps.
The new collage tool and external symmetrymanipulation programs can be combined to impose symmetry constraints on the fragments during multifragment docking. Fig. 5 provides an overview of the workflow. Standard volumetric map formats are converted to cubic lattices in Situs format with the map2map utility. Subsequently, the volume data are inspected and, if necessary, prepared for the fitting using a variety of map tools (§8). A datatype conversion (§3) is optional. The fitting of multiple PDB input files to the target map is handled by collage as described in §6. Symmetry constraints are outsourced to a separate program in the Unix bashshell script. The Situsnative pdbsymm tool allows the generation of multiple symmetryrelated copies following symmetryaxis conventions based on the target map. C, D and H (helical) symmetry options are currently supported. For other specialized cases (e.g. crystallographic symmetry), an alternative program (such as the CCP4 pdbset tool; http://www.ccp4.ac.uk/html/pdbset.html) may be substituted for pdbsymm in the workflow (Fig. 5). The entire Unix bashshell implementation then proceeds as follows.

The goal of the scripting approach is to keep Situs tools modular and to avoid having to write specialized tools for every possible symmetry scenario. This means that collage will technically still treat each fragment as independent. However, collage will take only a single conjugategradient step (step iv), after which the symmetry will again be enforced (in steps ii and vi). The net effect after several iterations of the loop is a symmetryenforced conjugategradient optimization. An example of this approach is available online in the multifragment tutorial (file run_tutorial.bash at http://situs.biomachina.org/tutorial_multi.html).
A new function was added to the voledit tool in Situs v.2.6 that conveniently allows the user to render twodimensional projections in addition to volume slices. Such projections are useful when comparing threedimensional maps or resolutionlowered atomic structures with twodimensional micrographs. Fig. 6 shows that discrepancies between two conformations can easily be detected in two dimensions. The projections (or slices) saved by voledit can be visualized by an external plotting program. Situs is mainly a threedimensional package and does not have a tool for computing differences between twodimensional images, but the difference of two twodimensional projections is identical to the projection of the threedimensional difference map, which can be computed with voldiff. The entire workflow used to create Fig. 6 consists of several steps that are most efficiently implemented in a shell script as follows.

8. Discussion
The Situs software consists of multiple standalone tools that can be combined in various creative ways. The modular design offers great advantages to the user who wishes to take advantage of shellscripting capabilities. To encourage experimentation, bash scripts are included with all online tutorials. This overview also summarizes important conventions that form the basis of Situs functionality.
One active research area that exceeds the scope of this article is flexible docking, which bring deviating features of multiresolution structures into register (Wriggers et al., 2004; Wriggers, 2010). Fig. 6 shows the differences between flexed structures of RNA polymerase in the form of twodimensional projections. Systematic tests and validations of flexible fitting with spatial interpolation have been carried out using the qplasty tool (Rusu et al., 2008) and experimental applications of flexible docking (using moleculardynamics refinement) were performed in collaboration with experimental laboratories on systems such as RNA polymerase (Darst et al., 2002) and the thick filament of tarantula muscle (Alamo et al., 2008).
Another noteworthy application that exceeds the scope of this paper is `volume algebra', in which map densities are modified according to simple algebraic operations (Wriggers et al., 2011). Volumealgebra operations are typically enabled by mapdensity registration (Fig. 1) and include map summation or averaging (volaver), difference mapping (voldiff), binary masking and map multiplication (voledit, volmult) and density matching using an affine transformation (volhist), as well as cropping, thresholding and segmentation (voledit). For application examples of volume algebra operations, see Fig. 1 in Wriggers et al. (2011) and the online tutorials.
The performance of multifragmentbased and 5) is the subject of ongoing research. Our empirical tests with experimental maps have shown that the is very robust. The normalization of the crosscorrelation coefficient (1) penalizes steric clashes between fragments, similar to features of a tabu search using genetic algorithms (Rusu & Birmanns, 2010). Also, the radius of convergence for the method appears to be rather large, reducing the need for exhaustive exploration. The performance will be evaluated further in future work.
(Figs. 4Situs has been ported to multiple platforms and the source code is freely available at http://situs.biomachina.org.
A comparison of various CCP4derived map formats and related Situs read and write conventions are provided in the Supplementary Material^{1}.
Supporting information
Supporting information file. DOI: https://doi.org/10.1107/S0907444911049791/ba5170sup1.pdf
Acknowledgements
I thank Michael Schatz, Mirabela Rusu, Stefan Birmanns, John Stone, Tom Goddard, Yao Cong and Steve Ludtke for helpful discussions. This work was supported in part by a grant from the National Institutes of Health (R01GM62968).
References
Alamo, L., Wriggers, W., Pinto, A., Bártoli, F., Salazar, L., Zhao, F.Q., Craig, R. & Padrón, R. (2008). J. Mol. Biol. 384, 780–797. CrossRef PubMed CAS Google Scholar
Alvarado, D., Klein, D. E. & Lemmon, M. A. (2009). Nature (London), 461, 287–291. CrossRef PubMed CAS Google Scholar
Birmanns, S., Rusu, M. & Wriggers, W. (2011). J. Struct. Biol. 173, 428–435. Web of Science CrossRef CAS PubMed Google Scholar
Birmanns, S. & Wriggers, W. (2007). J. Struct. Biol. 157, 271–280. CrossRef PubMed CAS Google Scholar
Chacón, P., Díaz, J. F., Morán, F. & Andreu, J. M. (2000). J. Mol. Biol. 299, 1289–1302. Web of Science CrossRef PubMed CAS Google Scholar
Chacón, P. & Wriggers, W. (2002). J. Mol. Biol. 317, 375–384. Web of Science CrossRef PubMed CAS Google Scholar
Darst, S., Opalka, N., Chacón, P., Polyakov, A., Richter, C., Zhang, G. & Wriggers, W. (2002). Proc. Natl. Acad. Sci. USA, 99, 4296–4301. CrossRef PubMed CAS Google Scholar
Fagan, R. P., AlbesaJové, D., Qazi, O., Svergun, D. I., Brown, K. A. & Fairweather, N. F. (2009). Mol. Microbiol. 71, 1308–1322. Web of Science CrossRef PubMed CAS Google Scholar
Frank, J. (2006). ThreeDimensional Electron Microscopy of Macromolecular Assemblies. Oxford University Press. Google Scholar
Harpaz, Y., Gerstein, M. & Chothia, C. (1994). Structure, 2, 641–649. CrossRef CAS PubMed Web of Science Google Scholar
Heymann, J. B., Chagoyen, M. & Belnap, D. M. (2005). J. Struct. Biol. 151, 196–207. Web of Science CrossRef PubMed Google Scholar
Humphrey, W. F., Dalke, A. & Schulten, K. (1996). J. Mol. Graph. 14, 33–38. CrossRef CAS PubMed Web of Science Google Scholar
Kobe, B. & Deisenhofer, J. (1996). J. Mol. Biol. 264, 1028–1043. CrossRef CAS PubMed Web of Science Google Scholar
Lindert, S., Stewart, P. L. & Meiler, J. (2009). Curr. Opin. Struct. Biol. 19, 218–225. CrossRef PubMed CAS Google Scholar
Lipfert, J., Das, R., Chu, V. B., Kudaravalli, M., Boyd, N., Herschlag, D. & Doniach, S. (2007). J. Mol. Biol. 365, 1393–1406. Web of Science CrossRef PubMed CAS Google Scholar
Ludtke, S. J., Baldwin, P. R. & Chiu, W. (1999). J. Struct. Biol. 128, 82–97. Web of Science CrossRef PubMed CAS Google Scholar
Mendelson, R. & Morris, E. P. (1997). Proc. Natl. Acad. Sci. USA, 94, 8533–8538. CrossRef CAS PubMed Google Scholar
Pettersen, E. F., Goddard, T. D., Huang, C. C., Couch, G. S., Greenblatt, D. M., Meng, E. C. & Ferrin, T. E. (2004). J. Comput. Chem. 25, 1605–1612. Web of Science CrossRef PubMed CAS Google Scholar
Radermacher, M. (1988). J. Electron Microsc. Tech. 9, 359–394. CrossRef CAS PubMed Web of Science Google Scholar
Rusu, M. & Birmanns, S. (2010). J. Struct. Biol. 170, 164–171. CrossRef CAS PubMed Google Scholar
Rusu, M., Birmanns, S. & Wriggers, W. (2008). Bioinformatics, 24, 2460–2466. Web of Science CrossRef PubMed CAS Google Scholar
Stenkamp, R. E. & Jensen, L. H. (1984). Acta Cryst. A40, 251–254. CrossRef CAS Web of Science IUCr Journals Google Scholar
Tagari, M., Newman, R., Chagoyen, M., Carazo, J. M. & Henrick, K. (2002). Trends Biochem. Sci. 27, 589. Web of Science CrossRef PubMed Google Scholar
Wriggers, W. (2010). Biophys. Rev. 2, 21–27. CrossRef PubMed Google Scholar
Wriggers, W., Alamo, L. & Padrón, R. (2011). J. Struct. Biol. 173, 445–450. CrossRef CAS PubMed Google Scholar
Wriggers, W. & Birmanns, S. (2001). J. Struct. Biol. 133, 193–202. Web of Science CrossRef PubMed CAS Google Scholar
Wriggers, W. & Chacón, P. (2001a). Structure, 9, 779–788. CrossRef PubMed CAS Google Scholar
Wriggers, W. & Chacón, P. (2001b). J. Appl. Cryst. 34, 773–776. CrossRef CAS IUCr Journals Google Scholar
Wriggers, W., Chacón, P., Kovacs, J., Tama, F. & Birmanns, S. (2004). Neurocomputing, 56, 365–379. CrossRef Google Scholar
Wriggers, W., Milligan, R. A. & McCammon, J. A. (1999). J. Struct. Biol. 125, 185–195. Web of Science CrossRef PubMed CAS Google Scholar
Wriggers, W., Milligan, R. A., Schulten, K. & McCammon, J. A. (1998). J. Mol. Biol. 284, 1247–1254. Web of Science CrossRef CAS PubMed Google Scholar
Yu, X. & Egelman, E. H. (1997). Nature Struct. Biol. 4, 101–104. CrossRef CAS PubMed Web of Science Google Scholar
This is an openaccess article distributed under the terms of the Creative Commons Attribution (CCBY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.