Parallel numerical solution to large-scale eigenvalue problem in master equation of protein folding kinetics

A master equation characterizes the time-evolution of trajectories, the transition of states in protein folding kinetics. Numerical solution of the master equation requires calculating eigenvalues for the corresponding large scale eigenvalue problem. In this paper, we present a parallel computing technique to compute the eigenvalues of the matrix with an N-dimensional vector of the instantaneous probability of the N conformations. Parallelization of the implicitly restarted Arnoldi method is successfully implemented on a PC-based Linux cluster. The parallelization scheme used in this work mainly partitions the operations of the matrix. For the Arnoldi factorization, we replicate the upper Hessenberg matrix H/sub m/ for each processor, and distribute the set of Arnoldi vectors V/sub m/ among processors. Each processor performs its own operations. This algorithm is implemented on a PC-based Linux cluster with message passing interface (MPI) libraries. Our preliminary numerical experiment performing on the 32-nodes PC-based Linux cluster has shown that the maximum difference among CPUs is within 10%. A 23 times speedup and 72% parallel efficiency are also attained for the tested cases. This approach enables us to explore large scale dynamics of protein folding.

[1]  S W Englander,et al.  Protein folding intermediates and pathways studied by hydrogen exchange. , 2000, Annual review of biophysics and biomolecular structure.

[2]  Guido Bologna,et al.  A comparison study on protein fold recognition , 2002, Proceedings of the 9th International Conference on Neural Information Processing, 2002. ICONIP '02..

[3]  Gerard L. G. Sleijpen,et al.  A Jacobi-Davidson Iteration Method for Linear Eigenvalue Problems , 1996, SIAM Rev..

[4]  Yiming Li.  Numerical Solutions of Master Equation for Protein Folding Kinetics , 2005, 11th International Conference on Parallel and Distributed Systems (ICPADS'05).

[5]  Yiming Li,et al.  A Practical Implementation of Parallel Dynamic Load Balancing for Adaptive Computing in VLSI Device Simulation , 2002, Engineering with Computers.

[6]  Gerard L. G. Sleijpen,et al.  A Jacobi-Davidson Iteration Method for Linear Eigenvalue Problems , 1996, SIAM J. Matrix Anal. Appl..

[7]  Kesheng Wu,et al.  A Block Orthogonalization Procedure with Constant Synchronization Requirements , 2000, SIAM J. Sci. Comput..

[8]  Eugene Santos,et al.  Local minima-based exploration for off-lattice protein folding , 2003, Computational Systems Bioinformatics. CSB2003. Proceedings of the 2003 IEEE Bioinformatics Conference. CSB2003.

[9]  John S. McCaskill,et al.  Folding Stabilizes the Evolution of Catalysts , 2004, Artificial Life.

[10]  K. Dill,et al.  From Levinthal to pathways to funnels , 1997, Nature Structural Biology.

[11]  David Baker,et al.  Simple physical models connect theory and experiment in protein folding kinetics. , 2002, Journal of molecular biology.

[12]  G. Stewart,et al.  Reorthogonalization and stable algorithms for updating the Gram-Schmidt QR factorization , 1976 .

[13]  Valerie Daggett,et al.  Molecular dynamics simulations of the protein unfolding/folding reaction. , 2002, Accounts of chemical research.

[14]  Masato Ishikawa,et al.  A three-dimensional animation system for protein folding simulation , 1994, 1994 Proceedings of the Twenty-Seventh Hawaii International Conference on System Sciences.

[15]  Yiming Li.  Numerical solutions of a master equation for protein folding kinetics , 2006 .

[16]  Danny C. Sorensen,et al.  Deflation Techniques for an Implicitly Restarted Arnoldi Iteration , 1996, SIAM J. Matrix Anal. Appl..

[17]  Amjad Umar.  Distributed And Parallel Computing , 1998, IEEE Concurrency.

[18]  D. Yee,et al.  Principles of protein folding — A perspective from simple exact models , 1995, Protein science : a publication of the Protein Society.

[19]  J. Ferkinghoff-Borg,et al.  Mean first-passage time analysis reveals rate-limiting steps, parallel pathways and dead ends in a simple model of protein folding , 2003 .

[20]  D. Sorensen.  IMPLICITLY RESTARTED ARNOLDI/LANCZOS METHODS FOR LARGE SCALE EIGENVALUE CALCULATIONS , 1996 .

[21]  Yiming Li,et al.  A Two-Dimensional Quantum Transport Simulation of Nanoscale Double-Gate MOSFETs Using Parallel Adaptive Technique , 2004, IEICE Trans. Inf. Syst..