Predicting zebrafish (Danio rerio) embryo developmental toxicity through a non-conformational QSAR approach

From molecular descriptors to the developmental toxicity prediction of pesticides/veterinary drugs/bio-pesticides against zebrafish embryo: Dual computational toxicological approaches for prioritization

The escalating introduction of pesticides/veterinary drugs into the environment has necessitated a rapid evaluation of their potential risks to ecosystems and human health. The developmental toxicity of pesticides/veterinary drugs was less explored, and much less the large-scale predictions for untested pesticides, veterinary drugs and bio-pesticides. Alternative methods like quantitative structure-activity relationship (QSAR) are promising because their potential to ensure the sustainable and safe use of these chemicals. We collected 133 pesticides and veterinary drugs with half-maximal active concentration (AC50) as the zebrafish embryo developmental toxicity endpoint. The QSAR model development adhered to rigorous OECD principles, ensuring that the model possessed good internal robustness (R2 > 0.6 and QLOO2 > 0.6) and external predictivity (Rtest2 > 0.7, QFn2 >0.7, and CCCtest > 0.85). To further enhance the predictive performance of the model, a quantitative read-across structure-activity relationship (q-RASAR) model was established using the combined set of RASAR and 2D descriptors. Mechanistic interpretation revealed that dipole moment, the presence of C-O fragment at 10 topological distance, molecular size, lipophilicity, and Euclidean distance (ED)-based RA function were main factors influencing toxicity. For the first time, the established QSAR and q-RASAR models were combined to prioritize the developmental toxicity of a vast array of true external compounds (pesticides/veterinary drugs/bio-pesticides) lacking experimental values. The prediction reliability of each query molecule was evaluated by leverage approach and prediction reliability indicator. Overall, the dual computational toxicology models can inform decision-making and guide the design of new pesticides/veterinary drugs with improved safety profiles.

Application of machine learning in the study of development, behavior, nerve, and genotoxicity of zebrafish

Machine learning (ML) as a novel model-based approach has been used in studying aquatic toxicology in the environmental field. Zebrafish, as an ideal model organism in aquatic toxicology research, has been widely used to study the toxic effects of various pollutants. However, toxicity testing on organisms may cause significant harm, consume considerable time and resources, and raise ethical concerns. Therefore, ML is used in related research to reduce animal experiments and assist researchers in conducting toxicological research. Although ML techniques have matured in various fields, research on ML-based aquatic toxicology is still in its infancy due to the lack of comprehensive large-scale toxicity databases for environmental pollutants and model organisms. Therefore, to better understand the recent research progress of ML in studying the development, behavior, nerve, and genotoxicity of zebrafish, this review mainly focuses on using ML modeling to assess and predict the toxic effects of zebrafish exposure to different toxic chemicals. Meanwhile, the opportunities and challenges faced by ML in the field of toxicology were analyzed. Finally, suggestions and perspectives were proposed for the toxicity studies of ML on zebrafish in future applications.

Ensemble multiclassification model for predicting developmental toxicity in zebrafish

In recent years, with the rapid development of society, organic compounds have been released into aquatic environments in various forms, posing a significant threat to the survival of aquatic organisms. The assessment of developmental toxicity is an important part of environmental safety risk systems, helping to identify the potential impacts of organic compounds on the embryonic development of aquatic organisms and enabling early detection and warning of potential ecological risks. Additionally, binary classification models cannot accurately classify organic compounds. Therefore, it is crucial to construct a multiclassification model for predicting the developmental toxicity of organic compounds. In this study, binary and multiclassification models were developed based on the ToxCast™ Phase I chemical library and literature data. The random forest, support vector machine, extreme gradient boosting, adaptive gradient boosting, and C5.0 decision tree algorithms, as well as 8 types of molecular fingerprint were used to establish a multiclassification base model for predicting developmental toxicity through 5-fold cross-validation and external validation. Ultimately, a multiclassification ensemble model was derived through a voting method. The performance of the binary ensemble model, as measured by the balanced accuracy, was 0.918, while that of the multiclassification model was 0.819. The developmental toxicity voting ensemble model (DT-VEM) achieved accuracies of 0.804, 0.834, and 0.855. Furthermore, by utilizing the XGBoost machine learning algorithm to construct separate models for molecular descriptors and substructure molecular fingerprints, we identified several substructures and physical properties related to developmental toxicity. Our research contributes to a more detailed classification of developmental toxicity, providing a new and valuable tool for predicting the developmental toxicity effects of unknown compounds. This supplement addresses the limitations of previous tools, as it offers an enhanced ability to predict potential developmental toxicity in novel compounds.

