Machine learning - based q-RASAR modeling to predict acute contact toxicity of binary organic pesticide mixtures in honey bees

The rat acute oral toxicity of trifluoromethyl compounds (TFMs): a computational toxicology study combining the 2D-QSTR, read-across and consensus modeling methods

All areas of the modern society are affected by fluorine chemistry. In particular, fluorine plays an important role in medical, pharmaceutical and agrochemical sciences. Amongst various fluoro-organic compounds, trifluoromethyl (CF3) group is valuable in applications such as pharmaceuticals, agrochemicals and industrial chemicals. In the present study, following the strict OECD modelling principles, a quantitative structure-toxicity relationship (QSTR) modelling for the rat acute oral toxicity of trifluoromethyl compounds (TFMs) was established by genetic algorithm-multiple linear regression (GA-MLR) approach. All developed models were evaluated by various state-of-the-art validation metrics and the OECD principles. The best QSTR model included nine easily interpretable 2D molecular descriptors with clear physical and chemical significance. The mechanistic interpretation showed that the atom-type electro-topological state indices, molecular connectivity, ionization potential, lipophilicity and some autocorrelation coefficients are the main factors contributing to the acute oral toxicity of TFMs against rats. To validate that the selected 2D descriptors can effectively characterize the toxicity, we performed the chemical read-across analysis. We also compared the best QSTR model with public OPERA tool to demonstrate the reliability of the predictions. To further improve the prediction range of the QSTR model, we performed the consensus modelling. Finally, the optimum QSTR model was utilized to predict a true external set containing many untested/unknown TFMs for the first time. Overall, the developed model contributes to a more comprehensive safety assessment approach for novel CF3-containing pharmaceuticals or chemicals, reducing unnecessary chemical synthesis whilst saving the development cost of new drugs.

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.

Molecular docking and MD simulation studies of 4-thiazol-N-(pyridin-2-yl)pyrimidin-2-amine derivatives as novel inhibitors targeted to CDK2/4/6

PURPOSE

Nowadays, cyclin-dependent kinase 4/6 (CDK4/6) inhibitors have been approved for treating metastatic breast cancer and have achieved inspiring curative effects. But some discoveries have indicated that CDK 4/6 are not the requisite factors in some cell types because CDK2 partly compensates for the inhibition of CDK4/6. Thus, it is urgent to design CDK2/4/6 inhibitors for significantly enhancing their potency. This study aims to explore the mechanism of the binding of CDK2/4/6 kinases and their inhibitors to design novel CDK2/4/6 inhibitors for significantly enhancing their potency in different kinds of cancers.

MATERIALS AND METHODS

A series of 72 disparately functionalized 4-substituted N-phenylpyrimidin-2-amine derivatives exhibiting potent inhibitor activities against CDK2, CDK4 and CDK6 were collected to apply to this research. The total set of these derivatives was divided into a training set (54 compounds) and a test set (18 compounds). The derivatives were constructed through the sketch molecule module in SYBYL 6.9 software. A Powell gradient algorithm and Tripos force field were used to calculate the minimal structural energy and the minimized structure was used as the initial conformation for molecular docking. By the means of 3D-QSAR models, partial least squares (PLS) analysis, molecular dynamics (MD) simulations and binding free energy calculations, we can find the relationship between structure and biological activity.

RESULTS

In this study, we used molecular docking, 3D-QSAR and molecular dynamics simulation methods to comprehensively analyze the interaction and structure-activity relationships of 72 new CDK2/4/6 inhibitors. We used detailed statistical data to reasonably verify the constructed 3D-QSAR models for three receptors (q2 of CDK2 = 0.714, R2pred = 0.764, q2 = 0.815; R2pred of CDK4 = 0.681, q2 = 0.757; R2pred of CDK6 = 0.674). MD simulations and decomposition energy analysis validated the reasonability of the docking results and identified polar interactions as crucial factors that influence the different bioactivities of the studied inhibitors of CDK2/4/6 receptors, especially the electrostatic interactions of Lys33/35/43 and Asp145/158/163. The nonpolar interaction with Ile10/12/19 was also critical for the differing potencies of the CDK2/4/6 inhibitors. We concluded that the following probably enhanced the bioactivity against CDK2/4/6 kinases: (1) electronegative groups at the N1-position and electropositive and moderate-sized groups at ring E; (2) electrogroups featured at R2; (3) carbon atoms at the X-position or ring C replaced by a benzene ring; and (4) an electrogroup as R4.

