View Article

  • In Silico Evaluation of Competitive Binding of Food Preservative Combinations at AChE, MAO-B and NMDA Receptors and Their Network Overlap with Neurodegeneration Pathways

  • * Department of Pharmaceutical Chemistry, Jagannath University, Chaksu, Jaipur.
    2 Institute of Pharmaceutical Sciences, University of Lucknow.
    2 Department of Zoology, Seth G.B. Podar College, Nawalgarh, Rajasthan.

Abstract

Background: Food preservatives such as sodium benzoate, potassium sorbate, sodium propionate, sodium nitrite, butylated hydroxyanisole (BHA), and butylated hydroxytoluene (BHT) are consumed together in packaged foods, yet regulatory limits and most neuropharmacological evaluations assess them as isolated chemicals.Objective: To evaluate, in silico, the individual and competitive pairwise binding affinities of these food preservatives at human acetylcholinesterase (AChE), monoamine oxidase B (MAO-B), and NMDA receptor (GluN2B), and identify their systems-level mechanistic overlap with Alzheimer's and Parkinson's disease pathways.Methods: Molecular docking (AutoDock Vina) was validated through native ligand redocking (RMSD < 2.0 Å) and compared against standard clinical references (donepezil, selegiline, memantine). Pairwise competitive binding was evaluated via sequential co-docking simulations. Biological target profiles (SwissTargetPrediction) were intersected with Alzheimer's and Parkinson's datasets (GeneCards, DisGeNET), mapped onto protein-protein interaction (STRING, Cytoscape cytoHubba), and characterized via KEGG and GO functional enrichment.Results: BHT and BHA displayed the strongest predicted solitary binding across all neuro-targets (AChE: -6.4 and -5.9 kcal/mol; MAO-B: -6.8 and -6.2 kcal/mol; NMDA: -6.3 and -5.8 kcal/mol), closely approaching standard reference drugs. Polar organic preservatives exhibited lower solitary binding (-2.6 to -5.5 kcal/mol). Sequential docking demonstrated prominent steric hindrance for bulky phenolic pairs, whereas small planar preservatives (benzoate) co-docked cooperatively with BHA (??G up to -0.8 kcal/mol). Network pharmacology highlighted 10 critical hub genes (AKT1, TP53, TNF, CASP3, EGFR, SRC, ESR1, PTGS2, MAPK3, MTOR) significantly converging on neuroactive ligand-receptor interactions, apoptosis, and neuroinflammation pathways (p < 10^-5).Conclusion: Dietary preservatives, especially synthetic phenolic antioxidants and multi-preservative combinations, demonstrate direct computational affinity for critical central neuroreceptors and converge on apoptotic and neuroinflammatory cascades, underscoring the urgent need for empirical multi-preservative neurotoxicological risk assessments.

Keywords

food preservatives, sodium benzoate, molecular docking, network pharmacology, acetylcholinesterase, MAO-B, NMDA receptor, neurodegeneration, synergistic binding

Introduction

× Popup Image

Synthetic food preservatives are widely utilized across modern food manufacturing to prevent microbial spoilage, delay enzymatic decay, and impede oxidative degradation. Among the most prevalent chemical additives are organic acid salts (sodium benzoate, potassium sorbate, sodium propionate), reactive nitrogen precursors (sodium nitrite), and synthetic sterically hindered phenolic antioxidants (butylated hydroxyanisole [BHA] and butylated hydroxytoluene [BHT]). Global regulatory bodies—including the Joint FAO/WHO Expert Committee on Food Additives (JECFA), the European Food Safety Authority (EFSA), and the US FDA—determine acceptable daily intakes (ADIs) based largely on single-substance toxicological evaluations. However, real-world human dietary habits entail concurrent exposure to multi-additive cocktails present in carbonated beverages, sauces, processed meats, and confectionery products.

Existing experimental literature reveals a complex and often contradictory picture regarding the neuropharmacological profiles of food additives. At supra-physiological dosages, rodent models demonstrate that sodium benzoate and sodium nitrite trigger systemic oxidative stress, hippocampal neuronal loss, and cognitive deficits. Conversely, clinical and pre-clinical trials have demonstrated that sodium benzoate acts as a D-amino acid oxidase (DAAO) inhibitor, enhancing synaptic D-serine availability and ameliorating cognitive and psychotic symptoms in clinical schizophrenia and Alzheimer's models. This striking duality highlights an urgent necessity to delineate how preservatives engage key central nervous system (CNS) enzymes and receptors at the atomic level.

Three macromolecular targets play pivotal, well-established roles in neurodegenerative cascades:

 (1) Acetylcholinesterase (AChE), which regulates cholinergic neurotransmission and whose peripheral anionic site accelerates amyloid-beta aggregation; (2) Monoamine Oxidase B (MAO-B), an outer mitochondrial membrane enzyme that governs dopamine metabolism and whose hyperactivation produces cytotoxic reactive oxygen species (ROS) in Parkinson's disease; and

(3) N-Methyl-D-Aspartate (NMDA) Receptors (GluN1/GluN2B), where channel disruption or unmodulated allosteric stimulation precipitates excitotoxic intracellular calcium overload. Most computational and toxicological evaluations assess single additives against single targets, leaving the questions of competitive active-site occupancy and synergistic pathway engagement unanswered.

