Tables for
Volume F
Crystallography of biological macromolecules
Edited by E. Arnold, D. M. Himmel and M. G. Rossmann

International Tables for Crystallography (2012). Vol. F, ch. 13.3, pp. 347-351   | 1 | 2 |

Chapter 13.3. Translation functions

L. Tonga*

aDepartment of Biological Sciences, Columbia University, New York, NY 10027, USA
Correspondence e-mail:

Translation functions are used to determine the positions of search models by the molecular-replacement method. The functions are generally based on the comparisons of observed and calculated structure-factor amplitudes (R factor and correlation coefficient) or the comparisons of the corresponding Patterson functions. Combined molecular replacement is a limited six-dimensional search where sampling of the rotational degrees of freedom is guided by the rotation function. Evaluating the packing arrangement of the model in the crystal is an important part of the molecular-replacement calculation. The locked translation function can be used in cases where the crystal contains point-group noncrystallographic symmetry.

13.3.1. Introduction

| top | pdf |

A structure determination by the molecular-replacement method traditionally proceeds in two steps (Rossmann, 1972[link], 1990[link], 2001[link]). The first step involves the determination of the orientation of the search model in the crystal unit cell by the rotation functions (see Chapter 13.2[link] ). Once the orientation of the search model is known, translation functions are employed in the second step to determine the location of the model in the crystal unit cell. This essentially reduces a six-dimensional problem (three rotational and three translational degrees of freedom) to two three-dimensional problems, which are computationally more manageable. With the speed of modern computers, a strict division between the rotational and the translational components of a molecular-replacement structure solution may no longer be necessary (see Section 13.3.7[link]).

Translation functions are normally formulated to achieve minimum or maximum values when the search model is at its correct position in the crystal unit cell. As with the rotation problem, the translation problem is frequently solved as a search. The positional parameters of the model are varied in the unit cell, generally on a grid. Translation functions are evaluated at these search grid points in order to identify those that minimize or maximize the functions.

Most translation functions involve a comparison between the observed structure-factor amplitudes (or squared amplitudes) and those calculated based on the search model. The R factor and the correlation coefficient can be used as indicators for translation searches (Section 13.3.2[link]). The correlation between the observed Patterson map and that based on the search model is the foundation of another translation function (Section 13.3.3[link]). If phase information is available from other sources, the correlation between the electron-density maps is the basis for phased translation functions (Section 13.3.4[link]). Proper packing of the search model in the crystal unit cell is an essential component of a solution to the translation problem (Section 13.3.5[link]). The power of the translation functions can be enhanced in the presence of noncrystallographic symmetry (Section 13.3.8[link]).

Besides R-factor, correlation-coefficient and Patterson-correlation functions, target functions based on the principle of maximum likelihood have been found to be powerful indicators for molecular-replacement calculations (Read, 2001[link]; McCoy et al., 2005[link], 2007[link]), and have been successful in helping determine structures where traditional molecular-replacement target functions have failed.

13.3.2. R-factor and correlation-coefficient translation functions

| top | pdf |

