278234 Direct Calculation of 1-Octanol–Water Partition Coefficients From Adaptive Biasing Force Molecular Dynamics Simulations

Wednesday, October 31, 2012: 10:00 AM
Conference A (Omni )
Navendu Bhatnagar, Chemical Engineering, Wayne State University, Detroit, MI, Ganesh Kamath, Chemistry, University of Missouri-Columbia, Columbia, MO, Isaac Chelst, Chemical Engineering and Materials Science, Wayne State University, Detroit, MI and Jeffrey J. Potoff, Chemcial Engineering and Materials Science, Wayne State University, Detroit, MI

The 1–octanol–water partition coefficient log Kow of a solute is a key parameter used in the prediction of a wide variety of complex phenomena such as drug availability and bioaccumulation potential of trace contaminants. In this work, the application of adaptive biasing force molecular dynamics simulations (ABF–MD)1-4 to the calculation of 1-octanol-water partition coefficients is presented.  Unlike traditional methods, such as thermodynamic integration and free energy perturbation, which produce relative free energy values, the ABF-MD method predicts absolute free energies of transfer and eliminates the need for reference solutes.  Two approaches are evaluated; the direct transfer of the solute from 1–octanol to water phase, and separate transfers of the solute from the water or 1–octanol phase to vacuum, with both methods yielding statistically indistinguishable results.   Simulations performed on a variety of system sizes are used to show that the calculated potential of mean force is remarkably insensitive to the microstructure of the solvent, a finding that suggests relatively small system sizes may be used to determine accurate free energies of transfer.  The effect of various simulation parameters on the sampling efficiency is discussed.

Calculations are performed on a variety of systems, including n-alkanes, alcohols, chemical warfare agents, energetic materials and ionic liquids. The effect of wet vs. dry octanol and choice of water model is assessed for each system.   The predictions of ABF-MD are found to be in excellent agreement with prior calculations performed with thermodynamic integration5 or Gibbs ensemble Monte Carlo6.

1.            E. Darve and A. Pohorille, J Chem Phys 115 (20), 9169-9183 (2001).

2.            E. Darve, M. A. Wilson and A. Pohorille, Mol Simulat 28 (1-2), 113-144 (2002).

3.            J. Henin and C. Chipot, J Chem Phys 121 (7), 2904-2914 (2004).

4.            D. Rodriguez-Gomez, E. Darve and A. Pohorille, J Chem Phys 120 (8), 3563-3578 (2004).

5.            N. M. Garrido, A. J. Queimada, M. Jorge, E. A. Macedo and I. G. Economou, J Chem Theory Comput 5 (9), 2436-2446 (2009).

6.            J. I. Siepmann and B. Chen, Journal of Physical Chemistry B 110 (8), 3555-3563 (2006).

Extended Abstract: File Not Uploaded
See more of this Session: Thermophysical Properties of Energetic Materials
See more of this Group/Topical: Particle Technology Forum