title: Comparison of thermodynamic integration and Bennett acceptance ratio for calculating relative protein-ligand binding free energies.
abstract: The performances of Bennett's acceptance ratio method and thermodynamic integration (TI) for the calculation of free energy differences in protein simulations are compared. For the latter, the standard trapezoidal rule, Simpson's rule, and Clenshaw-Curtis integration are used as numerical integration methods. We evaluate the influence of the number and definition of intermediate states on the precision, accuracy, and efficiency of the free energy calculations. Our results show that non-equidistantly spaced intermediate states are in some cases beneficial for the TI methods. Using several combinations of softness parameters and the λ power dependence, it is shown that these benefits are strongly dependent on the shape of the integrand. Although TI is more user-friendly due to its simplicity, it was found that Bennett's acceptance ratio method is the more efficient method. It is also the least dependent on the choice of the intermediate states, making it more robust than TI.
authors: Institute for Molecular Modeling and Simulation, BOKU-University of Natural Resources and Life Sciences, Muthgasse 18, 1190 Vienna, Austria.
authors: Anita
authors: A
authors: de Ruiter
authors: Stefan
authors: S
authors: Boresch
authors: Chris
authors: C
authors: Oostenbrink
mesh_terms: D001550:Benzamidines / Q000737:chemistry; D008956:Models, Chemical; D056004:Molecular Dynamics Simulation; D011485:Protein Binding; D013816:Thermodynamics; D014361:Trypsin Inhibitors / Q000737:chemistry
chemical_list: D001550:Benzamidines; D014361:Trypsin Inhibitors; C032157:benzamidine