Toward Physics-Based Solubility Computation for Pharmaceuticals to Rival Informatics
Autor: | Rui Guo, John B. O. Mitchell, Sarah L. Price, Daniel J. Fowles, David S. Palmer |
---|---|
Přispěvatelé: | University of St Andrews. EaSTCHEM, University of St Andrews. School of Chemistry, University of St Andrews. Biomedical Sciences Research Complex |
Rok vydání: | 2021 |
Předmět: |
Work (thermodynamics)
Phonon Computation NDAS Molecular Dynamics Simulation 01 natural sciences Article Machine Learning Molecular dynamics 0103 physical sciences QD Statistical physics Physical and Theoretical Chemistry Density Functional Theory Lattice energy 010304 chemical physics Water QD Chemistry Computer Science Applications Pharmaceutical Preparations Solubility Cheminformatics Thermodynamics Density functional theory Sublimation (phase transition) |
Zdroj: | Journal of Chemical Theory and Computation |
ISSN: | 1549-9626 1549-9618 |
Popis: | D.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. 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. Publisher PDF |
Databáze: | OpenAIRE |
Externí odkaz: |