Therefore, this investigation establishes a comprehensive in silico paradigm combining crystallographic molecular docking, sequential competitive co-docking simulations, pharmacokinetic ADMET profiling, and systems-level network pharmacology to evaluate individual and combined preservative interactions across AChE, MAO-B, and NMDA targets.

MATERIALS AND METHODS

2.1 Ligand Preparation

Three-dimensional chemical structures of six food preservatives—Sodium Benzoate (benzoic acid form, PubChem CID: 243), Potassium Sorbate (sorbic acid form, CID: 5281515), Sodium Propionate (propionic acid form, CID: 1032), Sodium Nitrite (CID: 24526), BHA (CID: 24667), and BHT (CID: 31404)—along with standard clinical neuro-therapeutics: Donepezil (CID: 3152, AChE inhibitor), Selegiline (CID: 26757, MAO-B inhibitor), and Memantine (CID: 4054, NMDA channel blocker) were obtained from NCBI PubChem in SDF format. Structures were imported into Open Babel within PyRx 0.8, assigned physiological ionization states at pH 7.4, energy-minimized using the MMFF94 force field (200 conjugate gradient steps), and converted into PDBQT format with atomic Gasteiger partial charges and defined rotatable bonds.

2.2 Protein Preparation

High-resolution X-ray crystal structures of human target proteins were retrieved from the RCSB Protein Data Bank (PDB): Human Recombinant AChE in complex with Donepezil (PDB ID: 4EY7, resolution 2.35 Å); Human MAO-B in complex with Selegiline/deprenyl analog (PDB ID: 2Z5X, resolution 1.65 Å); and Human NMDA Receptor GluN1/GluN2B ligand-binding/channel domain (PDB ID: 4PE5, resolution 2.50 Å). Protein cleanup was carried out using BIOVIA Discovery Studio Visualizer 2021 and PyMOL 2.5: water molecules, co-factors, and non-receptor crystallographic ligands were removed; polar hydrogens and Kollman charges were added; and structures were saved as rigid receptor PDBQT files.

2.3 Docking and Protocol Validation

Docking calculations were executed using AutoDock Vina v.1.2.0 via PyRx. Grid boxes were centered on the co-crystallized native inhibitor coordinates: AChE (4EY7: x = -13.98, y = -43.88, z = 27.91; size 24x24x24 Å); MAO-B (2Z5X: x = 52.82, y = 157.34, z = 25.11; size 22x22x22 Å); and NMDA (4PE5: x = -13.25, y = -18.72, z = -10.45; size 22x22x22 Å). An exhaustiveness level of 32 was maintained. Protocol validation was verified by extracting and re-docking native co-crystallized ligands, ensuring an RMSD < 2.0 Å relative to experimental crystallographic conformations.

2.4 Combination (Competitive) Docking

To computationally model simultaneous binary exposure within the catalytic binding pocket, a sequential co-docking workflow was implemented. For each preservative pair (A + B), Ligand A was first docked into the apoprotein. The lowest energy conformation (Pose 1) was merged into the receptor coordinate file, generating a modified 'Target+Ligand A' complex. Ligand B was subsequently docked into this modified pocket using identical grid boundaries. Binding energy shifts were computed as ΔΔG = ΔG (Ligand B in [Receptor+A]) - ΔG (Ligand B alone), where negative values denote cooperative pocket stabilization and positive values reflect competitive steric occlusion.

2.5 Pharmacokinetic and ADMET Prediction

Physicochemical descriptors, Lipinski's Rule of Five compliance, gastrointestinal absorption, and blood-brain barrier (BBB) penetration were determined using SwissADME. Quantitative organ toxicities, hepatotoxicity, and acute oral rat lethal dose 50 (LD50 in mg/kg) classifications were modeled via ProTox-II and pkCSM toxicology platforms.

2.6 Network Pharmacology, PPI, and Enrichment

Putative human macromolecular targets of the six preservatives were predicted using SwissTargetPrediction (Homo sapiens, probability score >= 0.1). Genes associated with 'Alzheimer's Disease' and 'Parkinson's Disease' were harvested from GeneCards (score >= 15.0) and DisGeNET (score >= 0.2). Target-disease intersections were identified using InteractiVenn and mapped into STRING 11.5 (high confidence score >= 0.70). Protein-protein interaction (PPI) networks were visualized in Cytoscape v3.9.1, and top hub bottleneck genes were extracted using cytoHubba. Gene Ontology (GO) and KEGG enrichment analyses were conducted via DAVID Bioinformatics Resources with Benjamini-Hochberg false discovery rate corrections (FDR p < 0.05).

RESULTS

3.1 Docking Validation

Redocking native co-crystallized inhibitors reproduced crystallographic binding poses with high fidelity across all three target receptors (Table 1), demonstrating RMSD values significantly below the 2.0 Å threshold and confirming the validity of the docking parameters.

Target

PDB ID

Native Ligand

Resolution (Å)

RMSD (Å)

Validated (< 2.0 Å)

AChE

4EY7

Donepezil (E20)

2.35

0.84

Yes

MAO-B

2Z5X

Safinamide analog

1.65

0.96

Yes

NMDA

4PE5

Ifenprodil analog

2.50

1.12

Yes

3.2 Binding of Individual Preservatives

Predicted binding free energies (Table 2) showed that sterically bulky, lipophilic preservatives (BHT and BHA) exhibited the strongest affinity across all three central targets, while low-molecular-weight polar organic acids displayed modest binding scores.

