Predicting the cytotoxicity of nanomaterials through explainable, extreme gradient boosting

Nanoparticles (NPs) are a wide class of materials currently used in several industrial and biomedical applications. Due to their small size (1-100 nm), NPs can easily enter the human body, inducing tissue damage. NP toxicity depends on physical and chemical NP properties (e.g., size, charge and surface area) in ways and magnitudes that are still unknown. We assess the average as well as the individual importance of NP atomic descriptors, along with chemical properties and experimental conditions, in determining cytotoxicity endpoints for several nanomaterials. We employ a multicenter cytotoxicity nanomaterial database (12 different materials with first and second dimensions ranging between 2.70 and 81.2 nm and between 4.10 and 4048 nm, respectively). We develop a regressor model based on extreme gradient boosting with hyperparameter optimization. We employ Shapley additive explanations to obtain good cytotoxicity prediction performance. Model performances are quantified as statistically significant Spearman correlations between the true and predicted values, ranging from 0.5 to 0.7. Our results show that i) size in situ and surface areas larger than 200 nm and 50 m2/g, respectively, ii) primary particles smaller than 20 nm; iii) irregular (i.e., not spherical) shapes and iv) positive Z-potentials contribute the most to the prediction of NP cytotoxicity, especially if lactate dehydrogenase (LDH) assays are employed for short experimental times. These results were moderately stable across toxicity endpoints, although some degree of variability emerged across dose quantification methods, confirming the complexity of nano-bio interactions and the need for large, systematic experimental characterization to reach a safer-by-design approach.