Journal
JOURNAL OF COMPUTATIONAL CHEMISTRY
Volume 44, Issue 20, Pages 1740-1749Publisher
WILEY
DOI: 10.1002/jcc.27124
Keywords
biomolecular simulations; enhanced conformational sampling; free energy calculation; high performance computing; MBAR; molecular dynamics; parallel computing; replica-exchange MD
Categories
Ask authors/readers for more resources
gREST is an enhanced sampling algorithm used for proteins or other systems with rugged energy landscapes. It differs from the REMD method as the solvent temperatures are the same in all replicas, while solute temperatures are different and frequently exchanged between replicas to explore various solute structures. The gREST scheme is applied to large biological systems using a large number of processors, reducing communication time and performing energy evaluations on-the-fly during simulations. These advanced schemes in the GENESIS software provide new possibilities for studying large biomolecular complex systems with slow conformational dynamics.
Generalized replica exchange with solute tempering (gREST) is one of the enhanced sampling algorithms for proteins or other systems with rugged energy landscapes. Unlike the replica-exchange molecular dynamics (REMD) method, solvent temperatures are the same in all replicas, while solute temperatures are different and are exchanged frequently between replicas for exploring various solute structures. Here, we apply the gREST scheme to large biological systems containing over one million atoms using a large number of processors in a supercomputer. First, communication time on a multi-dimensional torus network is reduced by matching each replica to MPI processors optimally. This is applicable not only to gREST but also to other multi-copy algorithms. Second, energy evaluations, which are necessary for the multistate bennet acceptance ratio (MBAR) method for free energy estimations, are performed on-the-fly during the gREST simulations. Using these two advanced schemes, we observed 57.72 ns/day performance in 128-replica gREST calculations with 1.5 million atoms system using 16,384 nodes in Fugaku. These schemes implemented in the latest version of GENESIS software could open new possibilities to answer unresolved questions on large biomolecular complex systems with slow conformational dynamics.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available