Compound

AChE (4EY7)

MAO-B (2Z5X)

NMDA (4PE5)

Strongest Target

Sodium Benzoate

-5.5

-5.2

-4.9

AChE

Potassium Sorbate

-4.8

-4.6

-4.2

AChE

Sodium Propionate

-3.8

-3.7

-3.4

AChE

Sodium Nitrite

-2.8

-2.6

-2.5

AChE

BHA

-5.9

-6.2

-5.8

MAO-B

BHT

-6.4

-6.8

-6.3

MAO-B

Donepezil (Ref)

-11.5

-8.1

-7.4

AChE

Selegiline (Ref)

-6.8

-7.8

-5.9

MAO-B

Memantine (Ref)

-5.7

-5.4

-7.2

NMDA

3.3 Key Molecular Interactions

Detailed active-site binding analysis (Table 3) showed that BHT and BHA engaged in extensive hydrophobic and pi-pi stacking contacts with aromatic residues within AChE and MAO-B, while sodium benzoate predominantly formed hydrogen bonds and electrostatic interactions.

Target

Ligand

dG (kcal/mol)

H-Bond Residues

Hydrophobic / Pi Residues

AChE

BHT

-6.4

Tyr124 (2.68 Å)

Trp86 (pi-alkyl), Tyr337 (pi-pi), Phe338

AChE

Benzoate

-5.5

Ser203 (2.12 Å), Gly121

Trp86 (pi-pi T-shaped), Tyr341, Phe295

MAO-B

BHT

-6.8

Cys172 (2.45 Å)

Phe343, Tyr398 (pi-pi), Leu171, Ile199

MAO-B

BHA

-6.2

Tyr435 (2.20 Å)

Phe168 (pi-alkyl), Trp119, Ile316, Tyr326

NMDA

BHT

-6.3

Glu236 (2.54 Å)

Phe176 (pi-alkyl), Trp223, Ile133, Val181

NMDA

Benzoate

-4.9

Arg115 (1.98 Å), Thr116

Phe114 (pi-pi), His113

3.4 Combination Docking

Sequential competitive docking of 15 preservative pairs (Table 4) revealed that combining two bulky phenolic preservatives (BHA + BHT) caused severe steric clash (positive delta-delta G), whereas pairing small planar molecules (benzoate) with phenolic antioxidants (BHA) resulted in cooperative affinity enhancement.

Pair (A + B)

AChE B alone

AChE B + A

MAO-B B alone

MAO-B B + A

NMDA B alone

NMDA B + A

Benzoate + Sorbate

-4.8

-4.9

-4.6

-4.5

-4.2

-4.1

Benzoate + Propionate

-3.8

-3.7

-3.7

-3.6

-3.4

-3.5

Benzoate + Nitrite

-2.8

-2.9

-2.6

-2.6

-2.5

-2.5

Benzoate + BHA

-5.9

-6.5

-6.2

-6.7

-5.8

-6.4

Benzoate + BHT

-6.4

-5.6

-6.8

-5.9

-6.3

-5.4

Sorbate + Propionate

-3.8

-3.8

-3.7

-3.7

-3.4

-3.3

Sorbate + Nitrite

-2.8

-2.7

-2.6

-2.6

-2.5

-2.4

Sorbate + BHA

-5.9

-6.2

-6.2

-6.5

-5.8

-6.1

Sorbate + BHT

-6.4

-5.8

-6.8

-6.1

-6.3

-5.7

Propionate + Nitrite

-2.8

-2.8

-2.6

-2.5

-2.5

-2.5

Propionate + BHA

-5.9

-6.1

-6.2

-6.3

-5.8

-6.0

Propionate + BHT

-6.4

-6.0

-6.8

-6.3

-6.3

-5.9

Nitrite + BHA

-5.9

-6.0

-6.2

-6.3

-5.8

-5.9

Nitrite + BHT

-6.4

-6.2

-6.8

-6.5

-6.3

-6.1

BHA + BHT

-6.4

-4.9

-6.8

-5.1

-6.3

-4.8

3.5 ADMET Prediction

Predicted pharmacokinetic profiles (Table 5) confirmed high gastrointestinal absorption for all preservatives. Sodium benzoate, BHA, and BHT were computationally predicted to penetrate the blood-brain barrier (BBB). Toxicological profiling identified sodium nitrite as Toxicity Class III (toxic) and BHT as Class IV (harmful).

Compound

MW (g/mol)

LogP

TPSA (Ų)

GI Abs.

BBB Perm.

LD50 (mg/kg) / Class

Sodium Benzoate

122.12

1.52

37.30

High

Yes

1700 (Class IV)

Potassium Sorbate

112.13

1.41

37.30

High

Yes

3800 (Class V)

Sodium Propionate

74.08

0.48

37.30

High

No

5100 (Class VI)

Sodium Nitrite

69.00

-0.65

50.14

High

No

180 (Class III)

BHA

180.24

2.81

29.46

High

Yes

2000 (Class IV)

BHT

220.35

4.17

20.23

High

Yes

890 (Class IV)

3.6 Network Analysis

Intersecting the predicted targets of the six preservatives with Alzheimer's and Parkinson's disease targetomes yielded 84 shared neuro-associated genes (Table 6). Topological analysis identified 10 key hub genes (Table 7), and pathway enrichment demonstrated significant convergence on neuroinflammatory and apoptotic signaling cascades (Table 8).

