This is an archive of papers published by the staff and faculty of Fox Chase Cancer Center. For questions about content, please contact Talbot Research Library
Last updated on
Xia J , Flynn WF , Gallicchio E , Zhang BW , He P , Tan Z , Levy RM
Large-scale asynchronous and distributed multidimensional replica exchange molecular simulations and efficiency analysis
J Comput Chem. 2015 Sep 5;36(23) :1772-85
PMID: 26149645 PMCID: PMC4512903
AbstractWe describe methods to perform replica exchange molecular dynamics (REMD) simulations asynchronously (ASyncRE). The methods are designed to facilitate large scale REMD simulations on grid computing networks consisting of heterogeneous and distributed computing environments as well as on homogeneous high-performance clusters. We have implemented these methods on NSF (National Science Foundation) XSEDE (Extreme Science and Engineering Discovery Environment) clusters and BOINC (Berkeley Open Infrastructure for Network Computing) distributed computing networks at Temple University and Brooklyn College at CUNY (the City University of New York). They are also being implemented on the IBM World Community Grid. To illustrate the methods, we have performed extensive (more than 60 ms in aggregate) simulations for the beta-cyclodextrin-heptanoate host-guest system in the context of one- and two-dimensional ASyncRE, and we used the results to estimate absolute binding free energies using the binding energy distribution analysis method. We propose ways to improve the efficiency of REMD simulations: these include increasing the number of exchanges attempted after a specified molecular dynamics (MD) period up to the fast exchange limit and/or adjusting the MD period to allow sufficient internal relaxation within each thermodynamic state. Although ASyncRE simulations generally require long MD periods (>picoseconds) per replica exchange cycle to minimize the overhead imposed by heterogeneous computing networks, we found that it is possible to reach an efficiency similar to conventional synchronous REMD, by optimizing the combination of the MD period and the number of exchanges attempted per cycle. (c) 2015 Wiley Periodicals, Inc.
NotesXia, Junchao Flynn, William F Gallicchio, Emilio Zhang, Bin W He, Peng Tan, Zhiqiang Levy, Ronald M P50 GM103368/GM/NIGMS NIH HHS/United States R01 GM030580/GM/NIGMS NIH HHS/United States United States J Comput Chem. 2015 Sep 5;36(23):1772-85. doi: 10.1002/jcc.23996. Epub 2015 Jul 7.