QSTR with extended topochemical atom (ETA) indices. VI. Acute toxicity of benzene derivatives to tadpoles (Rana japonica)

Elucidation of molecular mechanisms involved in tadpole toxicity employing QSTR and q-RASAR approach

Tadpoles, as early developmental stages of frogs, are vital indicators of toxicity and environmental health in ecosystems exposed to harmful organic compounds from industrial and runoff sources. Evaluating each compound individually is challenging, necessitating the use of in silico methods like Quantitative Structure Toxicity-Relationship (QSTR) and Quantitative Read-Across Structure-Activity Relationship (q-RASAR). Utilizing the comprehensive US EPA's ECOTOX database, which includes acute LC50 toxicity and chronic endpoints, we extracted crucial data such as study types, exposure routes, and chemical categories. Regression-based QSTR and q-RASAR models were developed from this dataset, emphasizing key chemical descriptors. Lipophilicity and unsaturation were significant for predicting acute toxicity, while electrophilicity, nucleophilicity, and molecular branching were crucial for chronic toxicity predictions. Additionally, q-RASAR models integrated with the "intelligent consensus" algorithm were employed to enhance predictive accuracy. The performance of these models was rigorously compared across various approaches. These refined models not only predict the toxicity of untested compounds but also reveal underlying structural influences. Validation through comparison with existing literature affirmed the relevance and robustness of our approach in ecotoxicology.

Does the accounting of the local symmetry fragments in quasi-SMILES improve the predictive potential of the QSAR models of toxicity toward tadpoles?

Models of toxicity to tadpoles have been developed as single parameters based on special descriptors which are sums of correlation weights, molecular features, and experimental conditions. This information is presented by quasi-SMILES. Fragments of local symmetry (FLS) are involved in the development of the model and the use of FLS correlation weights improves their predictive potential. In addition, the index of ideality correlation (IIC) and correlation intensity index (CII) are compared. These two potential predictive criteria were tested in models built through Monte Carlo optimization. The CII was more effective than IIC for the models considered here.

Quasi-SMILES: Self-consistent models for toxicity of organic chemicals to tadpoles

Simplified molecular input-line entry systems (SMILES) are the representation of the molecular structure that can be used to establish quantitative structure-property/activity relationships (QSPRs/QSARs) for various endpoints expressed as mathematical functions of the molecular architecture. Quasi-SMILES is extending the traditional SMILES by means of additional symbols that reflect experimental conditions. Using the quasi-SMILES models of toxicity to tadpoles gives the possibility to build up models by taking into account the time of exposure. Toxic effects of experimental situations expressed via 188 quasi-SMILES (the negative logarithm of molar concentrations which lead to lethal 50% tadpoles effected during 12 h, 24 h, 48 h, 72 h, and 96 h) were modelled with good results (the average determination coefficient for the validation sets is about 0.97). In this way, we developed new models for this amphibian endpoint, which is poorly studied.

A regression-based QSAR-model to predict acute toxicity of aromatic chemicals in tadpoles of the Japanese brown frog (Rana japonica): Calibration, validation, and future developments to support risk assessment of chemicals in amphibians

Amphibian populations are undergoing a global decline worldwide. Such decline has been attributed to their unique physiology, ecology, and exposure to multiple stressors including chemicals, temperature, and biological hazards such as fungi of the Batrachochytrium genus, viruses such as Ranavirus, and habitat reduction. There are limited toxicity data for chemicals available for amphibians and few quantitative structure-activity relationship (QSAR) models have been developed and are publicly available. Such QSARs provide important tools to assess the toxicity of chemicals particularly in a data poor context. QSARs provide important tools to assess the toxicity of chemicals particularly when no toxicological data are available. This manuscript provides a description and validation of a regression-based QSAR model to predict, in a quantitative manner, acute lethal toxicity of aromatic chemicals in tadpoles of the Japanese brown frog (Rana japonica). QSAR models for acute median lethal molar concentrations (LC50-12 h) of waterborne chemicals using the Monte Carlo method were developed. The statistical characteristics of the QSARs were described as average values obtained from five random distributions into training and validation sets. Predictions from the model gave satisfactory results for the overall training set (R² = 0.72 and RMSE = 0.33) and were even more robust for the validation set (R² = 0.96 and RMSE = 0.11). Further development of QSAR models in amphibians, particularly for other life stages and species, are discussed.

Chemometric modeling of acute toxicity of diverse aromatic compounds against Rana japonica

