Isomorphism extension theorem

The Kabsch algorithm, named after Wolfgang Kabsch, is a method for calculating the optimal rotation matrix that minimizes the RMSD (root mean squared deviation) between two paired sets of points. It is useful in graphics, cheminformatics to compare molecular structures, and also bioinformatics for comparing protein structures (in particular, see root-mean-square deviation (bioinformatics)).

The algorithm only computes the rotation matrix, but it also requires the computation of a translation vector. When both the translation and rotation are actually performed, the algorithm is sometimes called partial Procrustes superimposition (see also orthogonal Procrustes problem).

Description

The algorithm starts with two sets of paired points, P and Q. Each set of points can be represented as an N×3 matrix. The first row is the coordinates of the first point, the second row is the coordinates of the second point, the Nth row is the coordinates of the Nth point.

$\left(\begin{array}{ccc}{x}_1& y_1& z_1\\ x_2& y_2& z_2\\ ⋮& ⋮& ⋮\\ x_N& y_N& z_N\end{array}\right)$

The algorithm works in three steps: a translation, the computation of a covariance matrix, and the computation of the optimal rotation matrix.

Translation

Both sets of coordinates must be translated first, so that their centroid coincides with the origin of the coordinate system. This is done by subtracting from the point coordinates the coordinates of the respective centroid.

Computation of the covariance matrix

The second step consist of calculating a covariance matrix A. In matrix notation,

$A=P^{T}Q$

or, using summation notation,

$A_{ij}=\sum _{k=1}^N{P}_{ki}{Q}_{kj},$

Computation of the optimal rotation matrix

It is possible to calculate the optimal rotation U based on the matrix formula $U=(A^{T}A{)}^{1/2}{A}^{-1}$ but implementing a numerical solution to this formula becomes complicated when all special cases are accounted for (for example, the case of A not having an inverse).

If singular value decomposition (SVD) routines are available, the optimal rotation, U, can be calculated using the following simple algorithm.

First, calculate the SVD of the covariance matrix A.

$A=VS{W}^T$

Next, decide whether we need to correct our rotation matrix to ensure a right-handed coordinate system

$d=\mathrm{sign}(\det (W{V}^T))$

Finally, calculate our optimal rotation matrix, U, as

$U=W\left(\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 0& 0& d\end{array}\right)V^T$

Coutsias, Seok, and Dill have found an equivalent method that uses quaternions.

Generalizations

The algorithm was described for points in a three-dimensional space. The generalization to D dimensions is immediate.

External links

This SVD algorithm is described in more detail at http://cnx.org/content/m11608/latest/

A Matlab function is available at http://www.mathworks.com/matlabcentral/fileexchange/25746-kabsch-algorithm

A free PyMol plugin easily implementing Kabsch is Cealign. VMD uses the Kabsch algorithm for its alignment.

A Python script is available at http://github.com/charnley/rmsd

See also

Wahba's Problem

References

43 year old Petroleum Engineer Harry from Deep River, usually spends time with hobbies and interests like renting movies, property developers in singapore new condominium and vehicle racing. Constantly enjoys going to destinations like Camino Real de Tierra Adentro.

• Kabsch, Wolfgang, (1976) "A solution for the best rotation to relate two sets of vectors", Acta Crystallographica 32:922. 21 year-old Glazier James Grippo from Edam, enjoys hang gliding, industrial property developers in singapore developers in singapore and camping. Finds the entire world an motivating place we have spent 4 months at Alejandro de Humboldt National Park. with a correction in Kabsch, Wolfgang, (1978) "A discussion of the solution for the best rotation to relate two sets of vectors", "Acta Crystallographica", "A34", 827–828 21 year-old Glazier James Grippo from Edam, enjoys hang gliding, industrial property developers in singapore developers in singapore and camping. Finds the entire world an motivating place we have spent 4 months at Alejandro de Humboldt National Park.

• Lin Ying-Hung, Chang Hsun-Chang, Lin Yaw-Ling (2004) "A Study on Tools and Algorithms for 3-D Protein Structures Alignment and Comparison", International Computer Symposium, December 15–17, Taipei, Taiwan.