The strength of different hydrogen-bond acceptors is quite sensitive to the electronic properties of the acceptor, making it crucial to correctly account for the specific substitution of the molecule in question. For this reason, the efficient and accurate prediction of hydrogen-bond basicity has been a topic of active research in computational chemistry for decades. Since the 1990s, researchers have found that the minimum electrostatic potential V min in the region of lone pairs (Figure ) can be used to predict the strength of various hydrogenbond acceptors, with different scaling parameters employed for each functional group. Nevertheless, the computational cost and complexity of these workflows has limited their broader adoption in rational molecular optimization. Here, we report an efficient workflow for the V min -based automated prediction of hydrogen-bond basicity for unseen organic molecules.