International
Tables for
Crystallography
Volume B
Reciprocal space
Edited by U. Shmueli

International Tables for Crystallography (2010). Vol. B, ch. 2.3, pp. 269-272   | 1 | 2 |

## Section 2.3.7. Translation functions

L. Tong,c* M. G. Rossmanna and E. Arnoldb

aDepartment of Biological Sciences, Purdue University, West Lafayette, Indiana 47907, USA,bCABM & Rutgers University, 679 Hoes Lane, Piscataway, New Jersey 08854–5638, USA, and cDepartment of Biological Sciences, Columbia University, New York 10027, USA
Correspondence e-mail:  ltong@columbia.edu

### 2.3.7. Translation functions

| top | pdf |

#### 2.3.7.1. Introduction

| top | pdf |

The problem of determining the position of a noncrystallographic symmetry element in space, or the position of a molecule of known orientation in a unit cell, has been reviewed by Rossmann (1972), Colman et al. (1976), Karle (1976), Argos & Rossmann (1980), Harada et al. (1981) and Beurskens (1981). All methods depend on the prior knowledge of the object's orientation implied by the rotation matrix [C]. The various translation functions, T, derived below, can only be computed given this information.

The general translation function can be defined as where T is a six-variable function given by each of the three components that define and . Here and are equivalent reference positions of the objects, whose densities are and . The translation function searches for the optimal overlap of the two objects after they have been similarly oriented. Following the same procedure used for the rotation-function derivation, Fourier summations are substituted for and . It can then be shown that

Using the substitution and simplifying leads to The integral is the diffraction function (2.3.6.4). If the integration is taken over the volume U, centred at and , it follows that

#### 2.3.7.2. Position of a noncrystallographic element relating two unknown structures

| top | pdf |

The function (2.3.7.1) is quite general. For instance, the rotation function corresponds to a comparison of Patterson functions and at their origins. That is, the coefficients are , phases are zero and . However, the determination of the translation between two objects requires the comparison of cross-vectors away from the origin.

Consider, for instance, the determination of the precise translation vector parallel to a rotation axis between two identical molecules of unknown structure. For simplicity, let the noncrystallographic axis be a dyad (Fig. 2.3.7.1). Fig. 2.3.7.2 shows the corresponding Patterson of the hypothetical point-atom structure. Opposite sets of cross-Patterson vectors in Fig. 2.3.7.2 are related by a twofold rotation and a translation equal to twice the precise vector in the original structure. A suitable translation function would then compare a Patterson at S with the rotated Patterson at . Hence, substituting and in (2.3.7.1),

 Figure 2.3.7.1 | top | pdf |Crosses represent atoms in a two-dimensional model structure. The triangles are the points chosen as approximate centres of molecules A and B. has components t and s parallel and perpendicular, respectively, to the screw rotation axis. [Reprinted from Rossmann et al. (1964).]
 Figure 2.3.7.2 | top | pdf |Vectors arising from the structure in Fig. 2.3.7.1. The self-vectors of molecules A and B are represented by + and ·; the cross-vectors from molecules A to B and B to A by × and ○. Triangles mark the position of and . [Reprinted from Rossmann et al. (1964).]

The opposite cross-vectors can be superimposed only if an evenfold rotation between the unknown molecules exists. The translation function (2.3.7.2) is thus applicable only in this special situation. There is no published translation method to determine the interrelation of two unknown structures in a crystallographic asymmetric unit or in two different crystal forms. However, another special situation exists if a molecular evenfold axis is parallel to a crystallographic evenfold axis. In this case, the position of the noncrystallographic symmetry element can be easily determined from the large peak in the corresponding Harker section of the Patterson.

In general, it is difficult or impossible to determine the positions of noncrystallographic axes (or their intersection at a molecular centre). However, the position of heavy atoms in isomorphous derivatives, which usually obey the noncrystallographic symmetry, can often determine this information.

#### 2.3.7.3. Position of a known molecular structure in an unknown unit cell

| top | pdf |

