Multicanonical parallel simulations of proteins with continuous potentials

The determination of the three‐dimensional (3D) structure of a protein or peptide is a very important research problem in biological and medical sciences. Anfinsen's experiments (Science 1973, 181, 223) on renaturation of denatured proteins have shown that the native 3D structure of a (small) protein at low (room) temperatures is uniquely determined by its amino acid sequence, which suggests that it might be possible to determine the 3D structure of a protein from its amino acid sequence by pure computations. As a step toward that goal, in this article we present a simple approach for parallelization of multicanonical Monte Carlo simulations of proteins with continuous potentials. Our method is based on the parallel calculation of the protein energy function. The algorithm is tested by simulated annealing and multicanonical simulations of two small peptides, and known results are reproduced accurately. An acceptable degree of parallelization can be achieved in the simulation of Protein L using up to 30 PCs. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1287–1296, 2001

[1]  B. Berg.  Multicanonical recursions , 1995, hep-lat/9503019.

[2]  H. Scheraga,et al.  Energy parameters in polypeptides. 10. Improved geometrical parameters and nonbonded interactions for use in the ECEPP/3 algorithm, with application to proline-containing peptides , 1994 .

[3]  Yuko Okamoto,et al.  Comparative Study of Multicanonical and Simulated Annealing Algorithms in the Protein Folding Problem , 1994 .

[4]  Berg,et al.  Multicanonical ensemble: A new approach to simulate first-order phase transitions. , 1992, Physical review letters.

[5]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[6]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[7]  Y. Okamoto,et al.  GENERALIZED-ENSEMBLE MONTE CARLO METHOD FOR SYSTEMS WITH ROUGH ENERGY LANDSCAPE , 1997, cond-mat/9710306.

[8]  H. Scheraga,et al.  Intermolecular potentials from crystal data. 6. Determination of empirical potentials for O-H...O = C hydrogen bonds from packing configurations , 1984 .

[9]  Alan M. Ferrenberg,et al.  New Monte Carlo technique for studying phase transitions. , 1988, Physical review letters.

[10]  N. Wingreen,et al.  Emergence of Preferred Structures in a Simple Model of Protein Folding , 1996, Science.

[11]  Michael Griebel,et al.  Numerical Simulation in Fluid Dynamics: A Practical Introduction , 1997 .

[12]  W. Braun,et al.  Rapid calculation of first and second derivatives of conformational energy with respect to dihedral angles for proteins general recurrent equations , 1984, Comput. Chem..

[13]  B. Smit,et al.  Parallel Gibbs-ensemble simulations , 1995 .

[14]  L. Björck,et al.  Three-dimensional solution structure of an immunoglobulin light chain-binding domain of protein L. Comparison with the IgG-binding domains of protein G. , 1994, Biochemistry.

[15]  G. Parisi,et al.  Simulated tempering: a new Monte Carlo scheme , 1992, hep-lat/9205018.

[16]  Professor Dr. Dietrich Stauffer,et al.  Computer Simulation and Computer Algebra , 1993, Springer Berlin Heidelberg.

[17]  Wolfhard Janke,et al.  Multicanonical Monte Carlo simulations , 1998 .

[18]  Eli Gafni,et al.  A Distributed Implementation of Simulated Annealing , 1989, J. Parallel Distributed Comput..

[19]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[20]  Esselink,et al.  Parallel Monte Carlo simulations. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[21]  Bernd A. Berg.  Algorithmic aspects of multicanonical simulations , 1997 .

[22]  Hesselbo,et al.  Monte Carlo simulation and global optimization without parameters. , 1995, Physical review letters.

[23]  Y. Okamoto,et al.  Molecular dynamics, Langevin, and hybrid Monte Carlo simulations in multicanonical ensemble , 1996, physics/9710018.

[24]  B. Berg,et al.  Multicanonical algorithms for first order phase transitions , 1991 .

[25]  Yuko Okamoto,et al.  Prediction of peptide conformation by multicanonical algorithm: New approach to the multiple‐minima problem , 1993, J. Comput. Chem..

[26]  C. Tsallis.  Possible generalization of Boltzmann-Gibbs statistics , 1988 .

[27]  Yuko Okamoto,et al.  Tertiary Structure Prediction of C-Peptide of Ribonuclease A by Multicanonical Algorithm , 1998, physics/9806017.

[28]  F. Young.  Biochemistry , 1955, The Indian Medical Gazette.

[29]  Berend Smit,et al.  Understanding Molecular Simulation , 2001 .

[30]  C. Anfinsen.  Principles that govern the folding of protein chains. , 1973, Science.

[31]  R. Abagyan,et al.  New methodology for computer-aided modelling of biomolecular structure and dynamics. 1. Non-cyclic structures. , 1989, Journal of biomolecular structure & dynamics.

[32]  C. Epstein,et al.  The Genetic Control of Tertiary Protein Structure: Studies With Model Systems , 1963 .

[33]  H. Scheraga,et al.  Energy parameters in polypeptides. VII. Geometric parameters, partial atomic charges, nonbonded interactions, hydrogen bond interactions, and intrinsic torsional potentials for the naturally occurring amino acids , 1975 .

[34]  Chwan-Hwa Wu,et al.  A parallel Monte Carlo simulation for gaseous electronics on a dynamically reconfigurable multiprocessor system , 1991 .

[35]  N. Go,et al.  A Method of Rapid Calculation of a Second Derivative Matrix of Conformational Energy for Large Molecules , 1983 .

[36]  Y Okamoto,et al.  Dependence on the dielectric model and pH in a synthetic helical peptide studied by Monte Carlo simulated annealing , 1994, Biopolymers.

[37]  D. J. Wallace,et al.  Concurrency and parallelism in MC and MD simulations in physics , 1985 .

[38]  Julia M. Goodfellow,et al.  Parallelization strategies for molecular simulation using the Monte Carlo algorithm , 1993, J. Comput. Chem..

[39]  H. Scheraga,et al.  Energy parameters in polypeptides. 9. Updating of geometrical parameters, nonbonded interactions, and hydrogen bond interactions for the naturally occurring amino acids , 1983 .

[40]  A. Moatti,et al.  Parallel Monte Carlo calculations with many microcomputers , 1987 .

[41]  Y. Okamoto,et al.  Numerical comparisons of three recently proposed algorithms in the protein folding problem , 1997 .

[42]  Ian K. Snook,et al.  A massively parallel molecular dynamics algorithm for the MasPar supercomputer , 1998 .