English
 
Help Privacy Policy Disclaimer
  Advanced SearchBrowse

Item

ITEM ACTIONSEXPORT

Released

Journal Article

GROMACS implementation of free energy calculations with non-pairwise variationally derived intermediates

MPS-Authors
/persons/resource/persons231604

Reinhardt,  M.
Department of Theoretical and Computational Biophysics, MPI for Biophysical Chemistry, Max Planck Society;

/persons/resource/persons15155

Grubmüller,  H.
Department of Theoretical and Computational Biophysics, MPI for biophysical chemistry, Max Planck Society;

External Resource
No external resources are shared
Fulltext (public)

3264317.pdf
(Publisher version), 2MB

Supplementary Material (public)
There is no public supplementary material available
Citation

Reinhardt, M., & Grubmüller, H. (2021). GROMACS implementation of free energy calculations with non-pairwise variationally derived intermediates. Computer Physics Communications, 264: 107931. doi:10.1016/j.cpc.2021.107931.


Cite as: http://hdl.handle.net/21.11116/0000-0007-68CD-E
Abstract
Gradients in free energies are the driving forces of physical and biochemical systems. To predict free energy differences with high accuracy, Molecular Dynamics (MD) and other methods based on atomistic Hamiltonians conduct sampling simulations in intermediate thermodynamic states that bridge the configuration space densities between two states of interest ('alchemical transformations'). For uncorrelated sampling, the recent Variationally derived Intermediates (VI) method yields optimal accuracy. The form of the VI intermediates differs fundamentally from conventional ones in that they are non-pairwise, i.e., the total force on a particle in an intermediate states cannot be split into additive contributions from the surrounding particles. In this work, we describe the implementation of VI into the widely used GROMACS MD software package (2020, version 1). Furthermore, a variant of VI is developed that avoids numerical instabilities for vanishing particles. The implementation allows the use of previous non-pairwise potential forms in the literature, which have so far not been available in GROMACS. Example cases on the calculation of solvation free energies, and accuracy assessments thereof, are provided.