Compound

Predicted Targets

Overlap AD Genes

Overlap PD Genes

Total Unique Overlap

Sodium Benzoate

42

19

14

24

Potassium Sorbate

28

11

9

15

Sodium Propionate

35

16

13

20

Sodium Nitrite

19

9

8

12

BHA

67

34

26

41

BHT

58

29

22

36

Combined Unique

156

68

52

84

 

Rank

Hub Gene

Full Protein Name

Degree

Linked Preservatives

1

AKT1

RAC-alpha serine/threonine kinase

48

Benzoate, BHA, BHT

2

TP53

Cellular tumor antigen p53

44

Benzoate, Nitrite, BHA

3

TNF

Tumor necrosis factor alpha

41

Benzoate, Propionate, BHT

4

CASP3

Caspase-3 (Apoptosis executioner)

39

Benzoate, Nitrite, BHA, BHT

5

EGFR

Epidermal growth factor receptor

36

BHA, BHT

6

SRC

Proto-oncogene tyrosine kinase Src

32

Sorbate, BHA

7

ESR1

Estrogen receptor 1

30

BHA, BHT

8

PTGS2

Cyclooxygenase-2 (COX-2)

29

Benzoate, BHA, BHT

9

MAPK3

Mitogen-activated protein kinase 3

28

Propionate, BHA, BHT

10

MTOR

Mechanistic target of rapamycin

26

Benzoate, BHT

 

Pathway / Term ID

Description

Gene Count

FDR (p-value)

hsa04080

Neuroactive ligand-receptor interaction

22

1.42 x 10^-8

hsa05010

Alzheimer's disease pathway

19

3.84 x 10^-7

hsa04210

Apoptosis signaling pathway

16

8.21 x 10^-6

hsa05012

Parkinson's disease pathway

14

2.15 x 10^-5

hsa04151

PI3K-Akt signaling pathway

21

4.67 x 10^-5

GO:0006954

Inflammatory response (BP)

24

1.12 x 10^-7

GO:0008630

Intrinsic apoptotic signaling (BP)

13

6.54 x 10^-6

GO:0007268

Chemical synaptic transmission (BP)

17

1.88 x 10^-5

DISCUSSION

This investigation evaluated the direct molecular binding kinetics and systems-level network pharmacology of six ubiquitous dietary preservatives at human AChE, MAO-B, and NMDA receptor complexes. Our docking simulations demonstrated that synthetic phenolic antioxidants (BHT and BHA) exhibited moderate-to-high affinities across all three neuro-targets, approaching the binding energies of clinical reference standards such as Selegiline and Memantine. Structural interaction mapping confirmed that BHT and BHA insert stably into hydrophobic aromatic gorges (e.g., Trp86, Tyr337 in AChE; Phe343, Tyr398 in MAO-B), indicating the potential to interfere with physiological substrate binding.

Sequential co-docking simulations provided valuable insights into cocktail exposure dynamics. When bulky lipophilic compounds (BHA + BHT) competed for identical catalytic pockets, pronounced steric clash emerged, significantly destabilizing the binding energy of the secondary ligand. In contrast, small planar molecules (benzoate) co-docked alongside BHA exhibited cooperative stabilization (negative delta-delta G), suggesting that smaller additives may anchor to peripheral sites while lipophilic additives occupy catalytic cores.

Network analysis revealed that the targetome of these preservatives converges on essential regulatory hubs—notably AKT1, TP53, TNF, CASP3, and PTGS2 (COX-2)—which govern apoptosis, neuroinflammation, and synaptic plasticity. While low-dose sodium benzoate is recognized for its beneficial D-amino acid oxidase (DAAO) inhibitory effects, continuous multi-preservative co-exposure—particularly combinations involving nitrites and phenolic antioxidants—may trigger chronic neuroinflammatory signaling and cell death pathways.

4.1 Limitations

Molecular docking simulations utilize rigid receptor coordinates and cannot capture dynamic physiological conformational changes, metabolic transformations, or absolute in vivo concentrations. Sequential co-docking models steric feasibility rather than empirical non-equilibrium binding kinetics. Network pharmacology predictions are subject to database curation coverage and lack direct in vitro kinetic confirmation.

4.2 Future Directions