The most common type of translation function occurs when looking for the position of a known molecular structure in an unknown crystal. For instance, if the structure of an enzyme has previously been determined by the isomorphous replacement method, then the structure of the same enzyme from another species can often be solved by molecular replacement [e.g. Grau et al. (1981)]. However, there are some severe pitfalls when, for instance, there are gross conformational changes [e.g. Moras et al. (1980)]. This type of translation function could also be useful in the interpolation of E maps produced by direct methods. Here there may often be confusion as a consequence of a number of molecular images related by translations (Karle, 1976; Beurskens, 1981; Egert & Sheldrick, 1985).

Tollin's (1966) Q function and Crowther & Blow's (1967) translation function are essentially identical (Tollin, 1969) and depend on a prior knowledge of the search molecule as well as its orientation in the unknown cell. The derivation given here, however, is somewhat more general and follows the derivation of Argos & Rossmann (1980), and should be compared with the method of Harada et al. (1981).

If the known molecular structure is correctly oriented into a cell (p) of an unknown structure and placed at S with respect to a defined origin, then a suitable translation function is This definition is preferable to one based on an R-factor calculation as it is more amenable to computation and is independent of a relative scale factor.

The structure factor can be calculated by modifying expression (2.3.8.9) (see below). That is, where is the volume of cell (h) and is the position, in the nth crystallographic asymmetric unit, of cell (p) corresponding to S in known cell (h). Let which are the coefficients of the molecular transform for the known molecule placed into the nth asymmetric unit of the p cell. Thus or where and . Hence and then from (2.3.7.3) which is a Fourier summation with known coefficients such that T(S) will be a maximum at the correct molecular position.

Terms with in expression (2.3.7.4) can be omitted as they are independent of S and only contribute a constant to the value of T(S). For terms with , the indices take on special values. For instance, if the p cell is monoclinic with its unique axis parallel to b such that and , then would be (2p, 0, 2r). Hence, T(S) would be a two-dimensional function consistent with the physical requirement that the translation component, parallel to the twofold monoclinic axis, is arbitrary.

Crowther & Blow (1967) show that if are the structure factors of a known molecule correctly oriented within the cell of the unknown structure at an arbitrary molecular origin, then (altering the notation very slightly from above) where [C] is a crystallographic symmetry operator relative to which the molecular origin is to be determined. This is of the same form as (2.3.7.4) but concerns the special case where the h cell, into which the known molecule was placed, has the same dimensions as the p cell.

The translation function as defined by (2.3.7.4) is on an arbitrary scale, which makes it difficult to compare results from different calculations. Translation functions can also be defined based on the crystallographic R factor or a correlation coefficient (CC). In particular, CCs based on reflection intensities can be evaluated by Fourier methods (Navaza & Vernoslova, 1995), although it is still computationally more expensive than the evaluation of (2.3.7.4). Alternatively, the translation function can be calculated first with (2.3.7.4), and then the R factor and CC can be calculated for the resulting top solutions.

A correct solution should also produce satisfactory packing arrangements of the molecular models in the crystal. Packing functions have been derived that estimate the amount of overlap among the models (Harada et al., 1981; Hendrickson & Ward, 1976; Rabinovich & Shakked, 1984; Simpson et al., 2001), and such considerations can frequently limit the search volume very considerably. Alternatively, a simple enumeration of the actual close contacts among different molecules in the crystal (for example, Cα–Cα distances less than 3 Å) has also been found to be an effective way of eliminating those solutions that produce unreasonable crystal packing (Jogl et al., 2001; Tong, 1993). If conformational differences are expected between the search atomic model and the actual structure, care must be taken when applying this packing check.

#### 2.3.7.4. Position of a noncrystallographic symmetry element in a poorly defined electron-density map

| top | pdf |

If an initial set of poor phases, for example from an SIR derivative, are available and the rotation function has given the orientation of a noncrystallographic rotation axis, it is possible to search the electron-density map systematically to determine the translation axis position. The translation function must, therefore, measure the quality of superposition of the poor electron-density map on itself. Hence and the function (2.3.7.1) now becomes This real-space translation function has been used successfully to determine the intermolecular dyad axis for α-chymotrypsin (Blow et al., 1964) and to verify the position of immunoglobulin domains (Colman & Fehlhammer, 1976).

