A High-Throughput Steered Molecular Dynamics Study on the Free Energy Profile of Ion Permeation through Gramicidin A

Abstract

Steered molecular dynamics (SMD) simulations for the calculation of free energies are well suited for high-throughput molecular simulations on a distributed infrastructure due to the simplicity of the setup and parallel granularity of the runs. However, so far, the computational cost limited the estimation of the free energy typically over just a few pullings, thus impeding the evaluation of statistical uncertainties involved. In this work, we performed two thousand pulls for the permeation of a potassium ion in the gramicidin A pore by all-atom molecular dynamics in order to assess the bidirectional SMD protocol with a proper amount of sampling. The estimated free energy profile still shows a statistical error of several kcal/mol, while the work distributions are estimated to be non-Gaussian at pulling speeds of 10 Å/ns. We discuss the methodology and the confidence intervals in relation to increasing amounts of computed trajectories and how different permeation pathways for the potassium ion, knock-on and sideways, affect the sampling and the free energy estimation.

arXiv commentThis document is the Accepted Manuscript version of a Published Work that appeared in final form in J. Chem. Theory Comput., copyright 2011 after peer review and technical editing by the publisher. To access the final edited and published work see doi:10.1021/ct100707s

References in corpus (9)

• A nonequilibrium equality for free energy differences
• The Entropy Production Fluctuation Theorem and the Nonequilibrium Work Relation for Free Energy Differences
• Verification of the Crooks fluctuation theorem and recovery of RNA folding free energies
• Path ensembles averages in systems driven far-from-equilibrium
• Comparison of free energy methods for molecular systems
• Calculating potentials of mean force and diffusion coefficients from nonequilibirum processes without Jarzynski's equality
• Optimized free energies from bidirectional single-molecule force spectroscopy
• Calculating free energy profiles in biomolecular systems from fast non-equilibrium processes
• Determination of the chemical potential using energy-biased sampling

Cited by in corpus (5)

• METAGUI 3: a graphical user interface for choosing the collective variables in molecular dynamics simulations
• PLUMED-GUI: an environment for the interactive development of molecular dynamics analysis and biasing scripts
• Predicting ion channel conductance via dissipation-corrected targeted molecular dynamics and Langevin equation simulations
• Self-Organizing Maps of Unbiased Ligand-Target Binding Pathways and Kinetics
• Nonequilibrium path-ensemble averages for symmetric protocols