Closure for the Ornstein-Zernike equation with pressure and free energy consistency

Identifying Systematic Force Field Errors Using a 3D-RISM Element Counting Correction

Hydration free energies of small molecules are commonly used as benchmarks for solvation models. However, errors in predicting hydration free energies are partially due to the force fields used and not just the solvation model. To address this, we have used the 3D reference interaction site model (3D-RISM) of molecular solvation and existing benchmark explicit solvent calculations with a simple element count correction (ECC) to identify problems with the non-bond parameters in the general AMBER force field (GAFF). 3D-RISM was used to calculate hydration free energies of all 642 molecules in the FreeSolv database, and a partial molar volume correction (PMVC), ECC, and their combination (PMVECC) were applied to the results. The PMVECC produced a mean unsigned error of 1.01±0.04kcal/mol and root mean squared error of 1.44±0.07kcal/mol, better than the benchmark explicit solvent calculations from FreeSolv, and required less than 15 s of computing time per molecule on a single CPU core. Importantly, parameters for PMVECC showed systematic errors for molecules containing Cl, Br, I, and P. Applying ECC to the explicit solvent hydration free energies found the same systematic errors. The results strongly suggest that some small adjustments to the Lennard-Jones parameters for GAFF will lead to improved hydration free energy calculations for all solvent models.

Data-driven approximations to the bridge function yield improved closures for the Ornstein-Zernike equation

A key challenge for soft materials design and coarse-graining simulations is determining interaction potentials between components that give rise to desired condensed-phase structures. In theory, the Ornstein-Zernike equation provides an elegant framework for solving this inverse problem. Pioneering work in liquid state theory derived analytical closures for the framework. However, these analytical closures are approximations, valid only for specific classes of interaction potentials. In this work, we combine the physics of liquid state theory with machine learning to infer a closure directly from simulation data. The resulting closure is more accurate than commonly used closures across a broad range of interaction potentials.

Application of the 3D-RISM-KH molecular solvation theory for DMSO as solvent

The molecular solvation theory in the form of the Three-Dimensional Reference Interaction Site Model (3D-RISM) with Kovalenko-Hirata (KH) closure relation is benchmarked for use with dimethyl sulfoxide (DMSO) as solvent for (bio)-chemical simulation within the framework of integral equation formalism. Several force field parameters have been tested to correctly reproduce solvation free energy in DMSO, ion solvation in DMSO, and DMSO coordination prediction. Our findings establish a united atom (UA) type parameterization as the best model of DMSO for use in 3D-RISM-KH theory based calculations.

Predicting skin permeability using the 3D-RISM-KH theory based solvation energy descriptors for a diverse class of compounds

The state-of-the-art molecular solvation theory is used to predict skin permeability of a large set of compounds with available experimental skin permeability coefficient (logK_P). Encouraging results are obtained pointing to applicability of a novel quantitative structure activity model that uses statistical physics based 3D-RISM-KH theory for solvation free energy calculations as a primary descriptor for the prediction of logK_P with relative mean square error of 0.77 units.