Springer Verlag, Lecture Notes in Computer Science, p. 354-373
DOI: 10.1007/978-3-540-39763-2_26
Full text: Unavailable
Monte Carlo simulation (MCS) is a common methodology to compute pathways and thermodynamic properties of proteins. A sim- ulation run is a series of random steps in conformation space, each per- turbing some degrees of freedom of the molecule. A step is accepted with a probability that depends on the change in value of an energy func- tion. Typical energy functions sum many terms. The most costly ones to compute are contributed by atom pairs closer than some cuto distance. This paper introduces a new method that speeds up MCS by eciently computing the energy at each step. The method exploits the facts that proteins are long kinematic chains and that few degrees of freedom are changed at each step. A novel data structure, called the ChainTree, cap- tures both the kinematics and the shape of a protein at successive levels of detail. It is used to find all atom pairs contributing to the energy. It also makes it possible to identify partial energy sums left unchanged by a perturbation, thus allowing the energy value to be incrementally up- dated. Computational tests on four proteins of sizes ranging from 68 to 755 amino acids show that MCS with the ChainTree method is signifi- cantly faster (as much as 12 times faster for the largest protein) than with the widely used grid method. They also indicate that speed-up increases with larger proteins.