CONCLUSION

Previous studies, to our knowledge, only utilized a single approach of 3D-QSAR and did not integrate this method with other sophisticated techniques such as molecular dynamics simulations to discover new potential inhibitors of CDK2, CDK4, or CDK6. So we applied the intergenerational technology, such as 3D-QSAR technology, molecular docking simulation techniques, molecular dynamics simulations and MMPBSA19/MMGBSA20-binding free energy calculations to statistically explore the correlations between the structure with biological activities. The constructed 3D-QSAR models of the three receptors were reasonable and confirmed by the excellent statistical data. We hope the results obtained from this work will provide some useful references for the development of novel CDK2/4/6 inhibitors.

Combined interaction of fungicides binary mixtures: experimental study and machine learning-driven QSAR modeling

Fungicide mixtures are an effective strategy in delaying the development of fungicide resistance. In this research, a fixed ratio ray design method was used to generate fifty binary mixtures of five fungicides with diverse modes of action. The interaction of these mixtures was then analyzed using CA and IA models. QSAR modeling was conducted to assess their fungicidal activity through multiple linear regression (MLR), support vector machine (SVM), and artificial neural network (ANN). Most mixtures exhibited additive interaction, with the CA model proving more accurate than the IA model in predicting fungicidal activity. The MLR model showed a good linear correlation between selected theoretical descriptors by the genetic algorithm and fungicidal activity. However, both ML-based models demonstrated better predictive performance than the MLR model. The ANN model showed slightly better predictability than the SVM model, with R2 and R2cv at 0.91 and 0.81, respectively. For external validation, the R2test value was 0.845. In contrast, the SVM model had values of 0.91, 0.78, and 0.77 for the same metrics. In conclusion, the proposed ML-based model can be a valuable tool for developing potent fungicidal mixtures to delay fungicidal resistance emergence.

Toxicity of individual and mixture of organic compounds to P. Phosphoreum and S. Capricornutum using interpretable simple structural parameters

Human and environmental ecosystem beings are exposed to multicomponent compound mixtures but the toxicity nature of compound mixtures is not alike to the individual chemicals. This work introduces four models for the prediction of the negative logarithm of median effective concentration (pEC50) of individual chemicals to marine bacteria Photobacterium Phosphoreum (P. Phosphoreum) and algal test species Selenastrum Capricornutum (S. Capricornutum) as well as their mixtures to P. Phosphoreum, and S. Capricornutum. These models provide the simplest approaches for the forecast of pEC50 of some classes of organic compounds from their interpretable structural parameters. Due to the lack of adequate toxicity data for chemical mixtures, the largest available experimental data of individual chemicals (55 data) and their mixtures (99 data) are used to derive the new correlations. The models of individual chemicals are based on two simple structural parameters but chemical mixture models require further interaction terms. The new model's results are compared with the outputs of the best accessible quantitative structure-activity relationships (QSARs) models. Various statistical parameters are done on the new and comparative complex QSAR models, which confirm the higher reliability and simplicity of the new correlations.

Systematic approaches to machine learning models for predicting pesticide toxicity

Pesticides play an important role in modern agriculture by protecting crops from pests and diseases. However, the negative consequences of pesticides, such as environmental contamination and adverse effects on human and ecological health, underscore the importance of accurate toxicity predictions. To address this issue, artificial intelligence models have emerged as valuable methods for predicting the toxicity of organic compounds. In this review article, we explore the application of machine learning (ML) for pesticide toxicity prediction. This review provides a detailed summary of recent developments, prediction models, and datasets used for pesticide toxicity prediction. In this analysis, we compared the results of several algorithms that predict the harmfulness of various classes of pesticides. Furthermore, this review article identified emerging trends and areas for future direction, showcasing the transformative potential of machine learning in promoting safer pesticide usage and sustainable agriculture.

Prioritization of the ecotoxicological hazard of PAHs towards aquatic species spanning three trophic levels using 2D-QSTR, read-across and machine learning-driven modelling approaches

