Prerpints.org logo
Preprint
Article

Network Toxicology and Molecular Docking to Investigate the non-AChE Mechanisms of Organophosphate-Induced Neurodevelopmental Toxicity

Altmetrics

Downloads

149

Views

74

Comments

0

A peer-reviewed article of this preprint also exists.

supplementary.docx (2.83MB )

Submitted:

11 July 2023

Posted:

12 July 2023

You are already at the latest version

Alerts
Abstract
Organophosphate pesticides (OPs) are toxic substances that contaminate aquatic environments, interfere with the development of the nervous system, and induce Neurodevelopmental Toxicity (NDT) in animals and humans. The canonical mechanism of OP neurotoxicity involves the inhibi-tion of acetylcholinesterase (AChE), but other mechanisms non-AChE are also involved and not fully understood. We used network toxicology and molecular docking to identify molecular targets and toxicity mechanisms common to OPs. Targets related to diazinon-oxon, chlorpyrifos oxon, and paraoxon OPs were predicted using the Swiss Target Prediction and PharmMapper databases. Targets related to NDT were compiled from GeneCards and OMIM databases. In order to construct the protein-protein interaction (PPI) network, the common targets between OPs and NDT were imported into the STRING. Networks topological analyses identified EGFR, MET, HSP90AA1, and SRC as hub nodes common to the three OPs. Using Reactome pathway and gene ontology we found that signal transduction, axon guidance, cellular responses to stress, and glutamatergic signaling activation play key roles in OP-induced NDT.
Keywords: 
Subject: Biology and Life Sciences  -   Toxicology

1. Introduction

The widespread use of organophosphate pesticides (OPs) in agriculture still raises concerns about their ecological impact and human health implication [1]. Surface waters and sediments can be contaminated by OPs from rainfall, runoff, improper disposal, or leaching from groundwater. This contamination can affect a wide range of aquatic organisms, such as fish, invertebrates, and algae, as well as non-target terrestrial organisms [2].
There are numerous toxic effects associated with OPs on fish, including hematological disturbances, gill, kidney, and liver alterations [3], oxidative stress, immune disorders, alteration of the intestinal microbiota [4], and behavioral disorders including reduced predator escape response [5]. In humans, exposure to OPs during pregnancy and postnatal periods increases the risk of autism spectrum disorder [6], impaired IQ, and verbal comprehension [7].
Exposure to OPs poses a high risk of neurotoxicity, especially for the developing nervous system. The effects of this type of exposure can include cognitive, motor, and memory impairments, as well as changes in brain morphology, referred to as Neurodevelopmental Toxicity (NDT) [8].
Acute exposure to high doses of OPs, such as self-poisoning, can significantly inhibit the activity of the acetylcholinesterase (AChE) enzyme. As a result, acetylcholine accumulates in the synaptic cleft, and cholinergic receptors in the nervous system and neuromuscular junction are overactivated, leading to a cholinergic crisis. The symptoms of human poisoning include respiratory depression, cardiovascular complications, gastrointestinal problems, and mental confusion, and the mortality rate is 15-30% [9].
When environmental contamination occurs, OPs concentrations are typically low and may not significantly inhibit AChE activity. However, these low concentrations can still cause behavioral alterations and neurological damage, particularly during developmental stages [10]. Therefore, it is important to emphasize that non-AChE mechanisms are also involved in the NDT of OPs [8]. These non-AChE mechanisms may include glutamatergic excitotoxicity, apoptosis [11], disturbances in intracellular Ca2+ homeostasis, neuronal activity [12], modulation of BDNF expression, and effects on the GABAergic and serotonergic systems [13].
Even though AChE inhibition is a common mechanism of action for OPs, there is evidence that the non-AChE effects can vary among different OP types. For instance, chlorpyrifos (CP) and malathion exhibited contrasting effects on zebrafish larvae. While malathion treatment led to morphological abnormalities in the brains of larvae, no such effect was observed with CP. Furthermore, CP-exposed larvae exhibited decreased swimming speed and increased resting time, whereas malathion-exposed larvae showed increased swimming speed and decreased resting time [14].
A comparative analysis of seven OPs (acephate, CP, dichlorvos, diazinon, malathion, parathion, and profenofos) in planarians revealed shared toxicological endpoints among the OPs (except for acephate and parathion), such as abnormal body forms, increased viscosity, writhing behavior, and alterations in swimming speed during exposure to blue light. However, other toxicological endpoints differed between the OPs [15].
This study was motivated by the following question: Do OPs share non-AChE mechanisms that contribute to NDT? This information can be valuable in enhancing toxicity testing protocols, establishing more precise safety margins, and providing targets for more effective treatments in cases of exposure to these compounds.
Computational tools can enhance toxicity studies by offering fast and comprehensive results. A valuable tool among these methods is network toxicology, a branch of network pharmacology. Utilizing techniques like protein-protein interactions, network topology analysis, and virtual prediction, it becomes possible to investigate the toxic effects of compounds based on "toxic characterization - compound - gene - protein" relationship. Integrating network toxicology with molecular docking, it is possible to assesses the interactions between the investigated compounds and their predicted targets, thereby helping to determine their toxicological properties [16, 17].
The aim of this study was to predict non-AChE targets that are shared among OP pesticides that contribute to NDT using network toxicology and molecular docking.

2. Materials and Methods

2.1. Chemical Structures of OPs