Subsequent studies should prioritize 100-200 ns all-atom Molecular Dynamics (MD) simulations and MM-PBSA free energy calculations. Experimental validation via in vitro enzymatic assays (Ellman's AChE assay, fluorometric MAO-B inhibition assays) and in vivo multi-preservative dietary exposure models in rodents will be crucial to establish regulatory safety margins for additive mixtures.

CONCLUSION

This study demonstrates that common food preservatives, particularly synthetic phenolic antioxidants BHT and BHA, bind directly to human AChE, MAO-B, and NMDA receptor active sites. Sequential co-docking reveals distinct steric competition and cooperative binding modes depending on molecular structure. Furthermore, network pharmacology establishes significant overlap with apoptotic and neuroinflammatory cascades in Alzheimer's and Parkinson's pathologies. These computational insights highlight the necessity of updating regulatory toxicology paradigms to evaluate chronic dietary co-exposure to multi-preservative mixtures.

REFERENCES

  1. EFSA Panel on Food Additives and Nutrient Sources added to Food (ANS). Scientific Opinion on the re-evaluation of benzoic acid and its salts. EFSA J. 2016;14 (3):4433.
  2. World Health Organization, Joint FAO/WHO Expert Committee on Food Additives. Evaluation of certain food additives: eighty-fourth report of JECFA. WHO Tech Rep Ser. 2017;1007:1-84.
  3. Piper PW. Yeast safety and human health: the potential toxicity of the common food preservative sodium benzoate. Mol Biosyst. 2018;14 (2):292-300.
  4. Noorafshan A, Erfanizadeh M, Karbalay-Doust S. High dose of sodium benzoate induces neurotoxicity in the hippocampus of mice. Neurosci Lett. 2014;563:63-68.
  5. Khoshnood MJ, Hosseini A, Seydi E, Pourahmad J. Uncoupling of oxidative phosphorylation as a mechanism of sodium benzoate cytotoxicity. Toxicology. 2015;338:11-17.
  6. Lane HY, Lin CH, Green MF, et al. Add-on treatment of benzoate for schizophrenia: a randomized, double-blind, placebo-controlled trial. JAMA Psychiatry. 2013;70 (12):1267-1275.
  7. Lin CH, Chen PK, Wang SH, Lane HY. Sodium benzoate, a D-amino acid oxidase inhibitor, added to clozapine for the treatment of schizophrenia: a randomized, double-blind, placebo-controlled trial. Lancet Psychiatry. 2018;5 (5):406-415.
  8. Lane HY, Tsai GE. D-amino acid oxidase inhibitors in cognitive neuroscience and clinical psychiatry. Clin Psychopharmacol Neurosci. 2019;17 (3):341-352.
  9. Silman I, Sussman JL. Acetylcholinesterase: 'classical' and 'non-classical' functions and roles in neurodegenerative diseases. Physiol Rev. 2005;85 (4):1199-1222.
  10. Cheung J, Rudolph MJ, Burshteyn F, et al. Structures of human acetylcholinesterase in complex with pharmacologically important inhibitors. J Med Chem. 2012;55 (22):10282-10286.
  11. Binda C, Wang J, Pisani L, et al. Structures of human monoamine oxidase B complexes with selective noncovalent inhibitors. J Med Chem. 2007;50 (12):2775-2782.
  12. Edmondson DE, Binda C, Wang J, Upadhyay AK, Mattevi A. Molecular and structural basis of substrate and inhibitor selectivity of monoamine oxidases A and B. Curr Med Chem. 2009;16 (4):441-452.
  13. Karakas E, Furukawa H. Crystal structure of a heterotetrameric NMDA receptor ion channel. Science. 2014;344 (6187):992-997.
  14. Zhu S, Stein RA, Yoshioka C, et al. Mechanism of NMDA receptor inhibition and allosteric modulation. Nature. 2016;536 (7616):319-325.
  15. Trott O, Olson AJ. 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.
  16. Eberhardt J, Santos-Martins D, Tillack AF, Forli S. AutoDock Vina 1.2.0: New docking methods, expanded force field, and python bindings. J Chem Inf Model. 2021;61 (8):3891-3898.
  17. Dallakyan S, Olson AJ. Small-molecule library screening by docking with PyRx. Methods Mol Biol. 2015;1263:243-250.
  18. Daina A, Michielin O, Zoete V. SwissADME: a free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci Rep. 2017;7:42717.
  19. Pires DEV, Blundell TL, Ascher DB. pkCSM: predicting small-molecule pharmacokinetic and toxic properties using graph-based signatures. J Med Chem. 2015;58 (9):4066-4072.
  20. Banerjee P, Eckert AO, Schrey AK, Preissner R. ProTox-II: a webserver for the prediction of toxicity of chemicals. Nucleic Acids Res. 2018;46 (W1):W257-W263.
  21. Daina A, Michielin O, Zoete V. SwissTargetPrediction: updated data and new features for efficient prediction of protein targets of small molecules. Nucleic Acids Res. 2019;47 (W1):W357-W364.
  22. Stelzer G, Rosen N, Plaschkes I, et al. The GeneCards suite: from genomic information to robust knowledge of human disease. Curr Protoc Bioinformatics. 2016;54 (1):1.30.1-1.30.33.
  23. Piñero J, Ramírez-Anguita JM, Saüch-Pitarch J, et al. The DisGeNET knowledge platform for disease genomics: 2019 update. Nucleic Acids Res. 2020;48 (D1):D845-D855.
  24. Szklarczyk D, Gable AL, Nastou KC, et al. The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded datasets. Nucleic Acids Res. 2021;49 (D1):D605-D612.
  25. Shannon P, Markiel A, Ozier O, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13 (11):2498-2504.
  26. Chin CH, Chen SH, Wu HH, et al. cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst Biol. 2014;8 (Suppl 4):S11.
  27. Sherman BT, Hao M, Qiu J, et al. DAVID: a web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic Acids Res. 2022;50 (W1):W216-W221.
  28. Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28 (1):27-30.
  29. EFSA Panel on Food Additives and Nutrient Sources added to Food (ANS). Scientific Opinion on the re-evaluation of potassium sorbate (E 202). EFSA J. 2015;13 (6):4144.
  30. EFSA Panel on Food Additives and Nutrient Sources added to Food (ANS). Scientific Opinion on the re-evaluation of propionic acid and sodium propionate. EFSA J. 2014;12 (7):3779.
  31. EFSA Panel on Food Additives and Nutrient Sources added to Food (ANS). Re-evaluation of sodium nitrite (E 250) and potassium nitrite (E 249) as food additives. EFSA J. 2017;15 (6):4786.
  32. EFSA Panel on Food Additives and Nutrient Sources added to Food (ANS). Scientific Opinion on the re-evaluation of butylated hydroxyanisole - BHA (E 320) as a food additive. EFSA J. 2011;9 (10):2418.
  33. EFSA Panel on Food Additives and Nutrient Sources added to Food (ANS). Scientific Opinion on the re-evaluation of butylated hydroxytoluene - BHT (E 321) as a food additive. EFSA J. 2012;10 (3):2588.
  34. McCann D, Barrett A, Cooper A, et al. Food additives and hyperactive behaviour in 3-year-old and 8/9-year-old children in the community: a randomised, double-blinded, placebo-controlled trial. Lancet. 2007;370 (9598):1560-1567.
  35. Raposa B, Pónusz R, Gerencsér G, et al. Food additives: Sodium benzoate, potassium sorbate, and sodium nitrite induce behavioral and biochemical impairments in animal models. Nutrients. 2016;8 (6):334.
  36. Al-Ghamdi MS. The effect of long-term ingestion of sodium nitrite on neurobehavioral functions and brain histology in rats. Saudi Med J. 2003;24 (12):1302-1307.
  37. Park S, Choi Y, Kim S. Neurotoxic consequences of chronic dietary exposure to synthetic phenolic antioxidants BHA and BHT. Food Chem Toxicol. 2019;128:221-230.
  38. Lanfranchi M, Maesano N, De Caro M. Multiple chemical preservative interactions in ultra-processed matrices: an emerging neurological challenge. Crit Rev Food Sci Nutr. 2020;60 (14):2409-2422.
  39. O'Brien PJ. Peroxidase-catalyzed oxidation of phenolics and its toxicological implications. Free Radic Biol Med. 2000;29 (11):1089-1099.
  40. De la Monte SM, Tong M. Mechanisms of nitrosamine-mediated neurodegeneration: potential role of dietary sodium nitrite. J Alzheimers Dis. 2009;17 (1):81-99.
  41. Hardy J, Selkoe DJ. The amyloid hypothesis of Alzheimer's disease: progress and problems on the road to therapeutics. Science. 2002;297 (5580):353-356.
  42. Braak H, Del Tredici K, Bratzke H, et al. Staging of the intracerebral inclusion body pathology associated with idiopathic Parkinson's disease (Braak's staging). J Neural Transm. 2003;110 (5):517-536.
  43. Danysz W, Parsons CG. Alzheimer's disease, beta-amyloid, NMDA receptors and memantine—searching for a common ground. Br J Pharmacol. 2012;167 (2):324-352.
  44. Bolognesi ML, Cavalli A, Valgimigli L, et al. Multi-target-directed drug design strategy in Alzheimer's disease: from theory to clinical reality. Curr Opin Chem Biol. 2018;44:27-33.
  45. Hopkins AL. Network pharmacology: the next paradigm in drug discovery. Nat Chem Biol. 2008;4 (11):682-690.
  46. Lipinski CA, Lombardo F, Dominy BW, Feeney PJ. Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv Drug Deliv Rev. 2001;46 (1-3):3-26.
  47. Berman HM, Westbrook J, Feng Z, et al. The Protein Data Bank. Nucleic Acids Res. 2000;28 (1):235-242.
  48. Morris GM, Huey R, Lindstrom W, et al. AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J Comput Chem. 2009;30 (16):2785-2791.
  49. Wang R, Lu Y, Wang S. Comparative evaluation of 11 scoring functions for molecular docking. J Med Chem. 2003;46 (12):2287-2303.
  50. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B Stat Methodol. 1995;57 (1):289-300.
  51. Ellman GL, Courtney KD, Andres V, Featherstone RM. A new and rapid colorimetric determination of acetylcholinesterase activity. Biochem Pharmacol. 1961;7 (2):88-95.
  52. Tsopmo A, Awad AM. Food-derived additives and neurodegenerative risk: an integrated toxicology and structural perspective. Trends Food Sci Technol. 2021;112:352-364.

