Considering ionic state in modeling sorption of pharmaceuticals to sewage sludge

Information on the partitioning of chemicals between particulate matter and water in sewage treatment plants (STPs) can be used to predict their subsequent environmental fate. However, this information can be challenging to acquire, especially for pharmaceuticals that are frequently present in ionized forms. This study investigated the relationship between the ionization state of active pharmaceutical ingredients (APIs) and their partitioning between water and sludge in STPs. We also investigated the underlying mechanisms of sludge sorption by using chemical descriptors based on ionized structures, and evaluated the usefulness of these descriptors in quantitative structure-property relationship (QSPR) modeling. K_D values were collected for 110 APIs, which were classified as neutral, positive, or negative at pH 7. The models with the highest performance had the R²Y and Q² values of above 0.75 and 0.65, respectively. We found that the dominant intermolecular forces governing the interactions of neutral and positively charged APIs with sludge are hydrophobic, pi-pi, and dipole-dipole interactions, whereas the interactions of negatively charged APIs with sludge were mainly governed by covalent bonding as well as ion-ion, ion-dipole, and dipole-dipole interactions; hydrophobicity-driven interactions were rather unimportant. Including charge-related descriptors improved the models' performance by 5-10%, underlining the importance of electrostatic interactions. The use of descriptors calculated for ionized structures did not improve the model statistics for positive and negative APIs, but slightly increased model performance for neutral APIs. We attribute this to a better description of neutral zwitterions.