ChemRxiv
These are preliminary reports that have not been peer-reviewed. They should not be regarded as conclusive, guide clinical practice/health-related behavior, or be reported in news media as established information. For more information, please see our FAQs.
1/1
2 files
0/0

Alchemical Hydration Free Energy Calculations Using Molecular Dynamics with Explicit Polarization and Induced Polarity Decoupling: An On-the-Fly Approach

preprint
revised on 06.01.2020 and posted on 07.01.2020 by Braden Kelly, William Smith
We present a methodology using fixed charge force-fields for alchemical solvation free energy calculations which accounts for the change in polarity that the solute experiences as it transfers from the gas-phase to the condensed phase. We update partial charges using QM/MM snapshots, decoupling the electric field appropriately when updating the partial charges. We also show how to account for the cost of self-polarization. We test our methodology on 30 molecules ranging from small polar to large drug-like molecules. We use Minimum Basis Iterative Stockholder (MBIS), Restrained Electrostatic Potential (RESP) and AM1-BCC partial charge methodologies. Using our method with MP2/cc-pVTZ and MBIS partial charges yields an Average Absolute Deviation (AAD) of 6.3 kJ·mol−1 in comparison with the AM1-BCC result of 8.6 kJ·mol−1. AM1-BCC is within experimental uncertainty on 10% of the data compared to 30% with our method. We conjecture that results can be further improved by using Lennard-Jones and torsional parameters refitted to MBIS and RESP partial charge methods that use high levels of theory.

Funding

Natural Sciences and Engineering Research Council of Canada (Strategic Program grant no. STPGP 479466-15)

History

Email Address of Submitting Author

bkelly08@uoguelph.ca

Institution

University of Guelph

Country

Canada

ORCID For Submitting Author

0000-0003-3708-6517

Declaration of Conflict of Interest

no conflict of interest

Version Notes

Latest version (3) fixed bug in periodic boundary conditions of neighbor-list search of python script for updating partial charges which slightly changed some numerical results.

Exports