Fast intermolecular interaction energy calculation with the OPLS-AA force field

Abstract

A fast pipeline based on the OPLS-AA force field is presented that enables the automated calculation of non-bonding intermolecular (dimer) interaction energies derived from a set of small organic (monomer) molecules. To calculate the non-bonding contributions, optimized geometries of the monomer molecules as well as their OPLS-AA van der Waals and atomic partial charge parameters are required. The key advantage of the new pipeline lies in its fast, highly parallelized in-memory computations without slow I/O operations, which make it possible to thoroughly process comparatively large numbers of (more than a hundred) monomer molecules within acceptable time frames (hours and days): Compared to an analogous, more versatile, and comprehensive approach, the new computational scheme delivers comparable results while being more than two orders of magnitude faster. For locally estimating the required OPLS-AA force field parameters with the LigParGen and BOSS software packages, a user-friendly graphical user interface for the Windows operating system is provided. Scientific contribution Rapid calculation of mutual intermolecular energies for a set of monomer molecules based on the widely used OPLS-AA force field enables a considerable expansion of this type of calculation for practical purposes. The performance improvement can be used to achieve significant improvements of interaction energy averages as well as considerable expansions in the size of the monomer molecule set.