Multiple Topology Replica Exchange of Expanded Ensembles (MT-REXEE) for Multidimensional Alchemical Calculations
arXiv:2408.11038 · doi:10.1021/acs.jctc.4c01268
Abstract
Relative free energy calculations are now widely used in academia and industry, but the accuracy is often limited by poor sampling of the complexes conformational ensemble. To address this, we have developed a novel method termed Multi-Topology Replica Exchange of Expanded Ensembles (MT-REXEE). This method enables parallel expanded ensemble calculations, facilitating iterative relative free energy computations while allowing conformational exchange between parallel transformations. These iterative transformations are adaptable to any set of systems with a common backbone or central substructure. We demonstrate that the MT-REXEE method maintains thermodynamic cycle closure to the same extent as standard expanded ensemble for both solvation free energy and relative binding free energy. The transformations tested involve simple systems that incorporate diverse heavy atoms and multi-site perturbations of a small molecule core resembling multi-site dynamics, without necessitating modifications to the MD code, which in our initial implementation is GROMACS. We outline a systematic approach for topology set-up and provide instructions on how to perform inter-replicate coordinate modifications. This work shows that MT-REEXE can be used to perform accurate and reproducible free energy estimates and prompts expansion to more complex test systems and other molecular dynamics simulation infrastructures.
References in corpus (10)
- Canonical sampling through velocity-rescaling
- Statistically optimal analysis of samples from multiple equilibrium states
- Pressure control using stochastic cell rescaling
- Enhanced sampling methods for molecular dynamics simulations
- Best Practices for Alchemical Free Energy Calculations
- Replica exchange and expanded ensemble simulations as Gibbs sampling: Simple improvements for enhanced mixing
- A Unified Approach to Enhanced Sampling
- Wang-Landau Algorithm: a Theoretical Analysis of the Saturation of the Error
- Adding alchemical variables to metadynamics to enhance sampling in free energy calculations
- Lambda-ABF: Simplified, Portable, Accurate and Cost-effective Alchemical Free Energy Computations