Recovering position-dependent diffusion from biased molecular dynamics simulations
- Laboratory of Biophysics, Condensed Matter Physics Department, “Jožef Stefan” Institute, 1000 Ljubljana (Slovenia)
All atom molecular dynamics (MD) models provide valuable insight into the dynamics of biophysical systems, but are limited in size or length by the high computational demands. The latter can be reduced by simulating long term diffusive dynamics (also known as Langevin dynamics or Brownian motion) of the most interesting and important user-defined parts of the studied system, termed collective variables (colvars). A few hundred nanosecond-long biased MD trajectory can therefore be extended to millisecond lengths in the colvars subspace at a very small additional computational cost. In this work, we develop a method for determining multidimensional anisotropic position- and timescale-dependent diffusion coefficients (D) by analysing the changes of colvars in an existing MD trajectory. As a test case, we obtained D for dihedral angles of the alanine dipeptide. An open source Mathematica{sup ®} package, capable of determining and visualizing D in one or two dimensions, is available at https://github.com/lbf-ijs/DiffusiveDynamics . Given known free energy and D, the package can also generate diffusive trajectories.
- OSTI ID:
- 22255050
- Journal Information:
- Journal of Chemical Physics, Vol. 140, Issue 8; Other Information: (c) 2014 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9606
- Country of Publication:
- United States
- Language:
- English
Similar Records
Modeling solvation contributions to conformational free energy changes of biomolecules using a potential of mean force expansion
Tutorial: Langevin Dynamics methods for aerosol particle trajectory simulations and collision rate constant modeling