Abstract
Reliability of force fields is an essential aspect of protein-folding simulation. In this work, we introduced a newly developed on-the-fly charge-updating scheme called the polarized structure-specific backbone charge (PSBC) model. The PSBC model was designed with the purpose of building the polarizability of backbone hydrogen bonds into the force field by updating the partial charges of backbone hydrogen-bond donor and acceptor atoms during folding simulation. This implementation was intended to mimic the heterogeneity of the protein surrounding during folding. Multiple single-trajectory molecular dynamics simulations were performed to fold a polyalanine peptide, namely ER (Ac-A(EAAAR)(3)A-NH(2)), using both polarizable (PSBC) and nonpolarizable (Amber03) force fields. Through the PSBC model, ER was folded into a helical peptide with helix content that agrees well with experiments. Comparison between simulations performed using the aforementioned force fields demonstrably showed the importance of electrostatic polarization effect in the folding of the short α-helical peptide. The PSBC model was further validated by folding two other short peptides with different helicities.