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.

Atomevo: a web server combining protein modelling, docking, molecular dynamic simulation and MMPBSA analysis of Candida antarctica lipase B (CalB) fusion protein

Quantitative prediction of enantioselectivity of Candida antarctica lipase B by combining docking simulations and quantitative structure-activity relationship (QSAR) analysis

The HADDOCK2.4 web server for integrative modeling of biomolecular complexes

fastDRH: a webserver to predict and analyze protein-ligand complexes based on molecular docking and MM/PB(GB)SA computation

HawkDock: a Web Server to Predict and Analyze the Protein-Protein Complex Based on Computational Docking and MM/GBSA.

ZDOCK server: interactive docking prediction of protein-protein complexes and symmetric multimers

CovalentDock Cloud: a Web Server for Automated Covalent Docking.

Modeling of protein-peptide interactions using the CABS-dock web server for binding site search and flexible docking

farPPI: a webserver for accurate prediction of protein-ligand binding structures for small-molecule PPI inhibitors by MM/PB(GB)SA methods

MemBuilder: a web-based graphical interface to build heterogeneously mixed membrane bilayers for the GROMACS biomolecular simulation program

MMAEVe: a Versatile Program for Constructing Systems for Molecular Dynamics Simulations and Visualization

The ACPYPE web server for small-molecule MD topology generation

MoViES: Molecular Vibrations Evaluation Server for Analysis of Fluctuational Dynamics of Proteins and Nucleic Acids

A Combination of Computational and Experimental Approaches to Investigate the Binding Behavior of B.sub Lipase A Mutants with Substrate Pnpp

MDockPeP: An ab‐initio protein–peptide docking server

PIIMS Server: A Web Server for Mutation Hotspot Scanning at the Protein–Protein Interface

EnzyHTP: A High-Throughput Computational Platform for Enzyme Modeling

CHARMM-GUI PDB Reader and Manipulator: Covalent Ligand Modeling and Simulation

A resource for benchmarking the usefulness of protein structure models

MolDy: molecular dynamics simulation made easy

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