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.<h3 class="u-h4 u-margin-m-top u-margin-xs-bottom">Program summary</h3>Program Title: GROMACS-VI-ExtensionCPC Library link to program files: <a href="https://doi.org/10.17632/7yvc8mmnyv.1">https://doi.org/10.17632/7yvc8mmnyv.1</a>Developer's respository link: <a href="https://www.mpibpc.mpg.de/gromacs-vi-extension">https://www.mpibpc.mpg.de/gromacs-vi-extension</a> and <a href="https://www.gitlab.gwdg.de/martin.reinhardt/gromacs-vi-extension">https://www.gitlab.gwdg.de/martin.reinhardt/gromacs-vi-extension</a>Licensing provisions: LGPLProgramming language: C++14, CUDANature of problem: The free energy difference between two states of a thermodynamic system is calculated using samples generated by simulations based on atomistic Hamiltonians. Due to the high dimensionality of many applications as in, e.g., biophysics, only a small part of the configuration space can be sampled. The choice of the sampling scheme critically affects the accuracy of the final free energy estimate. The challenge is, therefore, to find the optimal sampling scheme that provides best accuracy for given computational effort.Solution method: Sampling is commonly conducted in intermediate states, whose Hamiltonians are defined based on the Hamiltonians of the two states of interest. Here, sampling is conducted in the variationally derived intermediates states that, under the assumption of uncorrelated sample points, yield optimal accuracy. These intermediates differ fundamentally from the common intermediates in that they are non-pairwise, i.e., the forces on a particle are only additive in the end state, whereas the total force in the intermediate states cannot be split into additive contributions from the surrounding particles.References[1] M. Reinhardt, H. Grubmüller, Determining Free-Energy Differences Through Variationally Derived Intermediates, Journal of Chemical Theory and Computation 16 (6) 3504-3512 (2020).

METAGUI 3: A graphical user interface for choosing the collective variables in molecular dynamics simulations

Molcontroller: a VMD Graphical User Interface for Manipulating Molecules

GROMACS implementation of free energy calculations with non-pairwise Variationally derived Intermediates

PLUMED-GUI: an environment for the interactive development of molecular dynamics analysis and biasing scripts

CHARMM-GUI Ligand Binder for Absolute Binding Free Energy Calculations and Its Application.

PACKMOL- GUI: An All-in-One VMD Interface for Efficient Molecular Packing

CHARMM-GUI PDB manipulator for advanced modeling and simulations of proteins containing nonstandard residues.

Geo-Measures: A PyMOL plugin for protein structure ensembles analysis

Grid-based state space exploration for molecular binding

CHARMM-GUI Free Energy Calculator for Absolute and Relative Ligand Solvation and Binding Free Energy Simulations

Metis - A Python-Based User Interface to Collect Expert Feedback for Generative Chemistry Models

ChemDyME: Kinetically Steered, Automated Mechanism Generation through Combined Molecular Dynamics and Master Equation Calculations

MolDy: molecular dynamics simulation made easy

StreaMD: the toolkit for high-throughput molecular dynamics simulations

MIST: A simple and efficient molecular dynamics abstraction library for integrator development

CHARMM-GUI Free Energy Calculator for Practical Ligand Binding Free Energy Simulations with AMBER

VMD as a Platform for Interactive Small Molecule Preparation and Visualization in Quantum and Classical Simulations

MEIGO: an open-source software suite based on metaheuristics for global optimization in systems biology and bioinformatics

Metadynamics: a method to simulate rare events and reconstruct the free energy in biophysics, chemistry and material science

Gaussian accelerated molecular dynamics: Principles and applications