Alchemical Free Energy Estimators and Molecular Dynamics Engines: Accuracy, Precision, and Reproducibility
Journal Article
·
· Journal of Chemical Theory and Computation
- Univ. College London (United Kingdom)
- Univ. College London (United Kingdom); Univ. of Amsterdam (Netherlands)
The binding free energy between a ligand and its target protein is an essential quantity to know at all stages of the drug discovery pipeline. Assessing this value computationally can offer insight into where efforts should be focused in the pursuit of effective therapeutics to treat a myriad of diseases. In this work, we examine the computation of alchemical relative binding free energies with an eye for assessing reproducibility across popular molecular dynamics packages and free energy estimators. The focus of this work is on 54 ligand transformations from a diverse set of protein targets: MCL1, PTP1B, TYK2, CDK2, and thrombin. These targets are studied with three popular molecular dynamics packages: OpenMM, NAMD2, and NAMD3 alpha. Trajectories collected with these packages are used to compare relative binding free energies calculated with thermodynamic integration and free energy perturbation methods. The resulting binding free energies show good agreement between molecular dynamics packages with an average mean unsigned error between them of 0.50 kcal/mol. The correlation between packages is very good, with the lowest Spearman’s, Pearson’s and Kendall’s tau correlation coefficients being 0.92, 0.91, and 0.76, respectively. Agreement between thermodynamic integration and free energy perturbation is shown to be very good when using ensemble averaging.
- Research Organization:
- Argonne National Laboratory (ANL), Argonne, IL (United States). Argonne Leadership Computing Facility (ALCF)
- Sponsoring Organization:
- National Science Foundation (NSF); USDOE
- Grant/Contract Number:
- AC02-06CH11357
- OSTI ID:
- 2470030
- Journal Information:
- Journal of Chemical Theory and Computation, Journal Name: Journal of Chemical Theory and Computation Journal Issue: 6 Vol. 18; ISSN 1549-9618
- Publisher:
- American Chemical SocietyCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Similar Records
Protein–Ligand Binding Free-Energy Calculations with ARROW–A Purely First-Principles Parameterized Polarizable Force Field
Studying the Collective Functional Response of a Receptor in Alchemical Ligand Binding Free Energy Simulations with Accelerated Solvation Layer Dynamics
Reduced Free Energy Perturbation/Hamiltonian Replica Exchange Molecular Dynamics Method with Unbiased Alchemical Thermodynamic Axis
Journal Article
·
Thu Dec 01 19:00:00 EST 2022
· Journal of Chemical Theory and Computation
·
OSTI ID:2439559
Studying the Collective Functional Response of a Receptor in Alchemical Ligand Binding Free Energy Simulations with Accelerated Solvation Layer Dynamics
Journal Article
·
Tue Apr 02 20:00:00 EDT 2024
· Journal of Chemical Theory and Computation
·
OSTI ID:2474889
Reduced Free Energy Perturbation/Hamiltonian Replica Exchange Molecular Dynamics Method with Unbiased Alchemical Thermodynamic Axis
Journal Article
·
Mon Sep 24 20:00:00 EDT 2018
· Journal of Physical Chemistry. B, Condensed Matter, Materials, Surfaces, Interfaces and Biophysical Chemistry
·
OSTI ID:1488538