“…(2) the 3D atomic coordinates of each molecular structure were calculated by structure embedding using the ETKDG method [28]; (3) the molecular geometry was optimized in the RDKit module [29] by the MMFF94 [30] force field with 200 iterations, selecting the minimum energy conformation; (4) the atomic partial charges were computed by means of the web-based tool Atomic Charge Calculator [31] using the B3LYP/6-311G/NPA approach [32]. The structural and physicochemical features of chemicals were then encoded through the following 18 molecular descriptors: molecular weight (MW), hydrophilic factor (Hy) [33], total topological polar surface area (TPSA) [34] and 15 charge descriptors (qpmax, qnmax, Qpos, Qneg, Qtot, Qmean, Q2, RPCG, RNCG, SPP, TE1, TE2, PCWTE1, PCWTE2, LDI) [35][36][37][38][39][40][41]. These descriptors were selected on the basis of the literature review and calculated by the software Dragon 7 [42].…”