A combinatorial algorithm for effective generation of long maximally compact lattice chains
- 1 November 1995
- journal article
- Published by AIP Publishing in The Journal of Chemical Physics
- Vol. 103 (17) , 7592-7604
- https://doi.org/10.1063/1.470277
Abstract
We investigate the problem of generation of maximally compact lattice chains which are useful in understanding folding of model proteins. The term, maximally compact chain, refers to a lattice self-avoiding walk that visits every lattice site. Generation of a representative sample of compact conformations is extremely difficult by conventional simulation methods such as static growth methods or dynamic Monte Carlo techniques. Growing a random walk is ineffective for generating long walks in a compact shape because a large number of walks are rejected due to overlap (attrition). In the interest of an unbiased sample, one needs to enumerate all possible compact conformations that are realizable or produce a representative sample, the former of which is intractable for long chains. In this paper a method is proposed for generation of compact chains on a lattice based on a mathematical programming approach. The method, which we refer to as the Hamiltonian path generation method, generates a random sample of lattice filling self-avoiding walks. A detailed description of a randomized generation algorithm is presented, which is effective for producing a static sample of compact lattice chains. There is a statistical evidence of fair generation of conformations from the conformational space using this scheme. This method generates a compact conformation on a 60×60×60 cubic lattice in forty minutes on a Sparc-2 workstation.Keywords
This publication has 16 references indexed in Scilit:
- Monte Carlo simulation of folding transitions of simple model proteins using a chain growth algorithmThe Journal of Chemical Physics, 1992
- Monte Carlo simulation of the collapse-coil transition in homopolymersThe Journal of Chemical Physics, 1992
- Enumeration of all compact conformations of copolymers with random sequence of linksThe Journal of Chemical Physics, 1990
- Origins of structure in globular proteins.Proceedings of the National Academy of Sciences, 1990
- The collapse transition of self-avoiding walks on a square lattice: A computer simulation studyThe Journal of Chemical Physics, 1989
- Monte Carlo simulation of lattice models for macromoleculesComputer Physics Reports, 1988
- The pivot algorithm: A highly efficient Monte Carlo method for the self-avoiding walkJournal of Statistical Physics, 1988
- Monte Carlo studies of polymer chain dimensions in the meltThe Journal of Chemical Physics, 1982
- Macromolecular dimensions obtained by an efficient Monte Carlo method without sample attritionThe Journal of Chemical Physics, 1975
- New Method for the Statistical Computation of Polymer DimensionsThe Journal of Chemical Physics, 1959