GAFF/IPolQ-Mod+LJ-Fit: Optimized Force Field Parameters for Solvation Free Energy Predictions

Overview

Journal ADMET DMPK

Specialty Pharmacology

Date 2022 Mar 18

PMID 35300308

Authors

Andreas Mecklenfeld

Gabriele Raabe

Affiliations

Soon will be listed here.

Abstract

Rational drug design featuring explicit solubility considerations can greatly benefit from molecular dynamics simulations, as they allow for the prediction of the Gibbs free energy of solvation and thus relative solubilities. In our previous work (A. Mecklenfeld, G. Raabe. J. Chem. Theory Comput. no. 12 (2017) 6266-6274), we have compared solvation free energy results obtained with the General Amber Force Field (GAFF) and its default restrained electrostatic potential (RESP) partial charges to those obtained by modified implicitly polarized charges (IPolQ-Mod) for an implicit representation of impactful polarization effects. In this work, we have adapted Lennard-Jones parameters for GAFF atom types in combination with IPolQ-Mod to further improve the accuracies of solvation free energy and liquid density predictions. We thereby focus on prominent atom types in common drugs. For the refitting, 357 respectively 384 systems were considered for free energies and densities and validation was performed for 142 free energies and 100 densities of binary mixtures. By the in-depth comparison of simulation results for default GAFF, GAFF with IPolQ-Mod and our new set of parameters, which we label GAFF/IPolQ-Mod+LJ-Fit, we can clearly highlight the improvements of our new model for the description of both relative solubilities and fluid phase behaviour.

Citing Articles

Flavones in pomelo peel resist fibril formation of human islet amyloid polypeptide.

Gao C, Wan Z, Liu Y, Meng Y, Chen X, Tang X Chin Herb Med. 2025; 17(1):166-177.

PMID: 39949806 PMC: 11814264. DOI: 10.1016/j.chmed.2024.02.002.

Prediction of toluene/water partition coefficients of SAMPL9 compounds: comparison of the molecular dynamics force fields GAFF/RESP and GAFF/IPolQ-Mod + LJ-fit.

Sprick M, Raabe G Phys Chem Chem Phys. 2024; 26(4):3126-3138.

PMID: 38189577 PMC: 10806620. DOI: 10.1039/d3cp04149b.

Strategies of solubility enhancement and perspectives in solubility measurements of pharmaceutical compounds.

Bergstrom C, Llinas A ADMET DMPK. 2022; 8(3):176-179.

PMID: 35300302 PMC: 8915603. DOI: 10.5599/admet.910.

References

van der Spoel D, Lindahl E, Hess B, Groenhof G, Mark A, Berendsen H . GROMACS: fast, flexible, and free. J Comput Chem. 2005; 26(16):1701-18. DOI: 10.1002/jcc.20291. View

Wang L, Martinez T, Pande V . Building Force Fields: An Automatic, Systematic, and Reproducible Approach. J Phys Chem Lett. 2015; 5(11):1885-91. PMC: 9649520. DOI: 10.1021/jz500737m. View

Klimovich P, Shirts M, Mobley D . Guidelines for the analysis of free energy calculations. J Comput Aided Mol Des. 2015; 29(5):397-411. PMC: 4420631. DOI: 10.1007/s10822-015-9840-9. View

Vanommeslaeghe K, Hatcher E, Acharya C, Kundu S, Zhong S, Shim J . CHARMM general force field: A force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fields. J Comput Chem. 2009; 31(4):671-90. PMC: 2888302. DOI: 10.1002/jcc.21367. View

Martinez L, Andrade R, Birgin E, Martinez J . PACKMOL: a package for building initial configurations for molecular dynamics simulations. J Comput Chem. 2009; 30(13):2157-64. DOI: 10.1002/jcc.21224. View

Muddana H, Sapra N, Fenley A, Gilson M . The SAMPL4 hydration challenge: evaluation of partial charge sets with explicit-water molecular dynamics simulations. J Comput Aided Mol Des. 2014; 28(3):277-87. PMC: 4006311. DOI: 10.1007/s10822-014-9714-6. View

Sousa da Silva A, Vranken W . ACPYPE - AnteChamber PYthon Parser interfacE. BMC Res Notes. 2012; 5:367. PMC: 3461484. DOI: 10.1186/1756-0500-5-367. View

Parameswaran S, Mobley D . Box size effects are negligible for solvation free energies of neutral solutes. J Comput Aided Mol Des. 2014; 28(8):825-9. DOI: 10.1007/s10822-014-9766-7. View

Peterson K, Shepler B, Figgen D, Stoll H . On the spectroscopic and thermochemical properties of ClO, BrO, IO, and their anions. J Phys Chem A. 2006; 110(51):13877-83. DOI: 10.1021/jp065887l. View

10.

Mecklenfeld A, Raabe G . Applicability of a thermodynamic cycle approach for a force field parametrization targeting non-aqueous solvation free energies. J Comput Aided Mol Des. 2019; 34(1):71-82. DOI: 10.1007/s10822-019-00261-5. View

11.

Vlcek L, Chialvo A . Rigorous force field optimization principles based on statistical distance minimization. J Chem Phys. 2015; 143(14):144110. DOI: 10.1063/1.4932360. View

12.

Mobley D, Wymer K, Lim N, Guthrie J . Blind prediction of solvation free energies from the SAMPL4 challenge. J Comput Aided Mol Des. 2014; 28(3):135-50. PMC: 4006301. DOI: 10.1007/s10822-014-9718-2. View

13.

Lemkul J, Huang J, Roux B, MacKerell Jr A . An Empirical Polarizable Force Field Based on the Classical Drude Oscillator Model: Development History and Recent Applications. Chem Rev. 2016; 116(9):4983-5013. PMC: 4865892. DOI: 10.1021/acs.chemrev.5b00505. View

14.

Pham T, Shirts M . Identifying low variance pathways for free energy calculations of molecular transformations in solution phase. J Chem Phys. 2011; 135(3):034114. DOI: 10.1063/1.3607597. View

15.

Hess B, Kutzner C, van der Spoel D, Lindahl E . GROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. J Chem Theory Comput. 2015; 4(3):435-47. DOI: 10.1021/ct700301q. View

16.

Wang J, Wolf R, Caldwell J, Kollman P, Case D . Development and testing of a general amber force field. J Comput Chem. 2004; 25(9):1157-74. DOI: 10.1002/jcc.20035. View

17.

Patel S, Brooks 3rd C . CHARMM fluctuating charge force field for proteins: I parameterization and application to bulk organic liquid simulations. J Comput Chem. 2003; 25(1):1-15. DOI: 10.1002/jcc.10355. View

18.

Shenfeld D, Xu H, Eastwood M, Dror R, Shaw D . Minimizing thermodynamic length to select intermediate states for free-energy calculations and replica-exchange simulations. Phys Rev E Stat Nonlin Soft Matter Phys. 2009; 80(4 Pt 2):046705. DOI: 10.1103/PhysRevE.80.046705. View

19.

Mayne C, Saam J, Schulten K, Tajkhorshid E, Gumbart J . Rapid parameterization of small molecules using the Force Field Toolkit. J Comput Chem. 2013; 34(32):2757-70. PMC: 3874408. DOI: 10.1002/jcc.23422. View

20.

Villa A, Mark A . Calculation of the free energy of solvation for neutral analogs of amino acid side chains. J Comput Chem. 2002; 23(5):548-53. DOI: 10.1002/jcc.10052. View