Consistent treatment of inter- and intramolecular polarization in molecular mechanics calculations

A protocol is described for the treatment of molecular polarization in force field calculations. The resulting model is consistent in that both inter- and intramolecular polarization are handled within a single scheme. An analytical formula for removing intramolecular polarization from a set of atomic multipoles for an arbitrary static structure or conformation is given. With the help of the intramolecular polarization, these permanent atomic multipoles can then be applied in modeling alternative conformations of a molecule. Equipped with this simple technique, one can derive transferable electrostatic parameters for peptides and proteins using flexible model compounds such as dipeptides. The proposed procedure is tested for its ability to describe the electrostatic potential around various configurations of the N-methylacetamide dimer. The effect of different intramolecular polarization schemes on the accuracy of a force field model of the electrostatic potential of alanine dipeptide is investigated. A group-based scheme for including direct intramolecular polarization is shown to be most successful in accounting for the conformational dependence of electrostatic potentials.