title: Calculating the sensitivity and robustness of binding free energy calculations to force field parameters.
abstract: Binding free energy calculations offer a thermodynamically rigorous method to compute protein-ligand binding, and they depend on empirical force fields with hundreds of parameters. We examined the sensitivity of computed binding free energies to the ligand's electrostatic and van der Waals parameters. Dielectric screening and cancellation of effects between ligand-protein and ligand-solvent interactions reduce the parameter sensitivity of binding affinity by 65%, compared with interaction strengths computed in the gas-phase. However, multiple changes to parameters combine additively on average, which can lead to large changes in overall affinity from many small changes to parameters. Using these results, we estimate that random, uncorrelated errors in force field nonbonded parameters must be smaller than 0.02 e per charge, 0.06 Å per radius, and 0.01 kcal/mol per well depth in order to obtain 68% (one standard deviation) confidence that a computed affinity for a moderately-sized lead compound will fall within 1 kcal/mol of the true affinity, if these are the only sources of error considered.
authors.affiliation: Department of Pharmaceutical Chemistry, University of California San Francisco, 1700 4 St, San Francisco California 94143-2550, USA ; Biophysics Graduate Program, University of California San Francisco, 1700 4 St, San Francisco California 94143-2550, USA.
authors.forename: Gabriel J
authors.initials: GJ
authors.lastname: Rocklin
authors.forename: David L
authors.initials: DL
authors.lastname: Mobley
authors.forename: Ken A
authors.initials: KA
authors.lastname: Dill
keywords: force fields; free energy calculations; ligand binding; molecular dynamics; parameter sensitivity