Identification of structural fingerprints for in vivo toxicity by using Monte Carlo based QSTR modeling of nitroaromatics

A novel procedure for selection of molecular descriptors: QSAR model for mutagenicity of nitroaromatic compounds

Nitroaromatic compounds (NACs) stand out as pervasive organic pollutants, prompting an imperative need to investigate their hazardous effects. Computational chemistry methods play a crucial role in this exploration, offering a safer and more time-efficient approach, mandated by various legislations. In this study, our focus lay on the development of transparent, interpretable, reproducible, and publicly available methodologies aimed at deriving quantitative structure-activity relationship models and testing them by modelling the mutagenicity of NACs against the Salmonella typhimurium TA100 strain. Descriptors were selected from Mordred and RDKit molecular descriptors, along with several quantum chemistry descriptors. For that purpose, the genetic algorithm (GA), as the most widely used method in the literature, and three alternative algorithms (Boruta, Featurewiz, and ForwardSelector) combined with the forward stepwise selection technique were used. The construction of models utilized the multiple linear regression method, with subsequent scrutiny of fitting and predictive performance, reliability, and robustness through various statistical validation criteria. The models were ranked by the Multi-Criteria Decision Making procedure. Findings have revealed that the proposed methodology for descriptor selection outperforms GA, with Featurewiz showing a slight advantage over Boruta and ForwardSelector. These constructed models can serve as valuable tools for the quick and reliable prediction of NACs mutagenicity.

In Silico Prediction of the Toxicity of Nitroaromatic Compounds: Application of Ensemble Learning QSAR Approach

In this work, a dataset of more than 200 nitroaromatic compounds is used to develop Quantitative Structure-Activity Relationship (QSAR) models for the estimation of in vivo toxicity based on 50% lethal dose to rats (LD₅₀). An initial set of 4885 molecular descriptors was generated and applied to build Support Vector Regression (SVR) models. The best two SVR models, SVR_A and SVR_B, were selected to build an Ensemble Model by means of Multiple Linear Regression (MLR). The obtained Ensemble Model showed improved performance over the base SVR models in the training set (R² = 0.88), validation set (R² = 0.95), and true external test set (R² = 0.92). The models were also internally validated by 5-fold cross-validation and Y-scrambling experiments, showing that the models have high levels of goodness-of-fit, robustness and predictivity. The contribution of descriptors to the toxicity in the models was assessed using the Accumulated Local Effect (ALE) technique. The proposed approach provides an important tool to assess toxicity of nitroaromatic compounds, based on the ensemble QSAR model and the structural relationship to toxicity by analyzed contribution of the involved descriptors.

Chemometric QSAR modeling of acute oral toxicity of Polycyclic Aromatic Hydrocarbons (PAHs) to rat using simple 2D descriptors and interspecies toxicity modeling with mouse

The information of the acute oral toxicity for most polycyclic aromatic hydrocarbons (PAHs) in mammals are lacking due to limited experimental resources, leading to a need to develop reliable in silico methods to evaluate the toxicity endpoint. In this study, we developed the quantitative structure-activity relationship (QSAR) models by genetic algorithm (GA) and multiple linear regression (MLR) for the rat acute oral toxicity (LD₅₀) of PAHs following the strict validation principles of QSAR modeling recommended by OECD. The best QSAR model comprised eight simple 2D descriptors with definite physicochemical meaning, which showed that maximum atom-type electrotopological state, van der Waals surface area, mean atomic van der Waals volume, and total number of bonds are main influencing factors for the toxicity endpoint. A true external set (554 compounds) without rat acute oral toxicity values, and 22 limit test compounds, were firstly predicted along with reliability assessment. We also compared our proposed model with the OPERA predictions and recently published literature to prove the prediction reliability. Furthermore, the interspecies toxicity (iST) models of PAHs between rat and mouse were also established, validated and employed for filling data gap. Overall, our developed models should be applicable to new or untested or not yet synthesized PAHs falling within the applicability domain (AD) of the models for rapid acute oral toxicity prediction, thus being important for environmental or personal exposure risk assessment under regulatory frameworks.

Quantitative Structure-Activity Relationship (QSAR) Studies on the Toxic Effects of Nitroaromatic Compounds (NACs): A Systematic Review

Nitroaromatic compounds (NACs) are ubiquitous in the environment due to their extensive industrial applications. The recalcitrance of NACs causes their arduous degradation, subsequently bringing about potential threats to human health and environmental safety. The problem of how to effectively predict the toxicity of NACs has drawn public concern over time. Quantitative structure-activity relationship (QSAR) is introduced as a cost-effective tool to quantitatively predict the toxicity of toxicants. Both OECD (Organization for Economic Co-operation and Development) and REACH (Registration, Evaluation and Authorization of Chemicals) legislation have promoted the use of QSAR as it can significantly reduce living animal testing. Although numerous QSAR studies have been conducted to evaluate the toxicity of NACs, systematic reviews related to the QSAR modeling of NACs toxicity are less reported. The purpose of this review is to provide a thorough summary of recent QSAR studies on the toxic effects of NACs according to the corresponding classes of toxic response endpoints.

Can the Monte Carlo method predict the toxicity of binary mixtures?

Risk assessment of toxicants mainly is a result of experiments with single substances. However, toxicity in natural ecosystems typically does not result from single toxicant exposure but is rather a result of exposure to mixtures of toxicants. It is not surprising a mixture of toxicity is a subject of eco-toxicological interest for several decades. A quantitative structure-activity relationships (QSAR)-based approach is an attractive approach to assessing the joint effects in the binary mixtures. The validity of the proposed approach was demonstrated by comparing the predicted values against the experimentally determined values. Simplified molecular input-line entry system (SMILES) is used for the representation of the molecular structures of components of two-component mixtures to build up QSAR. The SMILES-based models are improving if the Monte Carlo optimization aimed to define 2D-optimal descriptors apply the so-called index of ideality of correlation (IIC), which is a mathematical function of both the correlation coefficient and mean absolute error calculated for the positive and negative difference between observed and calculated values of toxicity. The average statistical quality of these models (for the validation set) is n=25, R²=0.95, and RMSE=0.375.

In vivo toxicity of nitroaromatic compounds to rats: QSTR modelling and interspecies toxicity relationship with mouse

Nitroaromatic compounds (NACs) in the environment can cause serious public health and environmental problems due to their potential toxicity. This study established quantitative structure-toxicity relationship (QSTR) models for the acute oral toxicity of NACs towards rats following the stringent OECD principles for QSTR modelling. All models were assessed by various internationally accepted validation metrics and the OECD criteria. The best QSTR model contains seven simple and interpretable 2D descriptors with defined physicochemical meaning. Mechanistic interpretation indicated that van der Waals surface area, presence of C-F at topological distance 6, heteroatom content and frequency of C-N at topological distance 9 are main factors responsible for the toxicity of NACs. This proposed model was successfully applied to a true external set (295 compounds), and prediction reliability was analysed and discussed. Moreover, the rat-mouse and mouse-rat interspecies quantitative toxicity-toxicity relationship (iQTTR) models were also constructed, validated and employed in toxicity prediction for true external sets consisting of 67 and 265 compounds, respectively. These models showed good external predictivity that can be used to rapidly predict the rat oral acute toxicity of new or untested NACs falling within the applicability domain of the models, thus being beneficial in environmental risk assessment and regulatory purposes.