Polycyclic aromatic hydrocarbons (PAHs) represent a common group of environmental pollutants that endanger various aquatic organisms via various pathways. To better prioritize the ecotoxicological hazard of PAHs to aquatic environment, we used 2D descriptors-based quantitative structure-toxicity relationship (QSTR) to assess the toxicity of PAHs toward six aquatic model organisms spanning three trophic levels. According to strict OECD guideline, six easily interpretable, transferable and reproducible 2D-QSTR models were constructed with high robustness and reliability. A mechanistic interpretation unveiled the key structural factors primarily responsible for controlling the aquatic ecotoxicity of PAHs. Furthermore, quantitative read-across and different machine learning approaches were employed to validate and optimize the modelling approach. Importantly, the optimum QSTR models were further applied for predicting the ecotoxicity of hundreds of untested/unknown PAHs gathered from Pesticide Properties Database (PPDB). Especially, we provided a priority list in terms of the toxicity of unknown PAHs to six aquatic species, along with the corresponding mechanistic interpretation. In summary, the models can serve as valuable tools for aquatic risk assessment and prioritization of untested or completely new PAHs chemicals, providing essential guidance for formulating regulatory policies.

Integrated predictive QSAR, Read Across, and q-RASAR analysis for diverse agrochemical phytotoxicity in oat and corn: A consensus-based approach for risk assessment and prioritization

In the modern fast-paced lifestyle, time-efficient and nutritionally rich foods like corn and oat have gained popularity for their amino acids and antioxidant contents. The increasing demand for these cereals necessitates higher production which leads to dependency on agrochemicals, which can pose health risks through residual present in the plant products. To first report the phytotoxicity for corn and oat, our study employs QSAR, quantitative Read-Across and quantitative RASAR (q-RASAR). All developed QSAR and q-RASAR models were equally robust (R2 = 0.680-0.762, Q2Loo = 0.593-0.693, Q2F1 = 0.680-0.860) and find their superiority in either oat or corn model, respectively, based on MAE criteria. AD and PRI had been performed which confirm the reliability and predictability of the models. The mechanistic interpretation reveals that the symmetrical arrangement of electronegative atoms and polar groups directly influences the toxicity of compounds. The final phytotoxicity and prioritization are performed by the consensus approach which results into selection of 15 most toxic compounds for both species.

Unveiling first report on in silico modeling of aquatic toxicity of organic chemicals to Labeo rohita (Rohu) employing QSAR and q-RASAR

Labeo rohita, a fish species within the Carp family, holds significant dietary and aquacultural importance in South Asian countries. However, the habitats of L. rohita often face exposure to various harmful pesticides and organic compounds originating from industrial and agricultural runoff. It is challenging to individually investigate the effects of each potentially harmful compound. In such cases, in silico techniques like Quantitative Structure-Activity Relationship (QSAR) and quantitative Read-Across Structure-Activity Relationship (q-RASAR) can be employed to construct algorithmic models capable of simultaneously assessing the toxicity of numerous compounds. We utilized the US EPA's ToxValDB database to curate data regarding acute median lethal concentration (LC50) toxicity for L. rohita. The experimental variables included study type (mortality), study duration (ranging from 0.25 h to 4 h), exposure route (static, flowthrough, and renewal), exposure method (drinking water), and types of chemicals (industrial chemicals and pharmaceuticals). Using this dataset, we developed regression-based QSAR and q-RASAR models to predict chemical toxicity to L. rohita based on chemical descriptors. The key descriptors for predicting the toxicity of L. rohita in the regression-based QSAR model include F05[S-Cl], SpMax_EA(ri), s4_relPathLength_2, and SpDiam_AEA(ed). These descriptors can be employed to estimate the toxicity of untested compounds and aid in the development of compounds with lower toxicity based on the presence or absence of these descriptors. Both the QSAR and q-RASAR models serve as valuable tools for understanding the chemicals' structural features responsible for toxicity and for filling gaps in aquatic toxicity data by predicting the toxicity of newly untested compounds in relation to L. rohita. Finally, the developed best model was employed to predict 297 external chemicals, the most toxic substances to L. rohita were identified as cyhalothrin, isobornyl thiocyanatoacetate, and paclobutrzol, while the least toxic ones included ethyl acetate, ethylthiourea, and n-butyric acid.

