Headed to ACS San Diego? Join us for Happy Hour!
← Back to glossary

MM/PBSA and MM/GBSA

Method
Method
Method

MM/PBSA (Molecular Mechanics/Poisson–Boltzmann Surface Area) and MM/GBSA (Molecular Mechanics/Generalized Born Surface Area) are computational methods used to estimate the binding free energy of biomolecular complexes, such as protein–ligand or protein–protein interactions. These approaches combine molecular mechanics energies for protein and ligand atoms with solvation terms calculated using either the Poisson–Boltzmann equation (PBSA) or the Generalized Born model (GBSA). PBSA and GBSA model solvent implicitly, thus are more computational efficient than energy perturbation methods that use explicit solvent models. Because protein and ligand atoms are modeled through MM, MM/PBSA and MM/GBSA provide a balance between accuracy and computational efficiency. The GB SA is an approximation of the PB SA, so it is the most efficient, but also the least accurate.

Importance in Computational Drug Discovery:

  • Enables rapid and cost-effective estimation of binding affinities for large sets of compounds after molecular docking or molecular dynamics simulations.
  • Facilitates ranking and prioritization of drug candidates based on predicted binding free energies.
  • Supports lead optimization by quantifying the energetic contributions of specific ligand modifications.
  • Complements experimental binding assays, guiding hypothesis-driven design and reducing experimental workload.
  • Integrates with molecular dynamics workflows to capture dynamic and entropic effects in binding.

Key Tools

  • AMBER: Widely used suite for molecular dynamics and MM/PBSA/MM/GBSA calculations.
  • GROMACS with g_mmpbsa: Supports MM/PBSA and MM/GBSA post-processing of MD trajectories.
  • Schrödinger Prime MM-GBSA: Provides user-friendly MM/GBSA calculations for protein–ligand complexes.

Literature

"The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities"

  • Publication Date: April 2, 2015
  • DOI: 10.1517/17460441.2015.1032936
  • Summary: This comprehensive review discusses the use of MM/PBSA and MM/GBSA methods for estimating ligand-binding affinities. It emphasizes the importance of calibration, testing, and validation, and highlights the methods' applications in reproducing experimental findings and improving virtual screening results.

"Assessing the performance of the MM/PBSA and MM/GBSA methods. 1. The accuracy of binding free energy calculations based on molecular dynamics simulations"

  • Publication Date: December 27, 2010
  • DOI: 10.1021/ci100275a
  • Summary: This study evaluates the accuracy of MM/PBSA and MM/GBSA methods in calculating binding free energies across various protein-ligand systems. It explores factors such as MD simulation length and solute dielectric constants, providing insights into optimizing these parameters for accurate predictions.

"Assessing the performance of MM/PBSA and MM/GBSA methods. 5. Improved docking performance using high solute dielectric constant MM/GBSA and MM/PBSA rescoring"

  • Publication Date: August 18, 2014
  • DOI: 10.1039/C4CP03179B
  • Summary: This paper investigates the impact of varying solute dielectric constants on the performance of MM/PBSA and MM/GBSA methods in docking studies. The findings suggest that higher dielectric constants can enhance the accuracy of binding affinity predictions, particularly in virtual screening applications.

"An Efficient Implementation of the Nwat-MMGBSA Method to Rescore Docking Results in Medium-Throughput Virtual Screenings"

  • Publication Date: February 13, 2018
  • DOI: 10.3389/fchem.2018.00043
  • Summary: This study introduces the Nwat-MMGBSA method, which incorporates explicit water molecules in MM/GBSA calculations to improve the accuracy of binding affinity predictions. The method demonstrates enhanced performance in both protein-protein interactions and ligand-receptor complexes.

"Assessing the performance of MM/PBSA and MM/GBSA methods. 3. The impact of force fields and ligand charge models"

  • Publication Date: July 18, 2013
  • DOI: 10.1021/jp404160y
  • Summary: This paper examines how different force fields and ligand charge models affect the performance of MM/PBSA and MM/GBSA methods. The study provides recommendations for selecting appropriate computational parameters to enhance the reliability of binding free energy calculations.