On the multiple simultaneous superposition of molecular structures by rigid body transformations

A method of optimally superimposing n coordinate sets on each other by rigid body transformations, which minimizes the sum of all n(n ‐ 1)/2 pairwise residuals, is presented. In the solution phase the work load is approximately linear on n, is independent of the size of the structures, is independent of their initial orientations, and terminates in one cycle if n = 2 or if the coordinate sets are exactly superposable, and otherwise takes a number of cycles dependent only on genuine shape differences. Enantiomorphism, if present, is detected, in which case the option exists to reverse or not to reverse the chirality of relevant coordinate sets. The method also offers a rational approach to the problem of multiple minima and has successfully identified four distinct minima in such a case. Source code, which is arranged to enable the study of the disposition of domains in multidomain structures, is available from the author.

[1]  James Hardy Wilkinson,et al.  The Calculation of the Eigenvectors of Codiagonal Matrices , 1958, Comput. J..

[2]  James Hardy Wilkinson,et al.  Householder's Method for the Solution of the Algebraic Eigenproblem , 1960, Comput. J..

[3]  James M. Ortega,et al.  On Sturm Sequences for Tridiagonal Matrices , 1960, JACM.

[4]  A. D. McLachlan,et al.  A mathematical procedure for superimposing atomic coordinates of proteins , 1972 .

[5]  R. Diamond On the comparison of conformations using linear and quadratic transformations , 1976 .

[6]  W. Kabsch A solution for the best rotation to relate two sets of vectors , 1976 .

[7]  W. Kabsch A discussion of the solution for the best rotation to relate two sets of vectors , 1978 .

[8]  A. Mclachlan Gene duplications in the structural evolution of chymotrypsin. , 1979, Journal of molecular biology.

[9]  A. D. McLachlan,et al.  Rapid comparison of protein structures , 1982 .

[10]  A. M. Lesk,et al.  A toolkit for computational molecular biology. II. On the optimal superposition of two sets of coordinates , 1986 .

[11]  R. Diamond A note on the rotational superposition problem , 1988 .

[12]  S. Kearsley On the orthogonal transformation used for structural comparisons , 1989 .

[13]  R. Diamond Chirality in rotational superposition , 1990 .

[14]  S. Kearsley An algorithm for the simultaneous superposition of a structural series , 1990 .

[15]  R J Williams,et al.  Topological mirror images in protein structure computation: An underestimated problem , 1991, Proteins.

[16]  A. Lesk,et al.  A method for multiple superposition of structures. , 1992, Acta crystallographica. Section A, Foundations of crystallography.

[17]  A Klug,et al.  Solution structures of two zinc-finger domains from SWI5 obtained using two-dimensional 1H nuclear magnetic resonance spectroscopy. A zinc-finger structure with a third strand of beta-sheet. , 1992, Journal of molecular biology.