Enhancing torsional sampling using fully adaptive simulated tempering
Enhancing torsional sampling using fully adaptive simulated tempering
Enhanced sampling algorithms are indispensable when working with highly disconnected multimodal distributions. An important application of these is the conformational exploration of particular internal degrees of freedom of molecular systems. However, despite the existence of many commonly used enhanced sampling algorithms to explore these internal motions, they often rely on system-dependent parameters, which negatively impact efficiency and reproducibility. Here, we present fully adaptive simulated tempering (FAST), a variation of the irreversible simulated tempering algorithm, which continuously optimizes the number, parameters, and weights of intermediate distributions to achieve maximally fast traversal over a space defined by the change in a predefined thermodynamic control variable such as temperature or an alchemical smoothing parameter. This work builds on a number of previously published methods, such as sequential Monte Carlo, and introduces a novel parameter optimization procedure that can, in principle, be used in any expanded ensemble algorithms. This method is validated by being applied on a number of different molecular systems with high torsional kinetic barriers. We also consider two different soft-core potentials during the interpolation procedure and compare their performance. We conclude that FAST is a highly efficient algorithm, which improves simulation reproducibility and can be successfully used in a variety of settings with the same initial hyperparameters.
Adaptively biased molecular dynamics,, Interpolation, Lysozyme, Monte Carlo methods, Proteins, Biomolecular dynamics, Replica-exchange molecular dynamics simulation, machine learning, molecular dynamics, optimization algorithms
Suruzhon, Miroslav
4ea4dd8b-0a98-4598-9eaa-756d943b5dca
Abdel-Maksoud, Khaled
836e55a9-b518-44b1-b423-deab8a0e69bb
Bodnarchuk, Michael S.
cb7c3390-a1e3-4e13-916c-200706d11f34
Ciancetta, Antonella
992e4150-4131-4adc-9e00-d7d415583ea5
Wall, Ian D.
059fa7bb-bde0-457f-8ba3-7e31cf72d3f1
Essex, Jonathan W.
1f409cfe-6ba4-42e2-a0ab-a931826314b5
19 April 2024
Suruzhon, Miroslav
4ea4dd8b-0a98-4598-9eaa-756d943b5dca
Abdel-Maksoud, Khaled
836e55a9-b518-44b1-b423-deab8a0e69bb
Bodnarchuk, Michael S.
cb7c3390-a1e3-4e13-916c-200706d11f34
Ciancetta, Antonella
992e4150-4131-4adc-9e00-d7d415583ea5
Wall, Ian D.
059fa7bb-bde0-457f-8ba3-7e31cf72d3f1
Essex, Jonathan W.
1f409cfe-6ba4-42e2-a0ab-a931826314b5
Suruzhon, Miroslav, Abdel-Maksoud, Khaled, Bodnarchuk, Michael S., Ciancetta, Antonella, Wall, Ian D. and Essex, Jonathan W.
(2024)
Enhancing torsional sampling using fully adaptive simulated tempering.
The Journal of Chemical Physics, 160 (15), [154110].
(doi:10.1063/5.0190659).
Abstract
Enhanced sampling algorithms are indispensable when working with highly disconnected multimodal distributions. An important application of these is the conformational exploration of particular internal degrees of freedom of molecular systems. However, despite the existence of many commonly used enhanced sampling algorithms to explore these internal motions, they often rely on system-dependent parameters, which negatively impact efficiency and reproducibility. Here, we present fully adaptive simulated tempering (FAST), a variation of the irreversible simulated tempering algorithm, which continuously optimizes the number, parameters, and weights of intermediate distributions to achieve maximally fast traversal over a space defined by the change in a predefined thermodynamic control variable such as temperature or an alchemical smoothing parameter. This work builds on a number of previously published methods, such as sequential Monte Carlo, and introduces a novel parameter optimization procedure that can, in principle, be used in any expanded ensemble algorithms. This method is validated by being applied on a number of different molecular systems with high torsional kinetic barriers. We also consider two different soft-core potentials during the interpolation procedure and compare their performance. We conclude that FAST is a highly efficient algorithm, which improves simulation reproducibility and can be successfully used in a variety of settings with the same initial hyperparameters.
Text
154110_1_5.0190659
- Version of Record
More information
Submitted date: 7 December 2023
Accepted/In Press date: 23 March 2024
Published date: 19 April 2024
Keywords:
Adaptively biased molecular dynamics,, Interpolation, Lysozyme, Monte Carlo methods, Proteins, Biomolecular dynamics, Replica-exchange molecular dynamics simulation, machine learning, molecular dynamics, optimization algorithms
Identifiers
Local EPrints ID: 491012
URI: http://eprints.soton.ac.uk/id/eprint/491012
ISSN: 0021-9606
PURE UUID: 6540378f-894a-47d4-8ce5-02dfceb47bb3
Catalogue record
Date deposited: 11 Jun 2024 16:37
Last modified: 12 Jun 2024 02:03
Export record
Altmetrics
Contributors
Author:
Miroslav Suruzhon
Author:
Khaled Abdel-Maksoud
Author:
Michael S. Bodnarchuk
Author:
Antonella Ciancetta
Author:
Ian D. Wall
Download statistics
Downloads from ePrints over the past year. Other digital versions may also be available to download e.g. from the publisher's website.
View more statistics