Study of charges transferability for use in force fields

Optimization of a molecular mechanics force field for type-II polyoxometalates focussing on electrostatic interactions: a case study

In this study, we have focussed on type-II polyanions such as [M(7)O(24)](6-), and we have developed and validated optimized force fields that include electrostatic and van der Waals interactions. These contributions to the total steric energy are described by the nonbonded term, which encompasses all interactions between atoms that are not transmitted through the bonds. A first validation of a stochastic technique based on genetic algorithms was previously made for the optimization of force fields dedicated to type-I polyoxometalates. To describe the new nonbonded term added in the functional, a fixed-charged model was chosen. Therefore, one of the main issues was to analyze that which partial atomic charges could be reliably used to describe these interactions in such inorganic compounds. Based on several computational strategies, molecular mechanics (MM) force field parameters were optimized using different types of atomic charges. Moreover, the influence of the electrostatic and van der Waals buffering constants and 1,4-interactions scaling factors used in the force field were also tested, either being optimized as well or fixed with respect to the values of CHARMM force field. Results show that some atomic charges are not well adapted to CHARMM parameters and lead to unrealistic MM-optimized structures or a MM divergence. As a result, a new scaling factor has been optimized for Quantum Theory of Atoms in Molecules charges and charges derived from the electrostatic potential such as ChelpG. The force fields optimized can be mixed with the CHARMM force field, without changing it, to study for the first time hepta-anions interacting with organic molecules.

Validating CHARMM parameters and exploring charge distribution rules in structure-based drug design

Using an extensive series of TIBO compounds that are non-nucleoside inhibitors of HIV-1 reverse transcriptase, we have systematically evaluated the quality of recently developed ligand parameters that are consistent with the CHARMM22 force field. Thermodynamic integration simulations for 44 pairs of TIBO compounds achieve a high level of success with an overall average unsigned error (AUE) in the relative binding affinities of 1.3 kcal/mol; however, the accuracy is strongly dependent on the size differential between the substituents sampled as well as the class of functional group. Low errors are observed among the alkyl, allyl, aldehyde, nitrile, trifluorinated methyl, and halide TIBO derivatives and large systematic errors among thioether derivatives. We have also investigated how different charge assignment schemes for small molecules impact the quality of computed binding affinities for a subset of this series. This study demonstrates the advantage of using model compounds to derive physically meaningful charge distributions and bond-charge increments for rapidly expanding fragment libraries for drug development applications. Specifically, in the absence of a bond-charge increment for a given pair of atom types, the strategy of adopting CHELPG charges from localized regions of model compounds provides reliable results when modeling with the CHARMM force field.

Charge parameterization of the metal centers in cytochrome c oxidase

Reliable atomic point charges are of key importance for a correct description of the electrostatic interactions when performing classical, force field based simulations. Here, we present a systematic procedure for point charge derivation, based on quantum mechanical methodology suited for the systems at hand. A notable difference to previous procedures is to include an outer region around the actual system of interest. At the cost of increasing the system sizes, here up to 265 atoms, including the surroundings achieves near-neutrality for the systems as well as structural stability, important factors for reliable charge distributions. In addition, the common problem of converting between C--H bonds and C--C bonds at the border vanishes. We apply the procedure to the four redox-active metal centers of cytochrome c oxidase: Cu(A), haem a, haem a(3), and Cu(B). Several relevant charge and ligand states are considered. Charges for two different force fields, CHARMM and AMBER, are presented.