Improvements to the confinement method for the calculation of conformational free energy differences are presented. By taking advantage of phase space overlap between simulations at different frequencies, significant gains in accuracy and speed are reached. The optimal frequency spacing for the simulations is obtained from extrapolations of the confinement energy, and relaxation time analysis is used to determine time steps, simulation lengths, and friction coefficients. At postprocessing, interpolation of confinement energies is used to significantly reduce discretization errors in the calculation of conformational free energies. The efficiency of this protocol is illustrated by applications to alanine n-peptides and lactoferricin. For the alanine-n-peptide, errors were reduced between 2- and 10-fold and sampling times between 8- and 67-fold, while for lactoferricin the long sampling times at low frequencies were reduced 10-100-fold.

Balancing Accuracy and Cost of Confinement Simulations by Interpolation and Extrapolation of Confinement Energies / F. Villemot, R. Capelli, G. Colombo, A. van der Vaart. - In: JOURNAL OF CHEMICAL THEORY AND COMPUTATION. - ISSN 1549-9618. - 12:6(2016 Jun), pp. 2779-2789. [10.1021/acs.jctc.5b01183]

Balancing Accuracy and Cost of Confinement Simulations by Interpolation and Extrapolation of Confinement Energies

R. Capelli
Secondo
;
2016

Abstract

Improvements to the confinement method for the calculation of conformational free energy differences are presented. By taking advantage of phase space overlap between simulations at different frequencies, significant gains in accuracy and speed are reached. The optimal frequency spacing for the simulations is obtained from extrapolations of the confinement energy, and relaxation time analysis is used to determine time steps, simulation lengths, and friction coefficients. At postprocessing, interpolation of confinement energies is used to significantly reduce discretization errors in the calculation of conformational free energies. The efficiency of this protocol is illustrated by applications to alanine n-peptides and lactoferricin. For the alanine-n-peptide, errors were reduced between 2- and 10-fold and sampling times between 8- and 67-fold, while for lactoferricin the long sampling times at low frequencies were reduced 10-100-fold.
free energy calculation; confinement method; molecular dynamics
Settore FIS/03 - Fisica della Materia
giu-2016
6-mag-2016
Article (author)
File in questo prodotto:
File Dimensione Formato  
s1-ln20681062012724402-1939656818Hwf1590306036IdV11907744622068106PDF_HI0001.pdf

Open Access dal 01/09/2017

Tipologia: Post-print, accepted manuscript ecc. (versione accettata dall'editore)
Dimensione 4.17 MB
Formato Adobe PDF
4.17 MB Adobe PDF Visualizza/Apri
Pubblicazioni consigliate

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/2434/399657
Citazioni
  • ???jsp.display-item.citation.pmc??? 1
  • Scopus 5
  • ???jsp.display-item.citation.isi??? 5
social impact