title: Accuracy and convergence of free energy differences calculated from nonequilibrium switching processes.
abstract: The molecular-dynamics-based calculation of accurate free energy differences for biomolecular systems is a challenging task. Accordingly, convergence and accuracy of established equilibrium methods has been subject of many studies, often focusing at small test systems. In contrast, the potential of more recently proposed nonequilibrium methods, derived from the Jarzynski and Crooks equalities, has not yet fully been explored. Here, we compare the performance of these methods by calculating free energy differences for test systems at different levels of complexity and varying extent of the involved perturbations. We consider the interconversion of ethane into methanol, the switching of a tryptophane-sidechain in a tripeptide, and the binding of two different ligands to the globular protein snurportin 1. On the basis of our results, we suggest and assess a new nonequilibrium free energy method, Crooks Gaussian Intersection (CGI), which combines the advantages of existing methods. CGI is highly parallelizable and, for the test systems considered here, is shown to outperform the other studied equilibrium and nonequilibrium methods.
authors.affiliation: Theoretical and Computational Biophysics Department, Max-Planck-Institute for Biophysical Chemistry, Am Fassberg 11, 37077 Göttingen, Germany.
authors.forename: Maik
authors.initials: M
authors.lastname: Goette
authors.forename: Helmut
authors.initials: H
authors.lastname: Grubmüller
mesh_terms: D003198:Computer Simulation; D004980:Ethane / Q000737:chemistry; D005998:Glycine / Q000737:chemistry; D006801:Humans; D008024:Ligands; D000432:Methanol / Q000737:chemistry; D008956:Models, Chemical; D009842:Oligopeptides / Q000737:chemistry; D039381:RNA Cap-Binding Proteins / Q000737:chemistry; D018160:Receptors, Cytoplasmic and Nuclear / Q000737:chemistry; D013816:Thermodynamics; D014364:Tryptophan / Q000737:chemistry
chemical_list: D008024:Ligands; D009842:Oligopeptides; D039381:RNA Cap-Binding Proteins; D018160:Receptors, Cytoplasmic and Nuclear; C114090:SNUPN protein, human; D014364:Tryptophan; D004980:Ethane; D005998:Glycine; D000432:Methanol