1-Octanol/water partition coefficients of n-alkanes from molecular simulations of absolute solvation free energies

Nuno M. Garrido, Antonio J. Queimada, Miguel Jorge, Eugenia A. Macedo, Ioannis G. Economou

Research output: Contribution to journalArticle

80 Citations (Scopus)

Abstract

The 1-octanol/water partition coefficient is an important thermodynamic variable usually employed to understand and quantify the partitioning of solutes between aqueous and organic phases. It finds widespread use in many empirical correlations to evaluate the environmental fate of pollutants as well as in the design of pharmaceuticals. The experimental evaluation of 1-octanol/water partition coefficients is an expensive and time-consuming procedure, and thus, theoretical estimation methods are needed, particularly when a physical sample of the solute may not yet be available, such as in pharmaceutical screening. 1-Octanol/water partition coefficients can be obtained from Gibbs free energies of solvation of the solute in both the aqueous and the octanol phases. The accurate evaluation of free energy differences remains today a challenging problem in computational chemistry. In order to study the absolute solvation Gibbs free energies in 1-octanol, a solvent that can mimic many properties of important biological systems, free energy calculations for n-alkanes in the range C-1-C-8 were performed using molecular simulation techniques, following the thermodynamic integration approach. In the first part of this paper, we test different force fields by evaluating their performance in reproducing pure 1-octanol properties. It is concluded that all-atom force fields can provide good accuracy but at the cost of a higher computational time compared to that of the united-atom force fields. Recent versions of united-atom force fields, such as Gromos and TraPPE, provide satisfactory results and are, thus, useful alternatives to the more expensive all-atom models. In the second part of the paper, the Gibbs free energy of solvation in 1-octanol is calculated for several n-alkanes using three force fields to describe the solutes, namely Gromos, TraPPE, and OPLS-AA. Generally, the results obtained are in excellent agreement with the available experimental data and are of similar accuracy to commonly used QSPR models. Moreover, we have estimated the Gibbs free energy of hydration for the different compounds with the three force fields, reaching average deviations from experimental data of less than 0.2 kcal/mol for the case of the Gromos force field. Finally, we systematically compare different strategies to obtain the 1-octanol/water partition coefficient from the simulations. It is shown that a fully predictive method combining the Gromos force field in the aqueous phase and the OPLS-AA/TraPPE force field for the organic phase can give excellent predictions for n-alkanes up to C-8 with an absolute average deviation of 0.1 log P units to the experimental data.

LanguageEnglish
Pages2436-2446
Number of pages11
JournalJournal of Chemical Theory and Computation
Volume5
Issue number9
DOIs
Publication statusPublished - Sep 2009

Fingerprint

1-Octanol
Alkanes
Solvation
Gibbs free energy
Paraffins
alkanes
field theory (physics)
Free energy
solvation
partitions
free energy
Atoms
Water
coefficients
Drug products
water
simulation
solutes
Thermodynamics
Computational chemistry

Keywords

  • thermodynamic integration
  • solvation
  • united-atom force fields
  • Gromos

Cite this

