Efficient analysis of the toxicity and mechanisms of Hexaconazole and two other triazole fungicides: insights from integrated network toxicology, molecular docking and bioinformatics data

Efficient analysis of the toxicity and mechanisms of Hexaconazole and two other triazole fungicides: insights from integrated network toxicology, molecular docking and bioinformatics data | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Efficient analysis of the toxicity and mechanisms of Hexaconazole and two other triazole fungicides: insights from integrated network toxicology, molecular docking and bioinformatics data Qing Lv, Xingke Zhu This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-5318182/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract Background: The safety of food grains is crucial for human health. Hexaconazole, Propiconazole, and Prothioconazole are widely used triazole fungicides primarily employed in agriculture for the control of plant diseases, aiming to enhance crop yield and quality. Objective: This research aims to elucidate the potential toxic targets and underlying molecular mechanisms of brain and liver damage induced by exposure to the three fungicides through network toxicology combined with molecular docking and bioinformatics data analysis. Method and results: Toxicity analysis through the ADMETlab database showed that the toxicity of the three fungicides was mainly associated with neurotoxicity and liver injury. Systematically utilizing GeneCards, OMIM and other databases and through Cytoscape tools, we identified potential and core targets (including EGFR , CASP3 , ESR1 , PPARG , TP53 , HSP90AA1 , and PTGS2 ) associated with fungicides and brain injury or liver injury. GO and KEGG enrichment analyses indicate that these targets are associated with pathways related to cancer, the PI3K-Akt signaling pathway, organic cyclic compounds, and organic nitrogen compounds. Molecular docking simulations conducted using AutoDock confirmed the high-affinity binding interactions between the fungicides and key target proteins. Conclusions: This study provides a theoretical foundation for understanding the molecular mechanisms underlying the neurotoxicity and hepatotoxicity induced by Hexaconazole, Propiconazole, and Prothioconazole, while establishing a foundational framework for the development of prevention and treatment strategies related to triazole pesticide-associated brain and liver injuries. Our findings underscore the potential risks these three pesticides pose to brain and liver health, highlighting the need for further epidemiological and clinical research in the future. Hexaconazole Propiconazole Prothioconazole Network toxicology Molecular docking Brain injury Liver injury Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Figure 8 1. Introduction Triazole fungicides represent a significant class of agricultural chemicals widely utilized for the control of various fungal pathogens affecting crops [1]. Characterized by their five-membered ring structure containing three nitrogen atoms, triazoles function primarily by inhibiting the biosynthesis of ergosterol, an essential component of fungal cell membranes [2]. Among the most commonly used triazole fungicides are Hexaconazole (HA), Propiconazole (PP) and Prothioconazole (PT) which exhibit broad-spectrum efficacy against numerous fungal diseases, including rusts, powdery mildews, and blights [c]. These compounds are particularly valued for their systemic properties, allowing them to be absorbed and translocated within plant tissues, thereby providing effective protection against both foliar and soil-borne pathogens [4]. Research on the toxicity of triazole fungicides has garnered increasing attention due to their widespread use in agriculture and potential environmental and health impacts. Triazoles, including well-known compounds such as Hexaconazole, Tebuconazole, and Propiconazole, have been demonstrated to possess varying degrees of toxicity to non-target organisms, including aquatic life, terrestrial invertebrates, and mammals [5–7].Furthermore, significant efforts have been made to evaluate the ecological risks associated with triazole fungicides. Toxicological assessments, including acute and chronic exposure studies, have revealed adverse effects on aquatic organisms such as fish and invertebrates, as well as terrestrial species like bees and earthworms [8, 9]. Recent studies have highlighted the neurotoxic potential of triazole compounds, linking exposure to cognitive impairments and neurodegenerative conditions [10–12]. Moreover, the involvement of triazoles in disrupting hepatic metabolic pathways raises concerns about their long-term effects on liver health [13, 14]. However, the precise biochemical pathways through which triazole fungicides induce neurotoxicity and hepatotoxicity remain poorly characterized. Therefore, the mechanism of brain injury and liver injury induced by triazole fungicides still needs further study. Advancements in toxicology have highlighted the importance of an integrated approach to elucidate the complex interactions between chemical agents and biological systems. Network toxicology offers a promising framework for this analysis, enabling researchers to examine the relationships between chemical exposure, biological pathways, and toxicological outcomes [15]. By leveraging bioinformatics tools, researchers can integrate large datasets, including gene expression profiles, protein interactions, and metabolic pathways, to identify key molecular targets and pathways affected by these fungicides. Molecular docking techniques further enhance our understanding of the interactions between triazole fungicides and their biological targets [16]. This computational method allows for the prediction of binding affinities and interaction modes of fungicides with specific proteins, providing insights into their mechanistic action at the molecular level [17]. By applying molecular docking alongside network toxicology and bioinformatics, we can gain a comprehensive understanding of how Hexaconazole and its counterparts exert their toxic effects. This study aims to efficiently analyze the toxicity and mechanisms of Hexaconazole, Propiconazole and Prothioconazole through an integrated approach that combines network toxicology, molecular docking, and bioinformatics data. By elucidating the key pathways and molecular interactions involved, this research will contribute to a better understanding of the risks associated with triazole fungicide use and inform future regulatory decisions and risk assessment strategies. Ultimately, our findings aim to promote safer agricultural practices and minimize the impact on human health of these widely used agrochemicals. 2. Materials and methods 2.1 Initial network toxicology predictive analysis of HA, PP and PT The terms “Hexaconazole (HA)”, “Propiconazole (PP)” and “Prothioconazole (PT)” were queried in the PubChem database ( https://pubchem.ncbi.nlm.nih.gov/ ) to obtain the standardized SMILES codes for these compounds. The SMILES codes of the three triazole compounds were then input into the ADMETlab 3.0 analysis tool ( https://admetlab3.scbdd.com/ ). Utilizing the structural modeling capabilities of this platform, we predicted the toxicological properties of HA, PP, and PT, providing an initial understanding of their toxicological profiles. The predicted results were downloaded. 2.2 Targets collection of HA, PP and PT First, a search was conducted in the PubChem database for “Hexaconazole”, “Propiconazole” and “Prothioconazole” to identify the best matches and verify the accuracy of their names and molecular formulas. Subsequently, the 2D structure SDF files were downloaded. The 2D structures are shown in Supplementary Fig. S1 , with CAS numbers 103556-63-4, 60207-90-1, and 178928-70-6, respectively. Based on the SDF files, potential targets of "Homo sapiens" specific HA, PP,PT were searched from ChEMBL ( https://www.ebi.ac.uk/chembl/ ) and Swiss Target Prediction database ( http://www.swisstargetprediction.ch/ ) [18]. In addition, the Swiss-Target-Prediction tool is used to identify any ignored targets. The structural information throughout the search results is carefully cross-verified for consistency. Then open the Excel table and select the target with probability > 0 as the target of HA, PP and PT. Finally, integration, de-duplication and gene naming coordination were carried out by UniProt database ( https://www.uniprot.org/ ) to establish the target library of HA, PP and PT. 2.3 Collection of targets related to brain and liver injury We conducted a comprehensive search of both domestic and international literature, as well as utilizing the GeneCards ( https://www.genecards.org/ ) and OMIM ( https://omim.org/ ) databases, using the keywords “brain injury” and “liver injury” to identify relevant targets. To ensure that the genes obtained were highly associated with brain and liver injuries, we set the “score” threshold at the median value, selecting genes with scores above the median to establish a library of neurotoxic targets. We considered the intersection of targets related to HA, PP, and PT with brain and liver injury targets as potential targets for pesticide-induced brain or liver damage. Venn diagrams were employed to filter the common potential targets. 2.4 Protein-protein interaction network construction and core targets screening The intersection genes representing potential targets for brain or liver injury induced by the three pesticides were individually input into the STRING database ( https://cn.string-db.org/ ). We limited the species to “Homo sapiens”, set the “minimum required interaction score (FDR)” to “medium confidence > 0.4” and selected “FDR strict value ≤ 0.05” for analysis. These parameters ensured a robust analysis of the active target proteins corresponding to our target genes. Subsequently, the results generated by STRING were imported into the network biology visualization application Cytoscape (version 3.10.0), which computes parameters for each node in the network diagram and displays molecular connections [19]. This allowed for the calculation of the topological properties of the network nodes and edges, resulting in a protein-protein interaction (PPI) network diagram. The criteria for screening core targets were as follows: nodes corresponding to targets that simultaneously met the conditions of betweenness centrality (BC) > median, closeness centrality (CC) > median, and degree centrality (DC) > median were selected as core targets for brain or liver injury induced by HA, PP and PT. Additionally, we utilized the MCODE plugin to validate the core genes within these PPI networks. 2.5 Functional enrichment and pathway analysis of targets for diseases induced by HA, PP and PT 2.5.1 GO and KEGG analysis of potential targets To analyze the enrichment of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, we imported the targets associated with brain or liver damage induced by HA, PP, and PT into the Metascape database ( https://metascape.org/gp/index.html ). This database serves as a comprehensive resource for gene functions, encompassing extensive information on biological processes (BP), cellular components (CC), and molecular functions (MF) [20]. The selection criteria included a minimum overlap of 3 and a minimum enrichment factor of 1.5, with a significance level set at P < 0.01. The overlapping targets were uploaded to the Metascape database, where human species were selected, and the results were submitted. The top 10 GO terms and the top 20 KEGG pathways were selected based on ascending p-values. GO and KEGG diagrams were generated using the Microbiotics website ( http://www.bioinformatics.com.cn/ ) and R software. The results of the bioinformatics analyses produced online KEGG pathway enrichment maps and GO term enrichment maps to visualize the interactions among pathways related to brain or liver damage induced by the three pesticides. 2.5.2 KEGG analysis of core targets We performed KEGG pathway enrichment analysis of the core targets associated with brain and liver damage induced by HA, PP, and PT using the Metascape database. The enrichment results were visualized to effectively interpret and present our findings. This multifaceted approach aims to comprehensively investigate the signaling pathways and biological processes mediated by the core targets related to brain and liver damage induced by these three pesticides, thereby elucidating and emphasizing the key mechanisms involved in this context. 2.6 Drug-target-pathway network of HA, PP and PT The drug-target-pathway network provides a clear overview of how drugs affect pathways through specific targets and illustrates the relationships between these key targets and pathways, allowing for the identification of the most critical targets and pathways. First, an attribute table was constructed in Excel, containing the core targets associated with the three triazole pesticides and all pathways corresponding to each core target. Next, Cytoscape was opened, and the attribute table was uploaded in the Network file to load the locations, set the starting points, and visualize the results. 2.7 Molecular docking of key targets with HA, PP and PT To investigate the molecular interactions and binding patterns between HA, PP, PT, and the identified key target proteins, we conducted molecular docking simulations. This theoretical structure-based approach enabled the predictive modeling of receptor-ligand binding geometries and affinities. The crystal structures of the core target proteins were obtained from the RCSB PDB database ( https://www.rcsb.org/ ) and preprocessed using PyMOL 3.0 ( https://www.pymol.org/ ) by removing solvent water molecules and any pre-existing ligands bound to the target proteins. The processed protein structures were then imported into AutoDock Tools 1.5.6 ( https://autodock.scripps.edu/ ) for hydrogen addition, charge calculation, and non-polar hydrogen merging to generate receptor models. The structures of HA, PP, and PT compounds were downloaded from the PubChem database. The mechanical structures were optimized using Chem3D ( https://revvitysignals.com/products/research/chemdraw ) to minimize energy and designated as ligands. The docking range was determined by adjusting the center coordinates and grid box size parameters. Following this, a genetic algorithm was selected, and molecular docking simulations were executed via command line using AutoDock Vina. Finally, the predicted binding modes and affinities were visualized and analyzed using PyMOL 3.0. 2.8 Collection and analysis of bioinformatics data 2.8.1 Identification and validation of brain and liver cancer-related DEGs To further elucidate the mechanistic impact of central targets associated with brain and liver damage induced by three pesticides on cancer, we obtained transcriptomic datasets from the NCBI GEO database ( https://www.ncbi.nlm.nih.gov/geo/ ), selecting GSE30563 and GSE117361 for detailed analysis. The GEOexplorer tool ( https://geoexplorer.rosalind.kcl.ac.uk/ ) was employed to determine differentially expressed genes (DEGs) between brain and liver tumors with normal tissue samples. Multichip analysis was conducted using the SVA and limma packages in R software ( https://www.r-project.org/ ), allowing for batch correction and calculation of fold-change (FC) and adjusted p -values. The thresholds for identifying DEGs were set at |log (FC)| > 1 and P < 0.05. Furthermore, the overlapping DEGs filtered from both datasets were classified as cancer-related targets. 2.8.2 Analysis of brain and liver cancer-related target genes Enrichment analysis (GO and KEGG) of differentially expressed genes (DEGs), commonly upregulated genes, and commonly downregulated genes was performed using the Metascape database ( https://metascape.org/gp/index.html ). The selection criteria included a minimum overlap of 3, a minimum enrichment factor of 1.5, and a significance level of P < 0.01. Additionally, we employed the Gene Expression Profiling Interactive Analysis (GEPIA) tool ( http://gepia.cancer-pku.cn/ ) to perform survival analysis on five key expressed genes. Our objective was to investigate the genetic alterations affecting each target gene and to test hypotheses regarding the recurrence patterns of genetic alteration events and their genomic context in brain tumors and liver tumors [21]. 3. Results 3.1 Toxicity analysis of HA, PP and PT Predictions from the ADMETlab 3.0 database revealed that the co-toxicity of HA, PP, and PT is primarily associated with neurotoxicity and liver damage (Table 1 ). Additionally, the database provided predictions regarding the genotoxicity, respiratory toxicity, and skin irritation potential of these three pesticides. A comprehensive review of both the database and relevant literature was conducted to predict the molecular mechanisms underlying these two diseases. Table 1 The prediction toxicity of Hexaconazole, Propiconazole and Prothioconazole. Ingredient Property Probability Hexaconazole Drug-induced Neurotoxicity 0.875 (++) Hexaconazole Human Hepatotoxicity 0.600 (+) Hexaconazole Genotoxicity 0.564 (+) Hexaconazole Respiratory toxicity 0.749 (++) Hexaconazole Skin sensitization 0.710 (++) Propiconazole Drug-induced Neurotoxicity 0.907 (+++) Propiconazole Human Hepatotoxicity 0.685 (+) Propiconazole Genotoxicity 0.872 (++) Propiconazole Respiratory toxicity 0.672 (+) Propiconazole Skin sensitization 0.658 (+) Prothioconazole Drug-induced Neurotoxicity 0.821 (++) Prothioconazole Human Hepatotoxicity 0.477 (-) Prothioconazole Genotoxicity 0.999 (+++) Prothioconazole Respiratory toxicity 0.723 (++) Prothioconazole Skin sensitization 0.505 (+) Classification endpoints of the prediction probability values: 0-0.1 (---), 0.1–0.3 (--), 0.3–0.5 (-), 0.5–0.7 (+), 0.7–0.9 (++), 0.9–1.0 (+++). Larger probability represents a greater risk of toxicity. 3.2 Targets analysis of HA, PP and PT induce diseases We initially integrated data from the ChEMBL and Swiss-Target-Prediction databases, removing duplicate entries to identify 120 HA targets, 32 PP targets, and 118 PT targets. Following the removal of duplicates, we also identified 5,464 targets highly associated with brain injury and 5,071 targets highly associated with liver injury from the GeneCards and OMIM databases. The list of these target genes can be found in Appendix A. Using a comprehensive Venn diagram approach, we identified a total of 94 HA cross-targets related to brain injury (Fig. 1 A), 28 PP cross-targets related to brain injury (Fig. 1 B), 84 PT cross-targets related to brain injury (Fig. 1 C), 92 HA cross-targets related to liver injury (Fig. 1 D), 24 PP cross-targets related to liver injury (Fig. 1 E), and 84 PT cross-targets related to liver injury (Fig. 1 F). These targets are considered potential candidates associated with brain or liver injuries induced by the three triazole pesticides. 3.3 PPI network and core target acquisition Potential targets associated with brain or liver injury induced by the three triazole pesticides were mapped to the protein-protein interaction (PPI) network using the STRING database (Supplementary Fig. S2 ), resulting in a highly interconnected network. Metrics such as the number of nodes, number of edges, average node degree, and PPI enrichment p -value are detailed in Appendix B. The topological analysis of the network, including degree and betweenness centrality metrics, was performed using Cytoscape. Additionally, we constructed a PPI network diagram to visually represent the results (Fig. 2 ). Then, we employed multiple complementary algorithms, specifically the Betweenness Centrality (BC), Closeness Centrality (CC), and Degree Centrality (DC) algorithms from the CytoHubba plugin, to identify highly connected hub targets. The targets identified as hubs resulted in the selection of a subset of 35 target genes closely related to HA-induced brain injury, 34 target genes closely related to HA-induced liver injury, 9 target genes closely related to PP-induced brain injury, 8 target genes closely related to PP-induced liver injury, 26 target genes closely related to PT-induced brain injury, and 26 target genes closely related to PT-induced liver injury (Supplementary Tables S1-S6). A centralized PPI network diagram was constructed to illustrate the interactions among these hub targets (Fig. 3 ). 3.4 GO analysis Gene Ontology (GO) analysis was conducted on potential targets associated with brain or liver injuries induced by three types of fungicides using the Metascape database. Our findings revealed significant GO terms (Appendix C). The top ten most enriched terms across the Biological Process (BP), Cellular Component (CC), and Molecular Function (MF) categories were prioritized based on gene count and are visually represented in the enrichment analysis plot (Fig. 4 ). In terms of Biological Process (BP), the potential targets associated with brain and liver injuries induced by the three triazole fungicides exhibited significant enrichment in pathways related to cellular responses to organic cyclic compounds and organic nitrogen compounds. Regarding CC, the potential targets showed significant enrichment in protein kinase binding pathways for HA and PT, while for PP, they were significantly enriched in pathways related to DNA and organic compound binding. With respect to MF, the potential targets demonstrated significant enrichment in receptor complex pathways. 3.5 KEGG pathway analysis KEGG pathway analysis was conducted using the Metascape database on potential targets associated with brain or liver injuries induced by three types of fungicides, revealing enriched pathways related to these targets (Appendix C). Overall, these targets were significantly enriched in pathways associated with cancer, the PI3K-Akt signaling pathway, and pathways related to chemical carcinogenesis. We generated a bubble plot (Fig. 4 ) displaying the top 20 most enriched KEGG signaling pathways, prioritized by gene count, to succinctly illustrate the interactions among these pathways and to highlight their significance in the pathogenesis and progression of brain and liver injuries caused by triazole fungicides. Additionally, we performed KEGG analysis on the core targets associated with brain or liver injuries induced by the three fungicides using the Metascape database. This analysis produced a Sankey bubble plot (Fig. 5 ) illustrating the top 20 most enriched KEGG signaling pathways, also prioritized by gene count. These targets showed significant enrichment in pathways related to cancer, the PI3K-Akt signaling pathway, and protein processing. 3.6 Network analysis of diseases induced by HA, PP and PT Using Cytoscape software, we constructed a target network diagram for HA, PP, and PT-induced brain and liver injuries, resulting in a "drug-target pathway" network diagram. As shown in Fig. 6 , deep red or deep blue represents hub targets, while light red or light blue denotes pathways. This diagram provides insights into the core targets, pathways, and the extent of brain or liver injuries induced by the three pesticides. 3.7 Molecular docking analysis of HA, PP and PT and key targets To further elucidate the interactions and relationships between HA, PP, and PT with the key targets (the top three core targets), we conducted comprehensive molecular docking simulations. Docking of each target was performed using AutoDock software, and the results indicated that all binding energies were below − 5.0 kcal/mol [22], demonstrating a strong binding affinity of HA, PP, and PT for these targets (Table 2 ). Our findings suggest that HA, PP, and PT can spontaneously bind to each key target, underscoring their critical roles in the molecular mechanisms underlying brain and liver injuries induced by these three pesticides. Visual representations of the lowest binding energy conformations between these pesticides and each key target were generated using Discovery Studio and PyMOL (Fig. 7 ). These results indicate that these three triazole pesticides may induce brain or liver injury through binding with the aforementioned targets. However, further experimental validation is required. Table 2 The result of molecular docking. Ligand Target Binding energy(kcal/mol) Hexaconazole CASP3 -6.5 Hexaconazole EGFR -6.0 Hexaconazole ESR1 -7.2 Propiconazole PPARG -7.3 Propiconazole TP53 -7.1 Propiconazole HSP90AA1 -6.6 Prothioconazole EGFR -6.4 Prothioconazole ESR1 -7.4 Prothioconazole PTGS2 -7.3 3.8 Screening for brain and liver cancer-related DEGs and analysis of target genes Building on our previous KEGG analysis, we revealed that the targets associated with brain and liver toxicity induced by the three triazole fungicides are enriched in pathways related to cancer. Therefore, we aimed to further elucidate the mechanistic impact of core targets for triazole fungicides-elicited brain and liver injury on brain and liver cancer. To identify genes associated with brain tumors and liver tumors, we utilized two gene expression profiles: GSE30563 and GSE117361. To visually highlight the unique distribution and significance of these differentially expressed genes (DEGs), we employed volcano plots to visualize all DEGs across the two datasets (Figs. 8 A-B). Detailed information regarding these datasets is provided in Appendix D. Venn diagrams were constructed to illustrate the overlap between the two sets of DEGs, revealing a total of 287 overlapping genes, of which 283 exhibited common expression patterns. Among these, 143 genes were significantly upregulated, while 140 genes were significantly downregulated (Fig. 8 C). We conducted GO and KEGG enrichment analyses on the 287 overlapping gene targets to explore the differential expression patterns in brain tumors and liver tumors compared to normal tissues. Detailed information regarding the enrichment results is provided in Appendix E. The top ten most enriched terms in BP, CC, and MF categories were prioritized based on gene count and visually represented in a histogram (Fig. 8 D). Then we generated a histogram of the top 20 most enriched KEGG signaling pathways, prioritized by gene count (Fig. 8 E). To further investigate the functional roles of the DEGs that were commonly upregulated or downregulated, we performed GO and KEGG enrichment analyses on the 283 genes exhibiting common expression patterns. The results of the GO and KEGG enrichment analyses are visualized in Fig. 8 F-G. In GO enrichment analysis,DEGs was significantly related to immune response and kinase binding pathways. In KEGG enrichment analysis, DEGs was significantly associated with immune disorders and infection-related pathways. This suggests that DEGs is highly correlated with immune responses in brain and liver tumors. Additionally, we selected key DEGs with a p-value less than 0.005, consisting of 15 genes with a common expression pattern. Among these, 5 genes were commonly upregulated, and 10 genes were commonly downregulated (Supplementary Fig. S4 -S5). Kaplan-Meier analysis was performed on the 5 commonly upregulated genes, and the results were visualized (Supplementary Fig. S6 ). We observed a significant difference in survival times between the high-expression and low-expression groups for these 5 genes. 4. Discussion The adoption of triazole fungicides has significantly increased in modern agriculture due to their high efficacy, low application rates, and relatively favorable environmental profiles compared to traditional classes of fungicides. However, the triazole fungicides have been detected in a variety of food products, beverages, and human biomonitoring studies. Given the widespread human exposure to these compounds, comprehensive investigations are needed to evaluate their potential safety hazards [23, 24]. In our study, we integrated a comprehensive approach utilizing network toxicology, molecular docking, and multi-layered bioinformatics to explore the effects of exposure to the triazole fungicides HA, PP, and PT on the brain and liver, as well as their potential mechanisms of action. After applying network-based computational assessment tools, we identified potential targets associated with brain or liver damage induced by HA, PP, and PT using the ChEMBL, Swiss Target Prediction, GeneCards, OMIM, and DisGeNet databases. Utilizing the STRING database and Cytoscape software, we constructed an interaction network for these potential targets and extracted core targets, including EGFR , CASP3 , ESR1 , PPARG , TP53 , HSP90AA1 , and PTGS2 . These targets were designated as critical for the induction of brain or liver damage by the three triazole fungicides. Epidermal growth factor receptor (EGFR), a transmembrane receptor tyrosine kinase, plays an important role in maintaining normal tissue and cell signaling pathways [25]. Dysregulation of EGFR has been implicated in neurodegenerative diseases, where altered signaling pathways contribute to neuronal apoptosis and cognitive decline [26]. EGFR signaling plays a central role in the regenerative response after liver injury and is involved in cell transformation associated with chronic injury [27]. In addition, studies have shown that EGFR is also associated with apoptotic cell death signaling in various hepatocytes, mitochondrial dysfunction, and acute liver necrosis [28]. Caspase-3 (CASP3), a widely expressed member of the conserved protein family, plays a key role in regulating the growth and homeostasis maintenance of normal and malignant cells and tissues in multicellular organisms [29]. Activation of CASP3 is strongly associated with neurodegenerative diseases and traumatic brain injury, and its upregulation leads to neuronal apoptosis, exacerbating cognitive decline and dysfunction [30, 31]. CASP3 mediates hepatocyte apoptosis in response to hepatotoxic substances, leading to liver damage and fibrosis, cirrhosis and other conditions [32]. The study suggests that modulating CASP3 activity may provide a therapeutic pathway for alleviating cell death in brain and liver injury [33]. Estrogen receptor alpha (ESR1) is a key regulator that regulates the effects of estrogen on a variety of tissues, including the brain and liver [34]. Studies have shown that ESR1 signaling has a neuroprotective effect, as estrogen can enhance neuronal survival and reduce apoptosis after trauma or ischemic injury [35]. And ESR1 activation is associated with improved cognitive function and reduced neuroinflammation [36]. In liver injury, ESR1 plays a complex role in promoting liver cell regeneration and preventing oxidative stress [37]. However, dysregulated ESR1 signaling can contribute to the progression of liver diseases, including fatty liver and fibrosis [38]. Peroxisome proliferator activating receptor (PPARG) is a nuclear receptor that regulates genes involved in lipid metabolism, inflammation, and insulin sensitivity [39]. Activation of PPARG has been shown to have neuroprotective effects, enhancing neuronal survival and reducing neuroinflammation after trauma or ischemic injury [40]. Studies have shown that activating PPARG can improve cognitive outcomes and promote repair mechanisms in neurodegenerative diseases [41]. In addition, PPARG plays a key role in lipid homeostasis and inflammation resolution [42]. Activation of PPARG is associated with a reduction in steatosis and fibrosis in various liver diseases [43]. The TP53 gene encodes the tumor suppressor protein p53, which plays a crucial role in cellular stress responses, including apoptosis, cell cycle regulation, and DNA repair [44]. In brain injury, activation of TP53 is often associated with neurodegenerative diseases, and in response to oxidative stress and DNA damage, it mediates neuronal apoptosis, leading to cognitive decline and neuron loss [45]. In liver injury, p53 can promote apoptosis of damaged liver cells and prevent tumorigenesis, but excessive activation can lead to chronic inflammation and fibrosis [46]. The HSP90AA1 gene encodes heat shock protein 90α (Hsp90), an important companion involved in protein folding and cellular stress response [47]. Previous studies have shown that Hsp90 exhibits neuroprotective properties, promoting the stabilization of client proteins and preventing neuronal cell apoptosis during ischemic or traumatic events [48]. On the other hand, Hsp90 plays an important role in maintaining the integrity of hepatocytes under stress conditions such as oxidative damage and inflammation [49]. The enzyme protein prostaglandin-endoperoxide synthase 2 (PTGS2) is involved in inflammation and other physiological processes by catalyzing the conversion of arachidonic acid to prostaglandin [50]. For brain injury, COX-2 expression is upregulated in response to neuronal injury, leading to neuroinflammation and pain [51]. While COX-2 is involved in protective mechanisms that promote healing and regeneration, overactivation can exacerbate neuronal damage and lead to chronic neurodegenerative diseases [51]. For liver damage, PTGS2 is involved in the inflammatory response associated with liver disease, including non-alcoholic fatty liver disease and hepatitis [52]. Elevated COX-2 levels are associated with increased inflammation and apoptosis of liver cells, which worsens liver damage [53]. Pathway analysis highlights key insights into the molecular mechanisms of brain and liver injury induced by triazole fungicides. The significant enrichment of targets in cancer-related pathways suggests that exposure to these fungicides may promote carcinogenic processes in these organs, raising concerns about long-term health effects [54]. Specifically, the PI3K-Akt signaling pathway plays an important role in cell survival, proliferation, and metabolism [55]. Dysregulation of this pathway may lead to neuroinflammatory responses and hepatotoxicity, which may exacerbate tissue damage. In addition, the identification of pathways associated with chemical carcinogenesis highlights the potential of these fungicides to initiate or enhance malignant transformation. 5. Conclusion In summary, the present study combined applied network toxicology, molecular docking, and multilevel bioinformatics data to systematically explore the toxicological and molecular mechanisms of three triazole fungicides, including Hexaconazole, Propiconazole, and Prothioconazole, in brain and liver injury. Our findings demonstrate that the toxicity profiles of these compounds are closely linked to their interactions with specific biological pathways, highlighting the potential risks they pose to human health and the environment. This integrated approach not only enhances our understanding of fungicide-induced toxicity but also offers valuable insights for future research aimed at developing safer agricultural practices and regulatory policies. Abbreviations HA, Hexaconazole; PP, Propiconazole; PT, Prothioconazole; DEGs, differential gene expression; BP, biological processes; CC, cellular components; MF, molecular functions; GP: Genetic Information Processing; EI: Environmental Information Processing; CP: Cellular Processes; OR: Organismal Systems; HD: Human Diseases. Declarations Funding: This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors. Author Contribution Author contributions: Qing Lv: Methodology, Conceptualization, Data analysis, Visualization, Writing – original draft. Xingke Zhu: Methodology, Conceptualization, Data analysis, Visualization, Writing – original draft, Writing – review & editing.Author Statement: By submitting this manuscript, we, the authors, affirm that we have contributed significantly to the research presented and are in agreement with the content and conclusions of the work. We confirm that this manuscript has not been published elsewhere and is not under consideration by any other journal. All authors have reviewed and approved the final version of the manuscript. We declare that there are no conflicts of interest related to this submission, and we have adhered to ethical standards in conducting our research. We agree to abide by the policies and procedures of the journal regarding the publication process. Acknowledgments: The author declares no financial conflicts of interest or financial relationships with any organizations or individuals that could potentially bias the results or interpretation of this study. This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors. Ethics approval and consent to participate: Not applicable. Running Title Toxicity and mechanisms of three triazole fungicides. Availability of data and materials Data will be made available on request. Further inquiries can be directed to the corresponding authors. References Fernández-Vizcaíno E, Mateo R, Fernández de Mera IG, Mougeot F, Camarero PR, Ortiz-Santaliestra ME. Transgenerational effects of triazole fungicides on gene expression and egg compounds in non-exposed offspring: A case study using Red-Legged Partridges (Alectoris rufa). Sci Total Environ . 2024;926:171546. doi:10.1016/j.scitotenv.2024.171546 He Z, Zhang J, Shi D, et al. Deoxynivalenol in Fusarium graminearum : Evaluation of Cyproconazole Stereoisomers In Vitro and In Planta . J Agric Food Chem . 2021;69(34):9735–9742. doi:10.1021/acs.jafc.1c02555 Raina R, Smith E. Determination of azole fungicides in atmospheric samples collected in the Canadian prairies by LC/MS/MS. J AOAC Int . 2012;95(5):1350–1356. doi:10.5740/jaoacint.sge_raina Mao CX, Luo J, Zhang Y, Zhang CQ. Targeted deletion of three CYP51s in Fusarium fujikuroi and their different roles in determining sensitivity to 14α-demethylase inhibitor fungicides. Pest Manag Sci . 2023;79(4):1324–1330. doi:10.1002/ps.7304 Zou PC, Zhang Y, Bian Y, et al. Triazoles in the environment: An update on sample pretreatment and analysis methods. Ecotoxicol Environ Saf . Published online October 8, 2024. doi:10.1016/j.ecoenv.2024.117156 Yao X, Lv H, Wang Q, et al. Novel Insights into Stereoselective Reproductive Toxicity Induced by Mefentrifluconazole in Earthworms ( Eisenia fetida ): First Report of Estrogenic Effects. J Agric Food Chem . 2024;72(35):19304–19311. doi:10.1021/acs.jafc.4c04168 Zhang Y, Guo J, Tang C, Xu K, Li Z, Wang C. Early life stage exposure to fenbuconazole causes multigenerational cardiac developmental defects in zebrafish and potential reasons. Environ Pollut . 2024;349:123938. doi:10.1016/j.envpol.2024.123938 Cang T, Lou Y, Zhu YC, et al. Mixture toxicities of tetrachlorantraniliprole and tebuconazole to honey bees (Apis mellifera L.) and the potential mechanism. Environ Int . 2023;172:107764. doi:10.1016/j.envint.2023.107764 Wang X, Liu Y, Xue M, Wang Z, Yu J, Guo X. Enantioselective degradation of chiral fungicides triticonazole and prothioconazole in soils and their enantioselective accumulation in earthworms Eisenia fetida. Ecotoxicol Environ Saf . 2019;183:109491. doi:10.1016/j.ecoenv.2019.109491 Wang X, Ni H, Xu W, et al. Difenoconazole induces oxidative DNA damage and mitochondria mediated apoptosis in SH-SY5Y cells. Chemosphere . 2021;283:131160. doi:10.1016/j.chemosphere.2021.131160 Sanchez CL, Souders CL 2nd, Pena-Delgado CJ, et al. Neurotoxicity assessment of triazole fungicides on mitochondrial oxidative respiration and lipids in differentiated human SH-SY5Y neuroblastoma cells. Neurotoxicology . 2020;80:76–86. doi:10.1016/j.neuro.2020.06.009 Yang Q, Deng P, Xing D, et al. Developmental Neurotoxicity of Difenoconazole in Zebrafish Embryos. Toxics . 2023;11(4):353. Published 2023 Apr 8. doi:10.3390/toxics11040353 Zerdoug A, Le Vée M, Le Mentec H, et al. Induction of drug metabolizing enzyme and drug transporter expression by antifungal triazole pesticides in human HepaSH hepatocytes. Chemosphere . Published online October 5, 2024. doi:10.1016/j.chemosphere.2024.143474 Li C, Yuan S, Zhou Y, et al. Microplastics reduce the bioaccumulation and oxidative stress damage of triazole fungicides in fish. Sci Total Environ . 2022;806(Pt 4):151475. doi:10.1016/j.scitotenv.2021.151475 Zeng Z, Hu J, Xiao G, et al. Integrating network toxicology and molecular docking to explore the toxicity of the environmental pollutant butyl hydroxyanisole: An example of induction of chronic urticaria. Heliyon . 2024;10(15):e35409. Published 2024 Jul 30. doi:10.1016/j.heliyon.2024.e35409 Tong L, Zhou Y, Zhang Y, et al. Unveiling the protein-lipid interaction mechanism: How the sturgeon lipids diminish the surimi gel properties. Food Chem . 2025;462:140926. doi:10.1016/j.foodchem.2024.140926 Yang Z, Han S, Keller M, et al. Structural basis of ligand binding modes at the neuropeptide Y Y 1 receptor. Nature . 2018;556(7702):520–524. doi:10.1038/s41586-018-0046-x Mendez D, Gaulton A, Bento AP, et al. ChEMBL: towards direct deposition of bioassay data. Nucleic Acids Res . 2019;47(D1):D930-D940. doi:10.1093/nar/gky1075 Liu M, Wang Y, Deng W, et al. Combining network pharmacology, machine learning, molecular docking and molecular dynamic to explore the mechanism of Chufeng Qingpi decoction in treating schistosomiasis. Front Cell Infect Microbiol . 2024;14:1453529. Published 2024 Sep 6. doi:10.3389/fcimb.2024.1453529 Zhang M, Zhang X, Pei J, et al. Identification of phytochemical compounds of Fagopyrum dibotrys and their targets by metabolomics, network pharmacology and molecular docking studies. Heliyon . 2023;9(3):e14029. doi:10.1016/j.heliyon.2023.e14029 Shen SH, Chen SF, Guo JF, Wang ZX. The GGCT and REST positive feedback loop promotes tumor growth in Glioma. Transl Oncol . 2024;49:102083. doi:10.1016/j.tranon.2024.102083 Zhang Q, Liang J, Li X, et al. Exploring antithrombotic mechanisms and effective constituents of Lagopsis supina using an integrated strategy based on network pharmacology, molecular docking, metabolomics, and experimental verification in rats. J Ethnopharmacol . 2025;336:118717. doi:10.1016/j.jep.2024.118717 Liu J, Nie Y, Niu Y, Li L, Jing X. Lignin-based emulsive liquid-liquid microextraction for detecting triazole fungicides in water, juice, vinegar, and alcoholic beverages via UHPLC-MS/MS. Food Chem . 2024;459:140407. doi:10.1016/j.foodchem.2024.140407 Marciano LPA, Kleinstreuer N, Chang X, Costa LF, Silvério ACP, Martins I. A novel approach to triazole fungicides risk characterization: Bridging human biomonitoring and computational toxicology. Sci Total Environ . 2024;953:176003. doi:10.1016/j.scitotenv.2024.176003 Damare R, Engle K, Kumar G. Targeting epidermal growth factor receptor and its downstream signaling pathways by natural products: A mechanistic insight. Phytother Res . 2024;38(5):2406–2447. doi:10.1002/ptr.8166 Mansour HM, Fawzy HM, El-Khatib AS, Khattab MM. Repurposed anti-cancer epidermal growth factor receptor inhibitors: mechanisms of neuroprotective effects in Alzheimer's disease. Neural Regen Res. 2022;17(9):1913–1918. doi:10.4103/1673-5374.332132 González L, Díaz ME, Miquet JG, Sotelo AI, Dominici FP. Growth Hormone Modulation of Hepatic Epidermal Growth Factor Receptor Signaling. Trends Endocrinol Metab . 2021;32(6):403–414. doi:10.1016/j.tem.2021.03.004 Bhushan B, Michalopoulos GK. Role of epidermal growth factor receptor in liver injury and lipid metabolism: Emerging new roles for an old receptor. Chem Biol Interact . 2020;324:109090. doi:10.1016/j.cbi.2020.109090 Eskandari E, Eaves CJ. Paradoxical roles of caspase-3 in regulating cell survival, proliferation, and tumorigenesis. J Cell Biol . 2022;221(6):e202201159. doi:10.1083/jcb.202201159 Unnisa A, Greig NH, Kamal MA. Inhibition of Caspase 3 and Caspase 9 Mediated Apoptosis: A Multimodal Therapeutic Target in Traumatic Brain Injury. Curr Neuropharmacol . 2023;21(4):1001–1012. doi:10.2174/1570159X20666220327222921 Rajesh Y, Kanneganti TD. Innate Immune Cell Death in Neuroinflammation and Alzheimer's Disease. Cells . 2022;11(12):1885. Published 2022 Jun 10. doi:10.3390/cells11121885 Ma XY, Zhang M, Fang G, et al. Ursolic acid reduces hepatocellular apoptosis and alleviates alcohol-induced liver injury via irreversible inhibition of CASP3 in vivo. Acta Pharmacol Sin . 2021;42(7):1101–1110. doi:10.1038/s41401-020-00534-y Wu J, Ye X, Yang S, Yu H, Zhong L, Gong Q. Systems Pharmacology Study of the Anti-Liver Injury Mechanism of Citri Reticulatae Pericarpium. Front Pharmacol . 2021;12:618846. Published 2021 Apr 12. doi:10.3389/fphar.2021.618846 Bartkowiak-Wieczorek J, Jaros A, Gajdzińska A, et al. The Dual Faces of Oestrogen: The Impact of Exogenous Oestrogen on the Physiological and Pathophysiological Functions of Tissues and Organs. Int J Mol Sci . 2024;25(15):8167. Published 2024 Jul 26. doi:10.3390/ijms25158167 Rettberg JR, Yao J, Brinton RD. Estrogen: a master regulator of bioenergetic systems in the brain and body. Front Neuroendocrinol . 2014;35(1):8–30. doi:10.1016/j.yfrne.2013.08.001IF: 6.5 Q1 Landel V, Millet P, Baranger K, Loriod B, Féron F. Vitamin D interacts with Esr1 and Igf1 to regulate molecular pathways relevant to Alzheimer's disease. Mol Neurodegener . 2016;11:22. Published 2016 Mar 1. doi:10.1186/s13024-016-0087-2 Palmisano BT, Zhu L, Stafford JM. Role of Estrogens in the Regulation of Liver Lipid Metabolism. Adv Exp Med Biol . 2017;1043:227–256. doi:10.1007/978-3-319-70178-3_12 Yin G, Liang H, Sun W, et al. Shuangyu Tiaozhi decoction alleviates non-alcoholic fatty liver disease by improving lipid deposition, insulin resistance, and inflammation in vitro and in vivo . Front Pharmacol . 2022;13:1016745. Published 2022 Nov 23. doi:10.3389/fphar.2022.1016745 Semple RK, Chatterjee VK, O'Rahilly S. PPAR gamma and human metabolic disease. J Clin Invest . 2006;116(3):581–589. doi:10.1172/JCI28003 Rudko OI, Tretiakov AV, Naumova EA, Klimov EA. Role of PPARs in Progression of Anxiety: Literature Analysis and Signaling Pathways Reconstruction. PPAR Res . 2020;2020:8859017. Published 2020 Nov 29. doi:10.1155/2020/8859017 Gold PW. The PPARg System in Major Depression: Pathophysiologic and Therapeutic Implications. Int J Mol Sci . 2021;22(17):9248. Published 2021 Aug 26. doi:10.3390/ijms22179248 Galbraith L, Leung HY, Ahmad I. Lipid pathway deregulation in advanced prostate cancer. Pharmacol Res . 2018;131:177–184. doi:10.1016/j.phrs.2018.02.022 Filali-Mouncef Y, Hunter C, Roccio F, et al. The ménage à trois of autophagy, lipid droplets and liver disease. Autophagy . 2022;18(1):50–72. doi:10.1080/15548627.2021.1895658 Voskarides K, Giannopoulou N. The Role of TP53 in Adaptation and Evolution. Cells . 2023;12(3):512. Published 2023 Feb 3. doi:10.3390/cells12030512 Dai CQ, Luo TT, Luo SC, et al. p53 and mitochondrial dysfunction: novel insight of neurodegenerative diseases. J Bioenerg Biomembr . 2016;48(4):337–347. doi:10.1007/s10863-016-9669-5 Staib F, Hussain SP, Hofseth LJ, Wang XW, Harris CC. TP53 and liver carcinogenesis. Hum Mutat . 2003;21(3):201–216. doi:10.1002/humu.10176 Zuehlke AD, Beebe K, Neckers L, Prince T. Regulation and function of the human HSP90AA1 gene. Gene . 2015;570(1):8–16. doi:10.1016/j.gene.2015.06.018 Smajić S, Prada-Medina CA, Landoulsi Z, et al. Single-cell sequencing of human midbrain reveals glial activation and a Parkinson-specific neuronal state. Brain . 2022;145(3):964–978. doi:10.1093/brain/awab446 Wei D, Tian X, Zhu L, Wang H, Sun C. USP14 governs CYP2E1 to promote nonalcoholic fatty liver disease through deubiquitination and stabilization of HSP90AA1. Cell Death Dis . 2023;14(8):566. Published 2023 Aug 26. doi:10.1038/s41419-023-06091-6 Martín-Vázquez E, Cobo-Vuilleumier N, López-Noriega L, Lorenzo PI, Gauthier BR. The PTGS2/COX2-PGE 2 signaling cascade in inflammation: Pro or anti? A case study with type 1 diabetes mellitus. Int J Biol Sci . 2023;19(13):4157–4165. Published 2023 Aug 6. doi:10.7150/ijbs.86492 Li Q, Han X, Lan X, et al. Inhibition of neuronal ferroptosis protects hemorrhagic brain. JCI Insight . 2017;2(7):e90777. Published 2017 Apr 6. doi:10.1172/jci.insight.90777 Zhang CY, Tan XH, Yang HH, et al. COX-2/sEH Dual Inhibitor Alleviates Hepatocyte Senescence in NAFLD Mice by Restoring Autophagy through Sirt1/PI3K/AKT/mTOR. Int J Mol Sci . 2022;23(15):8267. Published 2022 Jul 27. doi:10.3390/ijms23158267 Martín-Sanz P, Casado M, Boscá L. Cyclooxygenase 2 in liver dysfunction and carcinogenesis: Facts and perspectives. World J Gastroenterol . 2017;23(20):3572–3580. doi:10.3748/wjg.v23.i20.3572 Perdichizzi S, Mascolo MG, Silingardi P, et al. Cancer-related genes transcriptionally induced by the fungicide penconazole. Toxicol In Vitro . 2014;28(1):125–130. doi:10.1016/j.tiv.2013.06.006 Graziosi A, Sita G, Corrieri C, et al. Effects of Subtoxic Concentrations of Atrazine, Cypermethrin, and Vinclozolin on microRNA-Mediated PI3K/Akt/mTOR Signaling in SH-SY5Y Cells. Int J Mol Sci . 2022;23(23):14538. Published 2022 Nov 22. doi:10.3390/ijms232314538 Additional Declarations No competing interests reported. Supplementary Files GraphicalAbstract.tif SupplementaryContents.docx AppendixA.xlsx AppendixB.xlsx AppendixC.xlsx AppendixD.xlsx AppendixE.xlsx Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-5318182","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":370796022,"identity":"12b94b3d-bb9d-437b-a810-a34ee739d230","order_by":0,"name":"Qing Lv","email":"","orcid":"","institution":"Hubei University of Chinese Medicine","correspondingAuthor":false,"prefix":"","firstName":"Qing","middleName":"","lastName":"Lv","suffix":""},{"id":370796023,"identity":"3a6572b3-e7ef-42e0-845c-d27bceae35e9","order_by":1,"name":"Xingke Zhu","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAAu0lEQVRIie3QPwrCMBTH8RcKusS6JhTEIzwJZMpNXOLSzRt0iBR0KXQteIkeofDASfAK8QZ27+CfzUGabg75zu8z/B5ALPaHCWCOOQBcJgn5aUSeZjkGEoAPwRtfiyAiz3Tom4GUIg4IhdmOkizdlbLluda06Dxc8r0bIyvOjswLYzSlFpmjUILCqJKjCCLZm7TWaExCiaxYKZsuV4JeT7YhW8R1fu+rgTZ1TeQfhRkn39lp57FYLBb71RM+RzmDlSzp8QAAAABJRU5ErkJggg==","orcid":"","institution":"Hubei University of Chinese Medicine","correspondingAuthor":true,"prefix":"","firstName":"Xingke","middleName":"","lastName":"Zhu","suffix":""}],"badges":[],"createdAt":"2024-10-23 10:38:09","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-5318182/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-5318182/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":67736093,"identity":"02ca97d5-92ae-4966-9c10-75a01119b904","added_by":"auto","created_at":"2024-10-29 08:04:02","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":4394282,"visible":true,"origin":"","legend":"\u003cp\u003eVeen diagram of the potential targets of brain injury induced by Hexaconazole (A), Propiconazole (B) and Prothioconazole (C); Veen diagram of the potential targets of liver injury induced by Hexaconazole (D), Propiconazole (E) and Prothioconazole (F).\u003c/p\u003e","description":"","filename":"Fig.1.tif.png","url":"https://assets-eu.researchsquare.com/files/rs-5318182/v1/8270e6e54401d15e8b751e7e.png"},{"id":67735779,"identity":"2e23ccb4-ce66-494d-aa2a-3a91e800bb5e","added_by":"auto","created_at":"2024-10-29 07:56:02","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":14047985,"visible":true,"origin":"","legend":"\u003cp\u003ePPI network of the potential targets of brain injury induced by Hexaconazole (A), Propiconazole (B) and Prothioconazole (C); PPI network of the potential targets of liver injury induced by Hexaconazole (D), Propiconazole (E) and Prothioconazole (F).\u003c/p\u003e","description":"","filename":"Fig.2.tif.png","url":"https://assets-eu.researchsquare.com/files/rs-5318182/v1/f7732e4f0539d46935c8532d.png"},{"id":67736092,"identity":"3b0c73b4-f75d-4b5c-81dd-6c61cdb0bcc1","added_by":"auto","created_at":"2024-10-29 08:04:02","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":8264435,"visible":true,"origin":"","legend":"\u003cp\u003ePPI network of the core targets of brain injury induced by Hexaconazole (A), Propiconazole (B) and Prothioconazole (C); PPI network of the core targets of liver injury induced by Hexaconazole (D), Propiconazole (E) and Prothioconazole (F).\u003c/p\u003e","description":"","filename":"Fig.3.tif.png","url":"https://assets-eu.researchsquare.com/files/rs-5318182/v1/d17a1e5e1a000545358f3175.png"},{"id":67735778,"identity":"3b1797bb-89b9-4bb4-8d36-073a83c34de2","added_by":"auto","created_at":"2024-10-29 07:56:02","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":4783368,"visible":true,"origin":"","legend":"\u003cp\u003eGO enrichment analysis of the potential targets of brain injury induced by Hexaconazole (A), Propiconazole (B) and Prothioconazole (C); GO enrichment analysis of the potential targets of liver injury induced by Hexaconazole (D), Propiconazole (E) and Prothioconazole (F).\u003c/p\u003e","description":"","filename":"Fig.4.tif.png","url":"https://assets-eu.researchsquare.com/files/rs-5318182/v1/e1eb588363721230305d811c.png"},{"id":67736099,"identity":"2ce9dab8-1c29-4df9-9f17-f51be94bf7b4","added_by":"auto","created_at":"2024-10-29 08:04:02","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":17194088,"visible":true,"origin":"","legend":"\u003cp\u003eKEGG enrichment analysis of the core targets (top 10) of brain injury induced by Hexaconazole (A), Propiconazole (B) and Prothioconazole (C); KEGG enrichment analysis of the core targets (top 10) of liver injury induced by Hexaconazole (D), Propiconazole (E) and Prothioconazole (F).\u003c/p\u003e","description":"","filename":"Fig.5.tif.png","url":"https://assets-eu.researchsquare.com/files/rs-5318182/v1/610c395efa32831800b03144.png"},{"id":67737244,"identity":"e2b2b458-d984-4fbf-8e2a-b18c5553a6ca","added_by":"auto","created_at":"2024-10-29 08:12:02","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":16734547,"visible":true,"origin":"","legend":"\u003cp\u003e“Drug-target-pathway” network diagram of the brain injury induced by Hexaconazole (A), Propiconazole (B) and Prothioconazole (C); “Drug-target-pathway” network diagram of the liver injury induced by Hexaconazole (D), Propiconazole (E) and Prothioconazole (F).\u003c/p\u003e","description":"","filename":"Fig.6.tif.png","url":"https://assets-eu.researchsquare.com/files/rs-5318182/v1/ceb359d0f422d9b49149c58f.png"},{"id":67737246,"identity":"2121de16-00bd-41a8-9d9f-1b1baac6e603","added_by":"auto","created_at":"2024-10-29 08:12:02","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":5780587,"visible":true,"origin":"","legend":"\u003cp\u003eMolecular docking results of the lowest binding energy in each target protein with Hexaconazole, Propiconazole and Prothioconazole.\u003c/p\u003e","description":"","filename":"Fig.7.tif.png","url":"https://assets-eu.researchsquare.com/files/rs-5318182/v1/2ecae560d6c1588f48a2ef5d.png"},{"id":67736096,"identity":"7736ba47-0b07-4606-a24d-677e052e3837","added_by":"auto","created_at":"2024-10-29 08:04:02","extension":"png","order_by":8,"title":"Figure 8","display":"","copyAsset":false,"role":"figure","size":6847307,"visible":true,"origin":"","legend":"\u003cp\u003eCharacteristics of brain and liver cancer-related DEGs. A: Volcano plots showing DEGs between tumoral brain tissues and non-tumoral brain tissues. B: Volcano plots showing DEGs between tumoral liver tissues and non-tumoral liver tissues. C: Venn map of differentially expressed genes (\u003cem\u003eP\u003c/em\u003e \u0026lt; 0.05) from the GSE30563 and GSE117361 data sets. D: GO enrichment analysis of common DEGs in brain and liver tumors. E: KEGG enrichment analysis of common DEGs in brain and liver tumors. F: GO enrichment analysis of jointly upregulated genes in brain and liver tumors. G: KEGG enrichment analysis of jointly upregulated genes in brain and liver tumors.\u003c/p\u003e","description":"","filename":"Fig.8.tif.png","url":"https://assets-eu.researchsquare.com/files/rs-5318182/v1/0eeed90c15e5bdaac1bb716b.png"},{"id":67737639,"identity":"8e9ffff5-e6d3-440b-b23e-d9a9f8b8bdcc","added_by":"auto","created_at":"2024-10-29 08:20:47","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":76662152,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-5318182/v1/4195a63a-17b1-4991-a667-a9795c9659a2.pdf"},{"id":67735773,"identity":"f2ffc6e4-b0d4-4bd5-9226-2055d05476b0","added_by":"auto","created_at":"2024-10-29 07:56:02","extension":"tif","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":1104202,"visible":true,"origin":"","legend":"","description":"","filename":"GraphicalAbstract.tif","url":"https://assets-eu.researchsquare.com/files/rs-5318182/v1/a922f54e4304bab76d80cec3.tif"},{"id":67737566,"identity":"361f869f-3b43-44f5-be7f-c54720021404","added_by":"auto","created_at":"2024-10-29 08:20:02","extension":"docx","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":5513111,"visible":true,"origin":"","legend":"","description":"","filename":"SupplementaryContents.docx","url":"https://assets-eu.researchsquare.com/files/rs-5318182/v1/3cabaef7b2db6348bf7ae6f2.docx"},{"id":67735774,"identity":"5321d7c6-cf04-4e30-b8c9-acc36be85fb0","added_by":"auto","created_at":"2024-10-29 07:56:02","extension":"xlsx","order_by":3,"title":"","display":"","copyAsset":false,"role":"supplement","size":143322,"visible":true,"origin":"","legend":"","description":"","filename":"AppendixA.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-5318182/v1/29d42b27d25e416737c5517d.xlsx"},{"id":67735775,"identity":"9bc39dce-c9a0-412b-9dee-9f9e46fda18f","added_by":"auto","created_at":"2024-10-29 07:56:02","extension":"xlsx","order_by":4,"title":"","display":"","copyAsset":false,"role":"supplement","size":9447,"visible":true,"origin":"","legend":"","description":"","filename":"AppendixB.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-5318182/v1/46de700c4f83299ae831f7d4.xlsx"},{"id":67736094,"identity":"4c8ed4a5-c1e7-4790-8a73-a6b539b953f5","added_by":"auto","created_at":"2024-10-29 08:04:02","extension":"xlsx","order_by":5,"title":"","display":"","copyAsset":false,"role":"supplement","size":308389,"visible":true,"origin":"","legend":"","description":"","filename":"AppendixC.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-5318182/v1/9525c0aaa2531bc3fe85c2f3.xlsx"},{"id":67737243,"identity":"567c9a92-9b77-41a1-a63c-6822cdcfa5a4","added_by":"auto","created_at":"2024-10-29 08:12:02","extension":"xlsx","order_by":6,"title":"","display":"","copyAsset":false,"role":"supplement","size":1908990,"visible":true,"origin":"","legend":"","description":"","filename":"AppendixD.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-5318182/v1/67f37d2fb6b034db162ca922.xlsx"},{"id":67735782,"identity":"a9fd8579-dbf5-4d04-beaf-341d6d058019","added_by":"auto","created_at":"2024-10-29 07:56:02","extension":"xlsx","order_by":7,"title":"","display":"","copyAsset":false,"role":"supplement","size":68801,"visible":true,"origin":"","legend":"","description":"","filename":"AppendixE.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-5318182/v1/e2414f288502890e6e471cb4.xlsx"}],"financialInterests":"No competing interests reported.","formattedTitle":"Efficient analysis of the toxicity and mechanisms of Hexaconazole and two other triazole fungicides: insights from integrated network toxicology, molecular docking and bioinformatics data","fulltext":[{"header":"1. Introduction","content":"\u003cp\u003eTriazole fungicides represent a significant class of agricultural chemicals widely utilized for the control of various fungal pathogens affecting crops [1]. Characterized by their five-membered ring structure containing three nitrogen atoms, triazoles function primarily by inhibiting the biosynthesis of ergosterol, an essential component of fungal cell membranes [2]. Among the most commonly used triazole fungicides are Hexaconazole (HA), Propiconazole (PP) and Prothioconazole (PT) which exhibit broad-spectrum efficacy against numerous fungal diseases, including rusts, powdery mildews, and blights [c]. These compounds are particularly valued for their systemic properties, allowing them to be absorbed and translocated within plant tissues, thereby providing effective protection against both foliar and soil-borne pathogens [4].\u003c/p\u003e \u003cp\u003eResearch on the toxicity of triazole fungicides has garnered increasing attention due to their widespread use in agriculture and potential environmental and health impacts. Triazoles, including well-known compounds such as Hexaconazole, Tebuconazole, and Propiconazole, have been demonstrated to possess varying degrees of toxicity to non-target organisms, including aquatic life, terrestrial invertebrates, and mammals [5\u0026ndash;7].Furthermore, significant efforts have been made to evaluate the ecological risks associated with triazole fungicides. Toxicological assessments, including acute and chronic exposure studies, have revealed adverse effects on aquatic organisms such as fish and invertebrates, as well as terrestrial species like bees and earthworms [8, 9]. Recent studies have highlighted the neurotoxic potential of triazole compounds, linking exposure to cognitive impairments and neurodegenerative conditions [10\u0026ndash;12]. Moreover, the involvement of triazoles in disrupting hepatic metabolic pathways raises concerns about their long-term effects on liver health [13, 14]. However, the precise biochemical pathways through which triazole fungicides induce neurotoxicity and hepatotoxicity remain poorly characterized. Therefore, the mechanism of brain injury and liver injury induced by triazole fungicides still needs further study.\u003c/p\u003e \u003cp\u003eAdvancements in toxicology have highlighted the importance of an integrated approach to elucidate the complex interactions between chemical agents and biological systems. Network toxicology offers a promising framework for this analysis, enabling researchers to examine the relationships between chemical exposure, biological pathways, and toxicological outcomes [15]. By leveraging bioinformatics tools, researchers can integrate large datasets, including gene expression profiles, protein interactions, and metabolic pathways, to identify key molecular targets and pathways affected by these fungicides.\u003c/p\u003e \u003cp\u003eMolecular docking techniques further enhance our understanding of the interactions between triazole fungicides and their biological targets [16]. This computational method allows for the prediction of binding affinities and interaction modes of fungicides with specific proteins, providing insights into their mechanistic action at the molecular level [17]. By applying molecular docking alongside network toxicology and bioinformatics, we can gain a comprehensive understanding of how Hexaconazole and its counterparts exert their toxic effects.\u003c/p\u003e \u003cp\u003eThis study aims to efficiently analyze the toxicity and mechanisms of Hexaconazole, Propiconazole and Prothioconazole through an integrated approach that combines network toxicology, molecular docking, and bioinformatics data. By elucidating the key pathways and molecular interactions involved, this research will contribute to a better understanding of the risks associated with triazole fungicide use and inform future regulatory decisions and risk assessment strategies. Ultimately, our findings aim to promote safer agricultural practices and minimize the impact on human health of these widely used agrochemicals.\u003c/p\u003e"},{"header":"2. Materials and methods","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003e2.1 Initial network toxicology predictive analysis of HA, PP and PT\u003c/h2\u003e \u003cp\u003eThe terms \u0026ldquo;Hexaconazole (HA)\u0026rdquo;, \u0026ldquo;Propiconazole (PP)\u0026rdquo; and \u0026ldquo;Prothioconazole (PT)\u0026rdquo; were queried in the PubChem database (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://pubchem.ncbi.nlm.nih.gov/\u003c/span\u003e\u003cspan address=\"https://pubchem.ncbi.nlm.nih.gov/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e) to obtain the standardized SMILES codes for these compounds. The SMILES codes of the three triazole compounds were then input into the ADMETlab 3.0 analysis tool (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://admetlab3.scbdd.com/\u003c/span\u003e\u003cspan address=\"https://admetlab3.scbdd.com/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e). Utilizing the structural modeling capabilities of this platform, we predicted the toxicological properties of HA, PP, and PT, providing an initial understanding of their toxicological profiles. The predicted results were downloaded.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec4\" class=\"Section2\"\u003e \u003ch2\u003e2.2 Targets collection of HA, PP and PT\u003c/h2\u003e \u003cp\u003eFirst, a search was conducted in the PubChem database for \u0026ldquo;Hexaconazole\u0026rdquo;, \u0026ldquo;Propiconazole\u0026rdquo; and \u0026ldquo;Prothioconazole\u0026rdquo; to identify the best matches and verify the accuracy of their names and molecular formulas. Subsequently, the 2D structure SDF files were downloaded. The 2D structures are shown in Supplementary Fig. \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e, with CAS numbers 103556-63-4, 60207-90-1, and 178928-70-6, respectively. Based on the SDF files, potential targets of \"Homo sapiens\" specific HA, PP,PT were searched from ChEMBL (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.ebi.ac.uk/chembl/\u003c/span\u003e\u003cspan address=\"https://www.ebi.ac.uk/chembl/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e) and Swiss Target Prediction database (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://www.swisstargetprediction.ch/\u003c/span\u003e\u003cspan address=\"http://www.swisstargetprediction.ch/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e) [18]. In addition, the Swiss-Target-Prediction tool is used to identify any ignored targets. The structural information throughout the search results is carefully cross-verified for consistency. Then open the Excel table and select the target with probability\u0026thinsp;\u0026gt;\u0026thinsp;0 as the target of HA, PP and PT. Finally, integration, de-duplication and gene naming coordination were carried out by UniProt database (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.uniprot.org/\u003c/span\u003e\u003cspan address=\"https://www.uniprot.org/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e) to establish the target library of HA, PP and PT.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec5\" class=\"Section2\"\u003e \u003ch2\u003e2.3 Collection of targets related to brain and liver injury\u003c/h2\u003e \u003cp\u003eWe conducted a comprehensive search of both domestic and international literature, as well as utilizing the GeneCards (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.genecards.org/\u003c/span\u003e\u003cspan address=\"https://www.genecards.org/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e) and OMIM (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://omim.org/\u003c/span\u003e\u003cspan address=\"https://omim.org/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e) databases, using the keywords \u0026ldquo;brain injury\u0026rdquo; and \u0026ldquo;liver injury\u0026rdquo; to identify relevant targets. To ensure that the genes obtained were highly associated with brain and liver injuries, we set the \u0026ldquo;score\u0026rdquo; threshold at the median value, selecting genes with scores above the median to establish a library of neurotoxic targets. We considered the intersection of targets related to HA, PP, and PT with brain and liver injury targets as potential targets for pesticide-induced brain or liver damage. Venn diagrams were employed to filter the common potential targets.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec6\" class=\"Section2\"\u003e \u003ch2\u003e2.4 Protein-protein interaction network construction and core targets screening\u003c/h2\u003e \u003cp\u003eThe intersection genes representing potential targets for brain or liver injury induced by the three pesticides were individually input into the STRING database (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://cn.string-db.org/\u003c/span\u003e\u003cspan address=\"https://cn.string-db.org/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e). We limited the species to \u0026ldquo;Homo sapiens\u0026rdquo;, set the \u0026ldquo;minimum required interaction score (FDR)\u0026rdquo; to \u0026ldquo;medium confidence\u0026thinsp;\u0026gt;\u0026thinsp;0.4\u0026rdquo; and selected \u0026ldquo;FDR strict value\u0026thinsp;\u0026le;\u0026thinsp;0.05\u0026rdquo; for analysis. These parameters ensured a robust analysis of the active target proteins corresponding to our target genes.\u003c/p\u003e \u003cp\u003eSubsequently, the results generated by STRING were imported into the network biology visualization application Cytoscape (version 3.10.0), which computes parameters for each node in the network diagram and displays molecular connections [19]. This allowed for the calculation of the topological properties of the network nodes and edges, resulting in a protein-protein interaction (PPI) network diagram. The criteria for screening core targets were as follows: nodes corresponding to targets that simultaneously met the conditions of betweenness centrality (BC)\u0026thinsp;\u0026gt;\u0026thinsp;median, closeness centrality (CC)\u0026thinsp;\u0026gt;\u0026thinsp;median, and degree centrality (DC)\u0026thinsp;\u0026gt;\u0026thinsp;median were selected as core targets for brain or liver injury induced by HA, PP and PT. Additionally, we utilized the MCODE plugin to validate the core genes within these PPI networks.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec7\" class=\"Section2\"\u003e \u003ch2\u003e\u003cem\u003e2.5 Functional enrichment and pathway analysis of targets for diseases induced by HA, PP and PT\u003c/em\u003e\u003c/h2\u003e \u003cdiv id=\"Sec8\" class=\"Section3\"\u003e \u003ch2\u003e2.5.1 GO and KEGG analysis of potential targets\u003c/h2\u003e \u003cp\u003eTo analyze the enrichment of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, we imported the targets associated with brain or liver damage induced by HA, PP, and PT into the Metascape database (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://metascape.org/gp/index.html\u003c/span\u003e\u003cspan address=\"https://metascape.org/gp/index.html\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e). This database serves as a comprehensive resource for gene functions, encompassing extensive information on biological processes (BP), cellular components (CC), and molecular functions (MF) [20]. The selection criteria included a minimum overlap of 3 and a minimum enrichment factor of 1.5, with a significance level set at \u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.01. The overlapping targets were uploaded to the Metascape database, where human species were selected, and the results were submitted. The top 10 GO terms and the top 20 KEGG pathways were selected based on ascending p-values. GO and KEGG diagrams were generated using the Microbiotics website (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://www.bioinformatics.com.cn/\u003c/span\u003e\u003cspan address=\"http://www.bioinformatics.com.cn/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e) and R software. The results of the bioinformatics analyses produced online KEGG pathway enrichment maps and GO term enrichment maps to visualize the interactions among pathways related to brain or liver damage induced by the three pesticides.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec9\" class=\"Section3\"\u003e \u003ch2\u003e2.5.2 KEGG analysis of core targets\u003c/h2\u003e \u003cp\u003eWe performed KEGG pathway enrichment analysis of the core targets associated with brain and liver damage induced by HA, PP, and PT using the Metascape database. The enrichment results were visualized to effectively interpret and present our findings. This multifaceted approach aims to comprehensively investigate the signaling pathways and biological processes mediated by the core targets related to brain and liver damage induced by these three pesticides, thereby elucidating and emphasizing the key mechanisms involved in this context.\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e \u003cdiv id=\"Sec10\" class=\"Section2\"\u003e \u003ch2\u003e2.6 Drug-target-pathway network of HA, PP and PT\u003c/h2\u003e \u003cp\u003eThe drug-target-pathway network provides a clear overview of how drugs affect pathways through specific targets and illustrates the relationships between these key targets and pathways, allowing for the identification of the most critical targets and pathways. First, an attribute table was constructed in Excel, containing the core targets associated with the three triazole pesticides and all pathways corresponding to each core target. Next, Cytoscape was opened, and the attribute table was uploaded in the Network file to load the locations, set the starting points, and visualize the results.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec11\" class=\"Section2\"\u003e \u003ch2\u003e2.7 Molecular docking of key targets with HA, PP and PT\u003c/h2\u003e \u003cp\u003eTo investigate the molecular interactions and binding patterns between HA, PP, PT, and the identified key target proteins, we conducted molecular docking simulations. This theoretical structure-based approach enabled the predictive modeling of receptor-ligand binding geometries and affinities. The crystal structures of the core target proteins were obtained from the RCSB PDB database (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.rcsb.org/\u003c/span\u003e\u003cspan address=\"https://www.rcsb.org/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e) and preprocessed using PyMOL 3.0 (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.pymol.org/\u003c/span\u003e\u003cspan address=\"https://www.pymol.org/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e) by removing solvent water molecules and any pre-existing ligands bound to the target proteins. The processed protein structures were then imported into AutoDock Tools 1.5.6 (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://autodock.scripps.edu/\u003c/span\u003e\u003cspan address=\"https://autodock.scripps.edu/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e) for hydrogen addition, charge calculation, and non-polar hydrogen merging to generate receptor models. The structures of HA, PP, and PT compounds were downloaded from the PubChem database. The mechanical structures were optimized using Chem3D (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://revvitysignals.com/products/research/chemdraw\u003c/span\u003e\u003cspan address=\"https://revvitysignals.com/products/research/chemdraw\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e) to minimize energy and designated as ligands. The docking range was determined by adjusting the center coordinates and grid box size parameters. Following this, a genetic algorithm was selected, and molecular docking simulations were executed via command line using AutoDock Vina. Finally, the predicted binding modes and affinities were visualized and analyzed using PyMOL 3.0.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec12\" class=\"Section2\"\u003e \u003ch2\u003e2.8 Collection and analysis of bioinformatics data\u003c/h2\u003e \u003cdiv id=\"Sec13\" class=\"Section3\"\u003e \u003ch2\u003e2.8.1 Identification and validation of brain and liver cancer-related DEGs\u003c/h2\u003e \u003cp\u003eTo further elucidate the mechanistic impact of central targets associated with brain and liver damage induced by three pesticides on cancer, we obtained transcriptomic datasets from the NCBI GEO database (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.ncbi.nlm.nih.gov/geo/\u003c/span\u003e\u003cspan address=\"https://www.ncbi.nlm.nih.gov/geo/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e), selecting GSE30563 and GSE117361 for detailed analysis. The GEOexplorer tool (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://geoexplorer.rosalind.kcl.ac.uk/\u003c/span\u003e\u003cspan address=\"https://geoexplorer.rosalind.kcl.ac.uk/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e) was employed to determine differentially expressed genes (DEGs) between brain and liver tumors with normal tissue samples. Multichip analysis was conducted using the SVA and limma packages in R software (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.r-project.org/\u003c/span\u003e\u003cspan address=\"https://www.r-project.org/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e), allowing for batch correction and calculation of fold-change (FC) and adjusted \u003cem\u003ep\u003c/em\u003e-values. The thresholds for identifying DEGs were set at |log (FC)| \u0026gt; 1 and \u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.05. Furthermore, the overlapping DEGs filtered from both datasets were classified as cancer-related targets.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec14\" class=\"Section3\"\u003e \u003ch2\u003e2.8.2 Analysis of brain and liver cancer-related target genes\u003c/h2\u003e \u003cp\u003eEnrichment analysis (GO and KEGG) of differentially expressed genes (DEGs), commonly upregulated genes, and commonly downregulated genes was performed using the Metascape database (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://metascape.org/gp/index.html\u003c/span\u003e\u003cspan address=\"https://metascape.org/gp/index.html\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e). The selection criteria included a minimum overlap of 3, a minimum enrichment factor of 1.5, and a significance level of \u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.01.\u003c/p\u003e \u003cp\u003eAdditionally, we employed the Gene Expression Profiling Interactive Analysis (GEPIA) tool (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://gepia.cancer-pku.cn/\u003c/span\u003e\u003cspan address=\"http://gepia.cancer-pku.cn/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e) to perform survival analysis on five key expressed genes. Our objective was to investigate the genetic alterations affecting each target gene and to test hypotheses regarding the recurrence patterns of genetic alteration events and their genomic context in brain tumors and liver tumors [21].\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e"},{"header":"3. Results","content":"\u003cdiv id=\"Sec16\" class=\"Section2\"\u003e \u003ch2\u003e3.1 Toxicity analysis of HA, PP and PT\u003c/h2\u003e \u003cp\u003ePredictions from the ADMETlab 3.0 database revealed that the co-toxicity of HA, PP, and PT is primarily associated with neurotoxicity and liver damage (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). Additionally, the database provided predictions regarding the genotoxicity, respiratory toxicity, and skin irritation potential of these three pesticides. A comprehensive review of both the database and relevant literature was conducted to predict the molecular mechanisms underlying these two diseases.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eThe prediction toxicity of Hexaconazole, Propiconazole and Prothioconazole.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"3\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eIngredient\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eProperty\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eProbability\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eHexaconazole\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eDrug-induced Neurotoxicity\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.875 (++)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eHexaconazole\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eHuman Hepatotoxicity\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.600 (+)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eHexaconazole\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eGenotoxicity\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.564 (+)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eHexaconazole\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eRespiratory toxicity\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.749 (++)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eHexaconazole\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSkin\u0026nbsp;sensitization\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.710 (++)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003ePropiconazole\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eDrug-induced Neurotoxicity\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.907 (+++)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003ePropiconazole\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eHuman Hepatotoxicity\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.685 (+)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003ePropiconazole\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eGenotoxicity\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.872 (++)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003ePropiconazole\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eRespiratory toxicity\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.672 (+)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003ePropiconazole\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSkin\u0026nbsp;sensitization\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.658 (+)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eProthioconazole\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eDrug-induced Neurotoxicity\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.821 (++)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eProthioconazole\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eHuman Hepatotoxicity\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.477 (-)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eProthioconazole\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eGenotoxicity\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.999 (+++)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eProthioconazole\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eRespiratory toxicity\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.723 (++)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eProthioconazole\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSkin\u0026nbsp;sensitization\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.505 (+)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003ctfoot\u003e \u003ctr\u003e\u003ctd colspan=\"3\"\u003eClassification endpoints of the prediction probability values: 0-0.1 (---), 0.1\u0026ndash;0.3 (--), 0.3\u0026ndash;0.5 (-), 0.5\u0026ndash;0.7 (+), 0.7\u0026ndash;0.9 (++), 0.9\u0026ndash;1.0 (+++). Larger probability represents a greater risk of toxicity.\u003c/td\u003e\u003c/tr\u003e \u003c/tfoot\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec17\" class=\"Section2\"\u003e \u003ch2\u003e3.2 Targets analysis of HA, PP and PT induce diseases\u003c/h2\u003e \u003cp\u003eWe initially integrated data from the ChEMBL and Swiss-Target-Prediction databases, removing duplicate entries to identify 120 HA targets, 32 PP targets, and 118 PT targets. Following the removal of duplicates, we also identified 5,464 targets highly associated with brain injury and 5,071 targets highly associated with liver injury from the GeneCards and OMIM databases. The list of these target genes can be found in Appendix A.\u003c/p\u003e \u003cp\u003eUsing a comprehensive Venn diagram approach, we identified a total of 94 HA cross-targets related to brain injury (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eA), 28 PP cross-targets related to brain injury (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eB), 84 PT cross-targets related to brain injury (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eC), 92 HA cross-targets related to liver injury (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eD), 24 PP cross-targets related to liver injury (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eE), and 84 PT cross-targets related to liver injury (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eF). These targets are considered potential candidates associated with brain or liver injuries induced by the three triazole pesticides.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec18\" class=\"Section2\"\u003e \u003ch2\u003e3.3 PPI network and core target acquisition\u003c/h2\u003e \u003cp\u003ePotential targets associated with brain or liver injury induced by the three triazole pesticides were mapped to the protein-protein interaction (PPI) network using the STRING database (Supplementary Fig. \u003cspan refid=\"MOESM2\" class=\"InternalRef\"\u003eS2\u003c/span\u003e), resulting in a highly interconnected network. Metrics such as the number of nodes, number of edges, average node degree, and PPI enrichment \u003cem\u003ep\u003c/em\u003e-value are detailed in Appendix B. The topological analysis of the network, including degree and betweenness centrality metrics, was performed using Cytoscape. Additionally, we constructed a PPI network diagram to visually represent the results (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eThen, we employed multiple complementary algorithms, specifically the Betweenness Centrality (BC), Closeness Centrality (CC), and Degree Centrality (DC) algorithms from the CytoHubba plugin, to identify highly connected hub targets. The targets identified as hubs resulted in the selection of a subset of 35 target genes closely related to HA-induced brain injury, 34 target genes closely related to HA-induced liver injury, 9 target genes closely related to PP-induced brain injury, 8 target genes closely related to PP-induced liver injury, 26 target genes closely related to PT-induced brain injury, and 26 target genes closely related to PT-induced liver injury (Supplementary Tables S1-S6). A centralized PPI network diagram was constructed to illustrate the interactions among these hub targets (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec19\" class=\"Section2\"\u003e \u003ch2\u003e3.4 GO analysis\u003c/h2\u003e \u003cp\u003eGene Ontology (GO) analysis was conducted on potential targets associated with brain or liver injuries induced by three types of fungicides using the Metascape database. Our findings revealed significant GO terms (Appendix C). The top ten most enriched terms across the Biological Process (BP), Cellular Component (CC), and Molecular Function (MF) categories were prioritized based on gene count and are visually represented in the enrichment analysis plot (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e). In terms of Biological Process (BP), the potential targets associated with brain and liver injuries induced by the three triazole fungicides exhibited significant enrichment in pathways related to cellular responses to organic cyclic compounds and organic nitrogen compounds. Regarding CC, the potential targets showed significant enrichment in protein kinase binding pathways for HA and PT, while for PP, they were significantly enriched in pathways related to DNA and organic compound binding. With respect to MF, the potential targets demonstrated significant enrichment in receptor complex pathways.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec20\" class=\"Section2\"\u003e \u003ch2\u003e3.5 KEGG pathway analysis\u003c/h2\u003e \u003cp\u003eKEGG pathway analysis was conducted using the Metascape database on potential targets associated with brain or liver injuries induced by three types of fungicides, revealing enriched pathways related to these targets (Appendix C). Overall, these targets were significantly enriched in pathways associated with cancer, the PI3K-Akt signaling pathway, and pathways related to chemical carcinogenesis. We generated a bubble plot (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e) displaying the top 20 most enriched KEGG signaling pathways, prioritized by gene count, to succinctly illustrate the interactions among these pathways and to highlight their significance in the pathogenesis and progression of brain and liver injuries caused by triazole fungicides.\u003c/p\u003e \u003cp\u003eAdditionally, we performed KEGG analysis on the core targets associated with brain or liver injuries induced by the three fungicides using the Metascape database. This analysis produced a Sankey bubble plot (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e) illustrating the top 20 most enriched KEGG signaling pathways, also prioritized by gene count. These targets showed significant enrichment in pathways related to cancer, the PI3K-Akt signaling pathway, and protein processing.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec21\" class=\"Section2\"\u003e \u003ch2\u003e3.6 Network analysis of diseases induced by HA, PP and PT\u003c/h2\u003e \u003cp\u003eUsing Cytoscape software, we constructed a target network diagram for HA, PP, and PT-induced brain and liver injuries, resulting in a \"drug-target pathway\" network diagram. As shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e, deep red or deep blue represents hub targets, while light red or light blue denotes pathways. This diagram provides insights into the core targets, pathways, and the extent of brain or liver injuries induced by the three pesticides.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec22\" class=\"Section2\"\u003e \u003ch2\u003e3.7 Molecular docking analysis of HA, PP and PT and key targets\u003c/h2\u003e \u003cp\u003eTo further elucidate the interactions and relationships between HA, PP, and PT with the key targets (the top three core targets), we conducted comprehensive molecular docking simulations. Docking of each target was performed using AutoDock software, and the results indicated that all binding energies were below \u0026minus;\u0026thinsp;5.0 kcal/mol [22], demonstrating a strong binding affinity of HA, PP, and PT for these targets (Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e). Our findings suggest that HA, PP, and PT can spontaneously bind to each key target, underscoring their critical roles in the molecular mechanisms underlying brain and liver injuries induced by these three pesticides. Visual representations of the lowest binding energy conformations between these pesticides and each key target were generated using Discovery Studio and PyMOL (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003e). These results indicate that these three triazole pesticides may induce brain or liver injury through binding with the aforementioned targets. However, further experimental validation is required.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab2\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 2\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eThe result of molecular docking.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"3\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eLigand\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eTarget\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eBinding energy(kcal/mol)\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eHexaconazole\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eCASP3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-6.5\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eHexaconazole\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eEGFR\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-6.0\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eHexaconazole\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eESR1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-7.2\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003ePropiconazole\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003ePPARG\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-7.3\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003ePropiconazole\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eTP53\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-7.1\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003ePropiconazole\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eHSP90AA1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-6.6\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eProthioconazole\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eEGFR\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-6.4\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eProthioconazole\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eESR1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-7.4\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eProthioconazole\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003ePTGS2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-7.3\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec23\" class=\"Section2\"\u003e \u003ch2\u003e\u003cem\u003e3.8 Screening for brain and liver cancer-related DEGs and analysis of target genes\u003c/em\u003e\u003c/h2\u003e \u003cp\u003eBuilding on our previous KEGG analysis, we revealed that the targets associated with brain and liver toxicity induced by the three triazole fungicides are enriched in pathways related to cancer. Therefore, we aimed to further elucidate the mechanistic impact of core targets for triazole fungicides-elicited brain and liver injury on brain and liver cancer.\u003c/p\u003e \u003cp\u003eTo identify genes associated with brain tumors and liver tumors, we utilized two gene expression profiles: GSE30563 and GSE117361. To visually highlight the unique distribution and significance of these differentially expressed genes (DEGs), we employed volcano plots to visualize all DEGs across the two datasets (Figs.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003eA-B). Detailed information regarding these datasets is provided in Appendix D. Venn diagrams were constructed to illustrate the overlap between the two sets of DEGs, revealing a total of 287 overlapping genes, of which 283 exhibited common expression patterns. Among these, 143 genes were significantly upregulated, while 140 genes were significantly downregulated (Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003eC). We conducted GO and KEGG enrichment analyses on the 287 overlapping gene targets to explore the differential expression patterns in brain tumors and liver tumors compared to normal tissues. Detailed information regarding the enrichment results is provided in Appendix E. The top ten most enriched terms in BP, CC, and MF categories were prioritized based on gene count and visually represented in a histogram (Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003eD). Then we generated a histogram of the top 20 most enriched KEGG signaling pathways, prioritized by gene count (Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003eE). To further investigate the functional roles of the DEGs that were commonly upregulated or downregulated, we performed GO and KEGG enrichment analyses on the 283 genes exhibiting common expression patterns. The results of the GO and KEGG enrichment analyses are visualized in Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003eF-G. In GO enrichment analysis,DEGs was significantly related to immune response and kinase binding pathways. In KEGG enrichment analysis, DEGs was significantly associated with immune disorders and infection-related pathways. This suggests that DEGs is highly correlated with immune responses in brain and liver tumors. Additionally, we selected key DEGs with a p-value less than 0.005, consisting of 15 genes with a common expression pattern. Among these, 5 genes were commonly upregulated, and 10 genes were commonly downregulated (Supplementary Fig. \u003cspan refid=\"MOESM4\" class=\"InternalRef\"\u003eS4\u003c/span\u003e-S5). Kaplan-Meier analysis was performed on the 5 commonly upregulated genes, and the results were visualized (Supplementary Fig. \u003cspan refid=\"MOESM6\" class=\"InternalRef\"\u003eS6\u003c/span\u003e). We observed a significant difference in survival times between the high-expression and low-expression groups for these 5 genes.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e"},{"header":"4. Discussion","content":"\u003cp\u003eThe adoption of triazole fungicides has significantly increased in modern agriculture due to their high efficacy, low application rates, and relatively favorable environmental profiles compared to traditional classes of fungicides. However, the triazole fungicides have been detected in a variety of food products, beverages, and human biomonitoring studies. Given the widespread human exposure to these compounds, comprehensive investigations are needed to evaluate their potential safety hazards [23, 24].\u003c/p\u003e \u003cp\u003eIn our study, we integrated a comprehensive approach utilizing network toxicology, molecular docking, and multi-layered bioinformatics to explore the effects of exposure to the triazole fungicides HA, PP, and PT on the brain and liver, as well as their potential mechanisms of action. After applying network-based computational assessment tools, we identified potential targets associated with brain or liver damage induced by HA, PP, and PT using the ChEMBL, Swiss Target Prediction, GeneCards, OMIM, and DisGeNet databases. Utilizing the STRING database and Cytoscape software, we constructed an interaction network for these potential targets and extracted core targets, including \u003cem\u003eEGFR\u003c/em\u003e, \u003cem\u003eCASP3\u003c/em\u003e, \u003cem\u003eESR1\u003c/em\u003e, \u003cem\u003ePPARG\u003c/em\u003e, \u003cem\u003eTP53\u003c/em\u003e, \u003cem\u003eHSP90AA1\u003c/em\u003e, and \u003cem\u003ePTGS2\u003c/em\u003e. These targets were designated as critical for the induction of brain or liver damage by the three triazole fungicides.\u003c/p\u003e \u003cp\u003eEpidermal growth factor receptor (EGFR), a transmembrane receptor tyrosine kinase, plays an important role in maintaining normal tissue and cell signaling pathways [25]. Dysregulation of \u003cem\u003eEGFR\u003c/em\u003e has been implicated in neurodegenerative diseases, where altered signaling pathways contribute to neuronal apoptosis and cognitive decline [26]. \u003cem\u003eEGFR\u003c/em\u003e signaling plays a central role in the regenerative response after liver injury and is involved in cell transformation associated with chronic injury [27]. In addition, studies have shown that \u003cem\u003eEGFR\u003c/em\u003e is also associated with apoptotic cell death signaling in various hepatocytes, mitochondrial dysfunction, and acute liver necrosis [28].\u003c/p\u003e \u003cp\u003eCaspase-3 (CASP3), a widely expressed member of the conserved protein family, plays a key role in regulating the growth and homeostasis maintenance of normal and malignant cells and tissues in multicellular organisms [29]. Activation of \u003cem\u003eCASP3\u003c/em\u003e is strongly associated with neurodegenerative diseases and traumatic brain injury, and its upregulation leads to neuronal apoptosis, exacerbating cognitive decline and dysfunction [30, 31]. \u003cem\u003eCASP3\u003c/em\u003e mediates hepatocyte apoptosis in response to hepatotoxic substances, leading to liver damage and fibrosis, cirrhosis and other conditions [32]. The study suggests that modulating \u003cem\u003eCASP3\u003c/em\u003e activity may provide a therapeutic pathway for alleviating cell death in brain and liver injury [33].\u003c/p\u003e \u003cp\u003eEstrogen receptor alpha (ESR1) is a key regulator that regulates the effects of estrogen on a variety of tissues, including the brain and liver [34]. Studies have shown that ESR1 signaling has a neuroprotective effect, as estrogen can enhance neuronal survival and reduce apoptosis after trauma or ischemic injury [35]. And \u003cem\u003eESR1\u003c/em\u003e activation is associated with improved cognitive function and reduced neuroinflammation [36]. In liver injury, \u003cem\u003eESR1\u003c/em\u003e plays a complex role in promoting liver cell regeneration and preventing oxidative stress [37]. However, dysregulated \u003cem\u003eESR1\u003c/em\u003e signaling can contribute to the progression of liver diseases, including fatty liver and fibrosis [38].\u003c/p\u003e \u003cp\u003ePeroxisome proliferator activating receptor (PPARG) is a nuclear receptor that regulates genes involved in lipid metabolism, inflammation, and insulin sensitivity [39]. Activation of \u003cem\u003ePPARG\u003c/em\u003e has been shown to have neuroprotective effects, enhancing neuronal survival and reducing neuroinflammation after trauma or ischemic injury [40]. Studies have shown that activating \u003cem\u003ePPARG\u003c/em\u003e can improve cognitive outcomes and promote repair mechanisms in neurodegenerative diseases [41]. In addition, \u003cem\u003ePPARG\u003c/em\u003e plays a key role in lipid homeostasis and inflammation resolution [42]. Activation of \u003cem\u003ePPARG\u003c/em\u003e is associated with a reduction in steatosis and fibrosis in various liver diseases [43].\u003c/p\u003e \u003cp\u003eThe \u003cem\u003eTP53\u003c/em\u003e gene encodes the tumor suppressor protein p53, which plays a crucial role in cellular stress responses, including apoptosis, cell cycle regulation, and DNA repair [44]. In brain injury, activation of \u003cem\u003eTP53\u003c/em\u003e is often associated with neurodegenerative diseases, and in response to oxidative stress and DNA damage, it mediates neuronal apoptosis, leading to cognitive decline and neuron loss [45]. In liver injury, p53 can promote apoptosis of damaged liver cells and prevent tumorigenesis, but excessive activation can lead to chronic inflammation and fibrosis [46].\u003c/p\u003e \u003cp\u003eThe \u003cem\u003eHSP90AA1\u003c/em\u003e gene encodes heat shock protein 90α (Hsp90), an important companion involved in protein folding and cellular stress response [47]. Previous studies have shown that Hsp90 exhibits neuroprotective properties, promoting the stabilization of client proteins and preventing neuronal cell apoptosis during ischemic or traumatic events [48]. On the other hand, Hsp90 plays an important role in maintaining the integrity of hepatocytes under stress conditions such as oxidative damage and inflammation [49].\u003c/p\u003e \u003cp\u003eThe enzyme protein prostaglandin-endoperoxide synthase 2 (PTGS2) is involved in inflammation and other physiological processes by catalyzing the conversion of arachidonic acid to prostaglandin [50]. For brain injury, COX-2 expression is upregulated in response to neuronal injury, leading to neuroinflammation and pain [51]. While COX-2 is involved in protective mechanisms that promote healing and regeneration, overactivation can exacerbate neuronal damage and lead to chronic neurodegenerative diseases [51]. For liver damage, PTGS2 is involved in the inflammatory response associated with liver disease, including non-alcoholic fatty liver disease and hepatitis [52]. Elevated COX-2 levels are associated with increased inflammation and apoptosis of liver cells, which worsens liver damage [53].\u003c/p\u003e \u003cp\u003ePathway analysis highlights key insights into the molecular mechanisms of brain and liver injury induced by triazole fungicides. The significant enrichment of targets in cancer-related pathways suggests that exposure to these fungicides may promote carcinogenic processes in these organs, raising concerns about long-term health effects [54]. Specifically, the PI3K-Akt signaling pathway plays an important role in cell survival, proliferation, and metabolism [55]. Dysregulation of this pathway may lead to neuroinflammatory responses and hepatotoxicity, which may exacerbate tissue damage. In addition, the identification of pathways associated with chemical carcinogenesis highlights the potential of these fungicides to initiate or enhance malignant transformation.\u003c/p\u003e"},{"header":"5. Conclusion","content":"\u003cp\u003eIn summary, the present study combined applied network toxicology, molecular docking, and multilevel bioinformatics data to systematically explore the toxicological and molecular mechanisms of three triazole fungicides, including Hexaconazole, Propiconazole, and Prothioconazole, in brain and liver injury. Our findings demonstrate that the toxicity profiles of these compounds are closely linked to their interactions with specific biological pathways, highlighting the potential risks they pose to human health and the environment. This integrated approach not only enhances our understanding of fungicide-induced toxicity but also offers valuable insights for future research aimed at developing safer agricultural practices and regulatory policies.\u003c/p\u003e"},{"header":"Abbreviations","content":"\u003cp\u003eHA, Hexaconazole; PP, Propiconazole; PT, Prothioconazole; DEGs, differential gene expression; BP, biological processes; CC, cellular components; MF, molecular functions; GP: Genetic Information Processing; EI: Environmental Information Processing; CP: Cellular Processes; OR: Organismal Systems; HD: Human Diseases.\u003c/p\u003e\n\u003cp\u003e\u003cbr\u003e\u003c/p\u003e"},{"header":"Declarations","content":"\u003ch2\u003eFunding:\u003c/h2\u003e \u003cp\u003eThis research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.\u003c/p\u003e\u003ch2\u003eAuthor Contribution\u003c/h2\u003e\u003cp\u003eAuthor contributions: Qing Lv: Methodology, Conceptualization, Data analysis, Visualization, Writing \u0026ndash; original draft. Xingke Zhu: Methodology, Conceptualization, Data analysis, Visualization, Writing \u0026ndash; original draft, Writing \u0026ndash; review \u0026amp; editing.Author Statement: By submitting this manuscript, we, the authors, affirm that we have contributed significantly to the research presented and are in agreement with the content and conclusions of the work. We confirm that this manuscript has not been published elsewhere and is not under consideration by any other journal. All authors have reviewed and approved the final version of the manuscript. We declare that there are no conflicts of interest related to this submission, and we have adhered to ethical standards in conducting our research. We agree to abide by the policies and procedures of the journal regarding the publication process.\u003c/p\u003e\u003ch2\u003eAcknowledgments:\u003c/h2\u003e \u003cp\u003eThe author declares no financial conflicts of interest or financial relationships with any organizations or individuals that could potentially bias the results or interpretation of this study. This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.\u003c/p\u003e \u003cp\u003e \u003cstrong\u003eEthics approval and consent to participate:\u003c/strong\u003e \u003cp\u003eNot applicable.\u003c/p\u003e \u003c/p\u003e \u003cp\u003e \u003cstrong\u003eRunning Title\u003c/strong\u003e \u003cp\u003eToxicity and mechanisms of three triazole fungicides.\u003c/p\u003e \u003c/p\u003e\u003ch2\u003eAvailability of data and materials\u003c/h2\u003e \u003cp\u003eData will be made available on request. Further inquiries can be directed to the corresponding authors.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003e Fern\u0026aacute;ndez-Vizca\u0026iacute;no E, Mateo R, Fern\u0026aacute;ndez de Mera IG, Mougeot F, Camarero PR, Ortiz-Santaliestra ME. Transgenerational effects of triazole fungicides on gene expression and egg compounds in non-exposed offspring: A case study using Red-Legged Partridges (Alectoris rufa). \u003cem\u003eSci Total Environ\u003c/em\u003e. 2024;926:171546. doi:10.1016/j.scitotenv.2024.171546\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e He Z, Zhang J, Shi D, et al. Deoxynivalenol in \u003cem\u003eFusarium graminearum\u003c/em\u003e: Evaluation of Cyproconazole Stereoisomers \u003cem\u003eIn Vitro\u003c/em\u003e and \u003cem\u003eIn Planta\u003c/em\u003e. \u003cem\u003eJ Agric Food Chem\u003c/em\u003e. 2021;69(34):9735\u0026ndash;9742. doi:10.1021/acs.jafc.1c02555\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Raina R, Smith E. Determination of azole fungicides in atmospheric samples collected in the Canadian prairies by LC/MS/MS. \u003cem\u003eJ AOAC Int\u003c/em\u003e. 2012;95(5):1350\u0026ndash;1356. doi:10.5740/jaoacint.sge_raina\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Mao CX, Luo J, Zhang Y, Zhang CQ. Targeted deletion of three CYP51s in Fusarium fujikuroi and their different roles in determining sensitivity to 14α-demethylase inhibitor fungicides. \u003cem\u003ePest Manag Sci\u003c/em\u003e. 2023;79(4):1324\u0026ndash;1330. doi:10.1002/ps.7304\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Zou PC, Zhang Y, Bian Y, et al. Triazoles in the environment: An update on sample pretreatment and analysis methods. \u003cem\u003eEcotoxicol Environ Saf\u003c/em\u003e. Published online October 8, 2024. doi:10.1016/j.ecoenv.2024.117156\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Yao X, Lv H, Wang Q, et al. Novel Insights into Stereoselective Reproductive Toxicity Induced by Mefentrifluconazole in Earthworms (\u003cem\u003eEisenia fetida\u003c/em\u003e): First Report of Estrogenic Effects. \u003cem\u003eJ Agric Food Chem\u003c/em\u003e. 2024;72(35):19304\u0026ndash;19311. doi:10.1021/acs.jafc.4c04168\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Zhang Y, Guo J, Tang C, Xu K, Li Z, Wang C. Early life stage exposure to fenbuconazole causes multigenerational cardiac developmental defects in zebrafish and potential reasons. \u003cem\u003eEnviron Pollut\u003c/em\u003e. 2024;349:123938. doi:10.1016/j.envpol.2024.123938\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Cang T, Lou Y, Zhu YC, et al. Mixture toxicities of tetrachlorantraniliprole and tebuconazole to honey bees (Apis mellifera L.) and the potential mechanism. \u003cem\u003eEnviron Int\u003c/em\u003e. 2023;172:107764. doi:10.1016/j.envint.2023.107764\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Wang X, Liu Y, Xue M, Wang Z, Yu J, Guo X. Enantioselective degradation of chiral fungicides triticonazole and prothioconazole in soils and their enantioselective accumulation in earthworms Eisenia fetida. \u003cem\u003eEcotoxicol Environ Saf\u003c/em\u003e. 2019;183:109491. doi:10.1016/j.ecoenv.2019.109491\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Wang X, Ni H, Xu W, et al. Difenoconazole induces oxidative DNA damage and mitochondria mediated apoptosis in SH-SY5Y cells. \u003cem\u003eChemosphere\u003c/em\u003e. 2021;283:131160. doi:10.1016/j.chemosphere.2021.131160\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Sanchez CL, Souders CL 2nd, Pena-Delgado CJ, et al. Neurotoxicity assessment of triazole fungicides on mitochondrial oxidative respiration and lipids in differentiated human SH-SY5Y neuroblastoma cells. \u003cem\u003eNeurotoxicology\u003c/em\u003e. 2020;80:76\u0026ndash;86. doi:10.1016/j.neuro.2020.06.009\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Yang Q, Deng P, Xing D, et al. Developmental Neurotoxicity of Difenoconazole in Zebrafish Embryos. \u003cem\u003eToxics\u003c/em\u003e. 2023;11(4):353. Published 2023 Apr 8. doi:10.3390/toxics11040353\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Zerdoug A, Le V\u0026eacute;e M, Le Mentec H, et al. Induction of drug metabolizing enzyme and drug transporter expression by antifungal triazole pesticides in human HepaSH hepatocytes. \u003cem\u003eChemosphere\u003c/em\u003e. Published online October 5, 2024. doi:10.1016/j.chemosphere.2024.143474\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Li C, Yuan S, Zhou Y, et al. Microplastics reduce the bioaccumulation and oxidative stress damage of triazole fungicides in fish. \u003cem\u003eSci Total Environ\u003c/em\u003e. 2022;806(Pt 4):151475. doi:10.1016/j.scitotenv.2021.151475\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Zeng Z, Hu J, Xiao G, et al. Integrating network toxicology and molecular docking to explore the toxicity of the environmental pollutant butyl hydroxyanisole: An example of induction of chronic urticaria. \u003cem\u003eHeliyon\u003c/em\u003e. 2024;10(15):e35409. Published 2024 Jul 30. doi:10.1016/j.heliyon.2024.e35409\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Tong L, Zhou Y, Zhang Y, et al. Unveiling the protein-lipid interaction mechanism: How the sturgeon lipids diminish the surimi gel properties. \u003cem\u003eFood Chem\u003c/em\u003e. 2025;462:140926. doi:10.1016/j.foodchem.2024.140926\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Yang Z, Han S, Keller M, et al. Structural basis of ligand binding modes at the neuropeptide Y Y\u003csub\u003e1\u003c/sub\u003e receptor. \u003cem\u003eNature\u003c/em\u003e. 2018;556(7702):520\u0026ndash;524. doi:10.1038/s41586-018-0046-x\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Mendez D, Gaulton A, Bento AP, et al. ChEMBL: towards direct deposition of bioassay data. \u003cem\u003eNucleic Acids Res\u003c/em\u003e. 2019;47(D1):D930-D940. doi:10.1093/nar/gky1075\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Liu M, Wang Y, Deng W, et al. Combining network pharmacology, machine learning, molecular docking and molecular dynamic to explore the mechanism of Chufeng Qingpi decoction in treating schistosomiasis. \u003cem\u003eFront Cell Infect Microbiol\u003c/em\u003e. 2024;14:1453529. Published 2024 Sep 6. doi:10.3389/fcimb.2024.1453529\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Zhang M, Zhang X, Pei J, et al. Identification of phytochemical compounds of \u003cem\u003eFagopyrum dibotrys\u003c/em\u003e and their targets by metabolomics, network pharmacology and molecular docking studies. \u003cem\u003eHeliyon\u003c/em\u003e. 2023;9(3):e14029. doi:10.1016/j.heliyon.2023.e14029\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Shen SH, Chen SF, Guo JF, Wang ZX. The GGCT and REST positive feedback loop promotes tumor growth in Glioma. \u003cem\u003eTransl Oncol\u003c/em\u003e. 2024;49:102083. doi:10.1016/j.tranon.2024.102083\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Zhang Q, Liang J, Li X, et al. Exploring antithrombotic mechanisms and effective constituents of Lagopsis supina using an integrated strategy based on network pharmacology, molecular docking, metabolomics, and experimental verification in rats. \u003cem\u003eJ Ethnopharmacol\u003c/em\u003e. 2025;336:118717. doi:10.1016/j.jep.2024.118717\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Liu J, Nie Y, Niu Y, Li L, Jing X. Lignin-based emulsive liquid-liquid microextraction for detecting triazole fungicides in water, juice, vinegar, and alcoholic beverages via UHPLC-MS/MS. \u003cem\u003eFood Chem\u003c/em\u003e. 2024;459:140407. doi:10.1016/j.foodchem.2024.140407\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Marciano LPA, Kleinstreuer N, Chang X, Costa LF, Silv\u0026eacute;rio ACP, Martins I. A novel approach to triazole fungicides risk characterization: Bridging human biomonitoring and computational toxicology. \u003cem\u003eSci Total Environ\u003c/em\u003e. 2024;953:176003. doi:10.1016/j.scitotenv.2024.176003\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Damare R, Engle K, Kumar G. Targeting epidermal growth factor receptor and its downstream signaling pathways by natural products: A mechanistic insight. \u003cem\u003ePhytother Res\u003c/em\u003e. 2024;38(5):2406\u0026ndash;2447. doi:10.1002/ptr.8166\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Mansour HM, Fawzy HM, El-Khatib AS, Khattab MM. Repurposed anti-cancer epidermal growth factor receptor inhibitors: mechanisms of neuroprotective effects in Alzheimer's disease. Neural Regen Res. 2022;17(9):1913\u0026ndash;1918. doi:10.4103/1673-5374.332132\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Gonz\u0026aacute;lez L, D\u0026iacute;az ME, Miquet JG, Sotelo AI, Dominici FP. Growth Hormone Modulation of Hepatic Epidermal Growth Factor Receptor Signaling. \u003cem\u003eTrends Endocrinol Metab\u003c/em\u003e. 2021;32(6):403\u0026ndash;414. doi:10.1016/j.tem.2021.03.004\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Bhushan B, Michalopoulos GK. Role of epidermal growth factor receptor in liver injury and lipid metabolism: Emerging new roles for an old receptor. \u003cem\u003eChem Biol Interact\u003c/em\u003e. 2020;324:109090. doi:10.1016/j.cbi.2020.109090\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Eskandari E, Eaves CJ. Paradoxical roles of caspase-3 in regulating cell survival, proliferation, and tumorigenesis. \u003cem\u003eJ Cell Biol\u003c/em\u003e. 2022;221(6):e202201159. doi:10.1083/jcb.202201159\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Unnisa A, Greig NH, Kamal MA. Inhibition of Caspase 3 and Caspase 9 Mediated Apoptosis: A Multimodal Therapeutic Target in Traumatic Brain Injury. \u003cem\u003eCurr Neuropharmacol\u003c/em\u003e. 2023;21(4):1001\u0026ndash;1012. doi:10.2174/1570159X20666220327222921\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Rajesh Y, Kanneganti TD. Innate Immune Cell Death in Neuroinflammation and Alzheimer's Disease. \u003cem\u003eCells\u003c/em\u003e. 2022;11(12):1885. Published 2022 Jun 10. doi:10.3390/cells11121885\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Ma XY, Zhang M, Fang G, et al. Ursolic acid reduces hepatocellular apoptosis and alleviates alcohol-induced liver injury via irreversible inhibition of CASP3 in vivo. \u003cem\u003eActa Pharmacol Sin\u003c/em\u003e. 2021;42(7):1101\u0026ndash;1110. doi:10.1038/s41401-020-00534-y\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Wu J, Ye X, Yang S, Yu H, Zhong L, Gong Q. Systems Pharmacology Study of the Anti-Liver Injury Mechanism of Citri Reticulatae Pericarpium. \u003cem\u003eFront Pharmacol\u003c/em\u003e. 2021;12:618846. Published 2021 Apr 12. doi:10.3389/fphar.2021.618846\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Bartkowiak-Wieczorek J, Jaros A, Gajdzińska A, et al. The Dual Faces of Oestrogen: The Impact of Exogenous Oestrogen on the Physiological and Pathophysiological Functions of Tissues and Organs. \u003cem\u003eInt J Mol Sci\u003c/em\u003e. 2024;25(15):8167. Published 2024 Jul 26. doi:10.3390/ijms25158167\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Rettberg JR, Yao J, Brinton RD. Estrogen: a master regulator of bioenergetic systems in the brain and body. \u003cem\u003eFront Neuroendocrinol\u003c/em\u003e. 2014;35(1):8\u0026ndash;30. doi:10.1016/j.yfrne.2013.08.001IF: 6.5 Q1\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Landel V, Millet P, Baranger K, Loriod B, F\u0026eacute;ron F. Vitamin D interacts with Esr1 and Igf1 to regulate molecular pathways relevant to Alzheimer's disease. \u003cem\u003eMol Neurodegener\u003c/em\u003e. 2016;11:22. Published 2016 Mar 1. doi:10.1186/s13024-016-0087-2\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Palmisano BT, Zhu L, Stafford JM. Role of Estrogens in the Regulation of Liver Lipid Metabolism. \u003cem\u003eAdv Exp Med Biol\u003c/em\u003e. 2017;1043:227\u0026ndash;256. doi:10.1007/978-3-319-70178-3_12\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Yin G, Liang H, Sun W, et al. Shuangyu Tiaozhi decoction alleviates non-alcoholic fatty liver disease by improving lipid deposition, insulin resistance, and inflammation \u003cem\u003ein vitro\u003c/em\u003e and \u003cem\u003ein vivo\u003c/em\u003e. \u003cem\u003eFront Pharmacol\u003c/em\u003e. 2022;13:1016745. Published 2022 Nov 23. doi:10.3389/fphar.2022.1016745\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Semple RK, Chatterjee VK, O'Rahilly S. PPAR gamma and human metabolic disease. \u003cem\u003eJ Clin Invest\u003c/em\u003e. 2006;116(3):581\u0026ndash;589. doi:10.1172/JCI28003\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Rudko OI, Tretiakov AV, Naumova EA, Klimov EA. Role of PPARs in Progression of Anxiety: Literature Analysis and Signaling Pathways Reconstruction. \u003cem\u003ePPAR Res\u003c/em\u003e. 2020;2020:8859017. Published 2020 Nov 29. doi:10.1155/2020/8859017\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Gold PW. The PPARg System in Major Depression: Pathophysiologic and Therapeutic Implications. \u003cem\u003eInt J Mol Sci\u003c/em\u003e. 2021;22(17):9248. Published 2021 Aug 26. doi:10.3390/ijms22179248\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Galbraith L, Leung HY, Ahmad I. Lipid pathway deregulation in advanced prostate cancer. \u003cem\u003ePharmacol Res\u003c/em\u003e. 2018;131:177\u0026ndash;184. doi:10.1016/j.phrs.2018.02.022\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Filali-Mouncef Y, Hunter C, Roccio F, et al. The m\u0026eacute;nage \u0026agrave; trois of autophagy, lipid droplets and liver disease. \u003cem\u003eAutophagy\u003c/em\u003e. 2022;18(1):50\u0026ndash;72. doi:10.1080/15548627.2021.1895658\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Voskarides K, Giannopoulou N. The Role of \u003cem\u003eTP53\u003c/em\u003e in Adaptation and Evolution. \u003cem\u003eCells\u003c/em\u003e. 2023;12(3):512. Published 2023 Feb 3. doi:10.3390/cells12030512\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Dai CQ, Luo TT, Luo SC, et al. p53 and mitochondrial dysfunction: novel insight of neurodegenerative diseases. \u003cem\u003eJ Bioenerg Biomembr\u003c/em\u003e. 2016;48(4):337\u0026ndash;347. doi:10.1007/s10863-016-9669-5\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Staib F, Hussain SP, Hofseth LJ, Wang XW, Harris CC. TP53 and liver carcinogenesis. \u003cem\u003eHum Mutat\u003c/em\u003e. 2003;21(3):201\u0026ndash;216. doi:10.1002/humu.10176\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Zuehlke AD, Beebe K, Neckers L, Prince T. Regulation and function of the human HSP90AA1 gene. \u003cem\u003eGene\u003c/em\u003e. 2015;570(1):8\u0026ndash;16. doi:10.1016/j.gene.2015.06.018\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Smajić S, Prada-Medina CA, Landoulsi Z, et al. Single-cell sequencing of human midbrain reveals glial activation and a Parkinson-specific neuronal state. \u003cem\u003eBrain\u003c/em\u003e. 2022;145(3):964\u0026ndash;978. doi:10.1093/brain/awab446\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Wei D, Tian X, Zhu L, Wang H, Sun C. USP14 governs CYP2E1 to promote nonalcoholic fatty liver disease through deubiquitination and stabilization of HSP90AA1. \u003cem\u003eCell Death Dis\u003c/em\u003e. 2023;14(8):566. Published 2023 Aug 26. doi:10.1038/s41419-023-06091-6\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Mart\u0026iacute;n-V\u0026aacute;zquez E, Cobo-Vuilleumier N, L\u0026oacute;pez-Noriega L, Lorenzo PI, Gauthier BR. The PTGS2/COX2-PGE\u003csub\u003e2\u003c/sub\u003e signaling cascade in inflammation: Pro or anti? A case study with type 1 diabetes mellitus. \u003cem\u003eInt J Biol Sci\u003c/em\u003e. 2023;19(13):4157\u0026ndash;4165. Published 2023 Aug 6. doi:10.7150/ijbs.86492\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Li Q, Han X, Lan X, et al. Inhibition of neuronal ferroptosis protects hemorrhagic brain. \u003cem\u003eJCI Insight\u003c/em\u003e. 2017;2(7):e90777. Published 2017 Apr 6. doi:10.1172/jci.insight.90777\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Zhang CY, Tan XH, Yang HH, et al. COX-2/sEH Dual Inhibitor Alleviates Hepatocyte Senescence in NAFLD Mice by Restoring Autophagy through Sirt1/PI3K/AKT/mTOR. \u003cem\u003eInt J Mol Sci\u003c/em\u003e. 2022;23(15):8267. Published 2022 Jul 27. doi:10.3390/ijms23158267\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Mart\u0026iacute;n-Sanz P, Casado M, Bosc\u0026aacute; L. Cyclooxygenase 2 in liver dysfunction and carcinogenesis: Facts and perspectives. \u003cem\u003eWorld J Gastroenterol\u003c/em\u003e. 2017;23(20):3572\u0026ndash;3580. doi:10.3748/wjg.v23.i20.3572\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Perdichizzi S, Mascolo MG, Silingardi P, et al. Cancer-related genes transcriptionally induced by the fungicide penconazole. \u003cem\u003eToxicol In Vitro\u003c/em\u003e. 2014;28(1):125\u0026ndash;130. doi:10.1016/j.tiv.2013.06.006\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e Graziosi A, Sita G, Corrieri C, et al. Effects of Subtoxic Concentrations of Atrazine, Cypermethrin, and Vinclozolin on microRNA-Mediated PI3K/Akt/mTOR Signaling in SH-SY5Y Cells. \u003cem\u003eInt J Mol Sci\u003c/em\u003e. 2022;23(23):14538. Published 2022 Nov 22. doi:10.3390/ijms232314538\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":true,"hideJournal":true,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"Hexaconazole, Propiconazole, Prothioconazole, Network toxicology, Molecular docking, Brain injury, Liver injury","lastPublishedDoi":"10.21203/rs.3.rs-5318182/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-5318182/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003e\u003cem\u003eBackground:\u003c/em\u003e The safety of food grains is crucial for human health. Hexaconazole, Propiconazole, and Prothioconazole are widely used triazole fungicides primarily employed in agriculture for the control of plant diseases, aiming to enhance crop yield and quality.\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eObjective:\u003c/em\u003e This research aims to elucidate the potential toxic targets and underlying molecular mechanisms of brain and liver damage induced by exposure to the three fungicides through network toxicology combined with molecular docking and bioinformatics data analysis.\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eMethod and results:\u003c/em\u003e Toxicity analysis through the ADMETlab database showed that the toxicity of the three fungicides was mainly associated with neurotoxicity and liver injury. Systematically utilizing GeneCards, OMIM and other databases and through Cytoscape tools, we identified potential and core targets (including \u003cem\u003eEGFR\u003c/em\u003e, \u003cem\u003eCASP3\u003c/em\u003e, \u003cem\u003eESR1\u003c/em\u003e, \u003cem\u003ePPARG\u003c/em\u003e, \u003cem\u003eTP53\u003c/em\u003e, \u003cem\u003eHSP90AA1\u003c/em\u003e, and \u003cem\u003ePTGS2\u003c/em\u003e) associated with fungicides and brain injury or liver injury. GO and KEGG enrichment analyses indicate that these targets are associated with pathways related to cancer, the PI3K-Akt signaling pathway, organic cyclic compounds, and organic nitrogen compounds. Molecular docking simulations conducted using AutoDock confirmed the high-affinity binding interactions between the fungicides and key target proteins.\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eConclusions:\u003c/em\u003e This study provides a theoretical foundation for understanding the molecular mechanisms underlying the neurotoxicity and hepatotoxicity induced by Hexaconazole, Propiconazole, and Prothioconazole, while establishing a foundational framework for the development of prevention and treatment strategies related to triazole pesticide-associated brain and liver injuries. Our findings underscore the potential risks these three pesticides pose to brain and liver health, highlighting the need for further epidemiological and clinical research in the future.\u003c/p\u003e","manuscriptTitle":"Efficient analysis of the toxicity and mechanisms of Hexaconazole and two other triazole fungicides: insights from integrated network toxicology, molecular docking and bioinformatics data","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-10-29 07:55:57","doi":"10.21203/rs.3.rs-5318182/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"ab8b4fbf-9222-4988-94f6-f39ccb4abcf4","owner":[],"postedDate":"October 29th, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[],"tags":[],"updatedAt":"2024-10-29T07:55:57+00:00","versionOfRecord":[],"versionCreatedAt":"2024-10-29 07:55:57","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-5318182","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-5318182","identity":"rs-5318182","version":["v1"]},"buildId":"qtupq5eGEP_6zYnWcrvyt","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}