Phillip S Hudson1, Justin K White1, Fiona L Kearns1, Milan Hodoscek2, Stefan Boresch3, H Lee Woodcock4. 1. Department of Chemistry, University of South Florida, 4202 E. Fowler Ave., CHE205, Tampa, FL 33620-5250, USA. 2. Center for Molecular Modeling, National Institute of Chemistry, Hajdrihova 19, SI-1000 Ljubljana, Slovenia. 3. Department of Computational Biological Chemistry, Faculty of Chemistry, University of Vienna, Währingerstraße 17, A-1090 Vienna, Austria. 4. Department of Chemistry, University of South Florida, 4202 E. Fowler Ave., CHE205, Tampa, FL 33620-5250, USA. Electronic address: hlw@usf.edu.
Abstract
BACKGROUND: Accurately modeling condensed phase processes is one of computation's most difficult challenges. Include the possibility that conformational dynamics may be coupled to chemical reactions, where multiscale (i.e., QM/MM) methods are needed, and this task becomes even more daunting. METHODS: Free energy simulations (i.e., molecular dynamics), multiscale modeling, and reweighting schemes. RESULTS: Herein, we present two new approaches for mitigating the aforementioned challenges. The first is a new chain-of-replica method (off-path simulations, OPS) for computing potentials of mean force (PMFs) along an easily defined reaction coordinate. This development is coupled with a new distributed, highly-parallel replica framework (REPDstr) within the CHARMM package. Validation of these new schemes is carried out on two processes that undergo conformational changes. First is the simple torsional rotation of butane, while a much more challenging glycosidic rotation (in vacuo and solvated) is the second. Additionally, a new approach that greatly improves (i.e., possibly an order of magnitude) the efficiency of computing QM/MM PMFs is introduced and compared to standard schemes. Our efforts are grounded in the recently developed method for efficiently computing QM-based free energies (i.e., QM-Non-Boltzmann Bennett, QM-NBB). Again, we validate this new technique by computing the QM/MM PMF of butane's torsional rotation. CONCLUSIONS: The OPS-REPDstr method is a promising new approach that overcomes many limitations of standard pathway simulations in CHARMM. The combination of QM-NBB with pathway techniques is very promising as it offers significant advantages over current procedures. GENERAL SIGNIFICANCE: Efficiently computing potentials of mean force is a major, unresolved, area of interest. This article is part of a Special Issue entitled Recent developments of molecular dynamics.
BACKGROUND: Accurately modeling condensed phase processes is one of computation's most difficult challenges. Include the possibility that conformational dynamics may be coupled to chemical reactions, where multiscale (i.e., QM/MM) methods are needed, and this task becomes even more daunting. METHODS: Free energy simulations (i.e., molecular dynamics), multiscale modeling, and reweighting schemes. RESULTS: Herein, we present two new approaches for mitigating the aforementioned challenges. The first is a new chain-of-replica method (off-path simulations, OPS) for computing potentials of mean force (PMFs) along an easily defined reaction coordinate. This development is coupled with a new distributed, highly-parallel replica framework (REPDstr) within the CHARMM package. Validation of these new schemes is carried out on two processes that undergo conformational changes. First is the simple torsional rotation of butane, while a much more challenging glycosidic rotation (in vacuo and solvated) is the second. Additionally, a new approach that greatly improves (i.e., possibly an order of magnitude) the efficiency of computing QM/MM PMFs is introduced and compared to standard schemes. Our efforts are grounded in the recently developed method for efficiently computing QM-based free energies (i.e., QM-Non-Boltzmann Bennett, QM-NBB). Again, we validate this new technique by computing the QM/MM PMF of butane's torsional rotation. CONCLUSIONS: The OPS-REPDstr method is a promising new approach that overcomes many limitations of standard pathway simulations in CHARMM. The combination of QM-NBB with pathway techniques is very promising as it offers significant advantages over current procedures. GENERAL SIGNIFICANCE: Efficiently computing potentials of mean force is a major, unresolved, area of interest. This article is part of a Special Issue entitled Recent developments of molecular dynamics.
Authors: Gerhard König; Frank C Pickard; Jing Huang; Andrew C Simmonett; Florentina Tofoleanu; Juyong Lee; Pavlo O Dral; Samarjeet Prasad; Michael Jones; Yihan Shao; Walter Thiel; Bernard R Brooks Journal: J Comput Aided Mol Des Date: 2016-08-30 Impact factor: 3.686
Authors: Frank C Pickard; Gerhard König; Florentina Tofoleanu; Juyong Lee; Andrew C Simmonett; Yihan Shao; Jay W Ponder; Bernard R Brooks Journal: J Comput Aided Mol Des Date: 2016-09-19 Impact factor: 3.686
Authors: David K Wolfe; Joseph R Persichetti; Ajeet K Sharma; Phillip S Hudson; H Lee Woodcock; Edward P O'Brien Journal: J Chem Theory Comput Date: 2020-02-17 Impact factor: 6.006
Authors: Frank C Pickard; Gerhard König; Andrew C Simmonett; Yihan Shao; Bernard R Brooks Journal: Bioorg Med Chem Date: 2016-08-22 Impact factor: 3.641