Computing conformational free energy differences in explicit solvent: An efficient thermodynamic cycle using an auxiliary potential and a free energy functional constructed from the end points. Issue 15 (23rd December 2016)

• Record Type: Journal Article
• Authors: Harris, Robert C.; Deng, Nanjie; Levy, Ronald M.; Ishizuka, Ryosuke; Matubayasi, Nobuyuki
• Other Names: Hirst Jonathan guestEditor.; Im Wonpil guestEditor.; Shea Joan‐Emma guestEditor.
• Is Part Of: Journal of computational chemistry. Volume 38:Issue 15(2017)
• Page Start: 1198
• Page End: 1208
• Publication Date: 2016-12-23
• Subjects: energy representation -- molecular dynamics simulations -- distribution function -- conformational changes -- solvation free energy; Chemistry -- Data processing -- Periodicals; 542.85
• DOI: 10.1002/jcc.24668
• Languages: English
• ISSNs: 0192-8651

Abstract

Abstract : Many biomolecules undergo conformational changes associated with allostery or ligand binding. Observing these changes in computer simulations is difficult if their timescales are long. These calculations can be accelerated by observing the transition on an auxiliary free energy surface with a simpler Hamiltonian and connecting this free energy surface to the target free energy surface with free energy calculations. Here, we show that the free energy legs of the cycle can be replaced with energy representation (ER) density functional approximations. We compute: (1) The conformational free energy changes for alanine dipeptide transitioning from the right‐handed free energy basin to the left‐handed basin and (2) the free energy difference between the open and closed conformations of β ‐cyclodextrin, a "host" molecule that serves as a model for molecular recognition in host‐guest binding. β ‐cyclodextrin contains 147 atoms compared to 22 atoms for alanine dipeptide, making β ‐cyclodextrin a large molecule for which to compute solvation free energies by free energy perturbation or integration methods and the largest system for which the ER method has been compared to exact free energy methods. The ER method replaced the 28 simulations to compute each coupling free energy with two endpoint simulations, reducing the computational time for the alanine dipeptide calculation by about 70% and for the β ‐cyclodextrin by > 95%. The method works even when the distribution of conformations on the auxiliary free energy surface differs substantially from that on the target free energy surface, although some degree of overlap between the two surfaces is required. © 2016 Wiley Periodicals, Inc. Abstract : The free energy differences associated with conformational changes are difficult to compute in explicit solvent. Instead, these free energy differences can be computed on an auxiliary free energy surface and the desired free energy difference obtained by adding the free energies of transferring the end states from the auxiliary surface to the target surface. Here, we show that computing these transfer free energies with the energy representation method substantially reduces the cost of these calculations.