Developing a Machine Learning Model for Hydrogen Bond Acceptance Based on Natural Bond Orbital Descriptors

This study employs machine learning (ML) to assess the predictive power of electronic descriptors derived from natural bond orbital (NBO) analysis for hydrogen bond acceptance. Using a data set of 979 hydrogen bond complexes, each formed by a hydrogen bond acceptor and 4-fluorophenol as the donor, we optimized geometries via GFN2-xTB, followed by DFT single-point calculations. From these, NBO analysis was used to extract intramolecular donor-acceptor interactions, particularly the orbital stabilization energies (E⁽²⁾), which reflect electron delocalization and relate to canonical resonance structures. The E⁽²⁾ values served as features to train seven ML models, based on different techniques: KNN, Decision Tree, SVM, RF, MLP, XGBoost, and CatBoost. To our knowledge, this is the first work that uses E⁽²⁾ as a standalone ML descriptor for hydrogen bond acceptance. Even with a small set of descriptors, we achieved high predictive performance, with errors below 0.4 kcal mol^-1, surpassing previous studies that used heterogeneous descriptors, including quantum-chemical data. Our results highlight the utility of NBO-based features in building accurate, physically meaningful, and generalizable ML models for pK_BHX prediction.