Brian K. Radak, Christophe Chipot, Donghyuk Suh, Sunhwan Jo, Wei Jiang,
James C. Phillips, Klaus Schulten, and Benoît Roux.
Constant-pH molecular dynamics simulations for large biomolecular
systems.
Journal of Chemical Theory and Computation, 13:5933-5944,
2017.
(PMC: PMC5726918)
RADA2017
An increasingly important endeavor is to develop computational strategies
that enable molecular dynamics (MD) simulations of biomolecular systems
with spontaneous changes in protonation states under conditions of
constant pH. The present work describes our efforts to implement the
powerful constant-pH MD simulation method, based on a hybrid
nonequilibrium MD/Monte Carlo (neMD/MC) technique within the highly
scalable program NAMD. The constant-pH hybrid neMD/MC method has
several appealing features; it samples the correct semigrand canonical
ensemble rigorously, the computational cost increases linearly with the
number of titratable sites, and it is applicable to explicit solvent
simulations. The present implementation of the constant-pH hybrid
neMD/MC in NAMD is designed to handle a wide range of biomolecular
systems with no constraints on the choice of force field. Furthermore, the
sampling efficiency can be adaptively improved on-the-fly by adjusting
algorithmic parameters during the simulation. Illustrative examples
emphasizing medium- and large-scale applications on next-generation
supercomputing architectures are provided.
Download Full Text
The manuscripts available on our site are provided for your personal
use only and may not be retransmitted or redistributed without written
permissions from the paper's publisher and author. You may not upload any
of this site's material to any public server, on-line service, network, or
bulletin board without prior written permission from the publisher and
author. You may not make copies for any commercial purpose. Reproduction
or storage of materials retrieved from this web site is subject to the
U.S. Copyright Act of 1976, Title 17 U.S.C.