Garrido, Nuno M. ; Queimada, Antonio J. ; Jorge, Miguel ; Macedo, Eugenia A. ; Economou, Ioannis G. / 1-Octanol/water partition coefficients of n-alkanes from molecular simulations of absolute solvation free energies. In: Journal of Chemical Theory and Computation . 2009 ; Vol. 5, No. 9. pp. 2436-2446.
@article{adb20401f3494aed94294ebdab3f2d30,
title = "1-Octanol/water partition coefficients of n-alkanes from molecular simulations of absolute solvation free energies",
abstract = "The 1-octanol/water partition coefficient is an important thermodynamic variable usually employed to understand and quantify the partitioning of solutes between aqueous and organic phases. It finds widespread use in many empirical correlations to evaluate the environmental fate of pollutants as well as in the design of pharmaceuticals. The experimental evaluation of 1-octanol/water partition coefficients is an expensive and time-consuming procedure, and thus, theoretical estimation methods are needed, particularly when a physical sample of the solute may not yet be available, such as in pharmaceutical screening. 1-Octanol/water partition coefficients can be obtained from Gibbs free energies of solvation of the solute in both the aqueous and the octanol phases. The accurate evaluation of free energy differences remains today a challenging problem in computational chemistry. In order to study the absolute solvation Gibbs free energies in 1-octanol, a solvent that can mimic many properties of important biological systems, free energy calculations for n-alkanes in the range C-1-C-8 were performed using molecular simulation techniques, following the thermodynamic integration approach. In the first part of this paper, we test different force fields by evaluating their performance in reproducing pure 1-octanol properties. It is concluded that all-atom force fields can provide good accuracy but at the cost of a higher computational time compared to that of the united-atom force fields. Recent versions of united-atom force fields, such as Gromos and TraPPE, provide satisfactory results and are, thus, useful alternatives to the more expensive all-atom models. In the second part of the paper, the Gibbs free energy of solvation in 1-octanol is calculated for several n-alkanes using three force fields to describe the solutes, namely Gromos, TraPPE, and OPLS-AA. Generally, the results obtained are in excellent agreement with the available experimental data and are of similar accuracy to commonly used QSPR models. Moreover, we have estimated the Gibbs free energy of hydration for the different compounds with the three force fields, reaching average deviations from experimental data of less than 0.2 kcal/mol for the case of the Gromos force field. Finally, we systematically compare different strategies to obtain the 1-octanol/water partition coefficient from the simulations. It is shown that a fully predictive method combining the Gromos force field in the aqueous phase and the OPLS-AA/TraPPE force field for the organic phase can give excellent predictions for n-alkanes up to C-8 with an absolute average deviation of 0.1 log P units to the experimental data.",
keywords = "thermodynamic integration, solvation , united-atom force fields, Gromos",
author = "Garrido, {Nuno M.} and Queimada, {Antonio J.} and Miguel Jorge and Macedo, {Eugenia A.} and Economou, {Ioannis G.}",
note = "(c) American Chemical Society",
year = "2009",
month = "9",
doi = "10.1021/ct900214y",
language = "English",
volume = "5",
pages = "2436--2446",
journal = "Journal of Chemical Theory and Computation",
issn = "1549-9618",
publisher = "American Chemical Society",
number = "9",

}

1-Octanol/water partition coefficients of n-alkanes from molecular simulations of absolute solvation free energies. / Garrido, Nuno M.; Queimada, Antonio J.; Jorge, Miguel; Macedo, Eugenia A.; Economou, Ioannis G.

In: Journal of Chemical Theory and Computation , Vol. 5, No. 9, 09.2009, p. 2436-2446.

Research output: Contribution to journalArticle

TY - JOUR

T1 - 1-Octanol/water partition coefficients of n-alkanes from molecular simulations of absolute solvation free energies

AU - Garrido, Nuno M.

AU - Queimada, Antonio J.

AU - Jorge, Miguel

AU - Macedo, Eugenia A.

AU - Economou, Ioannis G.

N1 - (c) American Chemical Society

PY - 2009/9

Y1 - 2009/9

