A key
advantage of polarizable force fields is their ability to
model the atomic polarization effects that play key roles in the atomic
many-body interactions. In this work, we assessed the accuracy of
the recently developed polarizable Gaussian Multipole (pGM) models
in reproducing quantum mechanical (QM) interaction energies, many-body
interaction energies, as well as the nonadditive and additive contributions
to the many-body interactions for peptide main-chain hydrogen-bonding
conformers, using glycine dipeptide oligomers as the model systems.
Two types of pGM models were considered, including that with (pGM-perm)
and without (pGM-ind) permanent atomic dipoles. The performances of
the pGM models were compared with several widely used force fields,
including two polarizable (Amoeba13 and ff12pol) and three additive
(ff19SB, ff15ipq, and ff03) force fields. Encouragingly, the pGM models
outperform all other force fields in terms of reproducing QM interaction
energies, many-body interaction energies, as well as the nonadditive
and additive contributions to the many-body interactions, as measured
by the root-mean-square errors (RMSEs) and mean absolute errors (MAEs).
Furthermore, we tested the robustness of the pGM models against polarizability
parameterization errors by employing alternative polarizabilities
that are either scaled or obtained from other force fields. The results
show that the pGM models with alternative polarizabilities exhibit
improved accuracy in reproducing QM many-body interaction energies
as well as the nonadditive and additive contributions compared with
other polarizable force fields, suggesting that the pGM models are
robust against the errors in polarizability parameterizations. This
work shows that the pGM models are capable of accurately modeling
polarization effects and have the potential to serve as templates
for developing next-generation polarizable force fields for modeling
various biological systems.