A QSAR study for predicting malformation in zebrafish embryo

BACKGROUND

Developmental toxicity tests are extremely expensive, require a large number of animals, and are time-consuming. It is necessary to develop a new approach to simplify the analysis of developmental endpoints. One of these endpoints is malformation, and one group of ongoing methods for simplifying is in silico models. In this study, we aim to develop a quantitative structure-activity relationship (QSAR) model and identify the best algorithm for predicting malformations, as well as the most important and effective physicochemical properties associated with malformation.

METHODS

The dataset was extracted from a reliable database called COMPTOX. Physicochemical properties (descriptors) were calculated using Mordred and RDKit chemoinformatics software. The data were cleaned, preprocessed, and then split into training and testing sets. Machine learning algorithms, such as gradient boosting model (GBM) and logistic regression (LR), as well as deep learning models, including multilayer perceptron (MLP) and neural networks (NNs) trained with train set data and different sets of descriptors. The models were then validated with test set and various statistical parameters, such as Matthew's correlation coefficient (MCC) and balanced accuracy (BAC) score, were used to compare the models.

RESULTS

A set of descriptors containing with 78% AUC was identified as the best set of descriptors. Gradient boosting was determined to be the best algorithm with 78% predictive power.

CONCLUSIONS

The descriptors that were the most effective for developing models directly impact the mechanism of malformation, and GBM is the best model due to its MCC and BAC.

MCPNET: Development of an interpretable deep learning model based on multiple conformations of the compound for predicting developmental toxicity

The development of deep learning models for predicting toxicological endpoints has shown great promise, but one of the challenges in the field is the accuracy and interpretability of these models. The bioactive conformation of a compound plays a critical role for it to bind in the target. It is a big issue to figure out the bioactive conformation in deep learning without the co-crystal structure or highly precise molecular simulations. In this study, we developed a deep learning framework of Multi-Conformation Point Network (MCPNET) to construct classification and regression models, respectively, based on electrostatic potential distributions on vdW surfaces around multiple conformations of the compound using a dataset of compounds with developmental toxicity in zebrafish embryo. MCPNET applied 3D multi-conformational surface point cloud to extract the molecular features for model training, which may be critical for capturing the structural diversity of compounds. The models achieved an accuracy of 85 % on the classification task and R2 of 0.66 on the regression task, outperforming traditional machine learning models and other deep learning models. The key feature of our model is its interpretability with the component visualization to identify the factors contributing to the prediction and to understand the compound action mechanism. MCPNET may predict the conformation quietly close to the bioactive conformation of a compound by attention-based multi-conformation pooling mechanism. Our results demonstrated the potential of deep learning based on 3D molecular representations in accurately predicting developmental toxicity. The source code is publicly available at https://github.com/Superlit-CC/MCPNET.

The uses of zebrafish (Danio rerio) as an in vivo model for toxicological studies: A review based on bibliometrics

An in vivo model is necessary for toxicology. This review analyzed the uses of zebrafish (Danio rerio) in toxicology based on bibliometrics. Totally 56,816 publications about zebrafish from 2002 to 2023 were found in Web of Science Core Collection, with Toxicology as the top 6 among all disciplines. Accordingly, the bibliometric map reveals that "toxicity" has become a hot keyword. It further reveals that the most common exposure types include acute, chronic, and combined exposure. The toxicological effects include behavioral, intestinal, cardiovascular, hepatic, endocrine toxicity, neurotoxicity, immunotoxicity, genotoxicity, and reproductive and transgenerational toxicity. The mechanisms include oxidative stress, inflammation, autophagy, and dysbiosis of gut microbiota. The toxicants commonly evaluated by using zebrafish model include nanomaterials, arsenic, metals, bisphenol, and dioxin. Overall, zebrafish provide a unique and well-accepted model to investigate the toxicological effects and mechanisms. We also discussed the possible ways to address some of the limitations of zebrafish model, such as the combination of human organoids to avoid species differences.

