Q: a molecular dynamics program for free energy calculations and empirical valence bond simulations in biomolecular systems.

A new molecular dynamics program for free energy calculations in biomolecular systems is presented. It is principally designed for free energy perturbation simulations, empirical valence bond calculations, and binding affinity estimation by linear interaction energy methods. Evaluation of ligand-binding selectivity and free energy profiles for nucleophile activation in two protein tyrosine phosphatases as well as absolute binding affinity estimation for a lysine-binding protein are given as examples.