Gaussian Process Regression Models for the Prediction of Hydrogen Bond Acceptor Strengths
Christoph Alexander Bauer, Gisbert Schneider, Andreas H. Göller · Molecular Informatics · 2018
Abstract We present two approaches for the computation of hydrogen bond acceptor strengths, one by machine‐learning and one by a composite quantum‐mechanical protocol, both based on the well‐established pKBHX scale and dataset. The QM calculations after a necessary linear fit reproduce the complexation free energies in solution with an RMSE of 2.6 kJ mol−1, not far off the expected error of 2 kJ mol−1 obtained from the comparison of experimental data from two different sources. The second approach is by Gaussian Process Regression (GPR) machine‐learning. We describe the hydrogen bond acceptor atoms by a radial atomic reactivity descriptor that encodes their electronic and steric environment. The performance of the GPR model on an external test set corresponds to 3.3 kJ mol−1, which is also close to the experimental error. We apply the GPR model built on experimental data to model the hydrogen bond acceptor strengths of a series of hydrogen bond acceptor sites of 10 phosphodiesterase 10 A inhibitors. The predicted values correlate well with the experimentally measured IC50 values.