Influence of various parameters in the replica-exchange molecular dynamics method: Number of replicas, replica-exchange frequency, and thermostat coupling time constant

Abstract

The replica-exchange molecular dynamics (REMD) method has been used for conformational sampling of various biomolecular systems. To maximize sampling efficiency, some adjustable parameters must be optimized. Although it is agreed that shorter intervals between the replica-exchange attempts enhance traversals in the temperature space, details regarding the artifacts caused by these short intervals are controversial. In this study, we revisit this problem by performing REMD simulations on an alanine octapeptide in an implicit solvent. Fifty different sets of conditions, which are a combination of five replica-exchange periods, five different numbers of replicas, and two thermostat coupling time constants, were investigated. As a result, although short replica-exchange intervals enhanced the traversals in the temperature space, they led to artifacts in the ensemble average of the temperature, potential energy, and helix content. With extremely short replica-exchange intervals, i.e., attempted at every time step, the ensemble average of the temperature deviated from the thermostat temperature by ca. 7 K. Differences in the ensembles were observed even for larger replica-exchange intervals (between 100 and 1,000 steps). In addition, the shorter thermostat coupling time constant reduced the artifacts found when short replica-exchange intervals were used, implying that these artifacts are caused by insufficient thermal relaxation between the replica-exchange events. Our results will be useful to reduce the artifacts found in REMD simulations by adjusting some key parameters.

When MELD Meets GaMD: Accelerating Biomolecular Landscape Exploration

“…For example, the MELD Best1Pop of 3GB1 with seven replicas and 1FME with five is substantially lower than the previous results with 30 replicas. This is in agreement with previous results showing how the number of replicas could impact differently depending on the nature of the process. − …”

myPresto/omegagene 2020: a molecular dynamics simulation engine for virtual-system coupled sampling

“…These approaches can efficiently sample conformational ensembles of molecular systems. Simultaneously, generalized ensemble methods require users to acquire experience and knowledge for performing simulations appropriately [14][15][16][17][18][19][20].…”

Studies on Molecular Dynamics of Intrinsically Disordered Proteins and Their Fuzzy Complexes: A Mini-Review

“…Although the force fields have been extensively developed, the current classical ones still fail to reproduce some of the details of molecular behavior [35,103,104]. Artifacts from various conditions including the PBC [40] and conditions for enhanced sampling [55] are not fully understood for IDR simulations. Because users must acquire some experience and technical knowledge to appropriately apply some kinds of simulation methods (e.g., enhanced sampling methods and coarse-grained simulations), development and maintenance of easy-to-use software are essential.…”

“…In general, these enhanced sampling methods involve some technical issues such as fine tuning of adjustable parameters. For example, the resultant ensemble of a REMD simulation depends on the setting of the total length of trajectory, the number of replicas, and interval times for replica-exchange trials [[50], [51], [52], [53], [54], [55]]. Because there is no gold standard for adjusting these parameters, techniques and experiences to setup the simulation conditions are required for users.…”