Development of a QSAR for worst case estimates of acute toxicity of chemically reactive compounds

Ecotoxicological Evaluation of Bisphenol A and Alternatives: A Comprehensive In Silico Modelling Approach

Bisphenol A (BPA), a compound widely used in industrial applications, has raised concerns due to its environmental impact. As a key component in the manufacture of polycarbonate plastics and epoxy resins used in many consumer products, concerns about potential harm to human health and the environment are unavoidable. This study seeks to address these concerns by evaluating a range of potential BPA alternatives, focusing on their ecotoxicological properties. The research examines 76 bisphenols, including BPA derivatives, using a variety of in silico ecotoxicological models, although it should be noted that these models were not developed exclusively for this particular class of compounds. Consequently, interpretations should be made with caution. The results of this study highlight specific compounds of potential environmental concern and underscore the need to develop more specific models for BPA alternatives that will allow for more accurate and reliable assessment.

Direct Prediction of Physicochemical Properties and Toxicities of Chemicals from Analytical Descriptors by GC-MS

With advances in machine learning (ML) techniques, the quantitative structure–activity relationship (QSAR) approach is becoming popular for evaluating chemicals. However, the QSAR approach requires that the chemical structure of the target compound is known and that it should be convertible to molecular descriptors. These requirements lead to limitations in predicting the properties and toxicities of chemicals distributed in the environment as in the PubChem database; the structural information on only 14% of compounds is available. This study proposes a new ML-based QSAR approach that can predict the properties and toxicities of compounds using analytical descriptors of mass spectrum and retention index obtained via gas chromatography–mass spectrometry without requiring exact structural information. The model was developed based on the XGBoost ML method. The root-mean-square errors (RMSEs) for log K_o-w, log (molecular weight), melting point, boiling point, log (vapor pressure), log (water solubility), log (LD₅₀) (rat, oral), and log (LD₅₀) (mouse, oral) are 0.97, 0.052, 51, 23, 0.74, 1.1, 0.74, and 0.6, respectively. The model performed well on a chemical standard mixture measurement, with similar results to those of model validation. It also performed well on a measurement of contaminated oil with spectral deconvolution. These results indicate that the model is suitable for investigating unknown-structured chemicals detected in measurements. Any online user can execute the model through a web application named Detective-QSAR (http://www.mixture-platform.net/Detective_QSAR_Med_Open/). The analytical descriptor-based approach is expected to create new opportunities for the evaluation of unknown chemicals around us.

Environmental risk assessments of multiclass pharmaceutical active compounds: selection of high priority concern pharmaceuticals using entropy-utility functions

This research aimed to identify high-risk pharmaceutically active compounds (PhACs) by analyzing occurrence (O), persistence (P), bioaccumulation (B), and toxicity (T) of 62 drugs which are widely used in Iran. A comprehensive approach was taken in risk assessment of the selected PhACs and in their prioritization using multiple-criteria decision analysis (MCDA) such as utility functions and principal component analysis (PCA). In practice, assigning weight to each criterion (i.e., O, P, B, and T) for risk assessment of PhACs is a challenge. In this research, the impact of giving both equal and unequal weight to each criterion by using a quantitative entropy method was studied. For risk assessment, two exposure approaches (consumption rate and occurrence of PhACs) and three MCDA approaches (PCA and utility functions with and without equal weights for each criterion) were compared. The utility function using equal weights for all O, P, B, and T criteria showed that thioridazine, pimozide, chlorpromazine, sertraline, clomipramine, and aripiprazole were at the highest level of risk, with concern score of 0.75, 0.75, 0.67, 0.58, 0.58, and 0.58, respectively. Unequal weight approach included additional compounds such as fluoxetine, citalopram, and methadone as a priority. All three MCDA approaches showed that sedatives and antidepressants were prevalent PhACs in the risk-based priority lists. However, the exposure-based approaches showed antibiotics and analgesics as the pharmaceutical of the highest priority. Overall, selection of the high priority concern pharmaceuticals depends on the prioritization approach employed. However, the utility function using unequal weights is a more conservative and effective approach for prioritization.

Prediction of Human Lethal Doses and Concentrations of MEIC Chemicals from Rodent LD50 Values: An Attempt to Make Some Reparation

The prediction of human toxicities from animal toxicity tests is often poor, and is now discouraged and in some cases banned, especially those involving the LD₅₀ test. However, there is a vast number of historical LD₅₀ data in both public and in-house repositories that are being put to little use. This study examined the correlations between human lethality (doses and concentrations) of 36 MEIC chemicals and the median values of a large number of mouse and rat LD₅₀ values obtained for four different routes of administration. The best correlations were found with mouse and rat intraperitoneal LD₅₀ values (r² = 0.838 and 0.810 for human lethal dose, and r² = 0.753 and 0.785 for human lethal concentration). The results show that excellent prediction of human lethal dose and concentration can be made, for this series of chemicals at least, by using uncurated rodent LD₅₀ values, thus offering some reparation for the millions of rodent lives sacrificed in LD₅₀ testing.

Ionic Liquids for Topical Delivery in Cancer

:

The unique properties of ionic liquids make them quite appealing for diverse applications, from “green” solvents (1st generation ionic liquids) to finely tuned materials (2nd generation ionic liquids). A decade ago, a 3rd generation of ionic liquids emerged which is focused on their prospective clinical applications, either as drugs per se or as adjuvants in drug formulations. In recent years, research focused on the use of ionic liquids for topical drug delivery has been increasing and holds great promise towards clinical application against skin cancers. This article highlights the growing relevance of ionic liquids in medicinal chemistry and pharmaceutical technology, which is opening new windows of opportunity.

Consensus models to predict oral rat acute toxicity and validation on a dataset coming from the industrial context

We report predictive models of acute oral systemic toxicity representing a follow-up of our previous work in the framework of the NICEATM project. It includes the update of original models through the addition of new data and an external validation of the models using a dataset relevant for the chemical industry context. A regression model for LD₅₀ and multi-class classification model for toxicity classes according to the Global Harmonized System categories were prepared. ISIDA descriptors were used to encode molecular structures. Machine learning algorithms included support vector machine (SVM), random forest (RF) and naïve Bayesian. Selected individual models were combined in consensus. The different datasets were compared using the generative topographic mapping approach. It appeared that the NICEATM datasets were lacking some relevant chemotypes for chemical industry. The new models trained on enlarged data sets have applicability domains (AD) sufficiently large to accommodate industrial compounds. The fraction of compounds inside the models' AD increased from 58% (NICEATM model) to 94% (new model). The increase of training sets improved models' prediction performance: RMSE values decreased from 0.56 to 0.47 and balanced accuracies increased from 0.69 to 0.71 for NICEATM and new models, respectively.

In silico toxicology protocols

The present publication surveys several applications of in silico (i.e., computational) toxicology approaches across different industries and institutions. It highlights the need to develop standardized protocols when conducting toxicity-related predictions. This contribution articulates the information needed for protocols to support in silico predictions for major toxicological endpoints of concern (e.g., genetic toxicity, carcinogenicity, acute toxicity, reproductive toxicity, developmental toxicity) across several industries and regulatory bodies. Such novel in silico toxicology (IST) protocols, when fully developed and implemented, will ensure in silico toxicological assessments are performed and evaluated in a consistent, reproducible, and well-documented manner across industries and regulatory bodies to support wider uptake and acceptance of the approaches. The development of IST protocols is an initiative developed through a collaboration among an international consortium to reflect the state-of-the-art in in silico toxicology for hazard identification and characterization. A general outline for describing the development of such protocols is included and it is based on in silico predictions and/or available experimental data for a defined series of relevant toxicological effects or mechanisms. The publication presents a novel approach for determining the reliability of in silico predictions alongside experimental data. In addition, we discuss how to determine the level of confidence in the assessment based on the relevance and reliability of the information.

A new computational model for the prediction of toxicity of phosphonate derivatives using QSPR

Structural and electronic properties of a series of 25 phosphonate derivatives were analyzed applying density functional theory, with the exchange-correlation functional PBEPBE in combination with the 6-311++G** basis set for all atoms. The chemical reactivity of these derivatives has been interpreted using quantum descriptors such as frontier molecular orbitals (HOMO, LUMO), Hirshfeld charges, molecular electrostatic potential, and the dual descriptor [[Formula: see text]]. These descriptors are directly related to experimental median lethal dose ([Formula: see text], expressed as its decimal logarithm [[Formula: see text]([Formula: see text]] through a multiple linear regression equation. The proposed model predicts the toxicity of phosphonates in function of the volume (V), the load of the most electronegative atom of the molecule (q), and the eigenvalue of the molecular orbital HOMO ([Formula: see text]. The obtained values in the internal validation of the model are: [Formula: see text]%, [Formula: see text]%, [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], and [Formula: see text]%. The toxicity of nine phosphonate derivatives used as test molecules was adequately predicted by the model. The theoretical results indicate that the oxygen atom of the O=P group plays an important role in the interaction mechanism between the phosphonate and the acetylcholinesterase enzyme, inhibiting the removal of the proton of the ser-200 residue by the his-440 residue.

Development of thresholds of excess toxicity for environmental species and their application to identification of modes of acute toxic action

The acute toxicity of organic pollutants to fish, Daphnia magna, Tetrahymena pyriformis, and Vibrio fischeri was investigated. The results indicated that the Toxicity Ratio (TR) threshold of log TR =1, which has been based on the distribution of toxicity data to fish, can also be used to discriminate reactive or specifically acting compounds from baseline narcotics for Daphnia magna and Vibrio fischeri. A log TR=0.84 is proposed for Tetrahymena pyriformis following investigation of the relationships between the species sensitivity and the absolute averaged residuals (AAR) between the predicted baseline toxicity and the experimental toxicity. Less inert compounds exhibit relatively higher toxicity to the lower species (Tetrahymena pyriformis and Vibrio fischeri) than the higher species (fish and Daphnia magna). A greater number of less inert compounds with log TR greater than the thresholds was observed for Tetrahymena pyriformis and Vibrio fischeri. This may be attributed to the hydrophilic compounds which may pass more easily through cell membranes than the skin or exoskeleton of organisms and have higher bioconcentration factors in the lower species, leading to higher toxicity. Most of classes of chemical associated with excess toxicity to one species also exhibited excess toxicity to other species, however, a few classes with excess toxicity to one species exhibiting narcotic toxicity to other species and thus may have different MOAs between species. Some ionizable compounds have log TR much lower than one because of the over-estimated log K_OW. The factors that influence the toxicity ratio calculated from baseline level are discussed in this paper.

Antitumor Activity of Ionic Liquids Based on Ampicillin

Significant antiproliferative effects against various tumor cell lines were observed with novel ampicillin salts as ionic liquids. The combination of anionic ampicillin with appropriate ammonium, imidazolium, phosphonium, and pyridinium cations yielded active pharmaceutical ingredient ionic liquids (API-ILs) that show potent antiproliferative activities against five different human cancer cell lines: T47D (breast), PC3 (prostate), HepG2 (liver), MG63 (osteosarcoma), and RKO (colon). Some API-ILs showed IC50 values between 5 and 42 nM, activities that stand in dramatic contrast to the negligible cytotoxic activity level shown by the ampicillin sodium salt. Moreover, very low cytotoxicity against two primary cell lines-skin (SF) and gingival fibroblasts (GF)-indicates that the majority of these API-ILs are nontoxic to normal human cell lines. The most promising combination of antitumor activity and low toxicity toward healthy cells was observed for the 1-hydroxyethyl-3-methylimidazolium-ampicillin pair ([C2 OHMIM][Amp]), making this the most suitable lead API-IL for future studies.

Towards global QSAR model building for acute toxicity: Munro database case study

A series of 436 Munro database chemicals were studied with respect to their corresponding experimental LD50 values to investigate the possibility of establishing a global QSAR model for acute toxicity. Dragon molecular descriptors were used for the QSAR model development and genetic algorithms were used to select descriptors better correlated with toxicity data. Toxic values were discretized in a qualitative class on the basis of the Globally Harmonized Scheme: the 436 chemicals were divided into 3 classes based on their experimental LD50 values: highly toxic, intermediate toxic and low to non-toxic. The k-nearest neighbor (k-NN) classification method was calibrated on 25 molecular descriptors and gave a non-error rate (NER) equal to 0.66 and 0.57 for internal and external prediction sets, respectively. Even if the classification performances are not optimal, the subsequent analysis of the selected descriptors and their relationship with toxicity levels constitute a step towards the development of a global QSAR model for acute toxicity.

In silico prediction of chemical acute oral toxicity using multi-classification methods

Chemical acute oral toxicity is an important end point in drug design and environmental risk assessment. However, it is difficult to determine by experiments, and in silico methods are hence developed as an alternative. In this study, a comprehensive data set containing 12, 204 diverse compounds with median lethal dose (LD₅₀) was compiled. These chemicals were classified into four categories, namely categories I, II, III and IV, based on the criterion of the U.S. Environmental Protection Agency (EPA). Then several multiclassification models were developed using five machine learning methods, including support vector machine (SVM), C4.5 decision tree (C4.5), random forest (RF), κ-nearest neighbor (kNN), and naïve Bayes (NB) algorithms, along with MACCS and FP4 fingerprints. One-against-one (OAO) and binary tree (BT) strategies were employed for SVM multiclassification. Performances were measured by two external validation sets containing 1678 and 375 chemicals, separately. The overall accuracy of the MACCS-SVM(OAO) model was 83.0% and 89.9% for external validation sets I and II, respectively, which showed reliable predictive accuracy for each class. In addition, some representative substructures responsible for acute oral toxicity were identified using information gain and substructure frequency analysis methods, which might be very helpful for further study to avoid the toxicity.

Ranking of concern, based on environmental indexes, for pharmaceutical and personal care products: an application to the Spanish case

A wide range of Pharmaceuticals and Personal Care Products (PPCPs) are present in the environment, and many of their adverse effects are unknown. The emergence of new compounds or changes in regulations have led to dynamical studies of occurrence, impact and treatment, which consider geographical areas and trends in consumption and innovation in the pharmaceutical industry. A Quantitative study of Structure-Activity Relationship ((Q)SAR) was performed to assess the possible adverse effects of ninety six PPCPs and metabolites with negligible experimental data and establish a ranking of concern, which was supported by the EPA EPI Suite™ interface. The environmental and toxicological indexes, the persistence (P), the bioaccumulation (B), the toxicity (T) (extensive) and the occurrence in Spanish aquatic environments (O) (intensive) were evaluated. The most hazardous characteristics in the largest number of compounds were generated by the P index, followed by the T and B indexes. A high number of metabolites has a concern score equal to or greater than their parent compounds. Three PBT and OPBT rankings of concern were proposed using the total and partial ranking method (supported by a Hasse diagram) by the Decision Analysis by Ranking Techniques (DART) tool, which was recently recommended by the European Commission. An analysis of the sensibility of the relative weights of these indexes has been conducted. Hormones, antidepressants (and their metabolites), blood lipid regulators and all of the personal care products considered in this study were at the highest levels of risk according to the PBT and OPBT total rankings. Furthermore, when the OPBT partial ranking was performed, X-ray contrast media, H2 blockers and some antibiotics were included at the highest level of concern. It is important to improve and incorporate useful indexes for the predicted environmental impact of PPCPs and metabolites and thus focus experimental analysis on the compounds that require urgent attention.

Accurate prediction of the toxicity of benzoic acid compounds in mice via oral without using any computer codes

Most of benzoic acid derivatives are toxic, which may cause serious public health and environmental problems. Two novel simple and reliable models are introduced for desk calculations of the toxicity of benzoic acid compounds in mice via oral LD(50) with more reliance on their answers as one could attach to the more complex outputs. They require only elemental composition and molecular fragments without using any computer codes. The first model is based on only the number of carbon and hydrogen atoms, which can be improved by several molecular fragments in the second model. For 57 benzoic compounds, where the computed results of quantitative structure-toxicity relationship (QSTR) were recently reported, the predicted results of two simple models of present method are more reliable than QSTR computations. The present simple method is also tested with further 324 benzoic acid compounds including complex molecular structures, which confirm good forecasting ability of the second model.

QSAR modelling toxicity toward rats of inorganic substances by means of CORAL

CORAL (‘CORrelation And Logic’) is freeware available on the Internet www.insilico.eu/coral The aim of this program is to establish a correlation between an endpoint and descriptors calculated with a simplified molecular input line entry system (SMILES). Three models calculated by CORAL for toxicity towards rat (-pLD50) of inorganic substances (three random splits) have shown that CORAL could be a good tool to model this endpoint. The average statistical characteristics for the CORAL models are the following: n=38, r2=0.8461, q2=0.8298, s=0.273, F=198 (subtraining set); n=37, r2=0.8144, s=0.322, F=154 (calibration set); and n=10, r2=0.8004, Rm (test)2 =0.7815, s=0.240, F=32 (validation set).

Estimation of reliability of predictions and model applicability domain evaluation in the analysis of acute toxicity (LD50)

This study presents a new type of acute toxicity (LD(50)) prediction that enables automated assessment of the reliability of predictions (which is synonymous with the assessment of the Model Applicability Domain as defined by the Organization for Economic Cooperation and Development). Analysis involved nearly 75,000 compounds from six animal systems (acute rat toxicity after oral and intraperitoneal administration; acute mouse toxicity after oral, intraperitoneal, intravenous, and subcutaneous administration). Fragmental Partial Least Squares (PLS) with 100 bootstraps yielded baseline predictions that were automatically corrected for non-linear effects in local chemical spaces--a combination called Global, Adjusted Locally According to Similarity (GALAS) modelling methodology. Each prediction obtained in this manner is provided with a reliability index value that depends on both compound's similarity to the training set (that accounts for similar trends in LD(50) variations within multiple bootstraps) and consistency of experimental results with regard to the baseline model in the local chemical environment. The actual performance of the Reliability Index (RI) was proven by its good (and uniform) correlations with Root Mean Square Error (RMSE) in all validation sets, thus providing quantitative assessment of the Model Applicability Domain. The obtained models can be used for compound screening in the early stages of drug development and prioritization for experimental in vitro testing or later in vivo animal acute toxicity studies.

Quantitative structure-activity relationship modeling of rat acute toxicity by oral exposure

Few quantitative structure-activity relationship (QSAR) studies have successfully modeled large, diverse rodent toxicity end points. In this study, a comprehensive data set of 7385 compounds with their most conservative lethal dose (LD(50)) values has been compiled. A combinatorial QSAR approach has been employed to develop robust and predictive models of acute toxicity in rats caused by oral exposure to chemicals. To enable fair comparison between the predictive power of models generated in this study versus a commercial toxicity predictor, TOPKAT (Toxicity Prediction by Komputer Assisted Technology), a modeling subset of the entire data set was selected that included all 3472 compounds used in TOPKAT's training set. The remaining 3913 compounds, which were not present in the TOPKAT training set, were used as the external validation set. QSAR models of five different types were developed for the modeling set. The prediction accuracy for the external validation set was estimated by determination coefficient R(2) of linear regression between actual and predicted LD(50) values. The use of the applicability domain threshold implemented in most models generally improved the external prediction accuracy but expectedly led to the decrease in chemical space coverage; depending on the applicability domain threshold, R(2) ranged from 0.24 to 0.70. Ultimately, several consensus models were developed by averaging the predicted LD(50) for every compound using all five models. The consensus models afforded higher prediction accuracy for the external validation data set with the higher coverage as compared to individual constituent models. The validated consensus LD(50) models developed in this study can be used as reliable computational predictors of in vivo acute toxicity.

A novel two-step hierarchical quantitative structure-activity relationship modeling work flow for predicting acute toxicity of chemicals in rodents

BACKGROUND

Accurate prediction of in vivo toxicity from in vitro testing is a challenging problem. Large public-private consortia have been formed with the goal of improving chemical safety assessment by the means of high-throughput screening.

OBJECTIVE

A wealth of available biological data requires new computational approaches to link chemical structure, in vitro data, and potential adverse health effects.

METHODS AND RESULTS

A database containing experimental cytotoxicity values for in vitro half-maximal inhibitory concentration (IC(50)) and in vivo rodent median lethal dose (LD(50)) for more than 300 chemicals was compiled by Zentralstelle zur Erfassung und Bewertung von Ersatz- und Ergaenzungsmethoden zum Tierversuch (ZEBET; National Center for Documentation and Evaluation of Alternative Methods to Animal Experiments). The application of conventional quantitative structure-activity relationship (QSAR) modeling approaches to predict mouse or rat acute LD(50) values from chemical descriptors of ZEBET compounds yielded no statistically significant models. The analysis of these data showed no significant correlation between IC(50) and LD(50). However, a linear IC(50) versus LD(50) correlation could be established for a fraction of compounds. To capitalize on this observation, we developed a novel two-step modeling approach as follows. First, all chemicals are partitioned into two groups based on the relationship between IC(50) and LD(50) values: One group comprises compounds with linear IC(50) versus LD(50) relationships, and another group comprises the remaining compounds. Second, we built conventional binary classification QSAR models to predict the group affiliation based on chemical descriptors only. Third, we developed k-nearest neighbor continuous QSAR models for each subclass to predict LD(50) values from chemical descriptors. All models were extensively validated using special protocols.

CONCLUSIONS

The novelty of this modeling approach is that it uses the relationships between in vivo and in vitro data only to inform the initial construction of the hierarchical two-step QSAR models. Models resulting from this approach employ chemical descriptors only for external prediction of acute rodent toxicity.