QSAR modeling with ETA indices for cytotoxicity and enzymatic activity of diverse chemicals

A novel approach to predict the comprehensive EROD potency: Mechanism-based curve fitting of CYP1A1 activity by PAHs

Cytochrome P450 1A1 (CYP1A1) plays critical roles in polycyclic aromatic hydrocarbon (PAH) toxicity, including DNA adduction and ROS generation. Therefore, CYP1A1 activity quantified by the 7-ethoxyresorufin-O-deethylase (EROD) assay (named EROD potency) has been considered a typical biomarker of PAH exposure and toxicity. The EROD dose-response curve always presents a biphasic style, increasing at low concentrations and decreasing at high concentrations of PAHs, but relative effect potency (REP) commonly used in PAH risk assessment is only involved in the increasing phase. In this study, a full bell-shaped EROD curve fitting formula Eq. (1) was obtained by considering both CYP1A1 mRNA induction and enzyme inhibition to completely assess the EROD potency of PAHs. Correspondingly, in silico models of QSAR and docking methods successfully predicted the full EROD curves of PAHs, and the structure-activity relationship indicated that PAHs with heavy molecular weight and large diameter showed stronger EROD potency. Further EROD potency with predicted curve parameters (EC_50,ind and area index) was confirmed by the reported REP (R² = 0.697-0.977) and experimental data from human and mouse cells (R² = 0.700-0.804). This study provides a novel curve fitting for the EROD dose-response relationship and a prediction model for PAH EROD potency.

Critical target identification and human health risk ranking of metal ions based on mechanism-driven modeling

Huge amounts of metals have been released into environment due to various anthropogenic activities, such as smelting and processing of metals and subsequent application in construction, automobiles, batteries, optoelectronic devices, and so on, resulting in widespread detection in environmental media. However, some metal ions are considered as "Environmental health hazards", leading to serious human health concerns through affecting critical targets. Hence, it is necessary to quickly and effectively recognize the key target of metal ions in living organisms. Fortunately, the development of high-throughput analysis and in silico approaches offer a promising tool for target identification. In this study, the key oncogenic target (tumor suppressor protein, p53) was screened by network analysis based on the comparative toxicogenomics database (CTD). Some metal ions could bind to p53 core domain, impair its function and induce the development of cancer risk, but its mechanisms were still unclear. Therefore, a quantitative structure-activity relationship (QSAR) model was constructed to characterize the binding constants (K_a) between DNA binding domain of p53 (p53 DBD) and nine metal ions (Mg²⁺, Ca²⁺, Cu²⁺, Zn²⁺, Co²⁺, Ni²⁺, Mn²⁺, Fe³⁺ and Ba²⁺). It had good robustness and predictive ability, which could be used to predict the K_a values of other six metal ions (Li⁺, Ag⁺, Cs⁺, Cd²⁺, Hg²⁺ and Pb²⁺) within application domain. The results showed strong binding affinity between Cd²⁺/Hg²⁺/Pb²⁺ and p53 DBD. Subsequent mechanism analyses revealed that first hydrolysis constant (|logK_OH|) and polarization force (Z²/r) were key metal ion-characteristic parameters. The metal ions with weak hydrolysis constants and strong polarization forces could readily interact with N-containing histidine and S-containing cysteine of p53 DBD, which resulted in high K_a values. This study identified p53 as potential target for metal ions, revealed the key characteristics affecting the actions and provide a basic understanding of metal ions-p53 DBD interaction.

Towards a better understanding of deep convolutional neural network processes for recognizing organic chemicals of environmental concern

Deep convolutional neural network (DCNN) has proved to be a promising tool for identifying organic chemicals of environmental concern. However, the uncertainty associated with DCNN predictions remains to be quantified. The training process contains many random configurations, including dataset segmentation, input sequences, and initial weight, etc. Moreover, the DCNN working mechanism is non-linear and opaque. To increase confidence to use this novel approach, persistent, bioaccumulative, and toxic substances (PBTs) were utilized as representative chemicals of environmental concern to estimate the prediction uncertainty under five distinguished datasets and ten different molecular descriptor (MD) arrangements with 111,852 chemicals and 2424 available MDs. An internal correlation coefficient test indicated that the prediction confidence reached 0.98 when a mean of 50 DCNNs' predictions was used instead of a sing DCNN prediction. A threshold for PBT categorization was determined by considering costs between false-negative and false-positive predictions. As revealed by the guided backpropagation-class activation mapping (GBP-CAM) saliency images, only 12% of all selected MDs were activated by DCNN and influenced decision-making process. However, the activated MDs not only varied among chemical classes but also shifted with different DCNNs. Principal component analysis indicated that 2424 MDs could transform into 370 orthogonal variables. Both results suggest that redundancy exists among selected MDs. Yet, DCNN was found to adapt to redundant data by focusing on the most important information for better prediction performance.

Discrimination of active and inactive substances in cytotoxicity based on Tox21 10K compound library: Structure alert and mode of action

