Diagonalization in a mixed basis: A method to compute low-frequency normal modes for large macromolecules

Low-frequency collective motions in proteins are generally very important for their biological functions. To study such motions, harmonic dynamics proved most useful since it is a straightforward method; it consists of the diagonalization of the Hessian matrix of the potential energy, yielding the vibrational spectrum and the directions of internal motions. Unfortunately, the diagonalization of this matrix requires a large computer memory, which is a limiting factor when the protein contains several thousand atoms. To circumvent this limitation we have developed three methods that enable us to diagonalize large matrices using much less computer memory than the usual harmonic dynamics. The first method is approximate; it consists of diagonalizing small blocks of the Hessian matrix, followed by the coupling of the low-frequency modes obtained for each block. It yields the low-frequency vibrational spectrum with a maximum error of 20%. The second method consists, after diagonalizing small blocks, of coupling the high- and low-frequency modes using an iterative procedure. It yields the exact low-frequency normal modes, but requires a long computational time with convergence problems. The third method, DIMB (Diagonalization in a Mixed Basis), which has the best performance, consists of coupling the approximate low-frequency modes with the mass-weighted cartesian coordinates, also using an iterative procedure. It reduces significantly the required computer memory and converges rapidly. The eigenvalues and eigenvectors obtained by this method are without significant error in the chosen frequency range. Moreover, it is a general method applicable to any problem of diagonalization of a large matrix. We report the application of these methods to a deca-alanine helix, trypsin inhibitor, a neurotoxin, and lysozyme. © 1993 John Wiley & Sons, Inc.