The crystallographic R factor is often used as an indicator in translation searches. It is a measure of the percentage difference between the observed ([F_{h}^{o}]) and the calculated ([F_{h}^{c}]) structure-factor amplitudes, [R_{F} = 100 \times {\textstyle\sum\limits_{h}} \left|F_{h}^{o} - k_{F} F_{h}^{c}\right|\big/{\textstyle\sum\limits_{h}} F_{h}^{o}. \eqno (]A similar R factor can be defined based on the square of the structure-factor amplitudes, i.e. an R factor based on intensity, [R_{I} = 100 \times {\textstyle\sum\limits_{h}} \left|I_{h}^{o} - k_{I} I_{h}^{c}\right|\big/2{\textstyle\sum\limits_{h}} I_{h}^{o}. \eqno (]A factor of 2 is introduced in equation ([link]) to make [R_{I}] values fall in the same range as [R_{F}]. In equations ([link]) and ([link]) [k_{F}] and [k_{I}] are scale factors that bring the observed and the calculated structure factors to the same level. If these scale factors are calculated in shells of equal reciprocal volume, they can compensate for differences in the displacement factors between the observed and the calculated structure-factor amplitudes.

In addition to the R factor, the correlation coefficient between the observed and the calculated structure factors is also used in translation searches. Like the R factors, correlation coefficients can be defined based on the amplitude or the intensity of the reflections, and that based on the amplitude is shown below. [\eqalign{cc_{F} &= {{\textstyle\sum_{h}} \left(F_{h}^{o} - \left\langle F_{h}^{o} \right\rangle \right) \left(F_{h}^{c} - \left\langle F_{h}^{c} \right\rangle \right) \over \left[{\textstyle\sum_{h}} \left(F_{h}^{o} - \left\langle F_{h}^{o} \right\rangle \right)^{2} \left(F_{h}^{c} - \left\langle F_{h}^{c} \right\rangle \right)^{2}\right]^{1/2}}\cr &= {{\textstyle\sum_{h}} F_{h}^{o} F_{h}^{c} - {\textstyle\sum_{h}} F_{h}^{o} {\textstyle\sum_{h}} F_{h}^{c} / N \over \left\{\left[{\textstyle\sum_{h}} \left(F_{h}^{o}\right)^{2} - \left({\textstyle\sum_{h}} F_{h}^{o}\right)^{2} \big / N \right] \left[{\textstyle\sum_{h}} \left(F_{h}^{c}\right)^{2} - \left({\textstyle\sum_{h}} F_{h}^{c}\right)^{2} \big / N \right]\right\}^{1/2}}.} \eqno (]In ([link]), N is the number of reflections that are used in the calculation and [\left\langle F_{h} \right\rangle] denotes the average structure-factor amplitude over the reflections. Unlike the R factors, the correlation coefficients do not depend on the overall scale factor between the observed and the calculated structure factors. However, they can be affected by large differences in the overall displacement factors between the observed and the calculated structure factors. The correlation coefficient based on intensity can be calculated by the fast Fourier transform (FFT) technique for translation-function calculations (Navaza & Vernoslova, 1995[link]).

In order to evaluate R factors and correlation coefficients for a translation search, structure factors need to be calculated for the search model, with a given orientation, at different positions in the unit cell. For this special case, where only the positional parameters of the search model are varied, the calculation of the structure factors can be simplified (Nixon & North, 1976[link]; Rae, 1977[link]). The structure-factor equation can be written as a double summation – first over the atoms in one asymmetric unit of the unit cell and then over all the asymmetric units, [\bar{F}_{h}^{c} = {\textstyle\sum\limits_{n}} {\textstyle\sum\limits_{j}}\ f_{j} \exp \left\{2\pi ih\left(\left[T_{n}\right]x_j + t_{n}\right)\right\}, \eqno (]where j goes over all the atoms in the asymmetric unit and n goes over all the crystallographic asymmetric units. The nth crystallographic symmetry operator is given by [x_{n} = \left[T_{n}\right]x + t_{n}, \eqno (]where [\left[T_{n}\right]] is the rotational component and [t_{n}] is the translational component of the symmetry operator.

For simplicity, first consider the case where there is only one molecule in the asymmetric unit. In the translation search, the model will be placed at different positions in the unit cell, [x_{j} = x_{j}^{0} + v_{0}, \eqno (]where [v_{0}] is a translation vector which is applied to move the model from its starting position ([x_{j}^{0}]). Substituting [x_{j}] in ([link]) into ([link]) gives [\bar{F}_{h}^{c} = {\textstyle\sum\limits_{n}}\ \overline{f}_{h,  n} \exp \left(2\pi ih \left[T_{n}\right]v_{0}\right), \eqno (]where [\overline{f}_{h, n}] is the structure factor calculated based only on the nth symmetry-related molecule, [\overline{f}_{h,  n} = {\textstyle\sum\limits_{j}}\ f_{j} \exp \left\{2\pi ih \left(\left[T_{n}\right]x_{j}^{0} + t_{n}\right)\right\}. \eqno (]It can be calculated by placing the search model at the starting position in a P1 unit cell having the same cell dimensions as the crystal unit cell. The structure factors calculated for this P1 cell are related to [\overline{f}_{h,  n}] by [\overline{f}_{h, n} = \exp (2\pi iht_{n}) \overline{f}_{h[T_{n}], \, 1}. \eqno (]

Therefore, the summation over the atoms in the structure-factor calculation, a rather time-consuming process, needs to be performed only once, for the search model at the starting position. Subsequent structure-factor calculations after translation of the model are no longer dependent on the number of atoms present in the unit cell [equation ([link])]. The starting position is usually chosen such that the centre of the search model is at (0, 0, 0). Then the vector that is determined from the translation searches will define the centre of the model in the unit cell.

Equation ([link]) can be generalized to allow for the presence of other molecules that are to remain stationary during the translation search: [\bar{F}_{h}^{c} = \bar{A}_{h} + {\textstyle\sum\limits_{n}}\ \overline{f}_{h,  n} \exp(2\pi ih[T_{n}]v_{0}), \eqno (]where [\bar{A}_{h}] is the contribution from the stationary molecules. This formulation is useful if there is more than one molecule in the asymmetric unit. The position of one of the molecules can be determined first, and the model is then included as a stationary molecule for the position search of the next molecule. Incorporation of the contribution from stationary molecules can also improve the quality of the translation function (Zhang & Matthews, 1994[link]).

Equation ([link]) can also be generalized to allow for the presence of two (or more) search models that are to move independently of each other during the translation search. However, this will generally lead to a six- (or more) dimensional problem and is expensive in computation time. With improvements in computer technology (especially parallel processing), it may be feasible to carry out such searches in special cases. However, this aspect of translation functions will not be discussed further here.

13.3.3. Patterson-correlation translation function

| top | pdf |

Besides the R factor and the correlation coefficient, another commonly used translation search indicator is based on the correlation between the observed and the calculated Patterson maps (Crowther & Blow, 1967[link]; Tong, 1993[link]; Navaza & Vernoslova, 1995[link]). Rotation functions are based on the overlap of only a subset of the interatomic vectors in the Patterson map, i.e. only those near the origin of the unit cell, which generally contain the self vectors within each molecule. The correct orientation and position of a search model in the crystal unit cell should lead to the maximal overlap of both the self and the cross vectors, i.e. maximal overlap between the observed ([P^o]) and the calculated ([P^c)] Patterson maps throughout the entire unit cell (Tong, 1993[link]), [\eqalign{PC &= {\textstyle\int\limits_{\Omega}} P^{o} \left(u\right)P^{c} \left(u\right)\ \hbox{d}u\cr &= {\textstyle\sum\limits_{h}} \left(F_{h}^{o}\right)^{2} |\bar{F}_{h}^{c}|^{2}.} \eqno (]

The calculated structure factor is a function of the translation vector [v_{0}] [equation ([link])]. Combining equations ([link]) and ([link]) gives[\eqalign{PCTF\left(v_{0}\right) &= {\textstyle\sum\limits_{h}} \left(F_{h}^{o}\right)^{2} \left|A_{h}\right|^{2} + {\textstyle\sum\limits_{h}} {\textstyle\sum\limits_{n}} \left(F_{h}^{o}\right)^{2}\left|\,\overline{f}_{h,  n}\right|^{2}\cr &\quad + 2{\textstyle\sum\limits_{h}} {\textstyle\sum\limits_{n}} \left(F_{h}^{o}\right)^{2} \bar{A}_{h}\, \overline{f}_{h,  n}^{*} \exp(-2\pi ih [T_{n}]v_{0})\cr &\quad + {\textstyle\sum\limits_{h}} {\textstyle\sum\limits_{n}} {\textstyle\sum\limits_{m \neq n}} \left(F_{h}^{o}\right)^{2} \overline{f}_{h,  n}\, \overline{f}_{h,  m}^{*} \exp\left\{- 2\pi ih \left(\left[T_{m}\right] - \left[T_{n}\right]\right)v_{0}\right\}.} \eqno (]The first two terms in equation ([link]) contribute a constant to the correlation and are generally ignored in this calculation (but see Section 13.3.8[link]). The Patterson-correlation translation function is therefore a Fourier transform, with [h[T_{n}]] (from the third term) and [h([T_{m}] - [T_{n}])] (from the fourth term) as the indices. This function can be evaluated quickly by the FFT technique. The term [h([T_{m}] - [T_{n}])] often leads to the doubling of the original reflection indices. For example, [h([T_{1}] - [T_{2}])] will give rise to (2h, 0, 2l) for monoclinic space groups (b-unique setting). Therefore, the Patterson-correlation Fourier transform should normally be sampled with a grid size about 1/6 of the maximum resolution of the reflection data used in the calculation.

A disadvantage of the Patterson-correlation translation function as formulated in equation ([link]) is that the results from the calculation are on an arbitrary scale. This makes it difficult to compare results from different calculations. The R factor or the correlation coefficient can be calculated for the top peaks in the Patterson-correlation translation function [equation ([link])] to place the results on an `absolute' scale. Alternatively, the correlation, as defined by equation ([link]), can be normalized to become a Patterson-correlation coefficient (Harada et al., 1981[link]), [PC = {{\textstyle\sum_{h}} \left(F_{h}^{o}\right)^{2} \left|\bar{F}_{h}^{c}\right|^{2} \over \left[{\textstyle\sum_{h}} \left(F_{h}^{o}\right)^{4} {\textstyle\sum_{h}} \left|\bar{F}_{h}^{c}\right|^4\right]^{1/2}}. \eqno (]This correlation coefficient is equivalent to that defined by equation ([link]) for reflection intensities (Fujinaga & Read, 1987[link]). The evaluation of the [({\textstyle\sum_{h}} |\bar{F}_{h}^{c}|^{4})^{1/2}] term in the denominator of equation ([link]) is, however, time consuming, since the FFT technique cannot be applied directly. This term can be approximated by an expression which also measures the packing of the search molecules in the crystal (see Section 13.3.5[link]).

The intramolecular vectors can be removed (Crowther & Blow, 1967[link]) from the Patterson maps by subtracting, after appropriate scaling, the structure factors calculated from individual molecules [equation ([link])].

The translation functions, as defined above, are based on the structure-factor amplitudes. Normalized structure factors (E's) may provide better results under certain circumstances, since they increase the weight of the high-resolution reflections in the translation function (Tollin, 1966[link]; Harada et al., 1981[link]; Tickle, 1985[link]). Since the Patterson-correlation translation function is also based on reflection intensities, the `large term' approach can be used to accelerate the calculation (Tollin & Rossmann, 1966[link]).

13.3.4. Phased translation function

| top | pdf |

If an atomic model needs to be placed in an electron-density map that has been obtained through other methods (e.g. the multiple-isomorphous-replacement method or partial model phases), the phased translation function can be used (Colman et al., 1976[link]; Read & Schierbeek, 1988[link]; Bentley & Houdusse, 1992[link]; Tong, 1993[link]; Zhang & Matthews, 1994[link]). It is essentially based on the correlation between observed and calculated electron-density values throughout the unit cell: [PTF\left(v_{0}\right) = {\textstyle\sum\limits_{h}} {\textstyle\sum\limits_{n}}\ \bar{F}_{h}^{o}\, \bar{f}_{h,  n}^{*} \exp\left(- 2\pi ih\left[T_{n}\right]v_{0}\right). \eqno (]Therefore, the phased translation function can also be evaluated by the FFT technique. As with the Patterson-correlation translation function, the phased translation function can be placed on an `absolute' scale by introducing appropriate normalizing factors or by converting the results to R factors or correlation coefficients. It should be noted that the prior phase could be in the wrong hand, so both [\bar{F}_{h}^{o}] and [(\bar{F}_{h}^{o})^{*}] may need to be tried in the phased translation function, especially if the phase information is derived from isomorphous replacement.

The stationary molecules contribute a constant to the phased translation function and are not shown in equation ([link]). However, the phase information from the stationary molecules can be applied to the observed structure-factor amplitudes, and the phased translation function, rather than the Patterson-correlation translation function, can be used in the search for additional molecules (Read & Schierbeek, 1988[link]; Driessen et al., 1991[link]; Bentley & Houdusse, 1992[link]). This could prove especially useful in locating the last few molecules in cases where there are several molecules in the asymmetric unit. Modified phased translation functions have been designed to locate fragments in the electron-density map (Cowtan, 1998[link]).

13.3.5. Packing check in translation functions

| top | pdf |

A correct molecular-replacement solution should lead to the placement of the search model at the correct orientation and position in the crystal unit cell. For this solution, there should be no or minimal steric clashes among the crystallographically related and noncrystallographically related molecules in the unit cell (unless there are large conformational changes for the model in the crystal unit cell). Therefore, proper packing of the search model in the crystal unit cell is an important component of the molecular-replacement structure solution.

The packing of the search model in the unit cell can be estimated by determining the electron-density overlap among the molecules. This overlap can be calculated numerically, given the molecular envelope (Hendrickson & Ward, 1976[link]). It can also be estimated by an analytical function (Harada et al., 1981[link]), [O\left(v_{0}\right) = {\textstyle\sum\limits_{h}}\left|\bar{F}_{h}^{c} \left(v_{0}\right)\right|^{2}\Big/N{\textstyle\sum\limits_{h}}\left|\,\bar{f}_{h,  1}\right|^{2}, \eqno (]where N is the number of crystallographic symmetry operators. This overlap function assumes a value of 1 when there is no overlap among the molecules, and higher values when there is overlap. This function has been used to replace the [({\textstyle\sum_{h}}|\bar{F}_{h}^{c}|^{4})^{1/2}] term in the denominator of equation ([link]) (Harada et al., 1981[link]). Consequently, those positions that lead to steric clashes among the molecules will be down-weighted, thereby increasing the signal for the correct solution.

The overlap functions provide an overall estimate for the packing of the search model in the unit cell. A more detailed packing analysis can be based on the checking of atomic contacts (Tong, 1993[link]). For example, the number of [\hbox{C}_{\alpha} \cdots \hbox{C}_{\alpha}] contacts below a pre-specified distance cutoff (normally between 2 to 3 Å) in a protein crystal can be determined. Too many such contacts would indicate significant overlap of the molecules. For nucleic acid structures, a set of representative atoms (for example, P, [{\rm C}'_{1}]) can be selected from each nucleotide for this packing analysis.

13.3.6. The unique region of a translation function (the Cheshire group)

| top | pdf |

The region of the unit cell that should be covered during a translation search does not generally correspond to the asymmetric unit of the space group. Since the search model has a defined orientation, it can only reside in one of the asymmetric units in the unit cell. Lacking knowledge as to which asymmetric unit the model occupies, the entire unit cell would need to be searched. However, most space groups possess alternative origins, which means the position of a molecule in the unit cell can only be determined to within certain sets of translations. For example, in space group [P2_{1}2_{1}2_{1}], there are eight alternative origins at [(0, 0, 0)], [({1 \over 2}, 0, 0)], [(0, {1 \over 2}, 0)], [(0, 0, {1 \over 2})], [({1 \over 2}, {1 \over 2}, 0)], [({1 \over 2}, 0, {1 \over 2})], [(0, {1 \over 2}, {1 \over 2})] and [({1 \over 2}, {1 \over 2}, {1 \over 2})]. This implies that the region that should be searched to locate a molecule need only be [{1 \over 8}] of the volume of the unit cell [for example, [({1 \over 2}, {1 \over 2}, {1 \over 2})]]. In addition, for polar space groups, the position of the molecule along the polar axis is arbitrary. The symmetry, as defined by these unique regions, is also known as the Cheshire group (Hirshfeld, 1968[link]), and has been defined for all the 230 space groups.

Once the first molecule is positioned, the origin of the unit cell is fixed as well. The search for subsequent molecules will need to cover the entire unit cell, unless the crystal lattice is centred.

13.3.7. Combined molecular replacement

| top | pdf |

The traditional division of the molecular-replacement problem into two steps is partly due to limited computer power. Such a division has placed more pressure on the rotation function, since generally only a few rotation angles are examined by translation functions. The correct orientation, therefore, needs to be among the top few peaks either directly in the rotation functions or after Patterson-correlation refinement (Brünger, 1990[link]).

With modern computers, it is no longer necessary to maintain the strict division between the rotational and the translational components. Even though a full six-dimensional search is still generally impractical, a limited six-dimensional search can certainly be performed (Navaza, 1994[link]; Tong, 1996a[link]; Sheriff et al., 1999[link]; Jogl et al., 2001[link]). Target functions that can be evaluated quickly with the FFT technique are preferred for this purpose. Otherwise, the exploration of rotational degrees of freedom (Fujinaga & Read, 1987[link]) or the reflection data that are used in the calculation (Rabinovich & Shakked, 1984[link]) will be limited.

An automated molecular-replacement protocol has been implemented which automatically examines the top peaks in the rotation function by translation functions (Navaza, 1994[link]). This protocol has proven to be remarkably powerful. It assumes that the correct rotation solution is near the top peaks in the rotation function. A more general assumption is that the correct rotation for the search model should produce high values in the rotation function, even though they may not be near peaks in the rotation function. An error of 6° between the correct rotation angles and the peak in the rotation function, which often occurs, can make it impossible to obtain the correct translation-function solution (Fujinaga & Read, 1987[link]).

The combined-molecular-replacement protocol (Tong, 1996a[link]; Jogl et al., 2001[link]) therefore consists of examining all the grid points in the rotation function with heights greater than a defined cutoff value using the Patterson-correlation translation function. The top peaks (usually 20 to 50) in each translation function are all examined as possible solutions. The results from these translation functions are converted to the R factor or the correlation coefficient, enabling comparisons among the various orientations. This protocol allows the automatic examination of not only the top peaks in the rotation functions, but also those angles that produce high rotation-function values. In addition, packing of the model in the crystal is examined automatically to eliminate those solutions that have severe clashes among the molecules (see Section 13.3.5[link]). This generalized protocol has proven more powerful than conventional methods in structure determinations (Tong, 1996a[link]; Wu et al., 1997[link]). An alternative approach, examining the neighbourhood surrounding the rotation-function peaks, is also possible (Urzhumtsev & Podjarny, 1995[link]).

Stochastic or genetic algorithms have also been used successfully to simultaneously sample the six degrees of freedom (Chang & Lewis, 1997[link]; Kissinger et al., 1999[link]; Glykos & Kokkinidis, 2000[link]).

13.3.8. The locked translation function

| top | pdf |

In the presence of noncrystallographic symmetry (NCS), locked self-rotation functions can be used to determine the orientation of the NCS elements in the crystal unit cell (Tong & Rossmann, 1990[link]). Often, an atomic model for the monomer of the NCS assembly is available, but not a model of the entire assembly. This atomic model can be used in ordinary cross-rotation-function calculations. A more powerful technique is to use the locked cross-rotation function, which can define the orientations of all the molecules within the assembly at the same time (Tong & Rossmann, 1997[link]). With the knowledge of the orientations, several translation searches are needed to locate the individual monomers of the assembly. For cases where the assembly has high NCS, the translation searches to locate the first few molecules may not be very successful, since the search model only represents a small portion of the diffracting power of the crystal.

A locked translation function takes into account contributions from all the monomers of the assembly at the same time (Tong, 1996b[link]). It can determine the position of the monomer search model relative to the centre of the NCS assembly. With this knowledge, the entire assembly can be generated and can then be used in an ordinary translation search to locate the centre of this NCS assembly in the crystal unit cell.

Given the atomic model, [X_{j}^{0}] (in Cartesian coordinates), for the monomer at a starting position and the rotation, [F], that brings it into the same orientation as that of a monomer in the standard orientation, the model of the entire assembly in the standard orientation is given by [X_{j,  m} = [I_{m}] ([F]X_{j}^{0} + V_{0}), \eqno (]where [V_{0}] is a translation vector and the centre of the assembly is placed at (0, 0, 0). [[I_{m}]\ (m = 1, \ldots, M)] is the set of rotation matrices for the NCS point group in the standard orientation. The correct translation vector should give rise to the maximal overlap between the self vectors within the NCS assembly and the observed Patterson map. This overlap is given by the second term of equation ([link]). The locked translation function is therefore defined as [a comparison with equation ([link]) showing that the summation over the crystallographic symmetry operators can be omitted][\eqalignno{LTF(V_{0}) &= {\textstyle\sum\limits_{h}} {\textstyle\sum\limits_{m}} (F_{h}^{o})^{2}\ |\ \overline{f}_{h,  m}|^2\cr &\quad + {\textstyle\sum\limits_{h}} {\textstyle\sum\limits_{m}} {\textstyle\sum\limits_{n \neq m}} (F_{h}^{o})^{2} \overline{f}_{h,  m}\ \overline{f}_{h,  n}^{*} \exp \ \{- 2\pi ih([\theta_{n}] - [\theta_{m}])V_{0}\},\cr& &(}]where [\overline{f} _{h,  m} = \textstyle\sum\limits_{j} f_{j} \exp \ (2\pi ih[\theta_{m}][F]X_{j}^{0}) \eqno (]and [[\theta_{m}] = [\alpha][E][I_{m}]. \eqno (][E] is the rotation matrix that brings the standard orientation to that of the assembly in the crystal and [α] is the de-orthogonalization matrix (Rossmann & Blow, 1962[link]). Equation ([link]) can be evaluated indirectly by the FFT technique (Tong, 1996b[link]). As with the Patterson-correlation translation function, equation ([link]) can be converted to a correlation coefficient, although the evaluation will become more time-consuming. It should be noted that equation ([link]) bears much resemblance to equation ([link]), with the interchange of the crystallographic quantities [([T_{n}],\overline{f}_{h,  n})] and the noncrystallographic quantities [([\theta_{m}],\overline{f}_{h,  m})].

13.3.9. Miscellaneous translation functions

| top | pdf |

A variety of translation functions have been developed for special purposes. For example, a function has been proposed that can determine the translational component along an (NCS) twofold axis (Rossmann & Blow, 1962[link]).

If initial phase information is available for a crystal that possesses NCS, averaging among the NCS monomers is a powerful technique for improving the phase information (Rossmann, 1990[link]). Accurate parameters for the NCS elements (orientation and position) are essential for this averaging process. The orientations of the NCS elements can often be determined by self-rotation functions. The positions of the NCS elements can be determined by a special translation function based on the electron-density overlap (in a spherical volume) of NCS-related monomers (Tong, 1993[link]): [T([C],s_{1},s_{2}) = {\textstyle\sum\limits_{h}} {\textstyle\sum\limits_{p}} \overline{F}_{h} \overline{F}_{p} G_{hp} \exp (- 2\pi ihs_{1}) \exp (- 2\pi ips_{2}), \eqno (]where [s_{1}] and [s_{2}] are the centres of two molecules related by NCS, and [C] is the rotation matrix for the NCS. This equation can be used to obtain an initial estimate for the position of the NCS axis (given its orientation) (Tong, 1993[link]), to identify positions related by the NCS (Blow et al., 1964[link]) and to refine the NCS parameters iteratively (Tong et al., 1992[link], 1997[link]; Tong & Rossmann, 1997[link]).


Bentley, G. A. & Houdusse, A. (1992). Some applications of the phased translation function in macromolecular structure determination. Acta Cryst. A48, 312–322.
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.
Chang, G. & Lewis, M. (1997). Molecular replacement using genetic algorithms. Acta Cryst. D53, 279–289.
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.
Cowtan, K. (1998). Modified phased translation functions and their application to molecular-fragment location. Acta Cryst. D54, 750–756.
Crowther, R. A. & Blow, D. M. (1967). A method of positioning a known molecule in an unknown crystal structure. Acta Cryst. 23, 544–548.
Driessen, H. P. C., Bax, B., Slingsby, C., Lindley, P. F., Mahadevan, D., Moss, D. S. & Tickle, I. J. (1991). Structure of oligomeric βB2-crystallin: an application of the T2 translation function to an asymmetric unit containing two dimers. Acta Cryst. B47, 987–997.
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.
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.
Hirshfeld, F. L. (1968). Symmetry in the generation of trial structures. Acta Cryst. A24, 301–311.
Jogl, G., Tao, X., Xu, Y. & Tong, L. (2001). COMO: a program for combined molecular replacement. Acta Cryst. D57, 1127–1134.
Kissinger, C. R., Gehlhaar, D. K. & Fogel, D. B. (1999). Rapid automated molecular replacement by evolutionary search. Acta Cryst. D55, 484–491.
McCoy, A. J., Grosse-Kunstleve, R. W., Adams, P. D., Winn, M. D., Storoni, L. C. & Read, R. J. (2007). Phaser crystallographic software. J. Appl. Cryst. 40, 658–674.
McCoy, A. J., Grosse-Kunstleve, R. W., Storoni, L. C. & Read, R. J. (2005). Likelihood-enhanced fast translation functions. Acta Cryst. D61, 458–464.
Navaza, J. (1994). AMoRe: an automated package for molecular replacement. Acta Cryst. A50, 157–163.
Navaza, J. & Vernoslova, E. (1995). On the fast translation functions for molecular replacement. Acta Cryst. A51, 445–449.
Nixon, P. E. & North, A. C. T. (1976). Crystallographic relationship between human and hen-egg lysozymes. I. Methods for the establishment of molecular orientational and positional parameters. Acta Cryst. A32, 320–325.
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.
Read, R. J. (2001). Pushing the boundaries of molecular replacement with maximum likelihood. Acta Cryst. D57, 1373–1382.
Read, R. J. & Schierbeek, A. J. (1988). A phased translation function. J. Appl. Cryst. 21, 490–495.
Rossmann, M. G. (1972). Editor. The Molecular Replacement Method. New York: Gordon & Breach.
Rossmann, M. G. (1990). The molecular replacement method. Acta Cryst. A46, 73–82.
Rossmann, M. G. (2001). Molecular replacement – historical background. Acta Cryst. D57, 1360–1366.
Rossmann, M. G. & Blow, D. M. (1962). The detection of sub-units within the crystallographic asymmetric unit. Acta Cryst. 15, 24–31.
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.
Tickle, I. J. (1985). Review of space group general translation functions that make use of known structure information and can be expanded as Fourier series. In Proceedings of the Daresbury Study Weekend. Molecular Replacement, edited by P. A. Machin, pp. 22–26. Warrington: Daresbury Laboratory.
Tollin, P. (1966). On the determination of molecular location. Acta Cryst. 21, 613–614.
Tollin, P. & Rossmann, M. G. (1966). A description of various rotation function programs. Acta Cryst. 21, 872–876.
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., Choi, H.-K., Minor, W. & Rossmann, M. G. (1992). The structure determination of Sindbis virus core protein using isomorphous replacement and molecular replacement averaging between two crystal forms. Acta Cryst. A48, 430–442.
Tong, L., Qian, C., Davidson, W., Massariol, M.-J., Bonneau, P. R., Cordingley, M. G. & Lagacé, L. (1997). Experiences from the structure determination of human cytomegalovirus protease. Acta Cryst. D53, 682–690.
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.
Wu, H., Kwong, P. D. & Hendrickson, W. A. (1997). Dimeric association and segmental variability in the structure of human CD4. Nature (London), 387, 527–530.
Zhang, X.-J. & Matthews, B. W. (1994). Enhancement of the method of molecular replacement by incorporation of known structural information. Acta Cryst. D50, 675–686.

to end of page
to top of page