Toward physics-based solubility computation for pharmaceuticals to rival informatics
MetadataShow full item record
We demonstrate that physics-based calculations of intrinsic aqueous solubility can rival cheminformatics-based machine learning predictions. A proof-of-concept was developed for a physics-based approach via a sublimation thermodynamic cycle, building upon previous work that relied upon several thermodynamic approximations, notably the 2RT approximation, and limited conformational sampling. Here, we apply improvements to our sublimation free-energy model with the use of crystal phonon mode calculations to capture the contributions of the vibrational modes of the crystal. Including these improvements with lattice energies computed using the model-potential-based Ψmol method leads to accurate estimates of sublimation free energy. Combining these with hydration free energies obtained from either molecular dynamics free-energy perturbation simulations or density functional theory calculations, solubilities comparable to both experiment and informatics predictions are obtained. The application to coronene, succinic acid, and the pharmaceutical desloratadine shows how the methods must be adapted for the adoption of different conformations in different phases. The approach has the flexibility to extend to applications that cannot be covered by informatics methods.
Fowles , D J , Palmer , D S , Guo , R , Price , S L & Mitchell , J B O 2021 , ' Toward physics-based solubility computation for pharmaceuticals to rival informatics ' , Journal of Chemical Theory and Computation , vol. 17 , no. 6 , pp. 3700-3709 . https://doi.org/10.1021/acs.jctc.1c00130
Journal of Chemical Theory and Computation
Copyright © 2021 The Authors. Published by American Chemical Society. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.
DescriptionD.S.P. and D.J.F. thank the EPSRC for funding via Prosperity Partnership EP/S035990/1. D.S.P. and D.J.F. thank the ARCHIE-WeSt High-Performance Computing Centre (www.archie-west.ac.uk) for computational resources. The UCL authors thank Prof. Keith Refson for guidance with the phonon calculations, which used the ARCHER U.K. National Supercomputing Service (http://www.archer.ac.uk) as part of the U.K. HEC Materials Chemistry Consortium, which is funded by the EPSRC (EP/L000202, EP/R029431). R.G. was funded by MagnaPharm, a project funded by the European Union’s Horizon 2020 Research and Innovation programme under grant agreement number 736899.
Items in the St Andrews Research Repository are protected by copyright, with all rights reserved, unless otherwise indicated.