Reference

  1. EFSA Panel on Food Additives and Nutrient Sources added to Food (ANS). Scientific Opinion on the re-evaluation of benzoic acid and its salts. EFSA J. 2016;14 (3):4433.
  2. World Health Organization, Joint FAO/WHO Expert Committee on Food Additives. Evaluation of certain food additives: eighty-fourth report of JECFA. WHO Tech Rep Ser. 2017;1007:1-84.
  3. Piper PW. Yeast safety and human health: the potential toxicity of the common food preservative sodium benzoate. Mol Biosyst. 2018;14 (2):292-300.
  4. Noorafshan A, Erfanizadeh M, Karbalay-Doust S. High dose of sodium benzoate induces neurotoxicity in the hippocampus of mice. Neurosci Lett. 2014;563:63-68.
  5. Khoshnood MJ, Hosseini A, Seydi E, Pourahmad J. Uncoupling of oxidative phosphorylation as a mechanism of sodium benzoate cytotoxicity. Toxicology. 2015;338:11-17.
  6. Lane HY, Lin CH, Green MF, et al. Add-on treatment of benzoate for schizophrenia: a randomized, double-blind, placebo-controlled trial. JAMA Psychiatry. 2013;70 (12):1267-1275.
  7. Lin CH, Chen PK, Wang SH, Lane HY. Sodium benzoate, a D-amino acid oxidase inhibitor, added to clozapine for the treatment of schizophrenia: a randomized, double-blind, placebo-controlled trial. Lancet Psychiatry. 2018;5 (5):406-415.
  8. Lane HY, Tsai GE. D-amino acid oxidase inhibitors in cognitive neuroscience and clinical psychiatry. Clin Psychopharmacol Neurosci. 2019;17 (3):341-352.
  9. Silman I, Sussman JL. Acetylcholinesterase: 'classical' and 'non-classical' functions and roles in neurodegenerative diseases. Physiol Rev. 2005;85 (4):1199-1222.
  10. Cheung J, Rudolph MJ, Burshteyn F, et al. Structures of human acetylcholinesterase in complex with pharmacologically important inhibitors. J Med Chem. 2012;55 (22):10282-10286.
  11. Binda C, Wang J, Pisani L, et al. Structures of human monoamine oxidase B complexes with selective noncovalent inhibitors. J Med Chem. 2007;50 (12):2775-2782.
  12. Edmondson DE, Binda C, Wang J, Upadhyay AK, Mattevi A. Molecular and structural basis of substrate and inhibitor selectivity of monoamine oxidases A and B. Curr Med Chem. 2009;16 (4):441-452.
  13. Karakas E, Furukawa H. Crystal structure of a heterotetrameric NMDA receptor ion channel. Science. 2014;344 (6187):992-997.
  14. Zhu S, Stein RA, Yoshioka C, et al. Mechanism of NMDA receptor inhibition and allosteric modulation. Nature. 2016;536 (7616):319-325.
  15. Trott O, Olson AJ. 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.
  16. Eberhardt J, Santos-Martins D, Tillack AF, Forli S. AutoDock Vina 1.2.0: New docking methods, expanded force field, and python bindings. J Chem Inf Model. 2021;61 (8):3891-3898.
  17. Dallakyan S, Olson AJ. Small-molecule library screening by docking with PyRx. Methods Mol Biol. 2015;1263:243-250.
  18. Daina A, Michielin O, Zoete V. SwissADME: a free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci Rep. 2017;7:42717.
  19. Pires DEV, Blundell TL, Ascher DB. pkCSM: predicting small-molecule pharmacokinetic and toxic properties using graph-based signatures. J Med Chem. 2015;58 (9):4066-4072.
  20. Banerjee P, Eckert AO, Schrey AK, Preissner R. ProTox-II: a webserver for the prediction of toxicity of chemicals. Nucleic Acids Res. 2018;46 (W1):W257-W263.
  21. Daina A, Michielin O, Zoete V. SwissTargetPrediction: updated data and new features for efficient prediction of protein targets of small molecules. Nucleic Acids Res. 2019;47 (W1):W357-W364.
  22. Stelzer G, Rosen N, Plaschkes I, et al. The GeneCards suite: from genomic information to robust knowledge of human disease. Curr Protoc Bioinformatics. 2016;54 (1):1.30.1-1.30.33.
  23. Piñero J, Ramírez-Anguita JM, Saüch-Pitarch J, et al. The DisGeNET knowledge platform for disease genomics: 2019 update. Nucleic Acids Res. 2020;48 (D1):D845-D855.
  24. Szklarczyk D, Gable AL, Nastou KC, et al. The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded datasets. Nucleic Acids Res. 2021;49 (D1):D605-D612.
  25. Shannon P, Markiel A, Ozier O, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13 (11):2498-2504.
  26. Chin CH, Chen SH, Wu HH, et al. cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst Biol. 2014;8 (Suppl 4):S11.
  27. Sherman BT, Hao M, Qiu J, et al. DAVID: a web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic Acids Res. 2022;50 (W1):W216-W221.
  28. Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28 (1):27-30.
  29. EFSA Panel on Food Additives and Nutrient Sources added to Food (ANS). Scientific Opinion on the re-evaluation of potassium sorbate (E 202). EFSA J. 2015;13 (6):4144.
  30. EFSA Panel on Food Additives and Nutrient Sources added to Food (ANS). Scientific Opinion on the re-evaluation of propionic acid and sodium propionate. EFSA J. 2014;12 (7):3779.
  31. EFSA Panel on Food Additives and Nutrient Sources added to Food (ANS). Re-evaluation of sodium nitrite (E 250) and potassium nitrite (E 249) as food additives. EFSA J. 2017;15 (6):4786.
  32. EFSA Panel on Food Additives and Nutrient Sources added to Food (ANS). Scientific Opinion on the re-evaluation of butylated hydroxyanisole - BHA (E 320) as a food additive. EFSA J. 2011;9 (10):2418.
  33. EFSA Panel on Food Additives and Nutrient Sources added to Food (ANS). Scientific Opinion on the re-evaluation of butylated hydroxytoluene - BHT (E 321) as a food additive. EFSA J. 2012;10 (3):2588.
  34. McCann D, Barrett A, Cooper A, et al. Food additives and hyperactive behaviour in 3-year-old and 8/9-year-old children in the community: a randomised, double-blinded, placebo-controlled trial. Lancet. 2007;370 (9598):1560-1567.
  35. Raposa B, Pónusz R, Gerencsér G, et al. Food additives: Sodium benzoate, potassium sorbate, and sodium nitrite induce behavioral and biochemical impairments in animal models. Nutrients. 2016;8 (6):334.
  36. Al-Ghamdi MS. The effect of long-term ingestion of sodium nitrite on neurobehavioral functions and brain histology in rats. Saudi Med J. 2003;24 (12):1302-1307.
  37. Park S, Choi Y, Kim S. Neurotoxic consequences of chronic dietary exposure to synthetic phenolic antioxidants BHA and BHT. Food Chem Toxicol. 2019;128:221-230.
  38. Lanfranchi M, Maesano N, De Caro M. Multiple chemical preservative interactions in ultra-processed matrices: an emerging neurological challenge. Crit Rev Food Sci Nutr. 2020;60 (14):2409-2422.
  39. O'Brien PJ. Peroxidase-catalyzed oxidation of phenolics and its toxicological implications. Free Radic Biol Med. 2000;29 (11):1089-1099.
  40. De la Monte SM, Tong M. Mechanisms of nitrosamine-mediated neurodegeneration: potential role of dietary sodium nitrite. J Alzheimers Dis. 2009;17 (1):81-99.
  41. Hardy J, Selkoe DJ. The amyloid hypothesis of Alzheimer's disease: progress and problems on the road to therapeutics. Science. 2002;297 (5580):353-356.
  42. Braak H, Del Tredici K, Bratzke H, et al. Staging of the intracerebral inclusion body pathology associated with idiopathic Parkinson's disease (Braak's staging). J Neural Transm. 2003;110 (5):517-536.
  43. Danysz W, Parsons CG. Alzheimer's disease, beta-amyloid, NMDA receptors and memantine—searching for a common ground. Br J Pharmacol. 2012;167 (2):324-352.
  44. Bolognesi ML, Cavalli A, Valgimigli L, et al. Multi-target-directed drug design strategy in Alzheimer's disease: from theory to clinical reality. Curr Opin Chem Biol. 2018;44:27-33.
  45. Hopkins AL. Network pharmacology: the next paradigm in drug discovery. Nat Chem Biol. 2008;4 (11):682-690.
  46. Lipinski CA, Lombardo F, Dominy BW, Feeney PJ. Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv Drug Deliv Rev. 2001;46 (1-3):3-26.
  47. Berman HM, Westbrook J, Feng Z, et al. The Protein Data Bank. Nucleic Acids Res. 2000;28 (1):235-242.
  48. Morris GM, Huey R, Lindstrom W, et al. AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J Comput Chem. 2009;30 (16):2785-2791.
  49. Wang R, Lu Y, Wang S. Comparative evaluation of 11 scoring functions for molecular docking. J Med Chem. 2003;46 (12):2287-2303.
  50. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B Stat Methodol. 1995;57 (1):289-300.
  51. Ellman GL, Courtney KD, Andres V, Featherstone RM. A new and rapid colorimetric determination of acetylcholinesterase activity. Biochem Pharmacol. 1961;7 (2):88-95.
  52. Tsopmo A, Awad AM. Food-derived additives and neurodegenerative risk: an integrated toxicology and structural perspective. Trends Food Sci Technol. 2021;112:352-364.

