Implementation of a protein reduced point charge model toward molecular dynamics applications.

A reduced point charge model was developed in a previous work from the study of extrema in smoothed charge density distribution functions generated from the Amber99 molecular electrostatic potential. In the present work, such a point charge distribution is coupled with the Amber99 force field and implemented in the program TINKER to allow molecular dynamics (MD) simulations of proteins. First applications to two polypeptides that involve α-helix and β-sheet motifs are analyzed and compared to all-atom MD simulations. Two types of coarse-grained (CG)-based trajectories are generated using, on one hand, harmonic bond stretching terms and, on the other hand, distance restraints. Results show that the use of the unrestrained CG conditions are sufficient to preserve most of the secondary structure characteristics but restraints lead to a better agreement between CG and all-atom simulation results such as rmsd, dipole moment, and time-dependent mean square deviation functions.