Application of the Approximate 3D-Reference Interaction Site Model (RISM) Molecular Solvation Theory to Acetonitrile as Solvent

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.

Biomolecular Simulations with the Three-Dimensional Reference Interaction Site Model with the Kovalenko-Hirata Closure Molecular Solvation Theory

The statistical mechanics-based 3-dimensional reference interaction site model with the Kovalenko-Hirata closure (3D-RISM-KH) molecular solvation theory has proven to be an essential part of a multiscale modeling framework, covering a vast region of molecular simulation techniques. The successful application ranges from the small molecule solvation energy to the bulk phase behavior of polymers, macromolecules, etc. The 3D-RISM-KH successfully predicts and explains the molecular mechanisms of self-assembly and aggregation of proteins and peptides related to neurodegeneration, protein-ligand binding, and structure-function related solvation properties. Upon coupling the 3D-RISM-KH theory with a novel multiple time-step molecular dynamic (MD) of the solute biomolecule stabilized by the optimized isokinetic Nosé-Hoover chain thermostat driven by effective solvation forces obtained from 3D-RISM-KH and extrapolated forward by generalized solvation force extrapolation (GSFE), gigantic outer time-steps up to picoseconds to accurately calculate equilibrium properties were obtained in this new quasidynamics protocol. The multiscale OIN/GSFE/3D-RISM-KH algorithm was implemented in the Amber package and well documented for fully flexible model of alanine dipeptide, miniprotein 1L2Y, and protein G in aqueous solution, with a solvent sampling rate ~150 times faster than a standard MD simulation in explicit water. Further acceleration in computation can be achieved by modifying the extent of solvation layers considered in the calculation, as well as by modifying existing closure relations. This enhanced simulation technique has proven applications in protein-ligand binding energy calculations, ligand/solvent binding site prediction, molecular solvation energy calculations, etc. Applications of the RISM-KH theory in molecular simulation are discussed in this work.

A Force Field for a Manganese-Vanadium Water Oxidation Catalyst: Redox Potentials in Solution as Showcase

We present a molecular mechanics force field in AMBER format for the mixed-valence manganese vanadium oxide cluster [Mn4V4O17(OAc)3]3−—a synthetic analogue of the oxygen-evolving complex that catalyzes the water oxidation reaction in photosystem II—with parameter sets for two different oxidation states. Most force field parameters involving metal atoms have been newly parametrized and the harmonic terms refined using hybrid quantum mechanics/molecular mechanics reference simulations, although some parameters were adapted from pre-existing force fields of vanadate cages and manganese oxo dimers. The characteristic Jahn–Teller distortions of d4 MnIII ions in octahedral environments are recovered by the force field. As an application, the developed parameters have been used to calculate the redox potential of the [MnIIIMn3IV] ⇌ [Mn4IV]+e− half-reaction in acetonitrile by means of Marcus theory.

Predicting PAMPA permeability using the 3D-RISM-KH theory: are we there yet?

The parallel artificial membrane permeability assay (PAMPA), a non-cellular lab-based assay, is extensively used to measure the permeability of pharmaceutical compounds. PAMPA experiments provide a working mimic of a molecule passing through cells and PAMPA values are widely used to estimate drug absorption parameters. There is an increased interest in developing computational methods to predict PAMPA permeability values. We developed an in silico model to predict the permeability of compounds based on the PAMPA assay. We used the three-dimensional reference interaction site model (3D-RISM) theory with the Kovalenko-Hirata (KH) closure to calculate the excess chemical potentials of a large set of compounds and predicted their apparent permeability with good accuracy (mean absolute error or MAE = 0.69 units) when compared to a published experimental data set. Furthermore, our in silico PAMPA protocol performed very well in the binary prediction of 288 compounds as being permeable or impermeable (precision = 94%, accuracy = 93%). This suggests that our in silico protocol can mimic the PAMPA assay and could aid in the rapid discovery or screening of potentially therapeutic drug leads that can be delivered to a desired tissue.