Photo
Ankita Thakur
Corresponding author

Department of Pharmaceutical Chemistry, Jagannath University, Chaksu, Jaipur

Photo
Neha
Co-author

Department of Zoology, Seth G.B. Podar College, Nawalgarh, Rajasthan

Photo
Shivansh Mishra
Co-author

Institute of Pharmaceutical Sciences, Lucknow University

Photo
Amrita Thakur
Co-author

Institute of Pharmaceutical Sciences, Lucknow University

Ankita Thakur, Amrita Thakur, Shivansh Mishra, Neha, In Silico Evaluation of Competitive Binding of Food Preservative Combinations at AChE, MAO-B and NMDA Receptors and Their Network Overlap with Neurodegeneration Pathways, Int. J. Sci. R. Tech., 2026, 3 (10), 644-651. https://doi.org/10.5281/zenodo.23276281

More related articles
Mobile App to Manage Food Distribution to Needy Pe...
Pavithra S. Ganapa, Dr. R. N. Kulkarni, Sindhu K., Pradivya G....
Beyond The Cage: Advancing Science In In-Vivo, In-Vitro & In-Silico Models...
Gayatri S. Jawanjal, Aishwarya M. Hiwase, Snehal N. Uike, Rajeshwar R. Hingane, Jogeshwari S Burunga...
In Silico Molecular Docking And ADMET Evaluation Of Coumarin Derivatives As Pote...
Vivekanand Ekanath Teli, Dipali Vijay Patil, Bhumika Bhairu Mohite, Rutuja Ashok Vyawahare...
Related Articles
Phytochemical Screening, GC-MS Profiling And Molecular Docking Studies On The Et...
Cholaraja K., Manoj Kumar S., Madhu Priya K.B.R., Madhu Bala S., Logesh A., Lavanya S....
In Silico Design, Molecular Docking And ADME-Toxicity Evaluation of Substituted ...
Rutuja Purushottam Bhojane, Manoj Gangadhar Shinde, Tejal Hari Bare, Janhavi Harish Bhalerao , Aksha...
Mobile App to Manage Food Distribution to Needy People...
Pavithra S. Ganapa, Dr. R. N. Kulkarni, Sindhu K., Pradivya G....