Abstract:Approaches for computing small molecule binding free energies based on molecular simulations are now regularly being employed by academic and industry practitioners to study receptor-ligand systems and prioritize the synthesis of small molecules for ligand design. Given the variety of methods and implementations available, it is natural to ask how the convergence rates and final predictions of these methods compare. In this study, we describe the concept and results for the SAMPL6 SAMPLing challenge, the first challenge from the SAMPL series focusing on the assessment of convergence properties and reproducibility of binding free energy methodologies. We provided parameter files, partial charges, and multiple initial geometries for two octa-acid (OA) and one cucurbit[8]uril (CB8) host-guest systems. Participants submitted binding free energy predictions as a function of the number of force and energy evaluations for seven different alchemical and physical-pathway (i.e., potential of mean force and weighted ensemble of trajectories) methodologies implemented with the GROMACS, AMBER, NAMD, or OpenMM simulation engines. To rank the methods, we developed an efficiency statistic based on bias and variance of the free energy estimates. For the two small OA binders, the free energy estimates computed with alchemical and potential of mean force approaches show relatively similar variance and bias as a function of the number of energy/force evaluations, with the attach-pull-release (APR), GROMACS expanded ensemble, and NAMD double decoupling submissions obtaining the greatest efficiency. The differences between the methods increase when analyzing the CB8-quinine system, where both the guest size and correlation times for system dynamics are greater. For this system, nonequilibrium switching (GROMACS/NS-DS/SB) obtained the overall highest efficiency. Surprisingly, the results suggest that specifying force field parameters and partial charges is insufficient to generally ensure reproducibility, and we observe differences between seemingly converged predictions ranging approximately from 0.3 to 1.0 kcal/mol, even with almost identical simulations parameters and system setup (e.g., Lennard-Jones cutoff, ionic composition). Further work will be required to completely identify the exact source of these discrepancies. Among the conclusions emerging from the data, we found that Hamiltonian replica exchange-while displaying very small variance-can be affected by a slowly-decaying bias that depends on the initial population of the replicas, that bidirectional estimators are significantly more efficient than unidirectional estimators for nonequilibrium free energy calculations for systems considered, and that the Berendsen barostat introduces non-negligible artifacts in expanded ensemble simulations.

SAMPL9 blind predictions using nonequilibrium alchemical approaches

The SAMPL9 Host-Guest Blind Challenge: An Overview of Binding Free Energy Predictive Accuracy

Overview of the SAMPL6 host–guest binding affinity prediction challenge

SAMPL7: Host–guest binding prediction by molecular dynamics and quantum mechanics

The SAMPL6 SAMPLing challenge: assessing the reliability and efficiency of binding free energy calculations

Overview of the SAMPL5 host–guest challenge: Are we doing better?

Taming Multiple Binding Poses in Alchemical Binding Free Energy Prediction: the $β$-cyclodextrin Host-Guest SAMPL9 Blinded Challenge

Free-energy perturbation and quantum mechanical study of SAMPL4 octa-acid host–guest binding energies

Detailed potential of mean force studies on host–guest systems from the SAMPL6 challenge

Integrating Multiscale and Machine Learning Approaches towards the SAMPL9 LogP Challenge

Evaluation of log P, pKa, and log D predictions from the SAMPL7 blind challenge

The Thermodynamics of Guest Complexation to Octa-Acid and Tetra-Endo-methyl Octa-Acid: Reference Data for the Sixth Statistical Assessment of Modeling of Proteins and Ligands (SAMPL6)

Large scale free energy calculations for blind predictions of protein-ligand binding: the D3R Grand Challenge 2015

The SAMPL6 challenge on predicting octanol–water partition coefficients from EC-RISM theory

The SAMPL6 challenge on predicting aqueous pKa values from EC-RISM theory

Impact of domain knowledge on blinded predictions of binding energies by alchemical free energy calculations

Blinded prediction of protein-ligand binding affinity using Amber thermodynamic integration for the 2018 D3R grand challenge 4

Extended solvent-contact model approach to SAMPL4 blind prediction challenge for hydration free energies

Blinded predictions of binding modes and energies of HSP90-α ligands for the 2015 D3R grand challenge

Fast and accurate predictions of binding free energies using MM-PBSA and MM-GBSA