Novel procedure for thermal equilibration in molecular dynamics simulation

Author: Gallo Marco   Grant Barry   Teodoro Miguel   Melton Julia   Cieplak Piotr   Phillips George   Stec Boguslaw  

Publisher: Taylor & Francis Ltd

ISSN: 0892-7022

Source: Molecular Simulation, Vol.35, Iss.5, 2009-04, pp. : 349-357

Disclaimer: Any content in publications that violate the sovereignty, the constitution or regulations of the PRC is not accepted or approved by CNPIEC.

Previous Menu Next

Abstract

We describe a simple novel procedure for achieving thermal equilibration between a protein and a surrounding solvent during molecular dynamics (MD) simulation. The method uniquely defines the length of simulation time required to achieve thermal equilibrium over a broad range of parameters, thus removing ambiguities associated with the traditional heuristic approaches. The proposed protocol saves simulation time and avoids bias introduced by the inclusion of non-equilibrium events. The key element of the procedure involves coupling only the solvent atoms to a standard heat bath. Measuring progress towards thermal equilibration involves simply monitoring the difference in temperature between the solvent and the protein. Here, we report that the results of MD simulations using the above procedure are measurably improved relative to the traditional approaches in terms of root-mean-square deviations and principal components analysis both indicating significantly less undesirable divergence.