Polarizable Empirical Force Field for Halogen-Containing Compounds Based on the Classical Drude Oscillator

The quality of the force field is crucial to ensure the accuracy of simulations used in molecular modeling, including computer-aided drug design (CADD). To perform more accurate modeling and simulations of halogenated molecules, in this study the polarizable force field based on the classical Drude oscillator model was extended to both aliphatic and aromatic systems using halogenated ethane and benzene model compounds for the halogens F, Cl, Br, and I. The force field parameters were optimized targeting quantum mechanical dipole moments, water interactions, and molecular polarizabilities as well as experimental observables, including enthalpies of vaporization, molecular volumes, hydration free energies, and dielectric constants. The developed halogenated polarizable force field is capable of reproducing QM relative energies and geometries of both halogen bonds and halogen-hydrogen bond donor interactions at an unprecedented level due to the inclusion of a virtual particle and anisotropic atomic polarizability on the halogen and, notably, the inclusion of Lennard-Jones parameters on the halogen Drude particle. The model was validated on the basis of its ability to accurately reproduce pure solvent properties for halogenated naphthalenes and alkanes, including species analogous to those used as refrigerants. Accordingly, it is anticipated that the model will be applicable for the study of halogenated derivatives in CADD as well as in other chemical and biophysical studies.