Which corners to cut? Guidelines on choosing optimal settings to maximise sampling with limited computational resources

05 August 2022, Version 1
This content is a preprint and has not undergone peer review at the time of posting.

Abstract

Despite the increasingly wide availability of computational resources, it is still challenging for researchers to perform comprehensive molecular dynamics (MD) simulations on an industrial scale. With ensemble approaches, great accuracy and precision are achievable at the cost of considerable computational efforts. There is a trade-off between accuracy, precision and the compute cost of performing ensembles. A question frequently raised is what is the most optimal way to perform MD based calculation of one or more properties of interest? Here, we use our findings from extensive ensemble molecular dynamics simulations of ligand-protein systems to underpin the recommendations we are making in cases where computational cost is an important consideration. We recommend performing a minimum of 5 replicas for the ESMACS-style protocol and 3 replicas (per λ window) for the TIES-like protocol, and the use of a stepwise procedure to reduce costs so as to minimise the loss of accuracy and precision of the results obtained.

Keywords

Ensemble simulation
Free energy calculation
Computational cost
Accuracy

Comments

Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.