Automated Force Field Developer and Optimizer Platform: Torsion Reparameterization
Authors:
- Alejandro Blanco-GonzalezABAlejandro Blanco-Gonzalez
- Shi ZhangSZShi Zhang
- Timothy J. GieseTGTimothy J. Giese
- Zeke A. PiskulichZPZeke A. Piskulich
- Andreas W. GötzAGAndreas W. Götz
- Kenneth M. MerzKMKenneth M. Merz
- Darrin M. YorkDYDarrin M. York
- Hasan Metin AktulgaHAHasan Metin Aktulga
- Madushanka ManathungaMMMadushanka Manathunga
Open Access
- 9 March 2026
- journal article
- Published by American Chemical Society (ACS) in Journal of Chemical Information and Modeling
- Vol. 66 (6) , 3206-3219
Abstract
General force fields such as General Amber Force Field (GAFF) have been designed for broad applicability and are widely used in protein-ligand binding simulations in structure-based drug discovery. However, the force field parameters are not always transferable across ligand molecules, and custom reparameterization is sometimes necessary for accurate binding free energy simulations. This is especially true for torsion parameters, which are highly dependent on stereoelectronic and steric effects. Here, we report a novel, flexible, and user-friendly computational tool called the Automated Force Field Developer and Optimizer (AFFDO) platform that allows generating accurate, tailored GAFF2 torsion parameters for drug-like molecules. For a given ligand, AFFDO selects the most important torsions, carries out GPU-accelerated density functional theory calculations to collect reference data and fits torsion terms using a fast gradient-based optimizer that leverages automated differentiation. We benchmark AFFDO by parametrizing a series of drug-like molecules and carrying out protein-ligand relative binding free energy (RBFE) simulations. The results show that AFFDO can significantly improve GAFF2 torsion parameters against QM reference data, which in some cases translates into better agreement with experimental RBFE values within a reasonable computational time.Keywords
All Related Versions
Funding Information
- National Institute of General Medical Sciences (R01GM130641)
- National Institute of General Medical Sciences (R44GM150347)
- National Science Foundation (2209717)
- National Science Foundation (2209718)
- National Science Foundation (2435622)
- East Lansing Technology Innovation Center
- Lansing Regional SmartZone
This publication has 66 references indexed in Scilit:
- ParaMol: A Package for Automatic Parameterization of Molecular Mechanics Force FieldsJournal of Chemical Information and Modeling, 2021
- Extended tight‐binding quantum chemistry methodsWIREs Computational Molecular Science, 2020
- Lessons learned from comparing molecular dynamics engines on the SAMPL5 datasetJournal of Computer-Aided Molecular Design, 2016
- Paramfit: Automated optimization of force field parameters for molecular dynamics simulationsJournal of Computational Chemistry, 2014
- Rapid parameterization of small molecules using the force field toolkitJournal of Computational Chemistry, 2013
- Automatic atom type and bond type perception in molecular mechanical calculationsJournal of Molecular Graphics and Modelling, 2006
- Junmei Wang, Romain M. Wolf, James W. Caldwell, Peter A. Kollman, and David A. Case, "Development and testing of a general amber force field"Journal of Computational Chemistry(2004) 25(9) 1157–1174Journal of Computational Chemistry, 2004
- Development and testing of a general amber force fieldJournal of Computational Chemistry, 2004
- Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: II. Parameterization and validationJournal of Computational Chemistry, 2002
- Fast, efficient generation of high-quality atomic charges. AM1-BCC model: I. MethodJournal of Computational Chemistry, 2000