English
 
User Manual Privacy Policy Disclaimer Contact us
  Advanced SearchBrowse

Item

ITEM ACTIONSEXPORT
  The SAMPL6 SAMPLing challenge: Assessing the reliability and efficiency of binding free energy calculations.

Rizzi, A., Jensen, T., Slochower, D. R., Aldeghi, M., Gapsys, V., Ntekoumes, D., et al. (2020). The SAMPL6 SAMPLing challenge: Assessing the reliability and efficiency of binding free energy calculations. Journal of Computer-Aided Molecular Design, (in press). doi:10.1007/s10822-020-00290-5.

Item is

Basic

show hide
Item Permalink: http://hdl.handle.net/21.11116/0000-0005-991C-0 Version Permalink: http://hdl.handle.net/21.11116/0000-0005-9922-8
Genre: Journal Article

Files

show Files
hide Files
:
3191106.pdf (Preprint), 5MB
Name:
3191106.pdf
Description:
-
Visibility:
Public
MIME-Type / Checksum:
application/pdf / [MD5]
Technical Metadata:
Copyright Date:
-
Copyright Info:
-
License:
-
:
3191106_Suppl.pdf (Supplementary material), 8MB
Name:
3191106_Suppl.pdf
Description:
-
Visibility:
Public
MIME-Type / Checksum:
application/pdf / [MD5]
Technical Metadata:
Copyright Date:
-
Copyright Info:
-
License:
-

Locators

show

Creators

show
hide
 Creators:
Rizzi, A., Author
Jensen, T., Author
Slochower, D. R., Author
Aldeghi, M.1, Author              
Gapsys, V.1, Author              
Ntekoumes, D., Author
Bosisio, S., Author
Papadourakis, M., Author
Henriksen, N. M., Author
de Groot, B. L.1, Author              
Cournia, Z., Author
Dickson, A., Author
Michel, J., Author
Gilson, M. K., Author
Shirts, M. R., Author
Mobley, D. L., Author
Chodera, J. D., Author
Affiliations:
1Research Group of Computational Biomolecular Dynamics, MPI for biophysical chemistry, Max Planck Society, ou_578573              

Content

show
hide
Free keywords: SAMPL6; Host-guest; Binding affinity; Free energy calculations; Cucurbit[8]uril; Octa-acid; Sampling
 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.

Details

show
hide
Language(s): eng - English
 Dates: 2020-01-27
 Publication Status: Published online
 Pages: -
 Publishing info: -
 Table of Contents: -
 Rev. Method: Peer
 Identifiers: DOI: 10.1007/s10822-020-00290-5
 Degree: -

Event

show

Legal Case

show

Project information

show

Source 1

show
hide
Title: Journal of Computer-Aided Molecular Design
Source Genre: Journal
 Creator(s):
Affiliations:
Publ. Info: -
Pages: - Volume / Issue: - Sequence Number: (in press) Start / End Page: - Identifier: -