Predicting octanol/water partition coefficients for the SAMPL6 challenge using the SM12, SM8, and SMD solvation models

Blind predictions of octanol/water partition coefficients at 298 K for 11 kinase inhibitor fragment like compounds were made for the SAMPL6 challenge. We used the conventional, “untrained”, free energy based approach wherein the octanol/water partition coefficient was computed directly as the difference in solvation free energy in water and 1-octanol. We additionally proposed and used two different forms of a “trained” approach. Physically, the goal of the trained approach is to relate the partition coefficient computed using pure 1-octanol to that using water-saturated 1-octanol. In the first case, we assumed the partition coefficient using water-saturated 1-octanol and pure 1-octanol are linearly correlated. In the second approach, we assume the solvation free energy in water-saturated 1-octanol can be written as a linear combination of the solvation free energy in pure water and 1-octanol. In all cases here, the solvation free energies were computed using electronic structure calculations in the SM12, SM8, and SMD universal solvent models. In the context of the present study, our results in general do not support the additional effort of the trained approach.

[1]  Jonathan W. Essex,et al.  Theoretical determination of partition coefficients , 1992 .

[2]  J. Sangster Octanol-Water Partition Coefficients: Fundamentals and Physical Chemistry , 1997 .

[3]  J. S. Rowlinson,et al.  Molecular Thermodynamics of Fluid-Phase Equilibria , 1969 .

[4]  David S. Wishart,et al.  DrugBank 4.0: shedding new light on drug metabolism , 2013, Nucleic Acids Res..

[5]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[6]  P. Kollman,et al.  Automatic atom type and bond type perception in molecular mechanical calculations. , 2006, Journal of molecular graphics & modelling.

[7]  Andrew S. Paluch,et al.  Predicting the Solubility of Solid Phenanthrene: A Combined Molecular Simulation and Group Contribution Approach , 2013 .

[8]  Li Yang,et al.  Comparison of two simulation methods to compute solvation free energies and partition coefficients , 2013, J. Comput. Chem..

[9]  Ganesh Kamath,et al.  Direct calculation of 1-octanol-water partition coefficients from adaptive biasing force molecular dynamics simulations. , 2012, The Journal of chemical physics.

[10]  Peter W Carr,et al.  The chemical interpretation and practice of linear solvation energy relationships in chromatography. , 2006, Journal of chromatography. A.

[11]  David S. Wishart,et al.  HMDB: a knowledgebase for the human metabolome , 2008, Nucleic Acids Res..

[12]  Haiyang Zhang,et al.  Force Field Benchmark of Amino Acids. 2. Partition Coefficients between Water and Organic Solvents , 2018, J. Chem. Inf. Model..

[13]  Andrew S. Paluch,et al.  Using limiting activity coefficients to efficiently evaluate the ability of fixed-charge force fields to model miscible water plus cosolvent mixtures , 2018, Molecular Simulation.

[14]  A. Pohorille,et al.  Free energy calculations : theory and applications in chemistry and biology , 2007 .

[15]  Alán Aspuru-Guzik,et al.  Advances in molecular quantum chemistry contained in the Q-Chem 4 program package , 2014, Molecular Physics.

[16]  David S. Wishart,et al.  HMDB 4.0: the human metabolome database for 2018 , 2017, Nucleic Acids Res..

[17]  David S. Wishart,et al.  HMDB 3.0—The Human Metabolome Database in 2013 , 2012, Nucleic Acids Res..

[18]  Shinichiro Nakamura,et al.  Effect of Atomic Charges on Octanol–Water Partition Coefficient Using Alchemical Free Energy Calculation , 2018, Molecules.

[19]  J. Sangster,et al.  Octanol‐Water Partition Coefficients of Simple Organic Compounds , 1989 .

[20]  David S. Wishart,et al.  DrugBank 3.0: a comprehensive resource for ‘Omics’ research on drugs , 2010, Nucleic Acids Res..

[21]  Andrew S. Paluch,et al.  Microscopic Structure and Solubility Predictions of Multifunctional Solids in Supercritical Carbon Dioxide: A Molecular Simulation Study. , 2017, The journal of physical chemistry. B.

[22]  Binne Zwanenburg,et al.  A two-step chirality transfer from (−)-endo- to (−)-exo-tricyclo[5.2.1.026]deca-4,8-dien-3-one , 1993 .

[23]  Bin Chen,et al.  Microscopic structure and solvation in dry and wet octanol. , 2006, The journal of physical chemistry. B.

[24]  Ioannis G. Economou,et al.  Prediction of the n‐hexane/water and 1‐octanol/water partition coefficients for environmentally relevant compounds using molecular simulation , 2012 .

[25]  Julien Michel,et al.  Prediction of partition coefficients by multiscale hybrid atomic-level/coarse-grain simulations. , 2008, The journal of physical chemistry. B.

[26]  Daisy Y. Kyu,et al.  Calculating Partition Coefficients of Small Molecules in Octanol/Water and Cyclohexane/Water. , 2016, Journal of chemical theory and computation.

