Abstract:Alchemical free energy methods, such as free energy perturbation (FEP) and thermodynamic integration (TI), become increasingly popular and crucial for drug design and discovery. However, the system preparation of alchemical free energy simulation is an error-prone, time-consuming, and tedious process for a large number of ligands. To address this issue, we have recently presented CHARMM-GUI Free Energy Calculator that can provide input and postprocessing scripts for NAMD and GENESIS FEP molecular dynamics systems. In this work, we extended three submodules of Free Energy Calculator to work with the full suite of GPU-accelerated alchemical free energy methods and tools in AMBER, including input and postprocessing scripts. The BACE1 (β-secretase 1) benchmark set was used to validate the AMBER-TI simulation systems and scripts generated by Free Energy Calculator. The overall results of relatively large and diverse systems are almost equivalent with different protocols (unified and split) and with different timesteps (1, 2, and 4 fs), with R2 > 0.9. More importantly, the average free energy differences between two protocols are small and reliable with four independent runs, with a mean unsigned error (MUE) below 0.4 kcal/mol. Running at least four independent runs for each pair with AMBER20 (and FF19SB/GAFF2.1/OPC force fields), we obtained a MUE of 0.99 kcal/mol and root-mean-square error of 1.31 kcal/mol for 58 alchemical transformations in comparison with experimental data. In addition, a set of ligands for T4-lysozyme was used to further validate our free energy calculation protocol whose results are close to experimental data (within 1 kcal/mol). In summary, Free Energy Calculator provides a user-friendly web-based tool to generate the AMBER-TI system and input files for high-throughput binding free energy calculations with access to the full set of GPU-accelerated alchemical free energy, enhanced sampling, and analysis methods in AMBER.The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jcim.1c00747.Table S1, 58 ΔΔGbind pathways for 36 ligands; detailed ΔΔGbind values with unified and split protocols and with three timesteps; MUE and RMSE calculations; and ΔGbind calculations (PDF)Experimental ΔΔGbind values (XLSX)This article has not yet been cited by other publications.

Free energy simulations with the AMOEBA polarizable force field and metadynamics on GPU platform

Amber Free Energy Tools: Interoperable Software for Free Energy Simulations Using Generalized Quantum Mechanical/molecular Mechanical and Machine Learning Potentials.

Tinker-OpenMM: Absolute and relative alchemical free energies using AMOEBA on GPUs

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

Fast and flexible gpu accelerated binding free energy calculations within the amber molecular dynamics package

Exploring Biomolecular Conformational Dynamics with Polarizable Force Field AMOEBA and Enhanced Sampling Method Milestoning

Accurate Free Energy Calculation via Multiscale Simulations Driven by Hybrid Machine Learning and Molecular Mechanics Potentials

Quantum Fragment Based ab Initio Molecular Dynamics for Proteins.

Assessing many-body contributions to intermolecular interactions of the AMOEBA force field using energy decomposition analysis of electronic structure calculations.

CHARMM-GUI Free Energy Calculator for Practical Ligand Binding Free Energy Simulations with AMBER

A Force Balanced Fragmentation Method for Ab Initio Molecular Dynamic Simulation of Protein

Using AMBER18 for Relative Free Energy Calculations

Development and Comprehensive Benchmark of a High Quality AMBER-Consistent Small Molecule Force Field with Broad Chemical Space Coverage for Molecular Modeling and Free Energy Calculation

Assessing Ion-Water Interactions in the AMOEBA Force Field Using Energy Decomposition Analysis of Electronic Structure Calculations.

Boosting Free-Energy Perturbation Calculations with GPU-Accelerated NAMD

Prediction of hydration free energies for the SAMPL4 data set with the AMOEBA polarizable force field

Comparing the Performance of Different AMBER Protein Forcefields, Partial Charge Assignments, and Water Models for Absolute Binding Free Energy Calculations

Performance of the AMOEBA Water Model in the Vicinity of QM Solutes: A Diagnosis Using Energy Decomposition Analysis.

Free Energy Calculations Using the Movable Type Method with Molecular Dynamics Driven Protein-Ligand Sampling

AMBER Drug Discovery Boost Tools: Automated Workflow for Production Free-Energy Simulation Setup and Analysis (ProFESSA)