A review on in silico prediction of the environmental risks posed by pharmaceutical emerging contaminants

Computer-aided (in silico) prediction has shown good potential to support the environmental risk assessment (ERA) of pharmaceutical emerging contaminants (PECs), allowing low-cost, animal-free, high-throughput screening of multiple potential risks posed by a wide variety of pharmaceuticals in the environment based on insufficient toxicity data. This review provided recent insights regarding the application of in silico approaches in prediction for environmental risks of PECs. Based on the review of 20 included articles from 8 countries published since 2018, we found that the researchers' interest and concern in this research topic were sharply aroused since 2021. Recently, in silico approaches have been widely used for the prediction of bioaccumulation and biodegradability, lethal endpoints, developmental toxicity, mutagenicity, other eco-toxicological effects such as ototoxicity and hematological toxicity, and human health hazards of exposure to PECs. Particular attention has been given to the simultaneous discernment of multiple environmental risks and health effects of PECs based on mechanistic data of pharmaceuticals using advanced bioinformatic methods such as transcriptomic analysis and network pharmacology prediction. In silico software platforms and databases used in the included studies were diversified, and there is currently no standardized and accepted in silico model for ERA of PECs. Date suggested that in silico prediction of the environmental risks posed by PECs is still in its infancy. Considerable critical challenges need to be addressed, including consideration of environmental exposure concentration for PECs, interactions among mixtures of PECs and other contaminants coexisting in environments, and development of in silico models specific to ERA of PECs.

Machine Learning-Based Hazard-Driven Prioritization of Features in Nontarget Screening of Environmental High-Resolution Mass Spectrometry Data

Nontarget high-resolution mass spectrometry screening (NTS HRMS/MS) can detect thousands of organic substances in environmental samples. However, new strategies are needed to focus time-intensive identification efforts on features with the highest potential to cause adverse effects instead of the most abundant ones. To address this challenge, we developed MLinvitroTox, a machine learning framework that uses molecular fingerprints derived from fragmentation spectra (MS2) for a rapid classification of thousands of unidentified HRMS/MS features as toxic/nontoxic based on nearly 400 target-specific and over 100 cytotoxic endpoints from ToxCast/Tox21. Model development results demonstrated that using customized molecular fingerprints and models, over a quarter of toxic endpoints and the majority of the associated mechanistic targets could be accurately predicted with sensitivities exceeding 0.95. Notably, SIRIUS molecular fingerprints and xboost (Extreme Gradient Boosting) models with SMOTE (Synthetic Minority Oversampling Technique) for handling data imbalance were a universally successful and robust modeling configuration. Validation of MLinvitroTox on MassBank spectra showed that toxicity could be predicted from molecular fingerprints derived from MS2 with an average balanced accuracy of 0.75. By applying MLinvitroTox to environmental HRMS/MS data, we confirmed the experimental results obtained with target analysis and narrowed the analytical focus from tens of thousands of detected signals to 783 features linked to potential toxicity, including 109 spectral matches and 30 compounds with confirmed toxic activity.

QSAR assessment of aquatic toxicity potential of diverse agrochemicals

The fast-increasing number of commercially produced chemicals challenges the experimental ecotoxicity assessment methods, which are costly, time-consuming, and dependent on the sacrifice of animals. In this regard, Quantitative Structure-Property/Activity Relationships (QSPR/QSAR) have led the way in developing ecotoxicity assessment models. In this study, QSAR models have been developed using the pEC50 values of 82 diverse agrochemicals or agro-molecules against a planktonic crustacean Daphnia magna with easily interpretable 2D descriptors. Moreover, a link among octanol-water partition coefficient (KOW), bio-concentration factor (BCF), and critical body residue (CBR) has been addressed, and their imputation for the prediction of the toxicity endpoint (EC50) has been done with an objective of the advanced exploration of several ecotoxicological parameters for toxic chemicals. The developed partial least squares (PLS) models were validated rigorously and proved to be robust, sound, and immensely well-predictive. The final Daphnia toxicity model derived from experimental derived properties along with computed descriptors emerged better in statistical quality and predictivity than those obtained solely from computed descriptors. Additionally, the pEC50 and other important properties (log KOW, log BCF, and log CBR) for a set of external agro-molecules, not employed in model development, were predicted to show the predictive ability of the models.