In vitro cytotoxicity assay is an ideal alternative method for the in vivo toxicity in the risk assessment of pollutants in environment. However, modes of action (MOAs) of cytotoxicity have not been investigated for a wide range of compounds. In this paper, binomial and recursive partitioning analysis were carried out between the cytotoxicity and molecular descriptors for 8981 compounds. The results showed that cytotoxicity is strongly related to the chemical hydrophobicity and excess molar refraction, indicating the bio-uptake and chemical-receptor interaction through π and n electron pair play important roles in the cytotoxicity. The decision tree derived from recursive partitioning analysis revealed that the studied compounds could be divided into 25 groups and their structural characteristics could be used as structure alert to identify active and inactive compounds in cytotoxicity. The descriptors used in the decision tree revealed that chemical ionization and bioavailability could affect the cytotoxicity for ionizable and highly hydrophobic compounds. Comparison of MOAs based on Verhaar's classification scheme showed that many inert or less inert compounds were inactive substance, and many reactive or specifically-acting compounds were active substances in the cytotoxicity. In vitro toxicity assay instead of in vivo toxicity assay can be used in the environmental hazard and risk assessment of organic pollutants. The descriptors used in the binomial equation and decision tree reveal that chemical hydrophobicity, ionization and solubility play very important roles for identification of active and inactive compounds. The results obtained in this paper are valuable for understanding the modes of action in cytotoxicity and in vivo-in vitro toxicity relationship.

Review of the toxic effects of ionic liquids

Interest in ionic liquids (ILs), called green or designer solvents, has been increasing because of their excellent properties such as thermal stability and low vapor pressure; thus, they can replace harmful organic chemicals and help several industrial fields e.g., energy-storage materials production and biomaterial pretreatment. However, the claim that ILs are green solvents should be carefully considered from an environmental perspective. ILs, given their minimal vapor pressure, may not directly cause atmospheric pollution. However, they have the potential to cause adverse effects if leaked into the environment, for instance if they are spilled due to human mistakes or technical errors. To estimate the risks of ILs, numerous ILs have had their toxicity assessed toward several micro- and macro-organisms over the past few decades. Since the toxic effects of ILs depend on the method of estimating toxicity, it is necessary to briefly summarize and comprehensively discuss the biological effects of ILs according to their structure and toxicity testing levels. This can help simplify our understanding of the toxicity of ILs. Therefore, in this review, we discuss the key findings of toxicological information of ILs, collect some toxicity data of ILs to different species, and explain the influence of IL structure on their toxic properties. In the discussion, we estimated two different sensitivity values of toxicity testing levels depending on the experiment condition, which are theoretical magnitudes of the inherent sensitivity of toxicity testing levels in various conditions and their changes in biological response according to the change in IL structure. Finally, some perspectives, future research directions, and limitations to toxicological research of ILs, presented so far, are discussed.

Prediction of aquatic toxicity of chemical mixtures by the QSAR approach using 2D structural descriptors

The rapid industrialization has led to the generation of various organic chemicals and multi-component mixtures which affect the environment adversely. Although organic chemicals are often exposed to the environment as a form of chemical mixtures rather than individual compounds, there is insufficient toxicity data available for the chemical mixtures due to the associated complexities. Most importantly, the nature of toxicity of mixtures is completely different from the individual chemicals, which makes the evaluation more difficult and challenging. In this paper, we have developed QSAR models for various individual and mixture data sets for the prediction of the aquatic toxicity. We have used Partial Least Squares (PLS) regression as a statistical tool to build the models. The various structural features of the individual chemicals and the mixture components have been modeled against the toxicity end point pEC₅₀ (negative logarithm of median effective concentration in molar scale) of the aquatic organisms Photobacterium phosphoreum (marine bacterium) and Selenastrum capricornutum (freshwater algae). The mixture descriptors have been calculated by the weighted descriptor generation approach. The models were developed in accordance with OECD guidelines, and the quality of each model has been adjudged by strict validation parameters. The final models are robust, extremely predictive and interpretable mechanistically which can be used for the prediction of toxicity of untested chemical mixtures under the domain of applicability of the developed models.

QSPR modeling of octanol-water partition coefficient and organic carbon normalized sorption coefficient of diverse organic chemicals using Extended Topochemical Atom (ETA) indices

Octanol-water partition coefficient (logK_ow) and soil organic carbon content normalized sorption coefficient (logK_oc) values are two important physicochemical properties in the context of bioaccumulation and environmental fate of organic compounds and their environmental risk assessment. Simple, interpretable and easy-to-derive extended topochemical atom (ETA) indices obtained from 2D structural representation of compounds were used for quantitative structure-property relationship (QSPR) modeling of these two endpoints. Linear regression based models developed using only ETA indices show encouraging statistical and validation results. Based on the information obtained from developed QSPR models, we may conclude that molecular volume, branching pattern, presence of hydrophobic Cl atoms, cyclicity/fusion, polar environment, electron density, unsaturation content, hydrogen bonding propensity or hydrogen bond donor atoms, local topology, presence of heteroatoms and aromaticity are crucial factors in controlling the logK_ow and logK_oc values of the compounds. The suggested explanatory features for different classes of chemicals or the whole diverse set can help in safer designing of chemicals, which is one of the primary agenda of the "Green Chemistry" program.