The 2D structures of the oxon metabolites of the OP pesticides diazinon (DZ), CP, and parathion were retrieved from the PubChem database (https://pubchem.ncbi.nlm.nih.gov). We chose the oxon metabolites over of the original compounds due to their direct association with NDT [18, 19].

2.2. Acquisition of OPs Targets

The 2D chemical structures of paraoxon (PO), chlorpyrifos oxon (CPO), and diazinon-oxon (DO) were submitted to the Swiss Target Prediction database (http://www.swisstargetprediction.ch) and PharmMapper (http://www.lilab-ecust.cn/pharmmapper/) to predict target genes. We then entered the predicted target genes into the UniProt database (https://www.uniprot.org/) to retrieve the standard gene names, and duplicate entries were eliminated.

2.3. Prediction of NDT Targets

The GeneCards database (https://www.genecards.org/) and Online Mendelian Inheritance in Man (OMIM, http://omim.org/) were used to predict potential NDT targets. In February 2023, the related targets were collected using the keywords "neurodevelopmental abnormalities," "neurotoxicity," and "neurodevelopmental disorders," along with the Homo sapiens species. The UniProt database was utilized to retrieve the standard gene names, and after merging the predicted targets for the three keywords from both databases, duplicate entries were removed.

2.4. Venn Analysis and Construction of Protein-Protein Interaction (PPI) Network

In order to identify the intersection genes as potential targets of each OP related to NDT, a Venn diagram was created by using Draw Venn Diagram tool (http://bioinformatics.psb.ugent.be/webtools/Venn/).
Protein-protein interaction (PPI) networks were constructed using the STRING database (https://string-db.org/). A minimum score of 0.7 was set, indicating a high level of confidence, and only experimental data were selected as the source of active interaction. The PPI network analyzes direct interactions between the inputted proteins and generates a network of interactions where proteins are represented as nodes and the interactions between them are represented as edges.

2.5. Topological Analysis of PPI Networks

The PPI networks were imported into Cytoscape 3.9.1. The cytoHubba plugin was utilized to compute the top 10 hub nodes using the Maximal Clique Centrality (MCC) method. The Analyze Network tool was used for topological analyses of the network, including degree (k), betweenness centrality (BC), clustering coefficient, closeness centrality (CC), and average shortest pathway (ASPL).
The measure k represents the number of connections a node has with other nodes in the network, i.e., it quantifies the number of edges that connect to it. Nodes with a high k value are important for information transmission in the network as they have more control over the flow of information. On the other hand, BC is a measure that indicates a node's ability to connect other nodes in the network. A high BC value indicates that the node is important for communication between different parts of the network. In turn, the CC measure indicates how many direct neighbors of a node are interconnected, and a high CC value for a node indicates its involvement in a strongly connected cluster. Lastly, the ASPL measure indicates the average shortest path length between a node and all other nodes in the network [20].

2.6. Enrichment Analysis

An enrichment analysis was performed using the Consensus PathDB database (http://cpdb.molgen.mpg.de/) to explore Gene Ontology (GO) and Reactome Pathway. The Biological Process (BP), Molecular Function (MF), and Cellular Component (CC) categories of GO were chosen. GO terms and Reactome Pathways with a p-value <0.01 were considered statistically significant. The top 15 GO terms for BP, MF, and CC, as well as the top 30 Reactome pathways, were selected for each OP. Additionally, terms common to all three OPs that were directly related to the Nervous System superpathway, as classified by Reactome.org, were examined. GO and Reactome pathway graphs were generated with the ggplot2 package in RStudio.

2.7. Molecular Docking:

A molecular docking analysis was conducted to explore the potential interactions between OPs and their hub nodes. The analysis followed the following steps:
a) Macromolecule preparation: The 3D structures of the proteins were obtained from the PDB database (https://www.rcsb.org/) in PDB format. Discovery Studio 2021 was used to remove water molecules and ligands, add polar hydrogen bonds, and obtain the x, y, z coordinates for constructing the grid box. Energy minimization was performed using Swiss-PDB Viewer software.
b) Ligand preparation: The 3D structures of the ligands were retrieved from the PubChem database (https://pubchem.ncbi.nlm.nih.gov/) in SDF format and converted to PDB format using Open Babel.
c) Docking simulation was performed using Autodock Vina [21], and the analysis of the binding and generation of images was carried out using Discovery Studio 2021.

3. Results

3.1. Candidate Targets for OP-Induced NDT

After removing the duplicates, 388 targets were found for DZO, 343 for CPO, and 387 for PO in the Swiss Target Prediction and PharmMapper databases. Additionally, 3060 NDT-related targets were obtained from the Gene Cards and OMIM databases after removing the duplicates. Using a Venn diagram to identify the intersections of targets for each OP and NDT, we found 187, 121, and 181 targets for DZO, CPO, and PO, respectively. These targets were used to construct a PPI network using the STRING database, and the data were imported into Cytoscape for visualization and analysis of the hub nodes.
The analysis of PPI networks revealed 81 nodes and 113 edges for the DZO-NDT network, 37 nodes and 50 edges for the CPO-NDT network, and 71 nodes and 103 edges for the PO-NDT network (Figure 1). The top 10 hub nodes for each network were identified using the MCC algorithm in the cytoHubba plugin of Cytoscape. Topological measurements were obtained using the Network Analyzer tool. HSP90AA1, SRC, MET, and EGFR were hub nodes that were common among the three OPs. All hub nodes are listed in Table 1, along with their topological measurements. The measurements for the entire network are available in the supplementary material (Tables S1–S3).

3.2. GO and Reactome Pathway analysis

3.2.1. GO analysis

To analyze the predicted NDT targets for each OP, we conducted a GO analysis for BP, CC and MF (p-value < 0.01). We found 169 GO terms for DZO, with 102 in BP, 20 in MF, and 47 in CC. For CPO, we identified 143 GO terms, with 89 in BP, 18 in MF, and 36 in CC. PO had 155 GO terms, with 98 in BP, 18 in MF, and 39 in CC. For each OP, we selected the Top 15 processes for each term based on their p-values (Figure 2).
We observed that the top 15 GO terms for BP, MF, and CC were similar for DZO, CPO, and PO. The results indicated that the BP terms were mainly related to cellular response to stress, which could be a result of chemical exposure, as well as cellular processes such as communication, signal transduction, and cell death, along with processes involved in the development of an anatomical structure or organism.
Regarding the MF terms, we found that most targets for the three OPs were associated with binding, which refers to the interaction between one molecule and specific sites of another molecule. The results of the CC terms suggested most targets were found within the intracellular milieu, organelles, plasma membrane, and neurons.

3.2.2. Reactome Pathway Analysis

Based on Reactome Pathway analysis, 419 pathways were identified for DZO, 254 for CPO, and 358 for PO (p-value < 0.01). Figure 3 shows the top 30 Reactome pathways for each OP based on the p-value.
Reactome Pathway is a database that encompasses 2,546 pathways organized into 28 superpathways (e.g., signal transduction, disease, immune system, and developmental biology) and their subdivisions [22]. According to the analyses, signal transduction was identified as the most important superpathway for DZO, CPO, and PO, which included pathways related to signaling by receptor tyrosine kinases, MAPK family signaling cascades, and intracellular signaling by second messengers. Additionally, all three OPs shared the axon guidance pathway, which is directly related to nervous system development, suggesting that OPs interfere with neurogenesis.
We also evaluated the Reactome pathways for DZO, CPO and PO that showed significant values (p-value < 0.01) and were classified under the Neuronal system superpathway to identify the key processes in the nervous system affected by OP exposure. Activation of NMDA receptors and postsynaptic events was the common pathway for all three OPs, indicating the activation of N-methyl-D-aspartate receptors (NMDAR), an ionotropic glutamate receptor.

3.3. Molecular Docking

We used molecular docking to predict the interactions between each OP and their respective hub nodes (presented in Table 1). A lower binding energy between the compound and the target indicates a higher affinity. The binding energies between the OPs and macromolecules ranged from -7.2 to -4.2 kJ/mol (Figure 4). Among the targets shared by the OPs, HSP90AA1 had the lowest binding energy with the OPs, followed by MET, EGFR, and SRC. Various types of bonds were formed between the OPs and the amino acid residues of the targets, including hydrogen bonding and various pi interactions (e.g., Pi-Pi, Pi-sigma, and Pi-alkyl). It is noteworthy thar EGFR and SRC formed hydrogen bonds with all three OPs, suggesting a stronger interaction. Detailed results of the molecular docking analysis can be found in the supplementary material (Table S4 and S5 and Figure S1).

4. Discussion

The classic mechanism of neurotoxicity of organophosphates is the inhibition of the AChE activity, leading to the accumulation of the neurotransmitter acetylcholine in the synaptic cleft and consequently, overstimulation of cholinergic receptors. However, numerous studies have provided evidence that multiple mechanisms can contribute to the neurotoxicity of organophosphates, particularly during neurodevelopment [12].
Reactome pathway and Gene ontology enrichment analyzes were performed to identify potential mechanisms in the NTD of OP. The results indicate that these compounds impair signal transduction and axon guidance, disrupt cellular responses to stress, and activate NMDAR. In the OP-induced NTD PPI network, HSP90AA1, EGFR, MET, and SRC were identified as common hub nodes for DZO, CPO, and PO. For better understanding, the discussion will be divided into several topics according to the mechanism of toxicity.

4.1. Signal Transduction

Signal transduction is the process by which extracellular messengers bind to transmembrane receptors and provide information that triggers cellular responses, such as biochemical and biological transformations, as well as gene expression alterations [23]. According to Reactome pathway enrichment analyses, the major signal transduction pathways affected by organophosphates include the receptor Tyrosine Kinase (RTK) pathway, the mitogen-activated protein kinase (MAPK) pathway and signaling by second messengers. It is important to note that these pathways are interconnected, since RTK activation triggers phosphorylation and recruitment of effector proteins, initiating signaling cascades involving MAPKs, PI3K/Akt, phospholipase C-PKC, and STAT [24].
In vitro studies have explored how OP exposure affects signaling pathways. For instance, exposure of cell cultures to CP and monocrotophos activated the MAPK ERK1/2, JNK, and p38 pathways, thereby inducing apoptosis. Activation of these pathways were associated with reactive oxygen species (ROS) generation and oxidative stress, which may contribute to mitochondrial damage and cell death [25,26,27]. It is plausible that these pathways may also be influenced by the interaction of OPs with RTKs.
Based on our literature review, no experimental data were found indicating that OPs cause neurodevelopmental abnormalities via EGFR, MET, or non-receptor tyrosine kinase protein SRC. However, CP has previously been shown to increase migration and invasion of breast cancer cells through c-SRC pathway activation, amplifying downstream AKT and p-38 signaling [28]. Additionally, CP has promoted the growth of human colorectal adenocarcinoma H508 cells through increased EGFR/ERK1/2 signaling [29]. Although activation of the EGFR/ERK1/2 pathway generally stimulates cell growth, it can also promote neuronal cell death [30, 31].

4.2. Axon Guidance

Axon guidance is the process by which neurons send their axons toward their correct targets for synaptic formation, and it is crucial during neural circuit development. Growing axons have a structure at their tips known as the growth cone, which contains receptors that respond to extracellular messengers that attract or repel the axon, guiding them to their destination [32]. The hub nodes EGFR, MET, and SRC play important roles in axon guidance.
The EGFR (or ErbB1) belongs to the ErbB family, which includes three other members: ErbB2, ErbB3, and ErbB4. Activation of EGFR initiates signaling pathways that are important for neuronal proliferation, differentiation, migration, as well as neural circuit development, including the Ras-Raf-MEK-ERK1/2, STAT3, and PI3K-Akt-mTOR pathways [33]. Additionally, activated EGFR interacts with SRC, forming a complex that enhances EGFR phosphorylation, consequently amplifying downstream signaling [34]. SRC also interacts with other signaling pathways involved in axon guidance, such as the Ephrin receptor (EphA) pathways [35], Sonic hedgehog (Shh) pathway [36], Netrin-1 pathway [37], and MET pathway [38]. Similar to EGFR, the activation of MET triggers downstream signaling pathways of PI3K/Akt, Ras-Raf-MEK-ERK, and STAT3, which are crucial for neurodevelopment and synaptogenesis [38]
Axonal growth is a fundamental process in the formation of neural circuits and the development of neuronal connectivity. Previous studies have reported impairment of axonal growth by OPs through noncholinergic mechanisms. For example, low concentrations of CP and CPO inhibited axonal outgrowth in cultures of rat sympathetic neurons [39]. Moreover, exposure to CPO during neurodevelopment disrupted axonal outgrowth in zebrafish [40]. These results showed that CP and CPO did not cause axon retraction or a decrease in the number of axons, but rather reduced growth rate compared to the untreated control group. Furthermore, no morphological changes in the axonal growth cone were observed. Experimental data suggest that OPs altered axonal growth signaling [39], similar to the findings in our study that OPs impair axon guidance.

4.3. Cellular Response to Stress

Exposure to environmental contaminants during prenatal development can trigger cellular response to stress that is closely related to heat shock and oxidative stress systems [41]. The HSP90AA1 gene encodes the Heat shock protein (Hsps) HSP90, which acts as a chaperone to refold denatured proteins in the presence of stressors such as heat, oxidative stress, hypoxia, and exposure to cytotoxic agents. In the absence of stress, chaperones facilitate protein refolding, stabilization, translocation, and degradation, preventing protein aggregation. In neurodevelopment, Hsps regulate pathways related to cell growth and migration, such as the PI3K/Akt signaling pathway, and they are crucial mediators for axon guidance [42].
Cellular stress increases the expression of HSP90 and diverts its functions to cope with stress responses. This functional diversion or inhibition of HSP90 can alter the activity of various proteins and result in increased mutations, leading to neurodevelopmental disorders [43]. Increased expression of HSP90 has been observed in the liver, muscle, kidneys, and spleen of fish exposed to CP, indicating that OPs exposure can induce cellular response to stress [44, 45]. Although these data were not obtained from nervous system cells, common carp exposed to CP showed increased expression of HSP90 and hypoactivity, which may indicate neurotoxicity [44].
There are several ways cells respond to stress, including activating of survival pathways or apoptosis [46]. There is still much to be clarified about the pathways by which OPs induce apoptosis, but they are known to occur during vulnerable windows of neurodevelopment and is likely not related to AChE inhibition [47]. However, it has been observed that the AKT pathway and oxidative stress are mechanisms by which OPs induce neural apoptosis, but independently of each other [48]. Therefore, the relationship between OPs and HSP90 may be a pathway to understanding how these processes occur.

4.5. Activation of NMDAR

The ionotropic glutamate receptors N-Methyl-D-Aspartate receptors (NMDAR) are calcium-permeable channels. The synaptic activity of NMDAR contributes to axonal and dendritic growth, as well as the maturation of glutamatergic synapses [49]. These processes are crucial for the development of cognitive functions, neuronal plasticity, memory, and learning [50]. Excessive activation of NMDAR can increase calcium influx, leading to cellular damage and neuronal death, likely mediated by the PKC/ERK pathway [51, 52].
During exposure to high concentrations of OPs, ionotropic AMPA and NMDA glutamate receptors are activated due to increased excitatory signaling following AChE inhibition and excessive stimulation of muscarinic receptors, leading to seizures and neuronal death [53]. However, the mechanisms related to low concentrations of these chemicals have not been fully elucidated yet.
Studies indicate that glutamate-mediated excitotoxicity seems to be more related to the toxicity effects of CP than DZ. An NMDAR antagonist, MK-801, attenuated CP toxicity, but not DZ toxicity, in cortical neuronal and glial cultures [11]. Furthermore, in neonatal rats exposed to CP and DZ, both substances upregulated the expression of NMDAR subunits, but CP induced more significant changes than DZ [54]
In our study, we did not assess the direct activation of NMDAR by OPs, nor the genes encoding the subunits of this receptor. However, the activity of NMDAR can be regulated by SRC [55], EGFR [56], and MET [57], which have been identified as hub nodes in the PPI network of OPs. Furthermore, studies with chronic exposure to low doses of OPs is a more realistic way to assess environmental exposure, and further studies are needed to evaluate the relationship between this exposure and NMDAR activity, considering other biomarkers beyond the genes encoding the subunits of this receptor.

5. Conclusions

Using network toxicology and molecular docking analyses, this study identified the hub nodes HSP90AA1, EGFR, MET, and SRC as potential targets for OPs. Additionally, signal transduction, axon guidance, cellular response to stress, and activation of NMDAR were found to be key pathways involved in OP-induced NDT. These findings are important for understanding the mechanisms of neurotoxicity of these substances at environmentally relevant concentrations that do not involve cholinergic pathways.

Supplementary Materials

The following supporting information can be downloaded at the website of this paper posted on Preprints.org, Table S1: Topological measurements of nodes in the PPI network of diazinon oxon-induced developmental neurotoxicity, Table S2: Topological measurements of nodes in the PPI network of chlorpyrifos oxon-induced developmental neurotoxicity, Table S3: Topological measurements of nodes in the PPI network of paraoxon-induced developmental neurotoxicity, Table S4: Binding energies of molecular docking between organophosphates and their respective hub nodes, Table S5 - Intermolecular interactions of complexes between HSP90AA1, EGFR, MET, and SRC and the organophosphates, Figure S1: Protein-ligand interactions between organophosphates and hub nodes generated using BIOVIA Discovery Studio Visualizer.

Author Contributions

Conceptualization, J.A.C.R.S, D.F and T.S.; methodology, J.A.C.R.S. and T.S.; data curation, J.A.C.R.S., T.S, I.L.A.C.R.Q and D.F; writing—original draft preparation, J.A.C.R.S. and I.L.A.C.R.Q; writing—review and editing, J.A.C.R.S, I.L.A.C.R.Q, T.S and D.F.; visualization, J.A.C.R.S; supervision, T.S and D.F.; project administration, D.F.; funding acquisition, D.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Public Call n. 03 Produtividade em Pesquisa PROPESQ/PRPG/UFPB, João Pessoa, Brazil [grant number PVA13245-2020]; and the Public Call Demanda Universal FAPESQ, Campina Grande, Brazil [grant number 3045/2021].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The authors confirm that all data underlying the findings are fully available without restriction. All relevant data are within the paper.

Acknowledgments

We thank to Universidade Federal da Paraíba (UFPB, Brazil) and Fundação de Apoio à Pesquisa do Estado da Paraíba (FAPESQ, Brazil) for supporting this research with grants and scholarships.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Montuori, P.; De Rosa, E.; Di Duca, F.; De Simone, B.; Scippa, S.; Russo, I.; Sorrentino, M.; Sarnacchiaro, P.; Triassi, M. Occurrence, Distribution, and Risk Assessment of Organophosphorus Pesticides in the Aquatic Environment of the Sele River Estuary, Southern Italy. Toxics, 2022, 10(7), 377. [Google Scholar] [CrossRef]
  2. Sumon, K. A.; Rashid, H.; Peeters, E. T. H. M.; Bosma, R. H.; Van den Brink, P. J. Environmental Monitoring and Risk Assessment of Organophosphate Pesticides in Aquatic Ecosystems of North-West Bangladesh. Chemosphere, 2018, 206, 92–100. [Google Scholar] [CrossRef]
  3. Akter, R.; Pervin, M. A.; Jahan, H.; Rakhi, S. F.; Reza, A. H. M. M.; Hossain, Z. Toxic Effects of an Organophosphate Pesticide, Envoy 50 SC on the Histopathological, Hematological, and Brain Acetylcholinesterase Activities in Stinging Catfish (Heteropneustes Fossilis). The Journal of Basic and Applied Zoology, 2020, 81(1), 1–14. [Google Scholar]
  4. Tang, J.; Wang, W.; Jiang, Y.; Chu, W. Diazinon Exposure Produces Histological Damage, Oxidative Stress, Immune Disorders and Gut Microbiota Dysbiosis in Crucian Carp (Carassius auratus gibelio). Environmental Pollution, 2021, 269, 116129. [Google Scholar]
  5. Sandoval-Herrera, N.; Mena, F.; Espinoza, M.; Romero, A. Neurotoxicity of Organophosphate Pesticides Could Reduce the Ability of Fish to Escape Predation under Low Doses of Exposure. Scientific Reports 2019, 9(1), 1–11. [Google Scholar]
  6. Ehrenstein, O. S. von; Ling, C.; Cui, X.; Cockburn, M.; Park, A. S.; Yu, F.; Wu, J.; Ritz, B. Prenatal and Infant Exposure to Ambient Pesticides and Autism Spectrum Disorder in Children: Population Based Case-Control Study. Br Med J, 2019, 364, 962. [Google Scholar]
  7. Gunier, R. B.; Bradman, A.; Harley, K. G.; Kogut, K.; Eskenazi, B. Prenatal Residential Proximity to Agricultural Pesticide Use and IQ in 7-Year-Old Children. Environ Health Perspect, 2017, 125(5), 057002. [Google Scholar] [CrossRef]
  8. Burke, R. D.; Todd, S. W.; Lumsden, E.; Mullins, R. J.; Mamczarz, J.; Fawcett, W. P.; Gullapalli, R. P.; Randall, W. R.; Pereira, E. F. R.; Albuquerque, E. X. Developmental Neurotoxicity of the Organophosphorus Insecticide Chlorpyrifos: From Clinical Findings to Preclinical Models and Potential Mechanisms. J Neurochem, 2017, 142 (Suppl 2), 162. [Google Scholar]
  9. Eddleston, M.; Buckley, N. A.; Eyer, P.; Dawson, A. H. Management of Acute Organophosphorus Pesticide Poisoning. Lancet, 2008, 371(9612), 597. [Google Scholar]
  10. Perez-Fernandez, C.; Morales-Navas, M.; Guardia-Escote, L.; Garrido-Cárdenas, J. A.; Colomina, M. T.; Giménez, E.; Sánchez-Santed, F. Long-Term Effects of Low Doses of Chlorpyrifos Exposure at the Preweaning Developmental Stage: A Locomotor, Pharmacological, Brain Gene Expression and Gut Microbiome Analysis. Food and Chemical Toxicology, 2020, 135, 110865. [Google Scholar]
  11. Rush, T.; Liu, X. Q.; Hjelmhaug, J.; Lobner, D. Mechanisms of Chlorpyrifos and Diazinon Induced Neurotoxicity in Cortical Culture. Neuroscience 2010, 166(3), 899–906. [Google Scholar] [CrossRef]
  12. van Melis, L. V. J.; Heusinkveld, H. J.; Langendoen, C.; Peters, A.; Westerink, R. H. S. Organophosphate Insecticides Disturb Neuronal Network Development and Function via Non-AChE Mediated Mechanisms. Neurotoxicology, 2023, 94, 35–45. [Google Scholar]
  13. Perez-Fernandez, C.; Morales-Navas, M.; Guardia-Escote, L.; Colomina, M. T.; Giménez, E.; Sánchez-Santed, F. Postnatal Exposure to Low Doses of Chlorpyrifos Induces Long-Term Effects on 5C-SRTT Learning and Performance, Cholinergic and GABAergic Systems and BDNF Expression. Exp Neurol, 2020, 330, 113356. [Google Scholar]
  14. Richendrfer, H.; Creton, R. Chlorpyrifos and Malathion Have Opposite Effects on Behaviors and Brain Size That Are Not Correlated to Changes in AChE Activity. Neurotoxicology, 2015, 49, 50–58. [Google Scholar]
  15. Ireland, D.; Zhang, S.; Bochenek, V.; Hsieh, J.-H.; Rabeler, C.; Meyer, Z.; Collins, E.-M. S. Differences in Neurotoxic Outcomes of Organophosphorus Pesticides Revealed via Multi-Dimensional Screening in Adult and Regenerating Planarians. Frontiers in Toxicology, 2022, 4, 948455. [Google Scholar]
  16. Chen, Q.; Zhang, K.; Jiao, M.; Jiao, J.; Chen, D.; Yin, Y.; Zhang, J.; Li, F. Study on the Mechanism of Mesaconitine-Induced Hepatotoxicity in Rats Based on Metabonomics and Toxicology Network. Toxins (Basel), 2022, 14(7), 486. [Google Scholar] [CrossRef]
  17. Liu, C.; Zhang, C.; Wang, W.; Yuan, F.; He, T.; Chen, Y.; Wang, Q.; Huang, J. Integrated Metabolomics and Network Toxicology to Reveal Molecular Mechanism of Celastrol Induced Cardiotoxicity. Toxicol Appl Pharmacol, 2019, 383, 114785. [Google Scholar] [CrossRef]
  18. Flaskos, J. The Developmental Neurotoxicity of Organophosphorus Insecticides: A Direct Role for the Oxon Metabolites. Toxicol Lett, 2012, 209(1), 86–93. [Google Scholar]
  19. Androutsopoulos, V. P.; Hernandez, A. F.; Liesivuori, J.; Tsatsakis, A. M. A Mechanistic Overview of Health Associated Effects of Low Levels of Organochlorine and Organophosphorous Pesticides. Toxicology, 2013, 307, 89–94. [Google Scholar]
  20. Barabási, A. L.; Oltvai, Z. N. Network Biology: Understanding the Cell’s Functional Organization. Nature Reviews Genetics, 2004, 5(2), 101–113. [Google Scholar]
  21. Trott, O.; Olson, A. J. AutoDock Vina: Improving the Speed and Accuracy of Docking with a New Scoring Function, Efficient Optimization, and Multithreading. J Comput Chem, 2010, 31(2), 455–461. [Google Scholar] [CrossRef]
  22. Gillespie, M.; Jassal, B.; Stephan, R.; Milacic, M.; Rothfels, K.; Senff-Ribeiro, A.; Griss, J.; Sevilla, C.; Matthews, L.; Gong, C.; et al. The Reactome Pathway Knowledgebase 2022. Nucleic Acids Res, 2022, 50(D1), D687. [Google Scholar] [CrossRef]
  23. van der Geer, P. Signal Transduction. In Brenner’s Encyclopedia of Genetics, 2nd ed, 2013, Stanley Maloy, A., Kelly Hughes, B., Elsevier, 436–439.
  24. Katz, M.; Amit, I.; Yarden, Y. Regulation of MAPKs by Growth Factors and Receptor Tyrosine Kinases. Biochimica et Biophysica Acta (BBA) - Molecular Cell Research 2007, 1773(8), 1161–1176. [Google Scholar]
  25. Lee, J. E.; Park, J. H.; Shin, I. C.; Koh, H. C. Reactive Oxygen Species Regulated Mitochondria-Mediated Apoptosis in PC12 Cells Exposed to Chlorpyrifos. Toxicol Appl Pharmacol, 2012, 263(2), 148–162. [Google Scholar]
  26. Ki, Y. W.; Park, J. H.; Lee, J. E.; Shin, I. C.; Koh, H. C. JNK and P38 MAPK Regulate Oxidative Stress and the Inflammatory Response in Chlorpyrifos-Induced Apoptosis. Toxicol Lett, 2013, 218(3), 235–245. [Google Scholar]
  27. Jahan, S.; Kumar, D.; Singh, S.; Kumar, V.; Srivastava, A.; Pandey, A.; Rajpurohit, C. S.; Khanna, V. K.; Pant, A. B. Resveratrol Prevents the Cellular Damages Induced by Monocrotophos via PI3K Signaling Pathway in Human Cord Blood Mesenchymal Stem Cells. Mol Neurobiol, 2018, 55(11), 8278–8292. [Google Scholar] [CrossRef]
  28. Lasagna, M.; Ventura, C.; Hielpos, M. S.; Mardirosian, M. N.; Martín, G.; Miret, N.; Randi, A.; Núñez, M.; Cocca, C. Endocrine Disruptor Chlorpyrifos Promotes Migration, Invasion, and Stemness Phenotype in 3D Cultures of Breast Cancer Cells and Induces a Wide Range of Pathways Involved in Cancer Progression. Environ Res, 2022, 204, 111989. [Google Scholar]
  29. Suriyo, T.; Tachachartvanich, P.; Visitnonthachai, D.; Watcharasit, P.; Satayavivad, J. Chlorpyrifos Promotes Colorectal Adenocarcinoma H508 Cell Growth through the Activation of EGFR/ERK1/2 Signaling Pathway but Not Cholinergic Pathway. Toxicology, 2015, 338, 117–129. [Google Scholar]
  30. Subramaniam, S.; Unsicker, K. ERK and Cell Death: ERK1/2 in Neuronal Death. FEBS J, 2010, 277(1), 22–29. [Google Scholar]
  31. Botella, J. A.; Kretzschmar, D.; Kiermayer, C.; Feldmann, P.; Hughes, D. A.; Schneuwly, S. Deregulation of the Egfr/Ras Signaling Pathway Induces Age-Related Brain Degeneration in the Drosophila Mutant Vap. Mol Biol Cell, 2003, 14(1), 241. [Google Scholar] [CrossRef]
  32. Stoeckli, E. T. Understanding Axon Guidance: Are We Nearly There Yet? Development, 2018, 145(10), dev151415. [Google Scholar]
  33. Romano, R.; Bucci, C. Role of EGFR in the Nervous System. Cells, 2020, 9(8), 1887. [Google Scholar]
  34. Chen, Z.; Oh, D.; Dubey, A. K.; Yao, M.; Yang, B.; Groves, J. T.; Sheetz, M. EGFR Family and Src Family Kinase Interactions: Mechanics Matters? Curr Opin Cell Biol, 2018, 51, 97–102. [Google Scholar] [CrossRef]
  35. Knöll, B.; Drescher, U. Src Family Kinases Are Involved in EphA Receptor-Mediated Retinal Axon Guidance. The Journal of Neuroscience, 2004, 24(28), 6248. [Google Scholar] [CrossRef]
  36. Yam, P. T.; Langlois, S. D.; Morin, S.; Charron, F. Sonic Hedgehog Guides Axons through a Noncanonical, Src-Family-Kinase-Dependent Signaling Pathway. Neuron, 2009, 62(3), 349–362. [Google Scholar] [CrossRef]
  37. Poliak, S.; Morales, D.; Croteau, L. P.; Krawchuk, D.; Palmesino, E.; Morton, S.; Cloutier, J. F.; Charron, F.; Dalva, M. B.; Ackerman, S. L.; et al. Synergistic Integration of Netrin and Ephrin Axon Guidance Signals by Spinal Motor Neurons. Elife, 2015, 4, e10841. [Google Scholar]
  38. Desole, C.; Gallo, S.; Vitacolonna, A.; Montarolo, F.; Bertolotto, A.; Vivien, D.; Comoglio, P.; Crepaldi, T. HGF and MET: From Brain Development to Neurological Disorders. Front Cell Dev Biol, 2021, 9, 683609. [Google Scholar]
  39. Howard, A. S.; Bucelli, R.; Jett, D. A.; Bruun, D.; Yang, D.; Lein, P. J. Chlorpyrifos Exerts Opposing Effects on Axonal and Dendritic Growth in Primary Neuronal Cultures. Toxicol Appl Pharmacol, 2005, 207(2), 112–124. [Google Scholar] [CrossRef]
  40. Yang, D.; Lauridsen, H.; Buels, K.; Chi, L. H.; La Du, J.; Bruun, D. A.; Olson, J. R.; Tanguay, R. L.; Lein, P. J. Chlorpyrifos-Oxon Disrupts Zebrafish Axonal Growth and Motor Behavior. Toxicological Sciences, 2011, 121(1), 146. [Google Scholar] [CrossRef]
  41. Dowell, J.; Elser, B. A.; Schroeder, R. E.; Stevens, H. E. Cellular Stress Mechanisms of Prenatal Maternal Stress: Heat Shock Factors and Oxidative Stress. Neurosci Lett, 2019, 709, 134368. [Google Scholar]
  42. Miller, D. J.; Fort, P. E. Heat Shock Proteins Regulatory Role in Neurodevelopment. Front Neurosci, 2018, 12. [Google Scholar]
  43. Zohn, I. E. Hsp90 and Complex Birth Defects: A Plausible Mechanism for the Interaction of Genes and Environment. Neurosci Lett, 2020, 716, 134680. [Google Scholar] [PubMed]
  44. Xing, H.; Liu, T.; Zhang, Z.; Wang, X.; Xu, S. Acute and Subchronic Toxic Effects of Atrazine and Chlorpyrifos on Common Carp (Cyprinus Carpio L.): Immunotoxicity Assessments. Fish Shellfish Immunol, 2015, 45(2), 327–333. [Google Scholar] [PubMed]
  45. Eder, K. J.; Leutenegger, C. M.; Köhler, H. R.; Werner, I. Effects of Neurotoxic Insecticides on Heat-Shock Proteins and Cytokine Transcription in Chinook Salmon (Oncorhynchus tshawytscha). Ecotoxicol Environ Saf, 2009, 72(1), 182–190. [Google Scholar] [CrossRef]
  46. Fulda, S.; Gorman, A. M.; Hori, O.; Samali, A. Cellular Stress Responses: Cell Survival and Cell Death. Int J Cell Biol, 2010, 214074. [Google Scholar]
  47. Slotkin, T. A.; Seidler, F. J. Developmental Neurotoxicity of Organophosphates Targets Cell Cycle and Apoptosis, Revealed by Transcriptional Profiles in Vivo and in Vitro. Neurotoxicol Teratol, 2012, 34(2), 232–241. [Google Scholar] [CrossRef]
  48. Lin, J. W.; Fu, S. C.; Liu, J. M.; Liu, S. H.; Lee, K. I.; Fang, K. M.; Hsu, R. J.; Huang, C. F.; Liu, K. M.; Chang, K. C.; et al. Chlorpyrifos Induces Neuronal Cell Death via Both Oxidative Stress and Akt Activation Downstream-Regulated CHOP-Triggered Apoptotic Pathways. Toxicology in Vitro, 2023, 86, 105483. [Google Scholar]
  49. Ewald, R. C.; Cline, H. T. NMDA Receptors and Brain Development. In Biology of the NMDA Receptor. CRC Press/Taylor & Francis: Florida, United States, 2009.
  50. Burnashev, N.; Szepetowski, P. NMDA Receptor Subunit Mutations in Neurodevelopmental Disorders. Curr Opin Pharmacol, 2015, 20, 73–82. [Google Scholar]
  51. Chen, Z.; Wang, X.; Ashraf, U.; Zheng, B.; Ye, J.; Zhou, D.; Zhang, H.; Song, Y.; Chen, H.; Zhao, S.; et al. Activation of Neuronal N-Methyl-d-Aspartate Receptor Plays a Pivotal Role in Japanese Encephalitis Virus-Induced Neuronal Cell Damage. J Neuroinflammation, 2018, 15(1), 238. [Google Scholar] [CrossRef]
  52. Jiang, S.; Li, X.; Jin, W.; Duan, X.; Bo, L.; Wu, J.; Zhang, R.; Wang, Y.; Kang, R.; Huang, L. Ketamine-Induced Neurotoxicity Blocked by N-Methyl-d-Aspartate Is Mediated through Activation of PKC/ERK Pathway in Developing Hippocampal Neurons. Neurosci Lett, 2018, 673, 122–131. [Google Scholar]
  53. Eisenkraft, A.; Falk, A.; Finkelstein, A. The Role of Glutamate and the Immune System in Organophosphate-Induced CNS Damage. Neurotoxicity Research, 2013, 24(2), 265–279. [Google Scholar] [PubMed]
  54. Slotkin, T. A.; Seidler, F. J. Comparative Developmental Neurotoxicity of Organophosphates in Vivo: Transcriptional Responses of Pathways for Brain Cell Development, Cell Signaling, Cytotoxicity and Neurotransmitter Systems. Brain Res Bull, 2007, 72(4–6), 232–274. [Google Scholar] [PubMed]
  55. Banerjee, A.; Wang, H. Y.; Borgmann-Winter, K. E.; MacDonald, M. L.; Kaprielian, H.; Stucky, A.; Kvasic, J.; Egbujo, C.; Ray, R.; Talbot, K.; et al. Src Kinase as a Mediator of Convergent Molecular Abnormalities Leading to NMDAR Hypoactivity in Schizophrenia. Mol Psychiatry, 2015, 20(9), 1091. [Google Scholar] [PubMed]
  56. Tang, Y.; Ye, M.; Du, Y.; Qiu, X.; Lv, X.; Yang, W.; Luo, J. EGFR Signaling Upregulates Surface Expression of the GluN2B-Containing NMDA Receptor and Contributes to Long-Term Potentiation in the Hippocampus. Neuroscience, 2015, 304, 109–121. [Google Scholar] [PubMed]
  57. Gallo, S.; Vitacolonna, A.; Comoglio, P.; Crepaldi, T. MET Oncogene Controls Invasive Growth by Coupling with NMDA Receptor. Cancers (Basel), 2022, 14(18), 4408. [Google Scholar] [CrossRef] [PubMed]
Preprints 79143 g001
Figure 1. PPI Networks of organophosphate-induced developmental neurotoxicity: (a) diazinon-oxon, (b) chlorpyrifos oxon, (c) paraoxon. Each node represents a gene, and edges represent the interactions between them. The size of the node is directly related to the degree, and the color intensity represents the betweenness centrality of the nodes.
Figure 1. PPI Networks of organophosphate-induced developmental neurotoxicity: (a) diazinon-oxon, (b) chlorpyrifos oxon, (c) paraoxon. Each node represents a gene, and edges represent the interactions between them. The size of the node is directly related to the degree, and the color intensity represents the betweenness centrality of the nodes.
Preprints 79143 g001
Preprints 79143 g002aPreprints 79143 g002b
Figure 2. Figure 2. Results of Gene Ontology analyses for targets of neurodevelopmental toxicity induced by organophosphates (OPs) (a) diazinon-oxon, (b) chlorpyrifos oxon, and (c) paraoxon. The top 15 processes for Biological Process, Molecular Function, and Cellular Component for each OP were selected based on the p-value. The x-axis indicates the number of genes related to each process. The colors change from red to blue, indicating a decrease in the p-value.
Figure 2. Figure 2. Results of Gene Ontology analyses for targets of neurodevelopmental toxicity induced by organophosphates (OPs) (a) diazinon-oxon, (b) chlorpyrifos oxon, and (c) paraoxon. The top 15 processes for Biological Process, Molecular Function, and Cellular Component for each OP were selected based on the p-value. The x-axis indicates the number of genes related to each process. The colors change from red to blue, indicating a decrease in the p-value.
Preprints 79143 g002aPreprints 79143 g002b
Preprints 79143 g003
Figure 3. Results of Reactome pathway analyses for targets of organophosphate-induced developmental neurotoxicity (a) diazinon-oxon, (b) chlorpyrifos oxon, and (c) paraoxon. The top 30 pathways based on p-value are shown. The x-axis represents the number of genes associated with each pathway. The colors transition from red to blue, indicating a decrease in p-value.
Figure 3. Results of Reactome pathway analyses for targets of organophosphate-induced developmental neurotoxicity (a) diazinon-oxon, (b) chlorpyrifos oxon, and (c) paraoxon. The top 30 pathways based on p-value are shown. The x-axis represents the number of genes associated with each pathway. The colors transition from red to blue, indicating a decrease in p-value.
Preprints 79143 g003
Preprints 79143 g004
Figure 4. Heatmap of molecular docking of organophosphates (CPO, DZO, and PO) with their respective hub nodes, obtained from the PPI network analysis. The y-axis lists the hub nodes, and the x-axis lists the organophosphates. The color scale of the heatmap cells represents the predicted binding energy for each interaction. CPO: chlorpyrifos oxon, DZO: diazinon-oxon, PO: paraoxon.
Figure 4. Heatmap of molecular docking of organophosphates (CPO, DZO, and PO) with their respective hub nodes, obtained from the PPI network analysis. The y-axis lists the hub nodes, and the x-axis lists the organophosphates. The color scale of the heatmap cells represents the predicted binding energy for each interaction. CPO: chlorpyrifos oxon, DZO: diazinon-oxon, PO: paraoxon.
Preprints 79143 g004
Preprints 79143 g005aPreprints 79143 g005b
Figure 5. Interactions among the organophosphates diazinon-oxon (DZO), chlorpyrifos oxon (CPO), and paraoxon (PO) with the common hub nodes between them. Letters A, B, and C denote the interactions between EGFR and DZO, CPO, and PO, respectively. Letters D, E, and F represent the interactions between HSP90AA1 and DZO, CPO, and PO, respectively. Letters G, H, and I indicate the interactions between MET and DZO, CPO, and PO, respectively. Letters J, K, and L indicate the interactions between SRC and DZO, CPO, and PO, respectively.
Figure 5. Interactions among the organophosphates diazinon-oxon (DZO), chlorpyrifos oxon (CPO), and paraoxon (PO) with the common hub nodes between them. Letters A, B, and C denote the interactions between EGFR and DZO, CPO, and PO, respectively. Letters D, E, and F represent the interactions between HSP90AA1 and DZO, CPO, and PO, respectively. Letters G, H, and I indicate the interactions between MET and DZO, CPO, and PO, respectively. Letters J, K, and L indicate the interactions between SRC and DZO, CPO, and PO, respectively.
Preprints 79143 g005aPreprints 79143 g005b
Table 1. Topological measurements of hub nodes in the PPI networks of organophosphate-induced developmental neurotoxicity, identified using the MCC algorithm in the cytoHubba plugin of Cytoscape. k: degree; BC: Betweenness centrality; ASPL: Average shortest path length; CC: Closeness centrality.
Table 1. Topological measurements of hub nodes in the PPI networks of organophosphate-induced developmental neurotoxicity, identified using the MCC algorithm in the cytoHubba plugin of Cytoscape. k: degree; BC: Betweenness centrality; ASPL: Average shortest path length; CC: Closeness centrality.
Protein Node k Clustering Coefficient BC ASPL CC
Diazinon-oxon
Proto-oncogene tyrosine kinase SRC SRC 13 0.218 0.154 2.286 0.438
Epidermal growth factor receptor EGFR 14 0.165 0.275 2.143 0.467
Phosphatidylinositol 3-kinase regulatory subunit alpha PIK3R1 9 0.250 0.238 2.408 0.415
Tyrosine-protein phosphatase non-receptor type 11 PTPN11 7 0.4286 0.0164 2.7143 0.368
Heat shock protein HSP 90-alpha HSP90AA1 11 0.109 0.154 2.429 0.412
Gamma-secretase subunit APH-1A APH1A 4 0.8333 0.0667 1.5000 0.667
Nicastrin NCSTN 4 0.833 0.067 1.500 0.667
Gamma-secretase subunit PEN-2 PSENEN 4 0.8333 0.0667 1.5000 0.667
Signal transducer and activator of transcription 3 STAT3 8 0.179 0.142 2.327 0.430
Hepatocyte growth factor receptor MET 4 0.8333 0.0002 2.9184 0.343
Chlorpyrifos oxon
Heat shock protein HSP 90-alpha HSP90AA1 11 0.109 0.382 2.034 0.492
Epidermal growth factor receptor EGFR 8 0.179 0.276 2.172 0.460
Proto-oncogene tyrosine kinase SRC SRC
Estrogen receptor ESR1 7 0.143 0.269 2.103 0.475
Tyrosine-protein phosphatase non-receptor type 1 PTPN1 4 0.500 0.017 2.759 0.363
Hepatocyte growth factor receptor MET 3 1.000 0.000 2.793 0.358
Mitogen-activated protein kinase1 MAPK1 5 0.100 0.095 3.172 0.315
Calmodulin 1 CALM1 5 0.000 0.158 2.414 0.414
Histone acetyltransferase p300 EP300 5 0.000 0.308 2.552 0.392
Heat shock cognate 71 kDa protein HSPA8 3 0.667 0.008 2.379 0.420
Paraoxon
Proto-oncogene tyrosine kinase SRC SRC 14 0.220 0.264 1.875 0.533
Epidermal growth factor receptor EGFR 14 0.198 0.312 1.800 0.556
Tyrosine-protein phosphatase non-receptor type 11 PTPN11 9 0.361 0.059 2.225 0.449
Phosphatidylinositol 3-kinase regulatory subunit alpha PIK3R1 10 0.289 0.180 2.150 0.465
Steroid hormone receptor ERR2 ERRB2 6 0.667 0.026 2.100 0.476
Heat shock protein HSP 90-alpha HSP90AA1 15 0.095 0.378 1.925 0.519
Gamma-secretase subunit APH-1A APH1A 4 0.833 0.067 1.200 0.833
Nicastrin NCSTN 4 0.833 0.067 1.200 0.833
Hepatocyte growth factor receptor MET 4 0.833 0.001 2.450 0.408
Gamma-secretase subunit PEN-2 PSENEN 4 0.833 0.067 1.200 0.833
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated