Single molecule pulling with large time steps
- Author(s)
- Harald Oberhofer, Christoph Dellago, Stefan Boresch
- Abstract
Recently, we presented a generalization of the Jarzynski nonequilibrium work theorem for phase space mappings. The formalism shows that one can determine free energy differences from approximate trajectories obtained from molecular dynamics simulations in which very large time steps are used. In this work we test the method by simulating the force-induced unfolding of a deca-alanine helix in vacuum. The excellent agreement between results obtained with a small, conservative time step of 0.5 fs and results obtained with a time step of 3.2 fs (i.e., close to the stability limit) indicates that the large-time-step approach is practical for such complex biomolecules. We further adapt the method of Hummer and Szabo for the simulation of single molecule force spectroscopy experiments to the large-time-step method. While trajectories generated with large steps are approximate and may be unphysical-in the simulations presented here we observe a violation of the equipartition theorem-the computed free energies are exact in principle. In terms of efficiency, the optimum time step for the unfolding simulations lies in the range 1-3 fs. © 2007 The American Physical Society.
- Organisation(s)
- Computational Materials Physics, Computational and Soft Matter Physics, Department of Computational Biological Chemistry
- Journal
- Physical Review E
- Volume
- 75
- No. of pages
- 7
- ISSN
- 1539-3755
- DOI
- https://doi.org/10.1103/PhysRevE.75.061106
- Publication date
- 2007
- Peer reviewed
- Yes
- Austrian Fields of Science 2012
- 103036 Theoretical physics, 104022 Theoretical chemistry, 103029 Statistical physics
- Portal url
- https://ucrisportal.univie.ac.at/en/publications/11a8aaec-b319-46b5-aff4-22febee7bfaf