The accelerated weight histogram method for alchemical free energy calculations

被引:23
作者
Lundborg, M. [1 ]
Lidmar, J. [2 ]
Hess, B. [3 ]
机构
[1] ERCO Pharma AB, S-11439 Stockholm, Sweden
[2] KTH Royal Inst Technol, Dept Phys, S-10691 Stockholm, Sweden
[3] KTH Royal Inst Technol, Sci Life Lab, Dept Appl Phys, S-10691 Stockholm, Sweden
基金
欧盟地平线“2020”;
关键词
MOLECULAR-DYNAMICS METHOD; CHARMM FORCE-FIELD; SIMULATION; EFFICIENT; WATER; THERMODYNAMICS; MECHANICS; ALGORITHM;
D O I
10.1063/5.0044352
中图分类号
O64 [物理化学(理论化学)、化学物理学];
学科分类号
070304 ; 081704 ;
摘要
The accelerated weight histogram method is an enhanced sampling technique used to explore free energy landscapes by applying an adaptive bias. The method is general and easy to extend. Herein, we show how it can be used to efficiently sample alchemical transformations, commonly used for, e.g., solvation and binding free energy calculations. We present calculations and convergence of the hydration free energy of testosterone, representing drug-like molecules. We also include methane and ethanol to validate the results. The protocol is easy to use, does not require a careful choice of parameters, and scales well to accessible resources, and the results converge at least as quickly as when using conventional methods. One benefit of the method is that it can easily be combined with other reaction coordinates, such as intermolecular distances. (c) 2021 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
引用
收藏
页数:14
相关论文
共 66 条
[51]   THERMODYNAMIC LENGTH AND DISSIPATED AVAILABILITY [J].
SALAMON, P ;
BERRY, RS .
PHYSICAL REVIEW LETTERS, 1983, 51 (13) :1127-1130
[52]   Minimizing thermodynamic length to select intermediate states for free-energy calculations and replica-exchange simulations [J].
Shenfeld, Daniel K. ;
Xu, Huafeng ;
Eastwood, Michael P. ;
Dror, Ron O. ;
Shaw, David E. .
PHYSICAL REVIEW E, 2009, 80 (04)
[53]  
Shirts Michael R, 2013, Methods Mol Biol, V924, P271, DOI 10.1007/978-1-62703-017-5_11
[54]  
Shirts MR, 2007, ANN REP COMP CHEM, V3, P41, DOI 10.1016/S1574-1400(07)03004-6
[55]   Statistically optimal analysis of samples from multiple equilibrium states [J].
Shirts, Michael R. ;
Chodera, John D. .
JOURNAL OF CHEMICAL PHYSICS, 2008, 129 (12)
[56]   Identification of Novel Inhibitors of Mycobacterium tuberculosis PknG Using Pharmacophore Based Virtual Screening, Docking, Molecular Dynamics Simulation, and Their Biological Evaluation [J].
Singh, Nidhi ;
Tiwari, Sameer ;
Srivastava, Kishore K. ;
Siddiqi, Mohammad Imran .
JOURNAL OF CHEMICAL INFORMATION AND MODELING, 2015, 55 (06) :1120-1129
[57]   Thermodynamic Metrics and Optimal Paths [J].
Sivak, David A. ;
Crooks, Gavin E. .
PHYSICAL REVIEW LETTERS, 2012, 108 (19)
[58]  
sklearn Supervised Learning, 2020, PYTHON IMPLEMENTATIO
[59]   Replica-exchange molecular dynamics method for protein folding [J].
Sugita, Y ;
Okamoto, Y .
CHEMICAL PHYSICS LETTERS, 1999, 314 (1-2) :141-151
[60]   A COMPUTER-SIMULATION METHOD FOR THE CALCULATION OF EQUILIBRIUM-CONSTANTS FOR THE FORMATION OF PHYSICAL CLUSTERS OF MOLECULES - APPLICATION TO SMALL WATER CLUSTERS [J].
SWOPE, WC ;
ANDERSEN, HC ;
BERENS, PH ;
WILSON, KR .
JOURNAL OF CHEMICAL PHYSICS, 1982, 76 (01) :637-649