Skip to main content
Log in

LogP prediction performance with the SMD solvation model and the M06 density functional family for SAMPL6 blind prediction challenge molecules

  • Published:
Journal of Computer-Aided Molecular Design Aims and scope Submit manuscript

Abstract

This work presents a quantum mechanical model for predicting octanol-water partition coefficients of small protein-kinase inhibitor fragments as part of the SAMPL6 LogP Prediction Challenge. The model calculates solvation free energy differences using the M06-2X functional with SMD implicit solvation and the def2-SVP basis set. This model was identified as dqxk4 in the SAMPL6 Challenge and was the third highest performing model in the physical methods category with 0.49 log Root Mean Squared Error (RMSE) for predicting the 11 compounds in SAMPL6 blind prediction set. We also collaboratively investigated the use of empirical models to address model deficiencies for halogenated compounds at minimal additional computational cost. A mixed model consisting of the dqxk4 physical and hdpuj empirical models found improved performance at 0.34 log RMSE on the SAMPL6 dataset. This collaborative mixed model approach shows how empirical models can be leveraged to expediently improve performance in chemical spaces that are difficult for ab initio methods to simulate.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Institutional subscriptions

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6

Similar content being viewed by others

References

  1. Vlahovic F et al (2017) Density functional theory calculation of lipophilicity for organophosphate type pesticides. J Serb Chem Soc 82:104–104

    Article  Google Scholar 

  2. Michalík M, Lukeš V (2016) The validation of quantum chemical lipophilicity prediction of alcohols. Acta Chim Slov 9(2):89

    Article  Google Scholar 

  3. Zhang J et al (2017) Comparison of implicit and explicit solvent models for the calculation of solvation free energy in organic solvents. J Chem Theory Comput 13(3):1034–1043

    Article  CAS  Google Scholar 

  4. Bayat Z, Movaffagh J (2010) The 1-octanol/water partition coefficient of nucleoside analogs via free energy estimated in quantum chemical calculations. Russ J Phys Chem A 84(13):2293–2299

    Article  CAS  Google Scholar 

  5. Marenich AV, Cramer CJ, Truhlar DG (2009) 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. J Phys Chem B 113(18):6378–6396

    Article  CAS  Google Scholar 

  6. Neese F (2012) The ORCA program system. Wiley Interdiscip Rev: Computat Mol Sci 2(1):73–78

    CAS  Google Scholar 

  7. Neese F (2018) Software update: the ORCA program system, version 4.0. Wiley Interdiscip Rev: Comput Mol Sci 8(1):e1327

    Google Scholar 

  8. Zhao Y, Truhlar DG (2008) 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. Theor Chem Acc 120(1):215–241

    Article  CAS  Google Scholar 

  9. Nedyalkova MA et al (2019) Calculating the partition coefficients of organic solvents in octanol/water and octanol/air. J Chem Inf Model 59(5):2257–2263

    Article  CAS  Google Scholar 

  10. Jones MR, Brooks BR, Wilson AK (2016) Partition coefficients for the SAMPL5 challenge using transfer free energies. J Comput Aided Mol Des 30(11):1129–1138

    Article  CAS  Google Scholar 

  11. Zhao Y, Truhlar DG (2008) 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. Theor Chem Acc 120(1–3):215–241

    Article  CAS  Google Scholar 

  12. Zhao Y, Truhlar DG (2006) A new local density functional for main-group thermochemistry, transition metal bonding, thermochemical kinetics, and noncovalent interactions. J Chem Phys 125(19):194101

    Article  Google Scholar 

  13. Işık M et al (2019) Octanol-water partition coefficient measurements for the SAMPL6 blind prediction challenge. bioRxiv. https://doi.org/10.1007/s10822-019-00271-3

    Article  Google Scholar 

  14. OECD (2004) Test No. 117: Partition coefficient (n-octanol/water), HPLC method, OECD guidelines for the testing of chemicals, Section 1, OECD Publishing, Paris. https://doi.org/10.1787/9789264069824-en

  15. Berthold MR et al (2008) KNIME: the Konstanz information miner. Springer, Berlin

    Google Scholar 

  16. O'Boyle NM et al (2011) Open babel: an open chemical toolbox. J Cheminform 3(1):33

    Article  CAS  Google Scholar 

  17. RDKit: Open-source cheminformatics. http://www.rdkit.org

  18. Halgren TA (1996) Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94. J Comput Chem 17(5-6):490–519

    Article  CAS  Google Scholar 

  19. Weigend F, Ahlrichs R (2005) Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: design and assessment of accuracy. Phys Chem Chem Phys 7(18):3297–3305

    Article  CAS  Google Scholar 

  20. Grimme S et al (2010) A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J Chem Phys 132(15):154104

    Article  Google Scholar 

  21. Weigend F, Kattannek M, Ahlrichs R (2009) Approximated electron repulsion integrals: Cholesky decomposition versus resolution of the identity methods. J Chem Phys 130(16):164106

    Article  Google Scholar 

  22. Hohenstein EG, Chill ST, Sherrill CD (2008) Assessment of the performance of the M05-2X and M06-2X exchange-correlation functionals for noncovalent interactions in biomolecules. J Chem Theory Comput 4(12):1996–2000

    Article  CAS  Google Scholar 

  23. Goerigk L, Grimme S (2011) A thorough benchmark of density functional methods for general main group thermochemistry, kinetics, and noncovalent interactions. Phys Chem Chem Phys 13(14):6670–6688

    Article  CAS  Google Scholar 

  24. D. Laikov, C. van Wüllen. Lebedev-Laikov-Grids. http://www.ccl.net/cca/software/SOURCES/FORTRAN/Lebedev-Laikov-Grids/index.shtml

  25. Kossmann S, Neese F (2009) Comparison of two efficient approximate Hartee-Fock approaches. Chem Phys Lett 481:240–243

    Article  CAS  Google Scholar 

  26. Kozuch S, Martin JML (2013) Halogen bonds: benchmarks and theoretical analysis. J Chem Theory Comput 9(4):1918

    Article  CAS  Google Scholar 

  27. Basdogan Y, Keith JA (2018) A paramedic treatment for modeling explicitly solvated chemical reaction mechanisms. Chem Sci 9(24):5341–5346

    Article  CAS  Google Scholar 

  28. Viswanadhan VN et al (1989) Atomic physicochemical parameters for three dimensional structure directed quantitative structure-activity relationships. 4. Additional parameters for hydrophobic and dispersive interactions and their application for an automated superposition of certain naturally occurring nucleoside antibiotics. J Chem Inf Comput Sci 29(3):163–172

    Article  CAS  Google Scholar 

  29. Li W et al (2019) Efficient corrections for DFT noncovalent interactions based on ensemble learning models. J Chem Inf Model 59(5):1849–1857

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We acknowledge the National Institutes of Health for supporting the experimental work carried out in the SAMPL6 logP Prediction Challenge.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Slade Matthews.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Guan, D., Lui, R. & Matthews, S. LogP prediction performance with the SMD solvation model and the M06 density functional family for SAMPL6 blind prediction challenge molecules. J Comput Aided Mol Des 34, 511–522 (2020). https://doi.org/10.1007/s10822-020-00278-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s10822-020-00278-1

Keywords

Navigation