Chemicals used in our daily life show different toxic effects to the aquatic and terrestrial species and thus hamper the ecological balance. In the present time, amphibians are one of them, which are threatened to be extinct. Quantitative structure-activity relationship (QSAR) is an useful tool for prediction involving less time, money and manpower without requiring any animal experiments to assess the unavailable acute toxicity data for the untested molecules. In this study, we have developed QSAR models for ecotoxicity of some waterborne diverse aromatic compounds on an amphibian species Rana japonica (Japanese brown frog) employing Genetic Algorithm (GA) for variable selection followed by Partial Least Squares (PLS) regression method following recommendations of the Organization for Economic Co-operation and Development (OECD) for QSAR model development. Double cross-validation (DCV) followed by Best Subset Selection (BSS) were employed to select suitable models. The models displayed promising statistical quality in terms of R² (= 0.837-0.841), Q²_LOO (= 0.782-0.787), R²_pred or Q²_F1 (= 0.802-0.82) and some other internal and external validation metrics for tadpoles of Rana japonica (N_Training = 44, N_Test = 14). These models can be applied for data gap filling for a new untested compound falling within the applicability domain (AD) of the models.

The “ETA” Indices in QSAR/QSPR/QSTR Research

Descriptors are one of the most essential components of predictive Quantitative Structure-Activity/Property/Toxicity Relationship (QSAR/QSPR/QSTR) modeling analysis, as they encode chemical information of molecules in the form of quantitative numbers, which are used to develop mathematical correlation models. The quality of a predictive model not only depends on good modeling statistics, but also on the extraction of chemical features. A significant amount of research since the beginning of QSAR analysis paradigm has led to the introduction of a large number of predictor variables or descriptors. The Extended Topochemical Atom (ETA) indices, developed by the authors' group, successfully address the aspects of molecular topology, electronic information, and different types of bonded interactions, and have been extensively employed for the modeling of different types of activity/property and toxicity endpoints. This chapter provides explicit information regarding the basis, algorithm, and applicability of the ETA indices for a predictive modeling paradigm.

The “ETA” Indices in QSAR/QSPR/QSTR Research

Descriptors are one of the most essential components of predictive Quantitative Structure-Activity/Property/Toxicity Relationship (QSAR/QSPR/QSTR) modeling analysis, as they encode chemical information of molecules in the form of quantitative numbers, which are used to develop mathematical correlation models. The quality of a predictive model not only depends on good modeling statistics, but also on the extraction of chemical features. A significant amount of research since the beginning of QSAR analysis paradigm has led to the introduction of a large number of predictor variables or descriptors. The Extended Topochemical Atom (ETA) indices, developed by the authors' group, successfully address the aspects of molecular topology, electronic information, and different types of bonded interactions, and have been extensively employed for the modeling of different types of activity/property and toxicity endpoints. This chapter provides explicit information regarding the basis, algorithm, and applicability of the ETA indices for a predictive modeling paradigm.

QSTR with extended topochemical atom (ETA) indices. 15. Development of predictive models for toxicity of organic chemicals against fathead minnow using second-generation ETA indices

Modern industrialisation has led to the production of millions of toxic chemicals having hazardous effects on the ecosystem. It is impracticable to determine the toxic potential of a large number of chemicals in animal models, making the use of quantitative structure-toxicity relationship (QSTR) models an alternative strategy for toxicity prediction. Recently we introduced a set of second-generation extended topochemical atom (ETA) indices for predictive modelling. Here we have developed predictive toxicity models on a large dataset of 459 diverse chemicals against fathead minnow (Pimephales promelas) using the second-generation ETA indices. These descriptors can be easily calculated from two-dimensional molecular representation without the need of time-consuming conformational analysis and alignment, making the developed models easily reproducible. Considering the importance of hydrophobicity for toxicity prediction, AlogP98 was used as an additional predictor in all the models, which were validated rigorously using multiple strategies. The ETA models were comparable in predictability to those involving various non-ETA topological parameters and those previously reported using various descriptors including computationally demanding quantum-chemical ones.

On some novel extended topochemical atom (ETA) parameters for effective encoding of chemical information and modelling of fundamental physicochemical properties

Extended topochemical atom (ETA) indices developed by our group have been extensively applied in our previous reports for toxicity and ecotoxicity modelling in the field of quantitative structure-activity relationships (QSARs). In the present study these indices have been further explored by defining additional novel parameters to model n-octanol-water partition coefficient (two data sets; n = 168 and 139), water solubility (n = 193), molar refractivity (n = 166), and aromatic substituent constants π, MR, σ (m), and σ (p) (n = 99). All the models developed in the present study have undergone rigorous internal and external validation tests and the models have high statistical significance and prediction potential. In terms of Q² and r² values the models developed for the datasets of whole molecules are better than those previously reported, with topochemically arrived unique (TAU) indices on the same datasets of chemicals. An attempt has also been made to develop models using non-ETA topological and information indices. Interestingly, ETA and non-ETA models have been found to have similar predictive capacity.

