Theory and Practice in Replica-Exchange Molecular Dynamics Simulation
Date
2008-11-26
Authors
Advisors
Journal Title
Journal ISSN
Volume Title
Repository Usage Stats
views
downloads
Abstract
We study the comparison of computational simulations of biomolecules to experimental data. We study the convergence of these simulations to equilibrium and determine measures of variance of the data using statistical methods. We run replica-exchange molecular dynamics (REMD) simulations of eight helical peptides and compare the simulation helicity to the experimentally measured helicity of the peptides. We use one-way sensitivity analysis to determine which parameter changes have a large effect on helicity measurements and use Bayesian updating for a parameter of the AMBER potential. We then consider the theoretical convergence behavior of the REMD algorithm itself by evaluating the properties of the isothermal numerical integrators used in the underlying MD. The underlying constant-temperature integrators explored in this thesis represent a majority of the deterministic isothermal methods used with REMD simulations and we show that these methods either fail to be measure-invariant or are not ergodic. For each of the non-ergodic integrators we show that REMD fails to be ergodic when run with the integrator. We give computational results from examples to demonstrate the practical implications of non-ergodicity and describe hybrid Monte Carlo, a method that leads to ergodicity. Finally, we consider the use of stochastic Langevin dynamics to simulate isothermal MD. We show geometric ergodicity of the Langevin diffusion over a simplified system with the eventual goal of determining geometric ergodicity for Langevin dynamics over the full AMBER potential.
Type
Department
Description
Provenance
Subjects
Citation
Permalink
Citation
Cooke, Benjamin (2008). Theory and Practice in Replica-Exchange Molecular Dynamics Simulation. Dissertation, Duke University. Retrieved from https://hdl.handle.net/10161/923.
Collections
Except where otherwise noted, student scholarship that was shared on DukeSpace after 2009 is made available to the public under a Creative Commons Attribution / Non-commercial / No derivatives (CC-BY-NC-ND) license. All rights in student work shared on DukeSpace before 2009 remain with the author and/or their designee, whose permission may be required for reuse.