PROFASI: A Monte Carlo simulation package for protein folding and aggregation

We present a flexible and efficient program package written in C++, PROFASI, for simulating protein folding and aggregation. The systems are modeled using an all‐atom description of the protein chains with only torsional degrees of freedom, and implicit water. The program package has a modular structure that makes the interaction potential easy to modify. The currently implemented potential is able to fold several peptides with about 20 residues, and has also been used to study aggregation and force‐induced unfolding. The simulation methods implemented in PROFASI are Monte Carlo‐based and include a semilocal move and simulated tempering. Adding new updates is easy. The code runs fast in both single‐ and multi‐chain applications, as is illustrated by several examples. © 2006 Wiley Periodicals, Inc. J Comput Chem 27: 1548–1555, 2006

[1]  R. Leapman,et al.  Amyloid Fibril Formation by Aβ16-22, a Seven-Residue Fragment of the Alzheimer's β-Amyloid Peptide, and Structural Characterization by Solid State NMR† , 2000 .

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

[3]  A. Irbäck,et al.  Folding thermodynamics of peptides. , 2004, Biophysical journal.

[4]  U. Hansmann.  Parallel tempering algorithm for conformational studies of biological molecules , 1997, physics/9710041.

[5]  A. Sokal,et al.  The pivot algorithm: A highly efficient Monte Carlo method for the self-avoiding walk , 1988 .

[6]  P S Kim,et al.  Electrostatic screening of charge and dipole interactions with the helix backbone. , 1993, Science.

[7]  Wang,et al.  Replica Monte Carlo simulation of spin glasses. , 1986, Physical review letters.

[8]  Robert A. Grothe,et al.  Structure of the cross-β spine of amyloid-like fibrils , 2005, Nature.

[9]  Anthony Skjellum,et al.  A High-Performance, Portable Implementation of the MPI Message Passing Interface Standard , 1996, Parallel Comput..

[10]  Berk Hess,et al.  GROMACS 3.0: a package for molecular simulation and trajectory analysis , 2001 .

[11]  S Gnanakaran,et al.  Peptide folding simulations. , 2003, Current opinion in structural biology.

[12]  G. Favrin,et al.  Oligomerization of amyloid Abeta16-22 peptides using hydrogen bonds and hydrophobicity forces. , 2004, Biophysical journal.

[13]  Simon Mitternacht,et al.  Dissecting the mechanical unfolding of ubiquitin. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[14]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[15]  Stefan Wallin,et al.  Thermodynamics of alpha- and beta-structure formation in proteins. , 2003, Biophysical journal.

[16]  M. Lal,et al.  ‘Monte Carlo’ computer simulation of chain molecules , 1969 .

[17]  Yuko Okamoto,et al.  Secondary-structure preferences of force fields for proteins evaluated by generalized-ensemble simulations , 2004 .

[18]  Peter A. Kollman,et al.  AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of molecules , 1995 .

[19]  J. W. Neidigh,et al.  Designing a 20-residue protein , 2002, Nature Structural Biology.

[20]  D. Theodorou,et al.  A concerted rotation algorithm for atomistic Monte Carlo simulation of polymer melts and glasses , 1993 .

[21]  Erik Sandelin,et al.  Monte Carlo study of the phase structure of compact polymer chains , 1998, cond-mat/9812017.

[22]  G. Favrin,et al.  Monte Carlo update for chain molecules: Biased Gaussian steps in torsional space , 2001, cond-mat/0103580.

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

[24]  Ulrich H. E. Hansmann,et al.  SMMP) A modern package for simulation of proteins , 2001 .

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

[26]  Anders Irbäck,et al.  Folding thermodynamics of three β‐sheet peptides: A model study , 2003, Proteins.

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

[28]  A. Lyubartsev,et al.  New approach to Monte Carlo calculation of the free energy: Method of expanded ensembles , 1992 .

[29]  K. Hukushima,et al.  Exchange Monte Carlo Method and Application to Spin Glass Simulations , 1995, cond-mat/9512035.

[30]  L. Serrano,et al.  A short linear peptide that folds into a native stable β-hairpin in aqueous solution , 1994, Nature Structural Biology.

[31]  N. Skelton,et al.  Tryptophan zippers: Stable, monomeric β-hairpins , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[32]  STUDIES OF AN OFF-LATTICE MODEL FOR PROTEIN FOLDING: SEQUENCE DEPENDENCE AND IMPROVED SAMPLING AT FINITE TEMPERATURE , 1995, chem-ph/9505003.