QSTR with extended topochemical atom (ETA) indices. 14. QSAR modeling of toxicity of aromatic aldehydes to Tetrahymena pyriformis

Aldehydes are a toxic class of chemicals causing severe health hazards. In this background, quantitative structure-toxicity relationship (QSTR) models have been developed in the present study using Extended Topochemical Atom (ETA) indices for a large group of 77 aromatic aldehydes for their acute toxicity against the protozoan ciliate Tetrahymena pyriformis. The ETA models have been compared with those developed using various non-ETA topological indices. Attempt was also made to include the n-octanol/water partition coefficient (logK(o/w)) as an additional descriptor considering the importance of hydrophobicity in toxicity prediction. Thirty different models were developed using different chemometric tools. All the models have been validated using internal validation and external validation techniques. The statistical quality of the ETA models was found to be comparable to that of the non-ETA models. The ETA models have shown the important effects of steric bulk, lipophilicity, presence of electronegative atom containing substituents and functionality of the aldehydic oxygen to the toxicity of the aldehydes. The best ETA model (without using logK(o/w)) shows encouraging statistical quality (Q(int)(2)=0.709,Q(ext)(2)=0.744). It is interesting to note that some of the topological models reported here are better in statistical quality than previously reported models using quantum chemical descriptors.

QSTR with extended topochemical atom (ETA) indices. 12. QSAR for the toxicity of diverse aromatic compounds to Tetrahymena pyriformis using chemometric tools

We have developed QSTR models for the toxicity of 384 diverse aromatic compounds to Tetrahymena pyriformis with recently introduced extended topochemical atom (ETA) indices and compared the ETA models with those derived from various non-ETA topological descriptors and also combined set of descriptors encompassing the ETA and non-ETA parameters. The data set was split into test (25% compounds of total data points) and training (remaining 75%) sets based on K-mean clustering technique. Different statistical analyses (factor analysis followed by multiple linear regression (FA-MLR), stepwise regression and partial least squares (PLS)) were performed with the training set compounds to develop QSTR models using the topological descriptors. All the developed models were cross-validated using leave-one-out (LOO) technique. The best models were selected on the basis of predicted R(2) values for test set compounds. The best models (based on external validation) developed from different techniques came from the combined set of descriptors. The above results indicate that the use of ETA descriptors with non-ETA descriptors improved the statistical quality of the non-ETA models. From the best models involving ETA parameters, it is observed that functionality of halogen atoms (hydrophobicity), volume parameter (bulk) and nitrogen containing functionalities (polarity) are important for developing QSTR models for the current data set. This study suggests that ETA parameters are sufficient power to encode chemical information contributing significantly to the toxicity of diverse aromatic compounds to T. pyriformis.

QSPR of the bioconcentration factors of non-ionic organic compounds in fish using extended topochemical atom (ETA) indices

Bioconcentration refers to the absorption or uptake of a chemical from the media to an organism's tissues leading to greater concentration in tissues than that in the surrounding environment. Considering the importance of bioconcentration from the viewpoint of ecological safety assessment, a QSPR study was conducted based upon log BCF of 122 non-ionic organic compounds in fish using the recently introduced extended topochemical atom (ETA) indices. In deriving the models, principal component factor analysis (FA) followed by multiple linear regression (MLR), stepwise regression, partial least squares (PLS) and principal component regression analysis (PCRA) were applied as statistical tools. This was repeated with non-ETA (topological and physicochemical) descriptors and a combination set including both the ETA and non-ETA descriptors. The ETA indices suggested negative contributions of functionalities of nitro, amino and hydroxy substructures and positive contributions of branching, volume and functionality of chloro substituents. Again, the predictive ability of the developed models was compared with the previously reported models. Finally the validation of all the QSAR models was discussed based on random division, sorted log BCF data and K-means clusters for the factor scores of the original variable (ETA) matrix without the response property values. The results suggest that ETA parameters are sufficiently rich in chemical information to encode the structural features contributing to the bioconcentration of the non-ionic organic compounds in fish and thus these merit further assessment to explore their potential in QSAR/QSPR/QSTR modelling.