title: Overcoming dissipation in the calculation of standard binding free energies by ligand extraction.
abstract: This article addresses calculations of the standard free energy of binding from molecular simulations in which a bound ligand is extracted from its binding site by steered molecular dynamics (MD) simulations or equilibrium umbrella sampling (US). Host-guest systems are used as test beds to examine the requirements for obtaining the reversible work of ligand extraction. We find that, for both steered MD and US, marked irreversibilities can occur when the guest molecule crosses an energy barrier and suddenly jumps to a new position, causing dissipation of energy stored in the stretched molecule(s). For flexible molecules, this occurs even when a stiff pulling spring is used, and it is difficult to suppress in calculations where the spring is attached to the molecules by single, fixed attachment points. We, therefore, introduce and test a method, fluctuation-guided pulling, which adaptively adjusts the spring's attachment points based on the guest's atomic fluctuations relative to the host. This adaptive approach is found to substantially improve the reversibility of both steered MD and US calculations for the present systems. The results are then used to estimate standard binding free energies within a comprehensive framework, termed attach-pull-release, which recognizes that the standard free energy of binding must include not only the pulling work itself, but also the work of attaching and then releasing the spring, where the release work includes an accounting of the standard concentration to which the ligand is discharged.
authors.affiliation: Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, 9500 Gilman Dr. MC 0736, La Jolla, California, 92093-0736.
authors.forename: Camilo
authors.initials: C
authors.lastname: Velez-Vega
authors.forename: Michael K
authors.initials: MK
authors.lastname: Gilson
mesh_terms: D001643:Bridged Bicyclo Compounds / Q000737:chemistry*; D001952:Bridged-Ring Compounds / Q000737:chemistry*; D007093:Imidazoles / Q000737:chemistry*; D008024:Ligands; D056004:Molecular Dynamics Simulation; D009795:Octanes / Q000737:chemistry*; D013096:Spermine / Q000737:chemistry*; D013816:Thermodynamics; D047028:Macrocyclic Compounds; D000099246:Cucurbit(n)urils
chemical_list: D001643:Bridged Bicyclo Compounds; D001952:Bridged-Ring Compounds; D007093:Imidazoles; D008024:Ligands; D009795:Octanes; D013096:Spermine; C456276:cucurbit(7)uril; D047028:Macrocyclic Compounds; D000099246:Cucurbit(n)urils; C454805:cucurbit(6)uril
keywords: binding free energy; cucurbituril; ligand; molecular mechanics; statistical mechanics; steered molecular dynamics; umbrella sampling
references.citation: Phys Rev Lett. 2008 Dec 12;101(24):246105
references.pmid: 19113638
references.citation: Biophys J. 2006 Feb 1;90(3):864-77
references.pmid: 16284269
references.citation: J Chem Phys. 2009 Sep 28;131(12):124902
references.pmid: 19791916
references.citation: J Phys Chem B. 2008 Aug 21;112(33):10194-201
references.pmid: 18665626
references.citation: J Am Chem Soc. 2011 Mar 16;133(10):3570-81
references.pmid: 21341773
references.citation: J Chem Theory Comput. 2012 Mar 13;8(3):966-976
references.pmid: 22754402
references.citation: J Phys Chem B. 2008 May 15;112(19):5968-76
references.pmid: 18251532
references.citation: Biophys J. 2008 Oct;95(8):3575-90
references.pmid: 18621828
references.citation: Proteins. 2006 Nov 15;65(3):712-25
references.pmid: 16981200
references.citation: Proc Natl Acad Sci U S A. 2005 May 10;102(19):6825-30
references.pmid: 15867154
references.citation: Proc Natl Acad Sci U S A. 2001 Mar 27;98(7):3658-61
references.pmid: 11274384
references.citation: Biophys J. 1997 Mar;72(3):1047-69
references.pmid: 9138555
references.citation: Mol Inform. 2011 May 16;30(5):459-471
references.pmid: 21738559
references.citation: J Chem Theory Comput. 2008 Mar;4(3):435-47
references.pmid: 26620784
references.citation: J Comput Chem. 2004 Jul 15;25(9):1157-74
references.pmid: 15116359
references.citation: J Comput Aided Mol Des. 2012 May;26(5):475-87
references.pmid: 22366955
references.citation: J Am Chem Soc. 2010 Jun 2;132(21):7361-71
references.pmid: 20462212
references.citation: Electrophoresis. 1997 Dec;18(15):2714-23
references.pmid: 9504803
references.citation: J Comput Aided Mol Des. 2012 May;26(5):505-16
references.pmid: 22354755
references.citation: J Am Chem Soc. 2004 Oct 13;126(40):13156-64
references.pmid: 15469315
references.citation: Phys Rev Lett. 2007 Aug 10;99(6):068101
references.pmid: 17930869
references.citation: Carbohydr Res. 2004 Aug 2;339(11):1961-6
references.pmid: 15261589
references.citation: Science. 2002 Jun 7;296(5574):1832-5
references.pmid: 12052949
references.citation: J Phys Chem B. 2009 Jun 4;113(22):7836-43
references.pmid: 19425557
references.citation: J Phys Chem B. 2006 Mar 30;110(12):6372-8
references.pmid: 16553456
references.citation: J Comput Aided Mol Des. 2011 Aug;25(8):709-16
references.pmid: 21735261
references.citation: J Chem Theory Comput. 2010 Mar 9;6(3):637-646
references.pmid: 23794959
references.citation: J Chem Inf Comput Sci. 2003 Nov-Dec;43(6):1982-97
references.pmid: 14632449
references.citation: J Am Chem Soc. 2009 Mar 25;131(11):4012-21
references.pmid: 19133781
references.citation: Biophys J. 1999 Jan;76(1 Pt 1):188-97
references.pmid: 9876133
references.citation: J Chem Phys. 2009 Apr 28;130(16):164906
references.pmid: 19405629