2D-Quantitative structure-activity relationship modeling for risk assessment of pharmacotherapy applied during pregnancy

The risk evaluation for pharmacological therapy during pregnancy is critical for maternal and fetal health. The initial risk assessment stage, the risk measurement, begins with pregnancy-labeling categories (A, B, C, D, and X) for pharmaceuticals defined by the US Food and Drug Administration (FDA). Recently, in silico methods have been preferred in toxicology studies to eliminate ethical issues before conducting clinical toxicology studies and animal experiments. Quantitative structure-activity relationship (QSAR) modeling is one of the in silico methodologies. The research focuses on creating a QSAR model that predicts the five FDA pregnancy categories of medications. Our dataset included 868 pharmaceuticals, containing nearly every pharmacological group collected from the FDA. 2D-molecular descriptors were calculated using PaDEL software. Twenty-four QSAR models were developed, and the best four models were discussed. The results of the models were compared according to sensitivity, accuracy, F-score, specificity, receiver operating characteristic (ROC) values, and Matthews correlation coefficient. Considering the statistical results, random forest is the best model for determining the pregnancy risk category of drugs. The accuracy of the model was 76.49% for internal and 93.58% for external validation. According to the kappa statistics, there is an average agreement of 0.583 for internal validation and a perfect agreement of 0.893 for external validation. Because the error rates of the model are very close to 0, the model is highly accurate. Consequently, our novel QSAR model gives guidance on the safe use of pharmaceuticals during pregnancy without requiring animal tests or clinical trials on pregnant women.

Application of ToxCast/Tox21 data for toxicity mechanism-based evaluation and prioritization of environmental chemicals: Perspective and limitations

In response to the need to minimize the use of experimental animals, new approach methodologies (NAMs) using advanced technology have emerged in the 21st century. ToxCast/Tox21 aims to evaluate the adverse effects of chemicals quickly and efficiently using a high-throughput screening and to transform the paradigm of toxicity assessment into mechanism-based toxicity prediction. The ToxCast/Tox21 database, which contains extensive data from over 1400 assays with numerous biological targets and activity data for over 9000 chemicals, can be used for various purposes in the field of chemical prioritization and toxicity prediction. In this study, an overview of the database was explored to aid mechanism-based chemical prioritization and toxicity prediction. Implications for the utilization of the ToxCast/Tox21 database in chemical prioritization and toxicity prediction were derived. The research trends in ToxCast/Tox21 assay data were reviewed in the context of toxicity mechanism identification, chemical priority, environmental monitoring, assay development, and toxicity prediction. Finally, the potential applications and limitations of using ToxCast/Tox21 assay data in chemical risk assessment were discussed. The analysis of the toxicity mechanism-based assays of ToxCast/Tox21 will help in chemical prioritization and regulatory applications without the use of laboratory animals.

Analysis of the Expression of Neurotrophins and Their Receptors in Adult Zebrafish Kidney

Neurotrophins and their receptors are involved in the development and maintenance of neuronal populations. Different reports have shown that all neurotrophin/receptor pathways can also play a role in several non-neuronal tissues in vertebrates, including the kidney. These signaling pathways are involved in different events to ensure the correct functioning of the kidney, such as growth, differentiation, and regulation of renal tubule transport. Previous studies in some fish species have identified the neurotrophins and receptors in the kidney. In this study, for the first time, we compare the expression profiles (mRNA and protein) of all neurotrophin/receptor pathways in the kidney of the adult zebrafish. We quantify the levels of mRNA by using qPCR and identify the expression pattern of each neurotrophin/receptor pathway by in situ hybridization. Next, we detect the proteins using Western blotting and immunohistochemistry. Our results show that among all neurotrophins analyzed, NT-3/TrkC is the most expressed in the glomerule and tubule and in the hematopoietic cells, similar to what has been reported in the mammalian kidney.