Development of a read-across-derived classification model for the predictions of mutagenicity data and its comparison with traditional QSAR models and expert systems

Mutagenicity is considered an important endpoint from the regulatory, environmental and medical points of view. Due to the wide number of compounds that may be of concern and the enormous expenses (in terms of time, money, and animals) associated with rodent mutagenicity bioassays, this endpoint is a major target for the development of alternative approaches for screening and prediction. The majority of old-aged expert systems and quantitative structure-activity relationship (QSAR) models may show reduced performance over time for their application on newer chemical candidates; thus, researchers constantly try to improve the modeling strategies. In our report, we initially performed traditional classification-based linear discriminant analysis (LDA) QSAR modeling using the benchmark Ames dataset of diverse chemicals (6512 compounds) to recognize the relationship between the molecules and their potential mutagenic behavior. The classical LDA QSAR model is developed from a selected set of 2D descriptors. The LDA QSAR model was developed by using a total of 31 descriptors identified from the analysis of the most discriminating features. Additionally, we have used similarity-derived features obtained from the read-across (RA) to develop an RA-based QSAR model. The developed RA-based LDA QSAR model has better predictivity, transferability, and interpretability compared to the LDA QSAR model, and it uses a very small number of descriptors compared to the classical QSAR model. Different machine learning (ML) models were also developed using the descriptors appearing in the read-across-based LDA QSAR model for comparative studies. We have checked the prediction quality of 216 true external set compounds using the novel similarity-derived RA model. The performance of the OECD toolbox is also compared with the RA-derived LDA QSAR model for a true external set. The current study aimed to explore the significance of the read-across-based algorithm and its application to the most current experimental mutagenicity data to complement already available expert systems.

Quantitative Read-across structure-activity relationship (q-RASAR): A new approach methodology to model aquatic toxicity of organic pesticides against different fish species

We have developed quantitative toxicity prediction models for organic pesticides of agricultural importance considering different fish species using a novel quantitative Read-across structure-activity relationship (q-RASAR) approach. The current study uses experimental (Log 1/LC50) data of organic pesticides to various fish species, including Rainbow trout (RT: Oncorhynchus mykiss: 715 data points), Lepomis (LP: Lepomis macrochirus: 136 data points), and Miscellaneous (Pimephales promelas, Brachydanio rerio: 226 data points). This study has also discussed the validation of the developed models and the analysis of structural features that are important for aquatic toxicity towards fishes. The read-across-derived similarity, error, and concordance measures (RASAR descriptors) have been extracted from the preliminary 0D-2D descriptors; the combined pool of RASAR and selected 0D-2D descriptors have been used to develop the final models by employing partial least squares algorithm. All the q-RASAR models are acceptable in terms of goodness of fit, robustness, and external predictivity, superseding the quality of the respective QSAR models, as seen from the computed validation metrics. The q-RASAR is an effective approach that has the potential to be used as a good alternative way to enhance external predictivity, interpretability, and transferability for aquatic toxicity prediction as well as ecotoxicity potential identification.

Read-across-based intelligent learning: development of a global q-RASAR model for the efficient quantitative predictions of skin sensitization potential of diverse organic chemicals

Environmental chemicals and contaminants cause a wide array of harmful implications to terrestrial and aquatic life which ranges from skin sensitization to acute oral toxicity. The current study aims to assess the quantitative skin sensitization potential of a large set of industrial and environmental chemicals acting through different mechanisms using the novel quantitative Read-Across Structure-Activity Relationship (q-RASAR) approach. Based on the identified important set of structural and physicochemical features, Read-Across-based hyperparameters were optimized using the training set compounds followed by the calculation of similarity and error-based RASAR descriptors. Data fusion, further feature selection, and removal of prediction confidence outliers were performed to generate a partial least squares (PLS) q-RASAR model, followed by the application of various Machine Learning (ML) tools to check the quality of predictions. The PLS model was found to be the best among different models. A simple user-friendly Java-based software tool was developed based on the PLS model, which efficiently predicts the toxicity value(s) of query compound(s) along with their status of Applicability Domain (AD) in terms of leverage values. This model has been developed using structurally diverse compounds and is expected to predict efficiently and quantitatively the skin sensitization potential of environmental chemicals to estimate their occupational and health hazards.