Random walk in orthogonal space to achieve efficient free-energy simulation of complex systems
- 23 December 2008
- journal article
- research article
- Published by Proceedings of the National Academy of Sciences in Proceedings of the National Academy of Sciences of the United States of America
- Vol. 105 (51), 20227-20232
- https://doi.org/10.1073/pnas.0810631106
Abstract
In the past few decades, many ingenious efforts have been made in the development of free-energy simulation methods. Because complex systems often undergo nontrivial structural transition during state switching, achieving efficient free-energy calculation can be challenging. As identified earlier, the “Hamiltonian” lagging, which reveals the fact that necessary structural relaxation falls behind the order parameter move, has been a primary problem for generally low free-energy simulation efficiency. Here, we propose an algorithm by achieving a random walk in both the order parameter space and its generalized force space; thereby, the order parameter move and the required conformational relaxation can be efficiently synchronized. As demonstrated in both the alchemical transition and the conformational transition, a leapfrog improvement in free-energy simulation efficiency can be obtained; for instance, (i) it allows us to solve a notoriously challenging problem: accurately predicting the pKa value of a buried titratable residue, Asp-66, in the interior of the V66E staphylococcal nuclease mutant, and (ii) it allows us to gain superior efficiency over the metadynamics algorithm.Keywords
This publication has 23 references indexed in Scilit:
- pKa of Residue 66 in Staphylococal nuclease. I. Insights from QM/MM Simulations with Conventional SamplingThe Journal of Physical Chemistry B, 2008
- Quantitative computer simulations of biomolecules: A snapshotJournal of Computational Chemistry, 2007
- Synergistic approach to improve “alchemical” free energy calculation in rugged energy surfaceThe Journal of Chemical Physics, 2007
- High Apparent Dielectric Constant Inside a Protein Reflects Structural Reorganization Coupled to the Ionization of an Internal AspBiophysical Journal, 2007
- Optimal Estimates of Free Energies from Multistate Nonequilibrium Work DataPhysical Review Letters, 2006
- Efficient Sampling of Protein Structures by Model HoppingPhysical Review Letters, 2005
- Protein folding simulations in a deformed energy landscapeZeitschrift für Physik B Condensed Matter, 1999
- λ-dynamics: A new approach to free energy calculationsThe Journal of Chemical Physics, 1996
- Computational AlchemyAnnual Review of Physical Chemistry, 1992
- CHARMM: A program for macromolecular energy, minimization, and dynamics calculationsJournal of Computational Chemistry, 1983