Choice of force fields and water models for sampling solution conformations of bacteriophage T4 lysozyme

Abstract

A protein may exist as an ensemble of different conformations in solution, which cannot be represented by a single static structure. Molecular dynamics (MD) simulation has become a useful tool for sampling protein conformations in solution, but force fields and water models are important issues. This work presents a case study of the bacteriophage T4 lysozyme (T4L). We have found that MD simulations using a classic AMBER99SB force field and TIP4P water model cannot well describe hinge-bending domain motion of the wild-type T4L at the timescale of one microsecond. Other combinations, such as a residue-specific force field called RSFF2+ and a dispersion-corrected water model TIP4P-D, are able to sample reasonable solution conformations of T4L, which are in good agreement with experimental data. This primary study may provide candidates of force fields and water models for further investigating conformational transition of T4L.