N2 - The 1-octanol/water partition coefficient is an important thermodynamic variable usually employed to understand and quantify the partitioning of solutes between aqueous and organic phases. It finds widespread use in many empirical correlations to evaluate the environmental fate of pollutants as well as in the design of pharmaceuticals. The experimental evaluation of 1-octanol/water partition coefficients is an expensive and time-consuming procedure, and thus, theoretical estimation methods are needed, particularly when a physical sample of the solute may not yet be available, such as in pharmaceutical screening. 1-Octanol/water partition coefficients can be obtained from Gibbs free energies of solvation of the solute in both the aqueous and the octanol phases. The accurate evaluation of free energy differences remains today a challenging problem in computational chemistry. In order to study the absolute solvation Gibbs free energies in 1-octanol, a solvent that can mimic many properties of important biological systems, free energy calculations for n-alkanes in the range C-1-C-8 were performed using molecular simulation techniques, following the thermodynamic integration approach. In the first part of this paper, we test different force fields by evaluating their performance in reproducing pure 1-octanol properties. It is concluded that all-atom force fields can provide good accuracy but at the cost of a higher computational time compared to that of the united-atom force fields. Recent versions of united-atom force fields, such as Gromos and TraPPE, provide satisfactory results and are, thus, useful alternatives to the more expensive all-atom models. In the second part of the paper, the Gibbs free energy of solvation in 1-octanol is calculated for several n-alkanes using three force fields to describe the solutes, namely Gromos, TraPPE, and OPLS-AA. Generally, the results obtained are in excellent agreement with the available experimental data and are of similar accuracy to commonly used QSPR models. Moreover, we have estimated the Gibbs free energy of hydration for the different compounds with the three force fields, reaching average deviations from experimental data of less than 0.2 kcal/mol for the case of the Gromos force field. Finally, we systematically compare different strategies to obtain the 1-octanol/water partition coefficient from the simulations. It is shown that a fully predictive method combining the Gromos force field in the aqueous phase and the OPLS-AA/TraPPE force field for the organic phase can give excellent predictions for n-alkanes up to C-8 with an absolute average deviation of 0.1 log P units to the experimental data.

AB - The 1-octanol/water partition coefficient is an important thermodynamic variable usually employed to understand and quantify the partitioning of solutes between aqueous and organic phases. It finds widespread use in many empirical correlations to evaluate the environmental fate of pollutants as well as in the design of pharmaceuticals. The experimental evaluation of 1-octanol/water partition coefficients is an expensive and time-consuming procedure, and thus, theoretical estimation methods are needed, particularly when a physical sample of the solute may not yet be available, such as in pharmaceutical screening. 1-Octanol/water partition coefficients can be obtained from Gibbs free energies of solvation of the solute in both the aqueous and the octanol phases. The accurate evaluation of free energy differences remains today a challenging problem in computational chemistry. In order to study the absolute solvation Gibbs free energies in 1-octanol, a solvent that can mimic many properties of important biological systems, free energy calculations for n-alkanes in the range C-1-C-8 were performed using molecular simulation techniques, following the thermodynamic integration approach. In the first part of this paper, we test different force fields by evaluating their performance in reproducing pure 1-octanol properties. It is concluded that all-atom force fields can provide good accuracy but at the cost of a higher computational time compared to that of the united-atom force fields. Recent versions of united-atom force fields, such as Gromos and TraPPE, provide satisfactory results and are, thus, useful alternatives to the more expensive all-atom models. In the second part of the paper, the Gibbs free energy of solvation in 1-octanol is calculated for several n-alkanes using three force fields to describe the solutes, namely Gromos, TraPPE, and OPLS-AA. Generally, the results obtained are in excellent agreement with the available experimental data and are of similar accuracy to commonly used QSPR models. Moreover, we have estimated the Gibbs free energy of hydration for the different compounds with the three force fields, reaching average deviations from experimental data of less than 0.2 kcal/mol for the case of the Gromos force field. Finally, we systematically compare different strategies to obtain the 1-octanol/water partition coefficient from the simulations. It is shown that a fully predictive method combining the Gromos force field in the aqueous phase and the OPLS-AA/TraPPE force field for the organic phase can give excellent predictions for n-alkanes up to C-8 with an absolute average deviation of 0.1 log P units to the experimental data.

KW - thermodynamic integration

KW - solvation

KW - united-atom force fields

KW - Gromos

UR - http://pubs.acs.org/journal/jctcce

U2 - 10.1021/ct900214y

DO - 10.1021/ct900214y

M3 - Article

VL - 5

SP - 2436

EP - 2446

JO - Journal of Chemical Theory and Computation

T2 - Journal of Chemical Theory and Computation

JF - Journal of Chemical Theory and Computation

SN - 1549-9618

IS - 9

ER -