research papers
Automatic multiple-zone rigid-body
with a large convergence radiusaLawrence Berkeley National Laboratory, One Cyclotron Road, BLDG 64R0121, Berkeley, CA 94720, USA, bIGBMC, CNRS-INSERM-UdS, 1 rue Laurent Fries, BP 10142, 67404 Illkirch, France, cDépartement de Physique, Université Henri Poincaré, Nancy 1, BP 239, Faculté des Sciences et des Technologies, 54506 Vandoeuvre-lès-Nancy, France, and dDepartment of Bioenegineering, University of California Berkeley, CA 94720, USA
*Correspondence e-mail: pafonine@lbl.gov
Rigid-body phenix.refine program.
is the constrained coordinate of one or more groups of atoms that each move (rotate and translate) as a single body. The goal of this work was to establish an automatic procedure for rigid-body which implements a practical compromise between runtime requirements and convergence radius. This has been achieved by analysis of a large number of trial refinements for 12 classes of random rigid-body displacements (that differ in magnitude of introduced errors), using both least-squares and target functions. The results of these tests led to a multiple-zone protocol. The final parameterization of this protocol was optimized empirically on the basis of a second large set of test refinements. This multiple-zone protocol is implemented as part of theKeywords: rigid-body refinement; multiple-zone protocols.
1. Introduction
The vast majority of macromolecular crystal structures are solved either with experimental phasing methods (see, for example, Blow & Crick, 1959; Hendrickson, 1991) or with the method (Rossmann & Arnold, 2001, and references therein). In the case of experimental phasing the model is built into an The resulting model may contain many local errors, but significant concerted displacements are not expected. In contrast, models obtained via or with difference Fourier methods can be systematically displaced. In this situation rigid-body (Booth, 1947a,b, 1949; Cochran, 1948; Scheringer, 1963; Sussman et al., 1977; Hoard & Nordman, 1979; Huber & Schneider, 1985; Yeates & Rees, 1988; Derewenda, 1989; Urzhumtsev et al., 1989; Driessen et al., 1989; Yeates & Rini, 1990; Brünger, 1990a,b, 1991; Castellano et al., 1992; Noble et al., 1993; Navaza, 2001; Tronrud, 2004; McCoy, 2007; Lebedev et al., 2008) is a powerful method for correcting potentially large systematic displacements. Outside the field of crystallography, rigid-body is also an important tool when fitting models into envelopes (see, for example, Navaza et al., 2002, and numerous references therein). Rigid-body may also be a way of performing coordinate when only very low resolution data are available.
Rigid-body ; Sheldrick & Schneider, 1997; Tronrud, 2004). Least-squares (LS) is expected to be more affected than (ML), since the latter is designed to automatically weight down terms with poor model-to-data correspondence (Lunin et al., 2002), i.e. data at high resolution at the beginning of refinement.
moves groups of atoms as a whole, leaving the internal configuration of each group unchanged. It is well understood that the information about the large-scale distribution of atoms is contained in the low-resolution diffraction data. The high-resolution data convey information about the finer details of the atomic structure. Since these details are invariant during rigid-body it is expected that high-resolution data will be less important for this procedure than the low-resolution data. Inclusion of high-resolution data is known to hamper the progress of (see, for example, Sheldrick, 2008When choosing the high-resolution cutoff for SHELX (Sheldrick & Schneider, 1997). At the initial stage the convergence radius is large. The model is most likely to be moved closer to the correct position and orientation, but the accuracy is relatively low. At the subsequent stages the convergence radius is less critical, but the accuracy is improved by the inclusion of higher-resolution data. This approach can be robust but computationally expensive and requires ad hoc decisions about high-resolution data cutoff and the number of higher-resolution reflections to be added as the progresses. Here we report the results of numerical experiments aimed at finding a computationally economical and automated multiple-zone protocol that still results in a large convergence radius.
a practical balance between convergence radius, accuracy of the results and computational cost has to be found. Generally, moving the cutoff to lower resolution is expected to increase the radius of convergence, but at the cost of decreased accuracy. This suggests a multiple resolution approach with several sequential refinements using data at increasingly higher resolution, for example, as implemented by the STIR option inThe multiple-zone phenix.refine – a macromolecular structure program (Afonine et al., 2005b) that is under active development as part of the PHENIX project (Adams et al., 2002). Major development goals are increased automation and fast exploration of new approaches based on a modular architecture. Available features, among others, include various targets (maximum likelihood, twinned least squares, phased maximum-likelihood), of individual coordinates and ADPs (isotropic, anisotropic, group, TLS or any combination), automatic water picking built in to the robust bulk-solvent correction (Afonine et al., 2005a), Cartesian dynamics, simulated annealing, NCS restraints, at ultra-high resolution (Afonine et al., 2004, 2007), and joint using X-ray and neutron data. Here we describe the systematic investigation of rigid-body based on a large number of trial refinements in phenix.refine.
protocol is implemented in2. Methods
2.1. Parameterization of rigid-body motions
A rigid body is a group of atoms subject to a concerted motion, leaving the atoms fixed relative to each other. In rigid-body ; Weisstein, 2006). The Euler angles are commonly referred to as α, β, γ. One commonly used convention [e.g. AMoRe (Navaza, 2001) and REFMAC (Collaborative Computational Project 4, Number 4, 1994; Murshudov et al., 1997)] is to first rotate around the Cartesian z axis by the angle γ, then around the y axis by the angle β, and finally around the z axis again by the angle α; in this paper we refer to this convention as the zyz convention. At the usual starting point for rigid-body α = β = γ = 0°, α and γ are perfectly correlated, which could potentially lead to numerical instabilities. Another convention in common use (e.g. Urzhumtsev et al., 1989; Brünger et al., 1998; Kronenburg, 2004) differs in this respect. The first two rotations are as before, but the third rotation is around the x axis. Here we refer to this convention as the xyz convention. α and γ are perfectly correlated only if γ = ±90°, values that are highly unlikely to be reached in the course of rigid-body as the final rotations from the starting position are typically less than 20°.
a macromolecule is split into one or more non-overlapping rigid groups. The position of each rigid body is characterized by six The body translation is universally parameterized as three Cartesian or fractional coordinates. The body orientation is usually defined by three Euler angles. A large number of Euler angle conventions are in use (Urzhumtseva & Urzhumtsev, 19972.2. procedure
In the tests reported below, a rigid-body ) (the same minimizer is used in the CNS program) with a maximum of 25 iterations per macro cycle. During minimization, an LS or an ML target function [implemented as defined by Lunin & Skovoroda (1995) and Afonine et al. (2005a)] is used. All geometry restraints including nonbonded interactions are disabled. Thus the is purely based on the experimental data.
run is a series of macro cycles in each of which a bulk-solvent correction is followed by L-BFGS minimization (Liu & Nocedal, 1989It has been shown that a bulk-solvent correction of the low-resolution data is very important to achieve optimal ; Kostrewa, 1997; Badger, 1997). In the tests reported below, we used the bulk-solvent correction algorithm as described by Afonine et al. (2005a). In the context of rigid-body the model shifts are expected to be large and hence invalidate the bulk solvent mask calculated from the initial model. Therefore the bulk-solvent correction is tightly integrated into the and recomputed between macro-cycles if the model has moved beyond a certain default threshold.
results (Jiang & Brünger, 19942.3. Test data and models
Test data and models were taken from an in-house library of 56 structures collected over the course of some time. Table 1 lists reference information for all test structures. The original experimental data were used in all trial refinements reported below. To better approach typical practical situations, the models from the library were modified by deleting all atoms that are not part of a protein, RNA or DNA molecule.
As a last manipulation, all structures were subject to rigid-body ) with the entire model as one body. The refined corrections were typically very small. These refined models were considered as the best possible results and used as the ideal (reference) model for subsequent comparisons.
using data up to 3 Å (or the high-resolution limit shown in Table 1For tests with multiple rigid-bodies, seven out of the 56 structures were split into two to six bodies as indicated in Table 1 (column NB). We included calmodulin and gene-5 models even though they both consist of only one chain. Calmodulin was chosen because of the (P1), gene-5 because of the small size. In each model, the chain was split into two parts and a few atoms were deleted to avoid clashes of the two artificially created bodies.
One of the multi-body structures (1071B) was subject to multi-body rigid-body
at 3 Å to obtain a specific multi-body reference model, since the displacements with respect to the one-body model were significant.2.4. Random displacements
Our goal was to systematically sample the behavior of rigid-body refinements. For this we investigated three main variables:
(i) Averaging out model shape-related effects by using a large number of models (see previous section).
(ii) Sampling a matrix of displacement magnitudes, using combinations of translations with a 0, 2, 4 and 6 Å shift along a random vector and a 0, 5, 10 and 15° rotation around a random axis.
(iii) Averaging out effects due to interactions of model shape and translation vectors or rotation axes by sampling a large number (we used 100) of random vectors and axes for a given pair of displacement magnitudes.
The combination (translation, rotation) = (0 Å, 0°), i.e. no change in the starting model position or orientation, was excluded. To reduce the runtime for the tests, we also chose to omit the (4 Å, 10°), (6 Å, 10°) and (6 Å, 15°) combinations since the success rate (see §2.5 for the definition) for such very large displacements was known to be near zero, on the basis of preliminary trials. This left 12 combinations to be sampled 100 times for each of the 56 test structures, i.e. a total of 67 200 rigid-body runs per set of trial parameters. When determining the random translations, continuous allowed origin shifts (Grosse-Kunstleve, 1999) (e.g. parallel to the twofold axis in P2) were specifically taken into account: the translation vectors for the first body were chosen perpendicular to the allowed origin shifts. In P1 translations of the first body have no effect on the magnitudes and were therefore not considered.
In the tests with multiple rigid bodies, the random translations and rotations can lead to serious clashes, which are unlikely to occur in most practical situations since most ), the clashes have no direct effect and we decided to ignore them.
programs generally suppress configurations with clashes. However, since nonbonded interactions are not included in our rigid-body procedure (§2.22.5. Success rates
After each rigid-body ) was determined, taking allowed origin shifts into account. For each set of 100 random displacement magnitudes (previous section), we counted the number of refined models with an r.m.s.d. less than or equal to 1.0, 0.5 and 0.25 Å. These numbers are the success rates in percent, given the chosen r.m.s.d. value. Fig. 1 shows an example plot of the success rates.
run, the root-mean-square deviation (r.m.s.d.) with respect to the reference pre-refined model (§2.3When evaluating the effects of parameter changes, we compared the success rates using a tolerance to eliminate noise. Success rates with differences less than or equal to 2% were considered insignificant. Larger differences were considered significant and used as a guide in the optimization of the
protocol.3. Results
3.1. Refinements with fixed high-resolution cutoffs
Our initial test series was a systematic sampling of high-resolution cutoffs dmin = 3, 4, 6, 8 and 10 Å. Some data sets had an insufficient number of reflections given the 8 or 10 Å cutoffs. The largest resolution cutoffs used in these tests are 6 Å for gene-5 and hipip, 8 Å for gpatase, lysozyme, oat-gabaculine and rnase-p, and 10 Å for all other structures. The total number of rigid-body runs over all five resolution ranges was 32 6400 = (50 structures × 5 resolution cutoffs + 2 structures × 3 resolution cutoffs + 4 structures × 4 resolution cutoffs) × 12 rotation–translation shifts × 100 trials. The complete test series was run twice: once using the LS target function, then again using the ML target function. We manually reviewed the resulting 2 × 3264 success rate plots, where each plot was an average over 100 trials.
While there are significant individual differences between the test structures, the results show a general trend. This observation led us to prepare plots averaging the success rates over all structures (Figs. 2 and 3) so that each point in these plots shows the success rate of 100 × 56 refinements (with a few refinements less at 8 and 10 Å as explained above). This leads to the following observations:
(i) The success rates reach a plateau after a certain number of macro cycles. The height of the plateau depends on both the displacement magnitude and the high-resolution cutoff. The larger the displacement magnitudes, the lower the plateau. The larger the high-resolution cutoff, the higher the plateau.
(ii) The macro cycle at which the plateau is reached depends on the high-resolution cutoff. The larger the high-resolution cutoff, the more macro cycles are needed to reach the plateau.
(iii) The difference in the plateau heights for the three success rate cutoffs (1.0, 0.5, 0.25 Å) strongly depends on the high-resolution cutoff. With a high-resolution cutoff of 3.0 Å, the three plateaus in each plot have virtually identical heights. This suggests that, if a
converges to the solution, it is highly likely to be accurate. As the high-resolution cutoff is increased, the plateaus are at increasingly different heights. This means on average the solutions are increasingly less accurate.(iv) Increasing the high-resolution cutoff leads to a larger convergence radius. This effect is most pronounced for translational displacements when going from a 3.0 to a 4.0 Å high-resolution cutoff, or from 4.0 to 6.0 Å. Larger high-resolution cutoffs do not significantly increase the convergence radius. Furthermore, the effect is weaker for rotational displacements.
(v) The difference in results obtained with the least-squares target and the and 3 we observed that the least-squares target leads to slightly better success rates using high-resolution cutoffs of 6.0 Å or larger, and the target is slightly better using high-resolution cutoffs of 3.0 and 4.0 Å.
target are subtle. Comparing Figs. 23.2. Multiple-zone protocol
The observations reported in the previous section lead to the use of a multiple-zone protocol. The goal is to take advantage of the larger convergence radius at larger high-resolution cutoffs and higher accuracy at smaller resolution cutoffs. The multiple-zone protocol automates dmin. A straightforward approach is to decrement the high-resolution cutoff dmin by a certain amount after each round of rigid-body similar to the SHELX STIR option. Under this scheme the zones increase in size with the cube of the number of reflections. We chose to use a similar but more tunable function with a parameterization that is designed to be independent of the structure to be refined:
starting with a small number, n_ref(1), of low-resolution reflections (first zone), and successive addition of reflections up to a user-defined high-resolution cutoffzone_exponent is a user-defined value. The zone_factor is computed from a user-defined number of zones, n_zones, and the number of reflections at the highest resolution cutoff, n_ref(n_zones):
zone_exponent = 3 corresponds to the SHELX STIR option. With a smaller value the function is more linear, adding more reflections more quickly. With a larger value fewer reflections are added initially and more reflections in the later steps.
n_ref(1) is determined using the formula
Here n_ref(1)1 is a user-supplied value, n_bodies is the number of user-supplied atom selections for the rigid bodies and multi_body_factor is a tunable parameter. With multi_body_factor = 1 the number of reflections for the first resolution zone is a linear function of the number of rigid bodies. The phenix.refine program provides two alternatives for determining n_ref(1)1. The user can simply specify the value directly or specify a low-resolution cutoff from which n_ref(1)1 is computed. In all cases, default values are automatically chosen by the program, which can be overridden by user-defined values if required.
3.2.1. Exploration of parameter space
The parameterization presented in the previous section is designed to be independent of the structure to be refined. The critical variables are the number of zones n_zones, n_ref(1)1, the multi_body_factor and the zone_exponent. These values were optimized with a series of tests, using starting values derived from the results of the tests described above. To decrease the runtime requirements for a test series, we reviewed the refinements with fixed high-resolution cutoffs. In addition to the seven multi-body refinements, three single-body refinements were chosen with the aim of covering the distribution of number of atoms versus high-resolution limit of the diffraction data. The three selected structures are marked in Table 1, column NB.
A preliminary set of test runs with only seven structures split into multiple bodies indicated that n_zones = 5, zone_exponent = 4, multi_body_factor = 1 and n_ref(1)1 = 100 is a good default parameterization. After adding the three single-body refinements, we then started a second set of tests exploring the parameter space around these values. The results are reported in some detail below. Unless noted otherwise, the target was used in all resolution zones.
To determine the best choice for n_ref(1)1 we ran a test series with trial values 60, 80, 100, 120, 200 and 400. The values were chosen on the basis of the behavior of the trial refinements. Table 2 shows a mutual comparison of the n_ref(1)1 values of each trial with the others. Inspection suggests that the value 100 is the best choice overall for n_ref(1)1.
|
After identifying n_ref(1)1 = 100 as the best value, we explored values for multi_body_factor = 0.5, 1 and 2, fixing n_zones = 5 and zone_exponent = 4 as before. The results in Table 3 confirm our expectation that a simple linear coupling (multi_body_factor = 1) of the number of observations and the number of refineable parameters is optimal.
|
Fixing n_ref(1)1 = 100, multi_body_factor = 1 and n_zones = 5, we tried the alternative values zone_exponent = 1, 2, 3, 5 and 6. The results in Table 4 show that zone_exponent = 3 is the best choice. This result validates the approach used in the SHELX STIR algorithm.
|
Fixing n_ref(1)1 = 100, multi_body_factor = 1 and zone_exponent = 3, we tried the alternative values n_zones = 3, 4, 6, 7, 8 and 9. The disadvantage of using more zones is increased runtime. However, the results are expected to be better if more zones are used. This is largely confirmed by the data in Table 5. The use of more than five zones does not greatly improve the outcome compared with the use of three to five zones. Table 6 shows runtime statistics as a function of the number of zones. For example, using seven zones instead of five zones increases the runtime, on average, by about 25%. Thus, using five zones is a practical compromise between runtime considerations and expected benefit. However, for difficult cases it could be worth increasing the number of zones in order to increase the expected success rate, at the cost of increased runtime. It should also be noted that as rigid-body is often typically only performed once at the start of structure investing a some additional computing time to obtain the best solution might be well justified.
|
|
3.3. Effect of rotation convention
To analyze the role of different ways to describe the rotation we made a comparative shows the success rate comparison using the two different Euler angle conventions introduced in §2.1. All other parameters were fixed at the defaults [n_ref(1)1 = 100, multi_body_factor = 1, zone_exponent = 3 and n_zones = 5]. The results are surprisingly clear: the xyz convention drastically outperforms the zyz convention. Even though the L-BFGS minimizer used in the refinements is designed to tolerate singularities, it is evidently advantageous to avoid them.
in similar conditions using the two different Euler angle conventions. Table 7
|
3.4. Automatic switching between least-squares and target functions
In §3.1 it was found that the best choice of target function depends on the high-resolution cutoff. To take advantage of this knowledge a target_auto_switch_resolution parameter was introduced. For zones with a high-resolution cutoff larger than the value of this parameter, the LS target is used, and the ML target otherwise. With optimal values for the other parameters as presented in the previous sections, a new series of four tests were performed, with target_auto_switch_resolution = 4, 5, 6 and 7 Å. The corresponding success rate comparisons are shown in Table 8. These data indicate that switching at a lower value for the resolution cutoff is better than switching at a higher value. On the basis of the fixed resolution cutoff results (Figs. 2 and 3), 6 Å was selected as the default parameter.
|
4. Conclusion
of an atomic model as a rigid body or several independent rigid bodies is an important and routine step in macromolecular By combining multi-zone rigid body robust bulk solvent and scaling, methods, and large-scale optimizations of key parameters of the multi-zone protocol, it has been possible to increase the radius of convergence of rigid-body and make the process highly automated.
We observe that although likelihood methods do provide for some degree for automated weighting of data the extent of this weighting is not sufficient to provide a large radius of convergence when high-resolution data are used in the rigid-body
Thus, explicit removal of higher-resolution data during rigid-body even when using a likelihood target, significantly increases the radius of convergence.The idea of gradually increasing the
resolution from low to high or simply truncating the high-resolution data at some point between 3 and 6 Å resolution has been used by many practitioners for some time. However, the systematic investigation undertaken here has provided two very important enhancements. First, use of the number of reflections to define the resolution zones, instead of specific resolution limits, makes the process model-independent. For example, for relatively small structures cutting the data at 6 Å or even higher may not leave enough low-resolution reflections for In addition, defining the zones by the number of reflections always assures an adequate amount of data and an appropriate resolution for the Second, since large model shifts are expected during rigid-body it is essential to update the bulk solvent model as often as a model shifts beyond a certain threshold or additional reflections are included.Having performed hundreds of thousands of rigid-body refinements using a set of 56 models different in size, shape and packing, as well as having different quality experimental data sets associated with them (resolution and completeness), we have empirically confirmed our hypothesis that the xyz rotation parameterization performs better than the zyz parameterization since it avoids a singularity near the typical values for the rotation parameters encountered in Our results also define approximate convergence radii for gradient-driven rigid-body (Figs. 2 and 3).
Our systematic exploration of the parameter space (§3.2.1) was feasible only because we had access to a computer cluster with 200 fast CPUs. The tools that we have developed for running and analyzing the many refinements are being re-used for evaluating other algorithms and parameterizations. By making use of the latest computing technology, we can replace the very slow and subjective process of tuning parameters based on anecdotal evidence with a more scientific approach. With modern tools, practical experience that may have taken many years to accumulate in the past can now be obtained in a matter of days.
The described algorithms and protocols are implemented in the phenix.refine, which is available as part of the PHENIX package. The program is available from https://www.phenix-online.org/ . The core rigid-body calculations are part of open source libraries (https://cctbx.sourceforge.net ).
programAcknowledgements
This work was supported in part by the US Department of Energy under contract No. DE-AC03-76SF00098 and NIH/NIGMS grant 1P01GM063210. AU thanks Pole `Intelligence Logicielle', Nancy, for financial support. We thank L. Urzhumtseva for sharing her previous experience on rigid-body
and V. Lunin for useful discussions.References
Adams, P. D., Grosse-Kunstleve, R. W., Hung, L.-W., Ioerger, T. R., McCoy, A. J., Moriarty, N. W., Read, R. J., Sacchettini, J. C., Sauter, N. K. & Terwilliger, T. C. (2002). Acta Cryst. D58, 1948–1954. Web of Science CrossRef CAS IUCr Journals Google Scholar
Afonine, P. V., Grosse-Kunstleve, R. W. & Adams, P. D. (2005a). Acta Cryst. D61, 850–855. Web of Science CrossRef CAS IUCr Journals Google Scholar
Afonine, P. V., Grosse-Kunstleve, R. W. & Adams, P. D. (2005b). CCP4 Newsletter No. 42, https://www.ccp4.ac.uk/newsletters/newsletter42/content.html . Google Scholar
Afonine, P. V., Grosse-Kunstleve, R. W., Adams, P. D., Lunin, V. Y. & Urzhumtsev, A. (2007). Acta Cryst. D63, 1194–1197. Web of Science CrossRef CAS IUCr Journals Google Scholar
Afonine, P. V., Lunin, V. Y., Muzet, N. & Urzhumtsev, A. (2004). Acta Cryst. D60, 260–274. Web of Science CrossRef CAS IUCr Journals Google Scholar
Badger, J. (1997). Methods Enzymol. 277, 344–352. CrossRef PubMed CAS Web of Science Google Scholar
Berman, H. M., Westbrook, J., Feng, Z., Gilliland, G., Bhat, T. N., Weissig, H., Shindyalow, I. N. & Bourne, P. E. (2000). Nucleic Acids Res. 28, 235–242. Web of Science CrossRef PubMed CAS Google Scholar
Blow, D. M. & Crick, F. H. C. (1959). Acta Cryst. 12, 794–802. CrossRef CAS IUCr Journals Web of Science Google Scholar
Booth, A. D. (1947a). Nature (London), 160, 196. CrossRef Web of Science Google Scholar
Booth, A. D. (1947b). J. Chem. Phys. 15, 415. CrossRef Web of Science Google Scholar
Booth, A. D. (1949). Proc. R. Soc. London Ser. A, 197, 336–355. CrossRef CAS Google Scholar
Brünger, A. T. (1990a). Acta Cryst. A46, 46–57. CrossRef Web of Science IUCr Journals Google Scholar
Brünger, A. T. (1990b). X-PLOR Manual. Version 2.1. Yale University, New Haven, USA. Google Scholar
Brünger, A. T. (1991). Acta Cryst. A47, 195–204. CrossRef Web of Science IUCr Journals Google Scholar
Brünger, A. T., Adams, P. D., Clore, G. M., DeLano, W. L., Gros, P., Grosse-Kunstleve, R. W., Jiang, J.-S., Kuszewski, J., Nilges, M., Pannu, N. S., Read, R. J., Rice, L. M., Simonson, T. & Warren, G. L. (1998). Acta Cryst. D54, 905–921. Web of Science CrossRef IUCr Journals Google Scholar
Castellano, E. E., Oliva, G. & Navaza, J. (1992). J. Appl. Cryst. 25, 281–284. CrossRef CAS Web of Science IUCr Journals Google Scholar
Cochran, W. (1948). Acta Cryst. 1, 273. CrossRef IUCr Journals Web of Science Google Scholar
Collaborative Computational Project, Number 4 (1994). Acta Cryst. D50, 760–763. CrossRef IUCr Journals Google Scholar
Derewenda, Z. S. (1989). Acta Cryst. A45, 227–234. CrossRef CAS Web of Science IUCr Journals Google Scholar
Driessen, H., Haneef, M. I. J., Harris, G. W., Howlin, B., Khan, G. & Moss, D. S. (1989). J. Appl. Cryst. 22, 510–516. CrossRef CAS Web of Science IUCr Journals Google Scholar
Grosse-Kunstleve, R. W. (1999). Acta Cryst. A55, 383–395. Web of Science CrossRef CAS IUCr Journals Google Scholar
Hendrickson, W. A. (1991). Science, 254, 51–58. CrossRef PubMed CAS Web of Science Google Scholar
Hoard, L. G. & Nordman, C. E. (1979). Acta Cryst. A35, 1010–1015. CrossRef CAS IUCr Journals Web of Science Google Scholar
Huber, R. & Schneider, M. (1985). J. Appl. Cryst. 18, 165–169. CrossRef CAS Web of Science IUCr Journals Google Scholar
Jiang, J.-S. & Brünger, A. T. (1994). J. Mol. Biol. 243, 100–155. CrossRef CAS PubMed Web of Science Google Scholar
Kostrewa, D. (1997). CCP4 Newsl. 34, 9–22. Google Scholar
Kronenburg, M. J. (2004). Acta Cryst. A60, 250–256. Web of Science CrossRef CAS IUCr Journals Google Scholar
Liu, D. C. & Nocedal, J. (1989). Math. Program. 45, 503–528. CrossRef Web of Science Google Scholar
Lebedev, A. A., Vagin, A. A. & Murshudov, G. N. (2008). Acta Cryst. D64, 33–39. Web of Science CrossRef IUCr Journals Google Scholar
Lunin, V. Y., Afonine, P. V. & Urzhumtsev, A. G. (2002). Acta Cryst. A58, 270–282. Web of Science CrossRef CAS IUCr Journals Google Scholar
Lunin, V. Yu. & Skovoroda, T. P. (1995). Acta Cryst. A51, 880–887. CrossRef CAS Web of Science IUCr Journals Google Scholar
McCoy, A. J. (2007). Acta Cryst. D63, 32–41. Web of Science CrossRef CAS IUCr Journals Google Scholar
Murshudov, G. N., Vagin, A. A. & Dodson, E. J. (1997). Acta Cryst. D53, 240–255. CrossRef CAS Web of Science IUCr Journals Google Scholar
Navaza, J. (2001). Acta Cryst. D57, 1367–1372. Web of Science CrossRef CAS IUCr Journals Google Scholar
Navaza, J., Lepault, J., Rey, F. A., Álvarez-Rúa, C. & Borge, J. (2002). Acta Cryst. D58, 1820–1825. Web of Science CrossRef CAS IUCr Journals Google Scholar
Noble, M. E. M., Zeelen, J. P., Wierenga, R. K., Mainfroid, V., Goraj, K., Gohimont, A.-C. & Martial, J. A. (1993). Acta Cryst. D49, 403–417. CrossRef CAS Web of Science IUCr Journals Google Scholar
Rossmann, M. G. & Arnold, E. (2001). Editors. International Tables for Crystallography, Vol. F, pp. 263–292. Dorderecht: Kluwer Academic Publishers. Google Scholar
Scheringer, C. (1963). Acta Cryst. 16, 546–550. CrossRef CAS IUCr Journals Web of Science Google Scholar
Sheldrick, G. M. (2008). Acta Cryst. A64, 112–122. Web of Science CrossRef CAS IUCr Journals Google Scholar
Sheldrick, G. M. & Schneider, T. R. (1997). Methods Enzymol. 277, 319–343. CrossRef PubMed CAS Web of Science Google Scholar
Sussman, J. L., Holbrook, S. R., Church, G. M. & Kim, S.-H. (1977). Acta Cryst. A33, 800–804. CrossRef CAS IUCr Journals Web of Science Google Scholar
Tronrud, D. E. (2004). Acta Cryst. D60, 2156–2168. Web of Science CrossRef CAS IUCr Journals Google Scholar
Urzhumtsev, A. G., Lunin, V. Yu. & Vernoslova, E. A. (1989). J. Appl. Cryst. 22, 500–506. CrossRef CAS Web of Science IUCr Journals Google Scholar
Urzhumtseva, L. M. & Urzhumtsev, A. G. (1997). J. Appl. Cryst. 30, 402–410. CrossRef CAS Web of Science IUCr Journals Google Scholar
Weisstein, E. W. (2006). Euler Angles, https://mathworld.wolfram.com/EulerAngles.html . Google Scholar
Yeates, T. O. & Rees, D. C. (1988). J. Appl. Cryst. 21, 925–928. CrossRef CAS Web of Science IUCr Journals Google Scholar
Yeates, T. O. & Rini, J. M. (1990). Acta Cryst. A46, 352–359. CrossRef CAS Web of Science IUCr Journals Google Scholar
This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.