[27]  David S. Wishart,et al.  DrugBank: a comprehensive resource for in silico drug discovery and exploration , 2005, Nucleic Acids Res..

[28]  T. H. Dunning Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[29]  C. Cramer,et al.  Self-Consistent Reaction Field Model for Aqueous and Nonaqueous Solutions Based on Accurate Polarized Partial Charges. , 2007, Journal of chemical theory and computation.

[30]  David Weininger,et al.  SMILES, a chemical language and information system. 1. Introduction to methodology and encoding rules , 1988, J. Chem. Inf. Comput. Sci..

[31]  D. Truhlar,et al.  The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals , 2008 .

[32]  David L. Mobley,et al.  FreeSolv: a database of experimental and calculated hydration free energies, with input files , 2014, Journal of Computer-Aided Molecular Design.

[33]  Andrew S. Paluch,et al.  Developing a Predictive Form of MOSCED for Nonelectrolyte Solids Using Molecular Simulation: Application to Acetanilide, Acetaminophen, and Phenacetin , 2016 .

[34]  J. Ilja Siepmann,et al.  Partitioning of Alkane and Alcohol Solutes between Water and (Dry or Wet) 1-Octanol , 2000 .

[35]  J. Rumble CRC Handbook of Chemistry and Physics , 2019 .

[36]  D. van der Spoel,et al.  Prediction of Partition Coefficients of Environmental Toxins Using Computational Chemistry Methods , 2019, ACS omega.

[37]  David L. Mobley,et al.  Octanol–water partition coefficient measurements for the SAMPL6 blind prediction challenge , 2019, Journal of Computer-Aided Molecular Design.

[38]  S. Sandler,et al.  Prediction of Octanol−Water Partition Coefficients Using a Group Contribution Solvation Model , 1999 .

[39]  William L. Jorgensen,et al.  Relative partition coefficients for organic solutes from fluid simulations , 1990 .

[40]  David S. Wishart,et al.  DrugBank 5.0: a major update to the DrugBank database for 2018 , 2017, Nucleic Acids Res..

[41]  A. Leo,et al.  Partition coefficients and their uses , 1971 .

[42]  C. Cramer,et al.  Universal solvation model based on solute electron density and on a continuum model of the solvent defined by the bulk dielectric constant and atomic surface tensions. , 2009, The journal of physical chemistry. B.

[43]  M. Abraham,et al.  Abraham model correlations for solute transfer into 2-methyl-2-butanol based on measured activity coefficient and solubility data at 298.15 K , 2019, Journal of Molecular Liquids.

[44]  S. Sandler,et al.  Determination of infinite dilution activity coefficients and 1-octanol/water partition coefficients of volatile organic pollutants , 1994 .

[45]  Jonathan W. Essex,et al.  Relative partition coefficients from partition functions: a theoretical approach to drug transport , 1989 .

[46]  Junming Ho,et al.  Predicting Octanol-Water Partition Coefficients: Are Quantum Mechanical Implicit Solvent Models Better than Empirical Fragment-Based Methods? , 2019, The journal of physical chemistry. B.

[47]  David L. Mobley,et al.  Octanol–water partition coefficient measurements for the SAMPL6 blind prediction challenge , 2019, bioRxiv.

[48]  Justin L MacCallum,et al.  Structures of neat and hydrated 1-octanol from computer simulations. , 2002, Journal of the American Chemical Society.

[49]  Ganesh Kamath,et al.  Prediction of 1-octanol-water and air-water partition coefficients for nitro-aromatic compounds from molecular dynamics simulations. , 2013, Physical chemistry chemical physics : PCCP.

[50]  Johann Gasteiger,et al.  A new model for calculating atomic charges in molecules , 1978 .

[51]  Vasil Simeonov,et al.  Calculating the Partition Coefficients of Organic Solvents in Octanol/Water and Octanol/Air , 2019, J. Chem. Inf. Model..

[52]  F. Lombardo,et al.  Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. , 2001, Advanced drug delivery reviews.

[53]  Miguel Jorge,et al.  1-Octanol/Water Partition Coefficients of n-Alkanes from Molecular Simulations of Absolute Solvation Free Energies. , 2009, Journal of chemical theory and computation.

[54]  Sanket A. Deshmukh,et al.  Thermodynamic considerations for solubility and conformational transitions of poly-N-isopropyl-acrylamide. , 2013, Physical chemistry chemical physics : PCCP.

[55]  Donald G Truhlar,et al.  Generalized Born Solvation Model SM12. , 2013, Journal of chemical theory and computation.

[56]  G. Vesnaver,et al.  Thermodynamic investigation of the effect of octanol–water mutual miscibility on the partitioning and solubility of some guanine derivatives , 1995 .

[57]  Chris Morley,et al.  Open Babel: An open chemical toolbox , 2011, J. Cheminformatics.

[58]  Ying Zhang,et al.  HMDB: the Human Metabolome Database , 2007, Nucleic Acids Res..