Gmx_MMPBSA: A New Tool to Perform End-State Free Energy Calculations with GROMACS

Mario S. Valdés-Tresanco, Mario E. Valdés-Tresanco, Pedro A. Valiente, Ernesto Moreno

Research output: Contribution to journalArticlepeer-review

148 Scopus citations


Molecular mechanics/Poisson-Boltzmann (Generalized-Born) surface area is one of the most popular methods to estimate binding free energies. This method has been proven to balance accuracy and computational efficiency, especially when dealing with large systems. As a result of its popularity, several programs have been developed for performing MM/PB(GB)SA calculations within the GROMACS community. These programs, however, present several limitations. Here we present gmx_MMPBSA, a new tool to perform end-state free energy calculations from GROMACS molecular dynamics trajectories. gmx_MMPBSA provides the user with several options, including binding free energy calculations with different solvation models (PB, GB, or 3D-RISM), stability calculations, computational alanine scanning, entropy corrections, and binding free energy decomposition. Noteworthy, several promising methodologies to calculate relative binding free energies such as alanine scanning with variable dielectric constant and interaction entropy have also been implemented in gmx_MMPBSA. Two additional tools - gmx_MMPBSA_test and gmx_MMPBSA_ana - have been integrated within gmx_MMPBSA to improve its usability. Multiple illustrating examples can be accessed through gmx_MMPBSA_test, while gmx_MMPBSA_ana provides fast, easy, and efficient access to different graphics plotted from gmx_MMPBSA output files. The latest version (v1.4.3, 26/05/2021) is available free of charge (documentation, test files, and tutorials included) at

Original languageEnglish
Pages (from-to)6281-6291
Number of pages11
JournalJournal of Chemical Theory and Computation
Issue number10
StatePublished - 12 Oct 2021


Dive into the research topics of 'Gmx_MMPBSA: A New Tool to Perform End-State Free Energy Calculations with GROMACS'. Together they form a unique fingerprint.

Cite this