CGI – Accurate Free Energy Differences from Non-Equilibrium Simulations
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 on small test systems. In contrast, the potential of more recently proposed nonequilibrium methods (based on Fast Growth Thermodynamic Integration, FGTI), derived from the Jarzynski and Crooks equalities, has not yet fully been explored.
We compare the performance of these methods by calculating free energy differences for test systems at different levels of complexity and varying the extent of the involved perturbations. We consider the interconversion of ethane into methanol, the switching of a tryptophane side chain into 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.