The two-site coarse-grained model for the interactions of charged side chains, to be used with our coarse-grained UNRES force field for protein simulations proposed in the accompanying paper, has been extended to pairs of oppositely charged side chains. The potentials of mean force of four pairs of molecules modeling charged amino-acid side chains, i.e., propionate-n-pentylamine cation (for aspartic acid-lysine), butyrate-n-pentylamine cation (for glutamic acid-lysine), propionate-1-butylguanidine (for aspartic acid-arginine), and butyrate-1-butylguanidine (for glutamic acid-arginine) pairs were determined by umbrella-sampling molecular dynamics simulations in explicit water as functions of distance and orientation, and the analytical expression was fitted to the potentials of mean force. Compared to pairs of like-charged side chains discussed in the accompanying paper, an average quadrupole-quadrupole interaction term had to be introduced to reproduce the Coulombic interactions, and a multistate model of charge distribution had to be introduced to fit the potentials of mean force of all oppositely charged pairs well. The model reproduces all salt-bridge minima and, consequently, is likely to improve the performance of the UNRES force field.
Read full abstract