#### 2.3.7.5. Locked translation function

| top | pdf |

In a translation search, an atomic model with a given orientation is moved systematically through the unit cell. In such a situation, the structure-factor equation takes on the special form (Harada et al., 1981; Rae, 1977; Tong, 1993)where S is the translation vector and the summation goes over the crystallographic symmetry operators. is the structure factor calculated based only on the nth symmetry-related molecule,where represents the atomic position of the model at the reference position and the summation goes over all the atoms.

Noting equation (2.3.7.3), the translation function is given bywhere the second term is the ordinary translation function, analogous to (2.3.7.4). The first term of (2.3.7.5) depends on the orientation of the model. Maximization of this term, or its correlation coefficient equivalent, is the basis behind the Patterson-correlation refinement (Brünger, 1990; Tong, 1996b) and the direct rotation function (DeLano & Brünger, 1995). It is also related to the intensity-based domain refinement (Yeates & Rini, 1990).

In the presence of noncrystallographic symmetry, the locked self rotation function can be used to define the orientation of the noncrystallographic symmetry point group in the crystal. If an atomic model is available for the monomer but not for the entire oligomer, the locked cross rotation function can be used to determine the orientation of this monomer in the oligomer. The locked translation function can then be used to determine the position of this monomer relative to the centre of the noncrystallographic symmetry point group (Tong, 1996b, 2001a), which will produce a model for the entire oligomer. The centre of this oligomer in the crystal can be defined by a simple translation search.

With the knowledge of the orientation of one monomer of the oligomer, the first term of (2.3.7.5) is dependent on the position of this monomer relative to the centre of the noncrystallographic symmetry oligomer (Tong, 1996b). The atomic positions of the entire noncrystallographic symmetry oligomer in the standard orientation are given bywhere are the atomic positions of the monomer model, centred at (0, 0, 0); [F] is the orientation of this model in the oligomer in the standard orientation; V0 is the position of this monomer relative to the centre of the oligomer; and [In] is the nth noncrystallographic symmetry rotation matrix in the standard orientation. The atomic positions of the noncrystallographic symmetry oligomer in the crystal unit cell, centred at the origin, are given bywhere [E] is the orientation of the noncrystallographic symmetry in the crystal unit cell and is the deorthogonalization matrix.

By incorporating the calculated structure factors based on this noncrystallographic symmetry oligomer into the first term of (2.3.7.5), the locked translation function is given bywhere and . A constant term has been omitted from this equation.

Conceptually, the locked translation function is based on the overlap of intermolecular vectors within the noncrystallographic symmetry oligomer and the observed Patterson map (Tong, 1996b). The equation for the locked translation function, (2.3.7.6), bears remarkable resemblance to that for the ordinary Patterson-correlation translation function, (2.3.7.5), with the interchange of the crystallographic ([Tn]) and noncrystallographic symmetry parameters.

#### 2.3.7.6. Computer programs for rotation and translation function calculations

| top | pdf |

Several programs are currently in popular use for the calculation of rotation and translation functions. These include AMoRe (Navaza, 1994, 2001a), BEAST (Read, 2001b), CCP4 (Collaborative Computational Project, Number 4, 1994), CNS (Brünger et al., 1998), COMO (Jogl et al., 2001), EPMR (Kissinger et al., 1999), GLRF (part of the Replace package) (Tong, 1993, 2001a; Tong & Rossmann, 1990, 1997), Molrep (Vagin & Teplyakov, 2000) and Phaser (Storoni et al., 2004).

The correct placement of an atomic model in a crystal unit cell is generally a six-dimensional problem, with three degrees of rotational freedom and three degrees of translational freedom. Systematic examination of all six degrees of freedom at the same time is computationally expensive and cannot be used routinely (Fujinaga & Read, 1987; Rabinovich & Shakked, 1984; Sheriff et al., 1999). On the other hand, directed sampling of the six degrees of freedom, driven by a stochastic or genetic algorithm (Chang & Lewis, 1997; Glykos & Kokkinidis, 2000; Kissinger et al., 1999), has been successful in solving structures.

Traditionally, the calculations are divided into a rotational component (the rotation function) and a translational component (the translation function). Only a few rotation angles (for example the top few peaks of the rotation function) are manually passed to the translation function for examination (Fitzgerald, 1988). With the power of modern computers, it is now possible to perform limited six-dimensional searches, with the sampling of the rotational degrees of freedom guided by the rotation function. For example, the top peaks of the rotation function (Navaza, 1994) and their neighbours (Urzhumtsev & Podjarny, 1995) can be automatically examined by the translation function. A more general approach is to examine all rotation-function grid points with values greater than a certain threshold (Tong, 1996a). Such combined molecular replacement protocols have been found to be very powerful in solving new structures.

### References

Collaborative Computational Project, Number 4 (1994). The CCP4 suite: programs for protein crystallography. Acta Cryst. D50, 760–763.
Read, R. J. (2001b). Pushing the boundaries of molecular replacement with maximum likelihood. Acta Cryst. D57, 1373–1382.
Tong, L. (2001a). How to take advantage of non-crystallographic symmetry in molecular replacement: `locked' rotation and translation functions. Acta Cryst. D57, 1383–1389.
Argos, P. & Rossmann, M. G. (1980). Molecular replacement methods. In Theory and Practice of Direct Methods in Crystallography, edited by M. F. C. Ladd & R. A. Palmer, pp. 361–417. New York: Plenum.
Beurskens, P. T. (1981). A statistical interpretation of rotation and translation functions in reciprocal space. Acta Cryst. A37, 426–430.
Blow, D. M., Rossmann, M. G. & Jeffery, B. A. (1964). The arrangement of α-chymotrypsin molecules in the monoclinic crystal form. J. Mol. Biol. 8, 65–78.
Brünger, A. T. (1990). Extension of molecular replacement: a new search strategy based on Patterson correlation refinement. Acta Cryst. A46, 46–57.
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). Crystallography & NMR System: a new software suite for macromolecular structure determination. Acta Cryst. D54, 905–921.
Chang, G. & Lewis, M. (1997). Molecular replacement using genetic algorithms. Acta Cryst. D53, 279–289.
Colman, P. M. & Fehlhammer, H. (1976). Appendix: the use of rotation and translation functions in the interpretation of low resolution electron density maps. J. Mol. Biol. 100, 278–282.
Colman, P. M., Fehlhammer, H. & Bartels, K. (1976). Patterson search methods in protein structure determination: β-trypsin and immunoglobulin fragments. In Crystallographic Computing Techniques, edited by F. R. Ahmed, K. Huml & B. Sedlacek, pp. 248–258. Copenhagen: Munksgaard.
Crowther, R. A. & Blow, D. M. (1967). A method of positioning a known molecule in an unknown crystal structure. Acta Cryst. 23, 544–548.
DeLano, W. L. & Brünger, A. T. (1995). The direct rotation function: rotational Patterson correlation search applied to molecular replacement. Acta Cryst. D51, 740–748.
Egert, E. & Sheldrick, G. M. (1985). Search for a fragment of known geometry by integrated Patterson and direct methods. Acta Cryst. A41, 262–268.
Fitzgerald, P. M. D. (1988). MERLOT, an integrated package of computer programs for the determination of crystal structures by molecular replacement. J. Appl. Cryst. 21, 273–278.
Fujinaga, M. & Read, R. J. (1987). Experiences with a new translation-function program. J. Appl. Cryst. 20, 517–521.
Glykos, N. M. & Kokkinidis, M. (2000). A stochastic approach to molecular replacement. Acta Cryst. D56, 169–174.
Grau, U. M., Rossmann, M. G. & Trommer, W. E. (1981). The crystallization and structure determination of an active ternary complex of pig heart lactate dehydrogenase. Acta Cryst. B37, 2019–2026.
Harada, Y., Lifchitz, A., Berthou, J. & Jolles, P. (1981). A translation function combining packing and diffraction information: an application to lysozyme (high-temperature form). Acta Cryst. A37, 398–406.
Hendrickson, W. A. & Ward, K. B. (1976). A packing function for delimiting the allowable locations of crystallized macromolecules. Acta Cryst. A32, 778–780.
Jogl, G., Tao, X., Xu, Y. & Tong, L. (2001). COMO: a program for combined molecular replacement. Acta Cryst. D57, 1127–1134.
Karle, J. (1976). Partial structures and use of the tangent formula and translation functions. In Crystallographic Computing Techniques, edited by F. R. Ahmed, K. Huml & B. Sedlacek, pp. 155–164. Copenhagen: Munksgaard.
Kissinger, C. R., Gehlhaar, D. K. & Fogel, D. B. (1999). Rapid automated molecular replacement by evolutionary search. Acta Cryst. D55, 484–491.
Moras, D., Comarmond, M. B., Fischer, J., Weiss, R., Thierry, J. C., Ebel, J. P. & Giegé, R. (1980). Crystal structure of yeast tRNAAsp. Nature (London), 288, 669–674.
Navaza, J. (1994). AMoRe: an automated package for molecular replacement. Acta Cryst. A50, 157–163.
Navaza, J. (2001a). Implementation of molecular replacement in AMoRe. Acta Cryst. D57, 1367–1372.
Navaza, J. & Vernoslova, E. (1995). On the fast translation functions for molecular replacement. Acta Cryst. A51, 445–449.
Rabinovich, D. & Shakked, Z. (1984). A new approach to structure determination of large molecules by multi-dimensional search methods. Acta Cryst. A40, 195–200.
Rae, A. D. (1977). The use of structure factors to find the origin of an oriented molecular fragment. Acta Cryst. A33, 423–425.
Rossmann, M. G. (1972). The Molecular Replacement Method. New York: Gordon & Breach.
Sheriff, S., Klei, H. E. & Davis, M. E. (1999). Implementation of a six-dimensional search using the AMoRe translation function for difficult molecular-replacement problems. J. Appl. Cryst. 32, 98–101.
Simpson, A. A., Leiman, P. G., Tao, Y., He, Y., Badasso, M. O., Jardine, P. J., Anderson, D. L. & Rossmann, M. G. (2001). Structure determination of the head-tail connector of bacteriophage ϕ29. Acta Cryst. D57, 1260–1269.
Storoni, L. C., McCoy, A. J. & Read, R. J. (2004). Likelihood-enhanced fast rotation functions. Acta Cryst. D60, 432–438.
Tollin, P. (1966). On the determination of molecular location. Acta Cryst. 21, 613–614.
Tollin, P. (1969). A comparison of the Q-functions and the translation function of Crowther and Blow. Acta Cryst. A25, 376–377.
Tong, L. (1993). REPLACE, a suite of computer programs for molecular-replacement calculations. J. Appl. Cryst. 26, 748–751.
Tong, L. (1996a). Combined molecular replacement. Acta Cryst. A52, 782–784.
Tong, L. (1996b). The locked translation function and other applications of a Patterson correlation function. Acta Cryst. A52, 476–479.
Tong, L. & Rossmann, M. G. (1990). The locked rotation function. Acta Cryst. A46, 783–792.
Tong, L. & Rossmann, M. G. (1997). Rotation function calculations with GLRF program. Methods Enzymol. 276, 594–611.
Urzhumtsev, A. & Podjarny, A. (1995). On the solution of the molecular-replacement problem at very low resolution: Application to large complexes. Acta Cryst. D51, 888–895.
Vagin, A. & Teplyakov, A. (2000). An approach to multi-copy search in molecular replacement. Acta Cryst. D56, 1622–1624.
Yeates, T. O. & Rini, J. M. (1990). Intensity-based domain refinement of oriented but unpositioned molecular replacement models. Acta Cryst. A46, 352–359.