ORIGINAL RESEARCH article

Front. Chem., 30 May 2025

Sec. Medicinal and Pharmaceutical Chemistry

Volume 13 - 2025 | https://doi.org/10.3389/fchem.2025.1574702

Neuraminidase as a novel therapeutic management strategy for Alzheimer’s disease: evidenced through molecular docking, molecular dynamic simulation and gene expression analysis

  • 1. Department of Pharmacology, College of Pharmacy, Jouf University, Sakaka, Al-Jouf, Saudi Arabia

  • 2. King Salman Centre for Disability Research, Riyadh, Saudi Arabia

  • 3. Department of Pharmaceutics, College of Pharmacy, Jouf University, Sakaka, Al-Jouf, Saudi Arabia

  • 4. Department of Clinical Laboratory Sciences, College of Applied Medical Sciences, Jouf University, Sakaka, Al-Jouf, Saudi Arabia

  • 5. Department of Clinical Pharmacy, College of Pharmacy, Jouf University, Sakaka, Al-Jouf, Saudi Arabia

  • 6. Department of Pharmacology and Toxicology, College of Pharmacy, King Saud University, Riyadh, Saudi Arabia

  • 7. Department of Pharmaceutical Sciences, Pharmacy Program, Batterjee Medical College, Jeddah, Saudi Arabia

Abstract

Introduction:

Neuraminidase in humans is studied to see how well repurposed oseltamivir works for treating Alzheimer’s disease (AD) using methods like molecular docking, molecular dynamic (MD) simulation, and gene expression analysis. Gene enrichment analysis was also studied to understand the behaviour of neuraminidases in humans.

Methods:

Molecular docking was done using oseltamivir and the neuraminidase proteins with the PyRx tool, and the results were analysed using BIOVIA Discovery Studio. MD simulation (50 ns) of the oseltamivir and neuraminidase complex was performed using GROMACS tools. The gene expression analysis and gene enrichment study were done using GEO2R, which showed the results as log FC and significant values. Enricher tool-based gene enrichment analysis was done to determine the gene behaviour related to the AD.

Results:

The molecular docking showed a strong connection between oseltamivir and neuraminidase (−6.5 kcal/mol), acetylcholinesterase (−7.9 kcal/mol), CDKs (−6.5 kcal/mol), and GSKs (−6.6 kcal/mol), interacting with different amino acids in the protein sequences. MD simulations showed a strong interaction between the ligand and neuraminidase, with stable measurements indicating that both the protein and ligand remained consistent in size and energy, which is better explained through the results of MM_PBSA and MM_GBSA analysis of the complex, resulting in the ΔE_vdW, ΔE_elec, ΔG_polar, ΔG_nonpolar, ΔG_gas, (ΔE_vdW + ΔEEL), ΔG_solvation: (ΔG_polar + ΔG_nonpolar) and ΔG_bind: total energies suggesting the complex stayed stable in conditions similar to those resembling natural cell. The gene expression analysis expressed TUBB3 (formation of beta-tubulin), FABP3 (regulates alpha-synuclein uptake in dopaminergic neurons), and CALM1 (calcium signal transduction pathway) to be highly upregulated in the given conditions with kinase binding (p = 0.0006541) and protein phosphatase regulatory activity (p = 0.001357) were highly upregulated, implicating their importance in the AD.

Discussion:

The study ends on a hopeful note for using oseltamivir to treat neurological diseases, but it suggests that future research should include a solid cell line study, an in vitro study, and a clinical study.

Introduction

Alois Alzheimer pioneered the study of memory loss. He observed the existence of amyloid plaques and an immense harm to neurons during brain examinations of his very first patient who was suffering from cognitive loss and personality change before the patient died. Alois Alzheimer defined this particular condition as a disease of the cortical neurons in the cerebrum. Emil Kraepelin named this condition of serious concern as Alzheimer’s disease (AD), taken after the name of Alois Alzheimer (; ). Advanced reduction in cognition can be brought about by cerebro-cortical disorder like AD (). Other suggested causes are boozing, microorganism infections, anomaly of cardiovascular and/or pulmonary systems that can reduce oxygen supply to the brain, nutritional deficit, cyanocobalamin deficiency and cancers (; ; Rathod et al., 2016). Currently, about 50 million people are affected globally by AD. The numbers are expected to double every 5 years and to almost triple by 2050 (; ). AD overburden the affected individuals, their families and countries’ economies. In fact, the total health expenditure of AD alone as per researchers is over US$1 trillion per annum (). Until now, AD is incurable. The existing pharmacological interventions can only manage AD’s signs and symptoms (). Researchers classified AD as multifactorial disease; two main etiologies were hypothesised: the acetyl cholinergic and amyloid beta (Singh et al., 2024). At present, drugs which are approved to treat AD are acetylcholinesterase inhibitors and N-methyl d-aspartate (NMDA) antagonists (). Globally, researchers are trying to explore the mechanism involved behind the progressive neurodegeneration of cerebral cholinergic neurons, abnormal tau protein metabolism, beta-amyloid, inflammation and oxidative free radical damage. The aim is to develop promising pharmacological entities that can stop or modify the neurodegeneration in AD ()

Oseltamivir is an antiviral, which acts by inhibiting viral neuraminidase that is required for the virus to be released from the host cells. Researchers reported that oseltamivir reduces inflammation. It is one of the mechanisms suggested in the reduction of influenza symptoms () and in the production of antipyretic effect (Treanor et al., 2000). Oseltamivir is primarily used for the treatment of influenza A and B infections. Oseltamivir’s target viral enzyme neuraminidase cleaves the sialic acid residues present terminally on carbohydrate moieties of host cell membrane and influenza virus envelop (Treanor et al., 2000). This process breaks down the membrane barriers and causes the release of progeny viruses from the host cells. Through the same process, new virion infects other cells. The neuraminidase inhibitors bring about the static infectious condition and make available these virions to the macrophages intracellularly ().

Sialidases are recognised as analogues to viral neuraminidase in humans, and they have been termed as human neuraminidases (hNEU). Exploration of hNEU chemistry and pharmacology revealed to exhibit the same structure and function as that of viral neuraminidase (; ; ). Four isomers of sialidase, which are now known as hNEU, have been recognised, and they hydrolyse sialosides, gangliosides and glycoproteins in humans (Seo et al., 2021). These enzymes were found to play pivotal roles in different pathological states such as diabetes, (), cancers (Reddy et al., 2015) and neurodegenerative disorders (). hNEU-1 is the most abundant isoform of enzyme among all hNEU, which is present in lysosomes and plasma membrane and is a part of membrane complex (). hNEU-1 cleaves sialic acid from oligosaccharides and glycoproteins, and it has no effect on gangliosides. Sialidase (hNEU-2) localises in cytosol and catalyses sialic acid from verities of glycans. Sialidase (hNEU-3) is predominantly associated with plasma membrane and it selectively causes desialylation of gangliosides. Fourth isoform of sialidase (hNEU-4) exists mainly on membranes of cell organelles and possesses broad choice of glycan selectivity ().

Sialidases mediate cleavage reaction of glycoproteins, oligosaccharides and gangliosides. Sialidases also release sialic acid, which has been found to be involved in verities of biological processes including cancers, diabetes, neurodegeneration, inflammation and many more (Seo et al., 2021). Sialic acid releases play homeostatic as well as pathological roles. Sialic acid’s influence on neurons has been extensively studied by the researchers, and it has been found that sialidase mutations and generations of sialic acid are negatively associated with AD (Zhu et al., 2024). Sialic acid is required to activate many Siglec genes (primarily downregulated inflammation), which are required in stressful microglial environment and to protect the microglial cells (). However, this balance of microglial activation and inflammation disrupts in chronic cases and becomes the basis of neurodegeneration due to inflammatory upregulation (). Sialic acid on gangliosides interacts with Ca+2 and facilitates the synaptic transmission, which in excess can lead to excitoxicity (). Studies with increased serum sialic acid concentration in AD exhibited sialic to mediate the amyloid β (Aβ) association with gangliosides, which in turn leads to the Aβ plaque formation (Xiao et al., 2022). Sialic acid has been found to downregulate Siglec-3 (CD33). The stimulation of Siglec-3 is very well established in microglial protection ().

By reviewing the abovementioned literature and other related studies about the sialidase (hNEU), we came to hypothesise that hNEU can be a good target to inhibit the neurodegeneration from multiple fronts which are opened because of neuraminidase stimulation and liberation of sialic acid. Considering the fact that oseltamivir is a viral neuraminidase (vNEU) inhibitor and sialidases in humans are the structural analogues to (vNEU), we aimed to study the impact of neuraminidase inhibitor, i.e., oseltamivir against neuraminidase and different analogues of human neuraminidase including secretases, glycogen synthase kinase, cyclin dependent kinase, protein phosphatase 2A, sialyltransferase, acetylcholinesterase, interleukin-1B and tumour necrosis factor through preliminary pharmacological and toxicity profiling, molecular docking and molecular dynamic simulations. We aim to understand the role of different genes through gene expression analysis to establish the strength of interaction of oseltamivir against various enzymes mentioned and to shed light on its behaviour resembling cellular nature and its expression in the influence of different genes which are pivotal in human biology.

Methods

Preparing the protein for insilico processing

Proteins of interest, which play a significant role in the dysregulated metabolism, were identified as follows: neuraminidase (PDB ID: 2HTY), responsible for the entry of viral particles inside the healthy mammalian cells; (); beta secretase (PDB ID: 5YGX), responsible for the proteolytic cleavage of AD amyloid precursor; (Sambamurti et al., 2007); gamma secretase (PDB ID: 6YHF), responsible for the breakdown of the amyloid precursor protein to amyloid beta; (); alpha secretase (PDB ID: 6BE6), responsible for the breakdown of amyloid precursor protein to soluble amyloid precursor protein alpha; (); glycogen synthase kinase – 3 beta (GSK3B) (PDB ID: 1PYX), responsible for the regulation of signalling post endocytic transport; (); cyclin dependent kinase – 5 (CDK5) (PDB ID: 1H4L), responsible for organising the cytoskeleton and cellular growth; (Tian et al., 2022); protein phosphatase 2A (PP2A) (PDB ID: 1B3U), responsible for the cellular processes like autophagy, apoptosis, cell proliferation and DNA repair; (); sialyltransferase (PDB ID: 5BO6), responsible for the hydrolysis of sialic acid; (Zhang et al., 2010); acetylcholinesterase (PDB ID: 1ACJ), responsible for the metabolism of acetylcholine at the synapse of two different neurons; (Volkow et al., 2001); interleukin 1B (PDB ID: 1IOB), responsible for the endocrine and reproductive dysfunction; (); and, tumour necrosis factor (PDB ID: 1TNF), responsible for the modulation of various gene expression cellularly. (Schütze et al., 1992). The protein pdb files were downloaded from the RCSB PDB database (https://www.rcsb.org/) and () visualised through BIOVIA discovery studio () application, where the hetatms, water and co-crystals were removed and polar hydrogens were added to induce the charge around the protein structure to facilitate perfect interaction with the ligand in case of molecular docking and MD simulations.

Preparing oseltamivir as ligand

The structure of oseltamivir was downloaded in the .sdf format, which will be formatted as ligand of the study. Furthermore, it was downloaded from structural information rich databases like ‘PubChem’ directory (https://pubchem.ncbi.nlm.nih.gov/), () particularly oseltamivir (ethyl (3R,4R,5S)-4-acetamido-5-amino-3-pentan-3-yloxycyclohexene-1-carboxylate) (). The processing of the molecule for the conversion into the pdbqt was done by employing open babel tool ().

Identification of pharmacological and toxicological activity employing PASS tool

The prediction of activity spectra of substances (PASS) (https://www.way2drug.com/passonline/), an online tool, was selected to predict the possible activity of oseltamivir. The molecular SMILES formula was given as input in the interface of PASS, which predicts the possible pharmacological actions and toxicological activities. The structure-activity relationship was denoted by Pa: probability of activity and Pi: probability of inactivity. This approach provided an idea concerning the possible targets and activity for approaching through preparing the biomolecules of interest (enzymes) for further processing, employing molecular docking and MD simulation studies to develop an understanding with the evidence ().

Performing the proteins - ligand molecular docking

The processed protein crystal structures were uploaded into the PyRx interface (https://pyrx.sourceforge.io/) and converted into pdbqt format as a macromolecule of the study (). The ligands.sdf file was uploaded through the ‘Open Babel’ and was processed into pdbqt format to make ligand molecule for molecular docking after energy minimisation. The molecular docking was performed by creating a grid dimension (Table 2), and selecting complete protein. Twenty core CPUs with effectiveness were selected, and final docked molecules were saved for further processing. BIOVIA and discovery studios were the preferred methods to analyse the results obtained after molecular docking, including the 2D structures of proteins and ligands interactions, protein-ligand interactions and pocket identification (; ).

Preparing the enzyme and oseltamivir for MD simulation

Pymol interface was selected for converting the.dsv files into.pdb files. The separately imported protein and ligand files after molecular docking were joined through pymol. It was then saved in the form of pdb file for performing the MD simulation by means of GROMACS for a time period of 50 nanoseconds (ns), temperature of 300K, solvent–water, ions–sodium and chlorine, with one random seed and selected seed of seed_1234 for duality (; ; Sasumana and Kaushik, 2018; Varghese et al., 2024; Yuan et al., 2017; ).

The molecular dynamic (MD) simulation of a oseltamivir-neuraminidase complex was performed using GROMACS 2024, installed on linux (ubuntu) operating system, supported by the NIVIDIA RTX 4060 graphical processing unit and 16 GB RAM. Firstly, the topology files of the protein was generated using GROMACS compatible forcefield such as CHARMM36 and the ligand topology files were generated using Swissparam tools (; Zoete et al., 2011; Yesselman et al., 2012). Secondly, the ligand-protein complex was assembled and a simulation box was defined maintaining a minimum distance of 1.0 nm from the edges. TIP3P water model was chosen to solvate the simulation box was used. Sodium and chlorine was used as counterions to neutralise the system. To reduce steric clashes, energy minimising of the system was performed employing the steepest descent algorithm, followed by two equilibration steps such as NVT and NPT for a time period of 100 ps (picoseconds) to stabilise the temperature and pressure of the system using the V-rescale and Parrinello–Rahman algorithms. The production MD simulation run was conducted for 50 ns (nanoseconds) with a time step of 2 fs under stable boundary conditions saving the results every 10 ps.

Build in GROMACS utilities were used to perform post simulation analysis. The structural stability of the complex was assessed by calculating the root mean square deviation (RMSD) of the ligand and protein over time. While flexibility was evaluated using root mean square fluctuation (RMSF) per residue and the radius of gyration was assessed to determine the compactness of the protein. Binding free energies were carried employing the gmx_MMPBSA tool, externally installed using anaconda3 packages. The MD trajectory and topology files were used for processing gmx_MMPBSA (Molecular Mechanics Generalized Born Surface Area) and gmx_MMGBSA (Molecular Mechanics Poisson–Boltzmann Surface Area) binding free energies over represented snapshots. (). The energies included van der Waals, electrostatic, polar solvation and non-polar solvation contributions. This protocol provided insights into the stability, dynamics, and binding energetics of the ligand–protein complex.

Performing the gene expression analysis of sialyltransferase (analogue of neuraminidase)

Gene expression analysis was executed by means of GEO2R (https://www.ncbi.nlm.nih.gov/geo/geo2r/), an online tool of national centre for biotechnological information (NCBI) database (https://www.ncbi.nlm.nih.gov/) for variance expression analysis of huge datasets from gene expression omnibus (GEO). Datasets were normalised and comparisons between case and control groups were achieved. Significant genes were identified based on p-values (<0.05) and log 2-fold change thresholds for the gene analysis (). The GEO series accessions involving sialyltransferase genes expression data, i.e., GSE5281 (; ; Readhead et al., 2018; ), were selected for the study. The gene expression data was examined using dataset GSE5281, bearing the genes responsible for the development of AD and bearing the data of 61 control samples and 44 case samples. Differential expressional analysis was performed to recognise significant gene regulated within the groups. The study was completed using Benjamini and Hochberg method with NCBI generated data following the statistically significant level cut-off of p-values (<0.05), log 2-fold change threshold ‘zero’ with volcano and mean difference plot for control group vs. case group. The results are expressed in Table 4, bearing the data of twenty highly significant genes associated with the development or that played a crucial role in the pathway related to AD.

Gene enriching analysis

To identify the biological impact of twenty significant genes derived out of GEO analysis associated with AD, the gene enrichment analysis was performed employing the enrichr tool. The significant gene codes were given as input into the enrichr interface, (; ), which integrates the data from various databases. The analysis involves selecting relevant categories like gene ontology (GO), (), terms for cellular components, molecular functions and biological processes. Furthermore, disease specific annotations were explored employing Jensen diseases and Jensen tissue databases to understand the link of the genes with the AD pathology. Adjusted p-values (Padj) were used to ensure statistical significance. Combined scores, which integrate p-values and z-scores, were used to rank the enriched outcomes so that the results provide insights into the role of synaptic function, protein phosphorylation, lipid metabolism and neuroinflammation, which are critical in understanding the pathology of AD (Xie et al., 2021). The findings were visualised using volcano plots to understand and recognise key pathways and phenotypes for further investigation.

Results

The PASS prediction analysis points out several important pharmacological activities with high probable activities (Pa) and low possibility of inactivity (Pi), which suggests the oseltamivir is highly active against the targets like Neuraminidase inhibitor (Pa = 0.931) and Alpha-N-acetylglucosaminidase inhibitor (Pa = 0.455). The aforementioned targets could be of prime interest as the enzyme may be involved in the transferase reaction and may result in the production of abnormal proteins. Three different variants were identified based on the mutated gene expression, which could result in the development of neuronal disorders (). Inhibition of enzymes, highlighting its role in neuronal protection, misrepresentation or misexpression of the gene related to this enzyme, may result in neurological disorders, which is the point of interest of the current study Deficiency of the said enzyme may result in various other neurological disorders (Zhao and Neufeld, 2000). Additionally, it may result in toxicities like metabolic acidosis (0.511), sneezing, twitching, weight loss, hypocalcemia, dyspnea, gastrointestinal disturbance, ototoxicity, sensitization and thrombophlebitis with lower intensity, which are needed to be taken care of during the evaluation of the study. These activities express oseltamivir’s resourcefulness and therapeutic potential in diverse biological contexts related to Alzheimer’s disease (Table 1). (Annexure 1).

TABLE 1

Pharmacological activity prediction
PaPiPredicted activity
0.9310.000Neuraminidase (Influenza B) inhibitor
0.9320.001Antiviral (Influenza)
0.9200.001Antiviral (Influenza A)
0.6850.000Neuraminidase (influenza) inhibitor
0.6520.005Macrophage stimulant
0.4550.013Alpha-N-acetylglucosaminidase inhibitor
0.4730.048Muramoyltetrapeptide carboxypeptidase inhibitor
0.4610.042Antiviral (Rhinovirus)
0.4470.034Peptide-N4-(N-acetyl-beta-glucosaminyl) asparagine amidase inhibitor
0.4340.034Antimyopathies
Toxicity Prediction
PaPiPredicted activity
0.5110.063Acidosis, metabolic
0.4630.082Sneezing
0.5170.190Twitching
0.4170.110Weight loss
0.3590.065Hypocalcaemic
0.3750.097Dyspnea
0.3330.127Gastrointestinal disturbance
0.3260.129Ototoxicity
0.2590.084Sensitization
0.3630.195Thrombophlebitis

The activity prediction of oseltamivir.

Pa: probably activity, Pi: possibly inactivity, data derived through an online tool, i.e., PASS program (http://www.way2drug.com/PASSOnline/).

Molecular docking

The oseltamivir interacted with neuraminidase at a hydrogen bond length of 2.06Å with GLU C:174, expressing the lowest binding score of −6.5 kcal/mol at given grid values mentioned in Table 2. Further, the interaction between protein and ligand was formed even through Vander wall’s interaction with ARG C:172 along with six other amino acids. The oseltamivir also interacted with beta secretase at a hydrogen bond length of 1.86Å with TYR A:71, expressing the lowest binding score of −6.2 kcal/mol. Further, the interaction between protein and ligand was formed even through Vander wall’s interaction with ILE 118 along with nine other amino acids. The oseltamivir interacted with gamma secretase as well at a hydrogen bond length of 3.02Å with ALA 42, expressing the lowest binding score of −3.5 kcal/mol, a bit weaker interaction in comparison with other proteins. Further, the interaction between protein and ligand was formed even through Vander wall’s interaction with ILE 45 along with three other amino acids. Similarly, the oseltamivir interacted with alpha secretase at a hydrogen bond length of 2.14Å with ARG A:239, expressing the lowest binding score of - 6.4 kcal/mol. Further, the interaction between protein and ligand was formed even through Vander wall’s interaction with THR A:238 along with eight other amino acids. The oseltamivir interacted with glycogen synthase kinase 3 beta (GSK-3B) too at a hydrogen bond length of 2.05Å with GLY B:202, expressing the lowest binding score of - 6.6 kcal/mol. Further, the interaction between protein and ligand was formed even through Vander wall’s interaction with LYS A:292 along with six other amino acids. Additionally, the oseltamivir interacted with cyclin dependant kinase 5 (CDK-5) at a hydrogen bond length of 2.87Å with ASP A:86, expressing the lowest binding score of - 6.3 kcal/mol. Further, the interaction between protein and ligand was formed even through Vander wall’s interaction with 10 other amino acids. There was also an interaction of oseltamivir with protein phosphatase 2A (PP2A) at a hydrogen bond length of 1.94Å with ARG A:104, expressing the lowest binding score of - 5.5 kcal/mol. Further, the interaction between protein and ligand was formed even through Vander wall’s interaction with 12 other amino acids. Correspondingly, the oseltamivir interacted with sialyltransferase at a hydrogen bond length of 2.12Å with THR B:194, expressing the lowest binding score of - 6.4 kcal/mol. Further, the interaction between protein and ligand was formed even through Vander wall’s interaction with six other amino acids. In the same manner, the oseltamivir interacted with acetylcholinesterase at a hydrogen bond length of 2.55Å with TRP A:84, expressing the lowest binding score of - 7.9 kcal/mol, which is the highest binding affinity of the current study. Further, the interaction between protein and ligand was formed even through Vander wall’s interaction with 19 other amino acids. Adding to the list, the oseltamivir interacted also with interleukin 1B at a hydrogen bond length of 2.02Å with LEU A:62, expressing the lowest binding score of - 5.3 kcal/mol. Further, the interaction between protein and ligand was formed even through Vander wall’s interaction with nine other amino acids. Lastly, the oseltamivir interacted with tumour necrosis factor - alpha at a hydrogen bond length of 2.37Å with SER C:99, expressing the lowest binding score of - 7.4 kcal/mol, which is also one of the best results of the current study. Further, the interaction between protein and ligand was formed even through Vander wall’s interaction with 14 other amino acids. (Table 2) (Figure 1).

TABLE 2

Name of ligand (Oseltamivir)Name of proteinVarious receptor enzymes associated with Alzheimer’s as target proteins and oseltamivir as ligands
Highest affinity score
Kcal/mol
RMSDIonic or wander wall’s interactions (Å units)Hydrogen bond interactions (Å units)Amino acidsGrid dimensions
Centre Dimensions
Oseltamivir Neuraminidase (PDB ID: 2HTY)
−6.50ARG C:172, LEU A:127, SER A:101, ASN C:208, GLU C:209, GLU A:128, THR C:1912.30
2.06
2.47
2.82
PHE C:173, GLU C:174, LYS C:206, TYR A:100x = 1.5994
y = 50.1814
z = 100.559
x = 103.0816
y = 102.4211
z = 53.8510
Beta secretase (PDB ID: 5YGX)
−6.20ILE 118, ARG 128, SER36, GLY 34, TYR 198, ILE 226, ASP 228, VAL 332, ASP 32, THR 2311.86
2.15
TYR A:71
GLY A:230
x = 14.3015
y = 41.0494
z = 0.1888
x = 56.1419
y = 64.6465
z = 45.9909
Gamma secretase (PDB ID: 6YHF)
−3.50ILE 45, GLY 38, LEU 34, MET 353.02ALA 42x = - 2.2973
y = −3.2008
z = 11.7049
x = 26.2145
y = 48.5674
z = 13.6442
Alpha secretase (PDB ID: 6BE6)
−6.40THR A:238, LYS A:518, GLU A:240, ILE A:243, SER A:506, SER A:504, CYS A:473, ASP A:481, LYS A:4802.14ARG A:239x = 31.7668
y = 33.1965
z = 18.0663
x = 94.5299
y = 109.4674
z = 123.8035
GSK-3B (PDB ID: 1PYX)
−6.60LYS A:292, ASN B:95, GLN B:89, VAL B:87, LYS B:85, SER B:203, PHE B:672.51
2.57
2.05
ARG B:96
GLU B:97
GLY B:202
x = 24.2110
y = −0.3042
z = 21.4266
x = 72.9735
y = 80.3676
z = 103.6807
CDK-5 (PDB ID: 1H4L)
−6.30GLU A:81, ALA A:143, ASP A:84, GLN A:85, ILE A:10, GLN A:130, ASN A:131, ASP A:144, LYS A:33, VAL A:182.87ASP A:86x = 35.4578
y = −32.768
z = 26.1504
x = 77.4488
y = 60.5521
z = 49.0611
PP2A (PDB ID: 1B3U)
−5.50THR A:141, ALA A:183, SER A:186, LYS A:187, GLY A:148, CYS A:143, VAL A:108, SER A:145, THR A:101, ARG A:182, MET A:179, PHE A:1402.54
2.37
1.94
ASP A:105
THR A:144
ARG A:104
x = 34.6630
y = 46.7673
z = 23.0447
x = 97.3458
y = 107.1759
z = 122.2308
Sialyltransferase (PDB ID: 5BO6)
−6.40ALA B:192, LEU B:238, GLU B:195, SER B:234, ASN B; 230, PRO B:1902.12THR B:194x = −3.5210
y = 9.8904
z = −11.0690
x = 61.7766
y = 59.4714
z = 102.3049
Acetylcholine-esterase (PDB ID: 1ACJ)
−7.90GLU A: 199, GLY A:441, TYR A:130, HIS A:440, GLY A:117, TYR A:442, LEU A:127, SER A:124, GLY A:118, GLY A:123, SER A:122, ASN A:85, GLN A:69, TYR A:121, ASP A:72, GLY A:80, TRP A:432, SER A:81, TYR A:3342.55TRP A:84x = 4.7576
y = 65.5208
z = 56.8290
x = 64.6127
y = 61.3047
z = 56.9408
Interleukin 1B (PDB ID: 1IOB)
−5.30VAL A:85, ASN A:66, TYR A:90, ASN A:7, SER A:5, SER A:43, GLY A:61, GLU A:64, LYS A:652.02
2.39
LEU A:62
LYS A:63
x = 15.9224
y = 13.5665
z = 1.1895
x = 44.2915
y = 37.7782
z = 44.89.24
TNF alpha (PDB ID: 1TNF)
−7.40LYS C:98, PRO B:110, LYS A:98, GLU C:116, GLU B:116, SER B:99, GLU A:116, PRO A:100, TYR A:115, TRP A:114, CYS A:110, PRO C:100, GLN C:102, GLU C:1042.37
3.02
SER C:99
GLN B:102
x = 19.9745
y = 49.6645
z = 39.9343
x = 62.7877
y = 63.9782
z = 61.4350

Protein (Neuraminidase) -Ligand (Oseltamivir) Docking Results of ‘PyRx Python Prescription 0.8’ Analysed Through ‘BIOVIA-Discovery Studios Visualiser v.21.1.0.20298’.

FIGURE 1

The summary of the dual MD simulation (random simulation and seed_1234) of the oseltamivir-neuraminidase complex is explained in Table 4, which consist a of binding free energies components as integral part of MM-PBSA analysis. Significant high magnitude of van der Waals interactions (ΔE_vdW) changed form the negligible phase in pre-dynamic MM-PBSA (−0.00 ± 1.26 for random and 0.00 ± 0.40 for seed_1234) to a considerable negative value post-simulation (−20.22 ± 0.73 and −19.05 ± 0.50), which indicates high binding confirmation and attractive ness between the ligand and protein. Similar observations were noticed for electrostatic interactions (ΔE_elec), the energy quotient shifts from negligible to highly promising post-dynamic i.e., −284.37 ± 0.95 for random and −255.35 ± 1.46 for seed_1234, these results confirm the strong bonding between the ligand and protein. However, in case of polar solvation energy (ΔG_polar) the post simulation results of random seed was increased and seed_1234 was decreased, showing a system determined correction was pointed to the seed_1234 simulation being more favourable. The total gas-phase energy (ΔG_gas), transforms in to significantly favourable value post simulation, indicating stronger intermolecular connections in the complex environment. Finally, the binding free energies (ΔG_bind) change from the positive in pre-MM-PBSA value to negative in post MM-PBSA value exhibits improved binding characteristics post molecular dynamics.

The results MM-GBSA analysis of pre and post dual MD simulation of random seed and seed_1234 for determining their binding free energies is expressed in Table 5, the van der Waals interaction energy (ΔE_vdW) demonstrated shift from undermined values in pre dynamic simulation (0.00 ± 0.00 random and −0.00 ± 0.00 seed_1234) to recognisable post dynamic simulation (−18.78 ± 0.01 for random and −19.18 ± 0.78 for seed_1234) values representing better packing between ligand and protein post simulation. Electrostatic interactions (ΔE_elec) showed similar results, shifting from insignificant range of values to desirable values (−266.20 ± 0.46 for random and −259.86 ± 1.81 for seed_1234) post simulation. Conclusively, the total binding free energies (ΔG_bind) remained constantly negative in both pre and post MM-GBSA states, suggesting the stable and promising binding between the ligand and protein.

The MM-GBSA decomposition analysis for the ligand and protein complex, pre and post MD simulations in random and fixed (seed_1234) conditions revealed that key amino acid fragments such as GLU:276 and GLU:277 developed strong interactions in pre-MD simulation analysis with values −10.36 and −11.98 kcal/mol under random conditions. A great shift in the energies of the amino acids was identified in post MD simulation, i.e., GLU:276 reaching −15.61 kcal/mol and GLU:277 displaying a positive value of 3.87 kcal/mol, highlighting important changes. Residues like ARG:224 and ARG:292 better binding energies post MD simulations and LEU:134, TRP:178, SER:179, and ILE:222 maintained constant negative energies throughout simulations.

The interaction of ligand, oseltamivir, and enzyme protein, neuraminidase, was analysed through molecular dynamic simulation employing GROMACS tool. Various crucial parameters have been identified such as RMSD, radius of gyration (Rg), Vander Waal’s surface area (VSA) and the potential energy of the complex (PE). The protein RMSD ranged for Protein: P: 0.128424–0.493372 nm, ligand 0.122013–0.838708 nm under random seed simulation and under seed_1234 simulation it ranged from Protein 0.161276–2.858013 nm and Ligand 0.208401–3.217241 nm, which explains that the ligand-protein complex has undergone minimal structural disorientation and exhibited stable binding throughout the test. The radius of gyration remained constant between for protein at a range of 3.391135–3.518412 nm and for the ligand at a range of 0.354091–0.37452 nm under random simulation conditions, however, the protein and ligand showed the radius of gyration of about 3.390434–4.271447 nm and 0.341893–0.371202 nm under seed_1234 simulation conditions which explains minimal disorientation in seed_1234 simulation. The RMSF of protein and ligand under random and seed_1234 MD simulation conditions exhibited a value of P: 0.0451–0.2189 nm, L: 0.0153–0.1925 nm for random simulation and P: 1.2986–2.2136 nm and ligand L: 0.0160–0.2210 nm under seed_1234 simulation conditions exhibiting vigorous molecular interactions between the ligand and protein in a complex, which justify the results expressed in Tables 2, 3. Through this current simulation results, it is evident that the oseltamivir forms a stable and favourable bond with the enzyme protein like neuraminidase, which underscores the oseltamivir as a suitable candidate for development of pharmacophore against the enzyme protein of neuraminidase (Figure 2).

TABLE 3

SeedProteinLigandRoot mean square deviation (RMSD) nmRoot mean square fluctuations (RMSF) nmRadius of gyration (Rg) nm
Random seedNeuraminidaseOseltamivirP: 0.128424–0.493372
L: 0.122013–0.838708
P: 0.0451–0.2189
L: 0.0153–0.1925
P: 3.391135–3.518412
L: 0.354091–0.37452
Seed_1234NeuraminidaseOseltamivirP: 0.161276–2.858013
L: 0.208401–3.217241
P: 1.2986–2.2136
L: 0.0160–0.2210
P: 3.390434–4.271447
L: 0.341893–0.371202

Molecular dynamic (MD) simulation of oseltamivir and neuraminidase.

nm = Nanometres, P = protein, L = ligand.

FIGURE 2

Through these results, it can be explained that targeting the neuraminidase can ensure effective modulation of pathological pathway related to AD (), as neuraminidase or similar enzymes like sialic acid modifying agents are involved in neurodegenerative pathways (Zhang et al., 2023).

The dataset validation performed through GEO employing the accession number GSE5281 and focusing on the gene expression related to the Alzheimer’s disease for sialyltransferase, which is an analogue of neuraminidase in humans, has shown significant differential expression results. Based on their low p-value and high log fold-change (Log FC) ranging between 6.25 × 10−10 to 3.92 × 10−8, significant statistical evidence for changes related to gene expression was found. The log FC values of 1 shows two folds, 1.5 three folds and 2 four folds in the gene upregulation in the given dataset. The significantly expressed genes like tubulin beta 3 class III (TUBB3) and synaptosome associated protein 25 (SNAP25) are important for normal neuronal structure and function even at the synaptic level (; ), the processes which are dysregulated in AD. The four folds upregulation of ubiquitin C-terminal hydrolase L1 (UCHL1) and PPP2R2D, a subunit of protein phosphatase 2A, suggests altered signalling pathways and protein turnover in human brain (), which are associated with neurodegenerative diseases. The neurexin 3 (NRXN3) and calmodulin 1 (CALM1) are associated with synaptic connectivity, calcium signalling and disrupting, which are highly associated with AD. (; ) Other genes like fatty acid binding protein 3 (FABP3), (), which represents the lipid metabolism, and somatostatin (SST) (Samson et al., 2008) have neuromodulatory function highly affecting the cognitive processes. Dysregulation of which may develop the AD. (Samson et al., 2008) Aggregation of certain protein and inflammation of developing genes such as reticulon 3 (RTN3) () and C1q and tumour necrosis factor-related protein 4 (C1QTNF4) (Vester et al., 2021) can act as central to the development and prognosis of AD. Through this analysis, the network of dysregulated genes involved in synaptic function, neuroinflammation, calcium function and protein homeostasis can be explained. The significant alteration in the gene expression and homeostasis emphasises potential targets for novel therapeutic intervention in relation to AD. (Table 4) The volcano plot expresses the differential expressing gene between Alzheimer’s disease samples and normal samples obtained from the dataset GSE5281 (Figure 3), which shows distinct upregulated and downregulated genes in the context of AD. In the log FC values, positive values (red) show upregulated genes while negative values (blue) show downregulated genes. The y-axis (p-value) shows the significance of the study samples, where the non-significant genes are shown in grey colour. These genes may be involved in inflammation, synaptic dysfunction, amyloid metabolism or impaired neuroprotection. These findings provide potential biomarkers for the development of novel therapeutic agents. (Annexure II).

TABLE 4

Energy componentsPre MM-PBSA (Mean ± SD); random seedPost MM-PBSA (Mean ± SD); random seedPre MM-PBSA (Mean ± SD); Seed_1234Post MM-PBSA (Mean ± SD); Seed_1234
ΔE_vdW (van der Waals)−0.00 ± 1.26−20.22 ± 0.730.00 ± 0.40−19.05 ± 0.50
ΔE_elec (Electrostatic)−0.00 ± 2.22−284.37 ± 0.950.00 ± 3.20−255.35 ± 1.46
ΔG_polar (Polar solvation)246.03 ± 17.02289.19 ± 1.33267.15 ± 11.36263.47 ± 0.31
ΔG_nonpolar (Non-polar solvation)−3.56 ± 0.03−3.25 ± 0.04−3.57 ± 0.03−3.31 ± 0.05
ΔG_gas (Energy in vacuum)−229.62 ± 20.44−304.59 ± 1.20−269.62 ± 6.13−274.40 ± 1.56
ΔG_solvation (Total solvation energy)242.47 ± 17.02285.95 ± 1.33263.58 ± 11.36260.16 ± 0.31
ΔG_bind (Total Binding Free Energy)12.84 ± 26.60−18.65 ± 1.79−6.04 ± 12.91−14.24 ± 1.59

Results of pre and post MM-PBSA binding free energy components in kcal/mol of the MD simulated molecular complex of Oseltamivir and Neuraminidase.

ΔE_vdW: van der Waals interaction energy between the ligand and the protein.

ΔE_elec: electrostatic interaction energy.

ΔG_polar: solvation free energy due to polar interactions, calculated using the Poisson–Boltzmann (PB) or Generalized Born (GB) models.

ΔG_nonpolar: non-polar solvation free energy, often estimated from the solvent-accessible surface area (SASA).

ΔG_gas: Total energy in vacuum (ΔE_vdW + ΔEEL).

ΔG_solvation: Total solvation energy (ΔG_polar + ΔG_nonpolar).

ΔG_bind: total binding free energy, sum of the above components.

FIGURE 3

The enricher-based analysis of 20 highly significant genes has revealed an enriching analysis of genes implicated in Alzheimer’s disease, highlighting significant association among the cellular components, molecular function, mammalian phenotypes, tissues and related diseases. Within the cellular components, vesicles (adjusted p-value = 0.04964, combined score = 115.27) and glutamatergic synapses (adjusted p-value = 0.04964, combined score = 200.25) emerged as critical, emphasising the known roles of synaptic dysfunction in AD. In molecular function, kinase binding (adjusted p-value = 0.03205, combined score = 89.70) and protein phosphatase regulator activity (adjusted p-value = 0.03324, combined score = 278.70) highlighted the disruption in protein phosphorylation as a tau-related AD precursor (). Long chain fatty acid binding and phosphoserine residue binding showed strong results (adjusted p-value = 0.04964 for both, combined scores = 828.65), suggesting metabolic and post-translational monitoring dysregulation. In the case of mammalian phenotype level, the abnormal peripheral nervous system synaptic transmission with values (adjusted p-value = 0.008775, combined score = 3208.63) and abnormal endplate potential with values (adjusted p-value = 0.01166, combined score = 1658.66) were significant, which is consistent with the impact of AD on the synaptic transmission. Decreased paired-pulse facilitation with values (adjusted p-value = 0.01781, combined score = 961.11) proves impaired synaptic plasticity, further assisting neurodegeneration. From the tissue specific analysis, the cortex of the brain (adjusted p-value = 0.0001872, combined score = 511.81) and the ganglia with values (adjusted p-value = 0.0007525, combined score = 280.10) showed significant involvement, which are shown to be critically affected in AD. Within the disease association, strong links were noticed related to brain disease, pineoblastoma, dumping syndrome and dysembryoplastic neuroepithelial tumour, hinting at overlapping pathways of neurodegeneration. The analysis has given a hint related to the metabolic causes of the development of AD, where neuraminidase analogues like sialyltransferase plays a crucial role in disrupting the metabolism of various proteins and lipids. Hence, the enricher-based analysis of highly significant genes has revealed various potential targets for the novel development of therapeutic agents (Table 5).

TABLE 5

Energy componentsPre MM-GBSA (Mean ± SD); random seedPost MM-GBSA (Mean ± SD); random seedPre MM-GBSA (Mean ± SD); Seed_1234Post MM-GBSA (Mean ± SD); Seed_1234
ΔE_vdW (van der Waals)0.00 ± 0.00−18.78 ± 0.01−0.00 ± 0.00−19.18 ± 0.78
ΔE_elec (Electrostatic)0.00 ± 0.00−266.20 ± 0.46−0.00 ± 0.00−259.86 ± 1.81
ΔG_polar (Polar solvation, EGB)242.01 ± 5.04252.84 ± 0.07249.12 ± 6.73250.85 ± 0.52
ΔG_nonpolar (Non-polar solvation, ESURF)−4.21 ± 0.51−3.94 ± 0.04−4.67 ± 0.12−4.51 ± 0.08
ΔG_gas (Gas-phase energy)−270.30 ± 3.88−284.98 ± 0.48−279.86 ± 3.25−279.04 ± 1.97
ΔG_solvation (Total solvation)237.80 ± 5.06248.91 ± 0.08244.45 ± 6.73246.34 ± 0.52
ΔG_bind (Total binding energy)−32.50 ± 6.38−36.07 ± 0.49−35.41 ± 7.47−32.70 ± 2.04

Results of pre and post MM-GBSA binding free energy components in kcal/mol of the MD simulated molecular complex of Oseltamivir and Neuraminidase.

ΔE_vdW: van der Waals interaction energy between the ligand and the protein.

ΔE_elec: electrostatic interaction energy.

ΔG_polar: solvation free energy due to polar interactions, calculated using the Poisson–Boltzmann (PB) or Generalized Born (GB) models.

ΔG_nonpolar: non-polar solvation free energy, often estimated from the solvent-accessible surface area (SASA).

ΔG_gas: Total energy in vacuum (ΔE_vdW + ΔEEL).

ΔG_solvation: Total solvation energy (ΔG_polar + ΔG_nonpolar).

ΔG_bind: total binding free energy, sum of the above components.

Discussion

Oseltamivir, which inhibits neuraminidase and facilitates the viral entry into the cells, disarranges sialic acid from the glycoproteins. The sialic acid is also produced by the action of various human enzymes, such as secretases, glycogen synthesis kinases, cycling dependent kinases, protein phosphatase 2A, sialyltransferase, acetylcholinesterase, interleukin 1B and tumour necrosis factor, which plays a vital role either in the formation of insoluble sialic acid or in the production of neurodegeneration through inflammation. Hence, the current study was taken up to analyse the inhibitory potential of oseltamivir against all of the enzymes or cytokines related to Alzheimer’s disorders.

The PASS online tool was used to predict the pharmacological activity of a compound. It was predicted that oseltamivir functions as neuraminidase inhibitor and that it was a preferred antiviral agent against influenza A and B viruses. The ‘Pa’ values exceeding 0.9 were well-aligned with documented role of oseltamivir. A study published by Patricia et al. highlighted the antiviral activity of oseltamivir against influenza (Schirmer and Holodniy, 2009). The findings of this study confirmed the mechanism of action, which helped us in choosing the required topic.

However, the secondary activities such as macrophages stimulation () and alpha-N-acetylglucosaminidase inhibition had to settle with lower Pa values, which could be considered as the off-target effects and related to oseltamivir. The secondary activities are not yet well-developed or published, which opens up a potential change of studying the molecular-related activities concerning other diseases for the patient’s safety ().

The toxicity prediction has expressed that the oseltamivir may develop metabolic acidosis (Pa = 0.511), twitching (Pa = 0.517) and dyspnoea (Pa = 0.375). These effects were also noted in the clinical setting with higher doses or with prolonged use of oseltamivir. A published meta-analysis highlights gastrointestinal disturbances and no neuropsychiatric effects of oseltamivir in patients. Such studies confirm the application of oseltamivir with proper dose management (). Certain predicted toxicities, like thrombophlebitis and ototoxicity, lack substantial clinical evidences. To overcome these limitations of the insilico studies, a necessary clinical validation is recommended. Further, a published study explains that the adverse events related to oseltamivir are milder in intensity, which can be managed clinically (Smith et al., 2011). Nevertheless, the concerns raised by the predictive models as a part of insilico studies cannot be counted less. Hence, a thorough discussion and drug repurposing must be considered (). In order to avoid the adverse drug reactions or drug interactions of the oseltamivir with those that are milder in action, the drug can be formulated into a lipid or polymer-based carriers as targeted drug. Delivery system implicates its action directly in the targets like neuraminidase analogues or acetylcholinesterase (; ; ; ) for the application of novel molecules clinically in order to improve the quality of life of those with neurological disorders (Srujana et al., 2017).

The molecular docking studies that employed the PyRx revealed a wide range of binding affinities of oseltamivir to various proteins, included those which are of importance in Alzheimer’s disease. One such target of oseltamivir was acetylcholinesterase (−7.9 kcal/mol), followed by TNF-alpha (−7.4 kcal/mol), GSK-3B (−6.6 kcal/mol) and neuraminidase (−6.5 kcal/mol). These findings suggest the potential off-target interactions of oseltamivir outside its main mechanism of action and antiviral target (Walczak-Nowicka and Herbet, 2021).

Binding affinity to acetylcholinesterase and neuraminidase is particularly important as it highlights the possible role of oseltamivir in controlling cholinergic pathway and neurodegeneration, which explains the involvement of acetylcholinesterase in the development of Alzheimer’s disease (Walczak-Nowicka and Herbet, 2021). The observations of the current study provide high binding affinity of the oseltamivir with acetylcholinesterase; the said action can be validated through the pre-clinical or clinical investigations. Furthermore, the binding affinity of oseltamivir with neuraminidase of −6.5 kcal/mol was consistent with the role of oseltamivir as neuraminidase inhibitor. The observed amino acid interaction with ARG C:172 and GLU A: 128 confirms stable ligand binding aligning with structural studies of Larisa et al () Comparatively, the binding energies with beta and alpha secretase suggest weaker but potentially significant interaction with amyloid precursors protein-processing enzymes. These findings align with the results provided by Kioke et al., that antiviral agents may provide neuroprotective effect through off-target interactions, i.e., sialidase ().

The stability of neuraminidase-oseltamivir complex was confirmed by molecular dynamic simulations studies, as evident RMSD values for Protein: 0.128424–0.493372 nm, ligand 0.122013–0.838708 nm under random seed simulation and for seed_1234 simulation it ranged from Protein 0.161276–2.858013 nm and Ligand 0.208401–3.217241 nm confirming stability of the complex. The result of radius of gyration (for protein at a range of 3.391135–3.518412 nm and for the ligand at a range of 0.354091–0.37452 nm and 3.390434–4.271447 nm and 0.341893–0.371202 nm for seed_1234 simulation) further supports the structural stability of the complex. These findings can be related to the results of a research study published by Putra et al. (2018), which demonstrated similar RMSD for the neuraminidase inhibition forming protein-ligand complex. These results are similar with other earlier published computational studies, affirming strong interaction of oseltamivir with the primary target, i.e., neuraminidase. Similarly, the RMSF values P: 0.0451–0.2189 nm, L: 0.0153–0.1925 nm for random simulation and P: 1.2986–2.2136 nm and ligand L: 0.0160–0.2210 nm under seed_1234 explains extensive molecular interactions, which further validates the oseltamivir’s potency as neuraminidase inhibitor ().

Moreover, with regards to the acetylcholinesterase, the binding energies (-7.9 kcal/mol) also stands as an add-on to the off-target effects of oseltamivir. Combination of inhibition of neuraminidase and acetylcholinesterase in humans may prove to be beneficial in the management of neurological disorders like Alzheimer’s disease ().

Analysis of the results obtained post MD simulation MM-PBSA and MM-GBSA revealed increased binding affinities of the oseltamivir with Neuraminidase complex. Important energy components such as van der Waals (ΔE_vdW), electrostatic (ΔE_elec) interactions have shown a shift form the negligible (0 ± 0.00) to a favourable value in both productions and trajectories. Which indicates a stronger ligand-protein interactions post MD simulations. On the other hand, polar solvent energy (ΔG_polar) was identified with variable changes, and non-polar energies remains consistently constant. The total binding energies (ΔG_bind) became significantly favourable in random as well as seed_1234 simulations, which confirms improved complex stability. Decomposition level MM-GBSA analysis highlighted critical moieties like GLU:376, GLU:277 and LIG:469 with great shifts in the residual energies of ARG224 and ARG292 post MD simulations. These finding emphasize the value of MD simulations in capturing realistic energetic and structural refinements in ligand binding (Tables 5, 6).

TABLE 6

ResidueRandom seed pre-simulation resultsRandom seed post-simulation resultsSeed_1234 pre-simulation resultsSeed_1234 post-simulation results
Total Avg.Total Std. Dev.Total Avg.Total Std. Dev.Total Avg.Total Std. Dev.Total Avg.Total Std. Dev.
ARG:1180.16070967320.69705330.07313636410.867111430.486262812.272165310.16070967320.6970533
GLU:1190.14521705516.972887710.31266414519.516036180.096205613.843254060.14521705516.97288771
LEU:134−0.1003636365.004673121−0.0500909095.228447377−0.09454.94543138−0.1003636365.004673121
GLN:136−0.0009545457.7632184470.00455.52189690.0183.367739887−0.0009545457.763218447
ASP:1510.2980938918.9505987190.2159144737.731790941.04443964.992619370.2980938918.950598719
ARG:152−1.30888447313.319519080.52826596414.93811591−1.76584244.183743335−1.30888447313.31951908
ARG:1560.17711934511.311635940.10020214517.29397933−0.05319810.65081420.17711934511.31163594
TRP:178−0.7819319276.321986742−0.3167817824.945035206−0.7796341.907511047−0.7819319276.321986742
SER:179−0.1779925.067683233−0.2527096737.109482835−0.20320086.822918681−0.1779925.067683233
ALA:180−0.1265454553.673033655−0.0584545454.181619527−0.1185.358094251−0.1265454553.673033655
ILE:194−0.2177272735.611584371−0.1281363644.297757474−0.1514.03092049−0.2177272735.611584371
ILE:222−0.8060777454.719285829−1.3339929825.225018381−0.72422564.429854294−0.8060777454.719285829
LEU:2230.0375909094.2163637290.0624545453.4453867510.0661.0010684290.0375909094.216363729
ARG:2241.98446745511.95442862.42031003610.852856350.97649164.9832512771.98446745511.9544286
THR:225−0.3150341825.214389646−0.0964692365.560143728−0.0393.686086622−0.3150341825.214389646
GLU:227−1.12966098216.401068780.25111501812.70363795−0.62259645.087889271−1.12966098216.40106878
SER:2460.01786498216.358766770.0868314557.2700880860.08870522.3338317960.01786498216.35876677
ASN:247−0.0074029096.693001681−0.0157272735.9824052790.007726413.05948338−0.0074029096.693001681
HIS:274−0.0369545456.725702928−11.951824659.375090448−0.05153.514541713−0.0369545456.725702928
GLU:276−10.3573434510.12830086−15.612463613.24161797−10.57005843.694668202−10.3573434510.12830086
GLU:277−11.978354259.9609640893.87411145513.02749532−11.203475612.73485215−11.978354259.960964089
ARG:2923.61485021812.26243563−0.470602810.756493732.503060814.283426553.61485021812.26243563
ASN:2940.0260629098.481015303−0.2265560367.736542115−1.1376423.210434790.0260629098.481015303
TYR:347−0.6177384736.642141102−0.1721511274.10205902−0.9540446.307583365−0.6177384736.642141102
GLY:348−0.1183181826.1185482820.1515454558.723824196−0.02553.935544276−0.1183181826.118548282
ARG:3710.16268181811.276445020.26539821817.655789630.66295485.7265311470.16268181811.27644502
TYR:406−0.0489089095.856747652−0.4242908366.0797617460.20297125.262203844−0.0489089095.856747652
PRO:431−0.00854.373965825−0.0149545454.81802992−0.0471.539013945−0.00854.373965825
LIG:469−12.3937535611.253468380.0190909099.532751049−14.18274884.869457648−12.3937535611.25346838

Results of pre and post MM-GBSA binding free energy components in kcal/mol of the MD simulated molecular complex of Oseltamivir and Neuraminidase.

The genes, which play a vital role in the TUBB3 (formation of beta-tubulin), FABP3 (regulates alpha-Synuclein uptake in dopaminergic neurons) and CALM1 (Calcium signal transduction pathway) in Alzheimer’s disease, were upregulated. These findings align with the earlier studies published by Shu et al., which reported similar gene expression patterns in AD. The upregulation of neuronal structure and calcium signalling pathways provides further insights into the molecular basics of AD (Wang et al., 2012). Oseltamivir’s potential to modulate the gene expression, particularly in sialyltransferase-related pathways, requires further investigation (). The current study does not directly demonstrate the ligands’ impact on the genes. The observed docking ability of ligand against sialyltransferase (−6.4 kcal/mol) suggests a possible role in controlling glycosylation process related to neurodegeneration. This aligns with the study published by Jannis et al., which highlights the importance of sialylation in brain function and pathological development (Wißfeld et al., 2024) (Tables 7, 8).

TABLE 7

Accession numberGroupsIDAdj.
P- value
Log FCGene symbolGene title
GSE5281
Focusing on the Alzheimer’s disease
Control groups: 61 samples
Case group: 44 samples
213476_x_at6.25E-101.999TUBB3Tubulin beta 3 class III
202257_s_at9.34E-102.022CD2BP2CD2 cytoplasmic tail binding protein 2
230656_s_at1.00E-092.163UTP4UTP4, small subunit processome component
207232_s_at1.00E-092.027DZIP3DAZ interacting zinc finger protein 3
205738_s_at5.07E-092.028FABP3Fatty acid binding protein 3
213921_at8.05E-092.759SSTSomatostatin
202785_at8.40E-092.172NDUFA7NADH: ubiquinone oxidoreductase subunit A7
223708_at1.30E-082.636C1QTNF4C1q and tumour necrosis factor related protein 4
221772_s_at1.44E-081.926PPP2R2DProtein phosphatase 2 regulatory subunit B delta
1568603_at1.56E-082.250CADPSCalcium dependent secretion activator
201387_s_at1.63E-081.956UCHL1Ubiquitin C-terminal hydrolase L1
202507_s_at1.72E-082.319SNAP25Synaptosome associated protein 25
215020_at1.94E-082.147NRXN3Neurexin 3
1568604_a_at2.08E-082.130CADPSCalcium dependent secretion activator
219549_s_at2.58E-082.027RTN3Reticulon 3
213710_s_at2.84E-082.127CALM1Calmodulin 1
200638_s_at3.00E-081.948YWHAZTyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein zeta
219408_at3.15E-081.989PRMT7Protein arginine methyltransferase 7
215021_s_at3.44E-082.031NRXN3Neurexin 3
209953_s_at3.92E-081.893CDC37Cell division cycle 37

Gene expression analysis of sialyltransferase gene with accession number GSE5281.

TABLE 8

IndexNameP-valueAdjusted p-valueOdds RatioCombined score
Cellular components
1Vesicle0.0011470.0496417.03115.27
2Glutamatergic Synapse0.0022060.0496432.74200.25
Molecular function
1Kinase Binding0.00065410.0320512.2389.70
2Protein Phosphatase Regulator Activity0.0013570.0332442.21278.70
3Adenylate Cyclase Regulator Activity0.0044920.04964293.791588.07
4Long-Chain Fatty Acid Binding0.0071790.04964167.86828.65
5Phosphoserine Residue Binding0.0071790.04964167.86828.65
6Omega Peptidase Activity0.0080720.04964146.87707.80
Mammalian phenotype level 4
1Abnormal PNS Synaptic Transmission0.000034280.008775312.093208.63
2Abnormal Endplate Potential0.000091120.01166178.291658.66
3Decreased Paired-Pulse Facilitation0.00020870.01781113.41961.11
Jensen tissue
1Brain cortex cell line9.549e-70.000187236.92511.81
2Ganglion0.0000076780.000752523.78280.10
3Interstitial cell of Cajal0.000023610.00154329.76317.08
4Temporal lobe0.00042180.015166.8853.46
5Parietal lobe0.00058330.015168.1460.61
Jensen diseases
1Brain disease0.0021510.0392933.18203.78
2Pineoblastoma0.0053880.03929235.021227.64
3Dumping syndrome0.0071790.03929167.86828.65
4Dysembryoplastic neuroepithelial tumour0.0071790.03929167.86828.65

The results of significant gene expression from the data of 20 highly upregulated and significantly expressed genes.

Functional enrichment analysis revealed significant association with molecular function, such as kinase binding (p = 0.0006541) and protein phosphatase regulatory activity (p = 0.001357). These findings explore the role of oseltamivir in modulating signal pathways indirectly through protein interactions. Additionally, the action on the cellular components, such as vesicle and glutaminergic synapse, underlines the potential impact of oseltamivir on synaptic transmission. Furthermore, the association of oseltamivir’s activity with brain cortex cell lines and disease such as dysembryoplastic neuroepithelial tumour (p = 0.007179) raises a point regarding the higher implications of the oseltamivir use in neurological conditions. These findings also align with the earlier published studies by Pamela et al., which reported off-target effects on the neural pathways. However, a thorough experimental validation is still required to justify these activities obtained through in silico studies ().

Conclusion

The current study aims to explore the impact of oseltamivir as an inhibitor on neuraminidase and its analogues in humans, along with the significant genes and pathways involved in the AD. The results highlight the binding affinities of oseltamivir with neuraminidase, acetylcholinesterase and sialyltransferase, along with the importance of synaptic functions, vesicle transport and kinase-regulating genes, suggesting their pivotal role in the progression of AD. The binding affinities, synaptic functions, vesicle transport and kinase regulating are prior expressed in various research publications linking to cognition decline and development of AD. The identified enzymes can be considered as important molecular targets for developing novel drugs to treat AD in humans. However, a small number of target proteins and insilico methods employed in the study needs to be validated experimentally. Future research anchoring on the same idea should focus on the understanding of molecular mechanisms and clinical validations to solidify the findings with the therapeutic potential of these targets in AD. The current study provides a deeper understanding of the various patterns related to AD, which could help in the future development of novel agents.

Limitations and future directions

This study provides comprehensive yet limited insights regarding oseltamivir’s pharmacology and toxicology profile. The insilico predictions and docking studies require experimental validations to confirm their applicability in clinical settings. Additionally, the observed gene expression analysis and pathway association are correlative and do not establish causality.

Future studies should focus on experimental validation of the pharmacological and toxicological activities predicted related to oseltamivir, along with its off-target interaction. For oseltamivir repurposing, both in vitro and in vivo studies are required to confirm its effect on the neurological pathways and to correct the neurological diseases. Additionally, integrating multi-omics data such as proteomics and metabolomics could possibly elucidate the molecular mechanisms highlighting oseltamivir’s pharmacological effects. However, through this study, the potential pharmacological, toxicological and gene regulatory activities of the oseltamivir are uncovered as neuraminidase inhibitor. Likewise, the potential secondary activities and off-target interactions are also uncovered. These findings provide a broader therapeutic applications and safety profile of oseltamivir. Moreover, experimental validation is being considered in collaboration with renowned research institutes and we hope to explore this aspect in the upcoming series works.

Statements

Data availability statement

The data can be found in the given link below: https://figshare.com/s/e2b28ab2503735cd6302.

Author contributions

SA: Conceptualization, Formal Analysis, Funding acquisition, Writing – original draft. OA: Conceptualization, Investigation, Resources, Writing – original draft. HA: Writing – review and editing. AA: Methodology, Project administration, Writing – review and editing. TA: Conceptualization, Formal Analysis, Investigation, Writing – review and editing. MeA: Writing – original draft. MuA: Methodology, Visualization, Writing – original draft, Writing – review and editing. MJ: Formal Analysis, Methodology, Software, Writing – original draft.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. The authors extend their appreciation to the King Salman center For Disability Research for funding this work through Research Group no KSRG-2024-340.

Acknowledgments

The authors extend their appreciation to the King Salman center For Disability Research for funding this work through Research Group no KSRG-2024-340.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.

Generative AI statement

The author(s) declare that Generative AI was used in the creation of this manuscript. During the preparation of this work the authors used SciSpace (https://typeset.io/) in order to improve the quality of english language, and to align the names of drugs with the text. After using this tool/service, the authors reviewed and edited the content as needed and takes full responsibility for the content of the publication.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2025.1574702/full#supplementary-material

References

  • 1

    AbdullahiM.UzairuA.ShallangwaG. A.MamzaP. A.IbrahimM. T.ChandraA.et al (2023). Unveiling 1,3-thiazine derivative as a potential neuraminidase inhibitor: molecular docking, molecular dynamics, ADMET and DFT studies. Chem. Afr.6 (6), 29572967. 10.1007/s42250-023-00713-4

  • 2

    AlroyJ.García-MolinerM. L.LeeR. E. (2014). “The pathology of the skeleton in lysosomal storage diseases,” in Pathobiology of human disease. Editors McManusL. M.MitchellR. N. (San Diego: Academic Press), 874892.

  • 3

    AshburnerM.BallC. A.BlakeJ. A.BotsteinD.ButlerH.CherryJ. M.et al (2000). Gene Ontology: tool for the unification of biology. Nat. Genet.25 (1), 2529. 10.1038/75556

  • 4

    BaazaouiN.IqbalK. (2022). Alzheimer’s disease: challenges and a therapeutic opportunity to treat it with a neurotrophic compound. Biomolecules12 (10), 1409. 10.3390/biom12101409

  • 5

    BaraskarK.ThakurP.ShrivastavaR.ShrivastavaV. K. (2021). Female obesity: association with endocrine disruption and reproductive dysfunction. Obes. Med.28, 100375. 10.1016/j.obmed.2021.100375

  • 6

    BarorohS. S. M. B. U.MuscifaZ. S.DestiaraniW.RohmatullahF. G.YusufM. (2023). Molecular interaction analysis and visualization of protein-ligand docking using Biovia Discovery Studio Visualizer. Indonesian J. Comput. Biol. (IJCB)2 (1), 22. 10.24198/ijcb.v2i1.46322

  • 7

    BartholdD.JoyceG.FeridoP.DraboE. F.MarcumZ. A.GrayS. L.et al (2020). Pharmaceutical treatment for alzheimer's disease and related dementias: utilization and disparities. J. Alzheimers Dis.76 (2), 579589. 10.3233/jad-200133

  • 8

    BekkerH.BerendsenH. J. C.DijkstraE. J.AchteropS.VondrumenR.VanderspoelD.et al (1993). “Gromacs - a parallel computer for molecular-dynamics simulations,” in 4th international conference on computational physics (PC 92). Editors DeGrootR. A.NadrchalJ. (World Scientific Publishing: SINGAPORE), 252256.

  • 9

    BermanH. M.WestbrookJ.FengZ.GillilandG.BhatT. N.WeissigH.et al (2000). The protein data bank. Nucleic Acids Res.28 (1), 235242. 10.1093/nar/28.1.235

  • 10

    BirdN. L.OlsonM. R.HurtA. C.OshanskyC. M.OhD. Y.ReadingP. C.et al (2015). Oseltamivir prophylaxis reduces inflammation and facilitates establishment of cross-strain protective T cell memory to influenza viruses. PLOS ONE10 (6), e0129768. 10.1371/journal.pone.0129768

  • 11

    BishopP.RoccaD.HenleyJ. M. (2016). Ubiquitin C-terminal hydrolase L1 (UCH-L1): structure, distribution and roles in brain function and dysfunction. Biochem. J.473 (16), 24532462. 10.1042/bcj20160082

  • 12

    BreijyehZ.KaramanR. (2020). Comprehensive review on Alzheimer’s disease: causes and treatment. Molecules25 (24), 5789. 10.3390/molecules25245789

  • 13

    BugnonM.GoullieuxM.RöhrigU. F.PerezM. A. S.DainaA.MichielinO.et al (2023). SwissParam 2023: a modern web-based tool for efficient small molecule parametrization. J. Chem. Inf. Model.63 (21), 64696475. 10.1021/acs.jcim.3c01053

  • 14

    BythewayI.DarleyM. G.PopelierP. L. A. (2008). The calculation of polar surface area from first principles: an application of quantum chemical topology to drug design. ChemMedChem3 (3), 445453. 10.1002/cmdc.200700262

  • 15

    CarrollC. M.LiY.-M. (2016). Physiological and pathological roles of the γ-secretase complex. Brain Res. Bull.126, 199206. 10.1016/j.brainresbull.2016.04.019

  • 16

    ChahinianH.YahiN.FantiniJ. (2024). Glutamate, gangliosides, and the synapse: electrostatics at work in the brain. Int. J. Mol. Sci.25 (16), 8583. 10.3390/ijms25168583

  • 17

    ChenE. Y.TanC. M.KouY.DuanQ.WangZ.MeirellesG. V.et al (2013). Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinforma.14 (1), 128. 10.1186/1471-2105-14-128

  • 18

    ChenS.CaoZ.NandiA.CountsN.JiaoL.PrettnerK.et al (2024). The global macroeconomic burden of Alzheimer's disease and other dementias: estimates and projections for 152 countries or territories. Lancet Glob. Health12 (9), e1534e1543. 10.1016/s2214-109x(24)00264-x

  • 19

    ChinD.MeansA. R. (2000). Calmodulin: a prototypical calcium sensor. Trends Cell Biol.10 (8), 322328. 10.1016/s0962-8924(00)01800-6

  • 20

    CorboC. P.AlonsoA. d.C. (2011). “Chapter 2 - therapeutic targets in alzheimer's disease and related tauopathies,” in Progress in molecular biology and translational science. Editor RahmanS. (Academic Press), 4783.

  • 21

    CoutsiasE. A.WesterM. J. (2019). RMSD and symmetry. J. Comput. Chem.40 (15), 14961508. 10.1002/jcc.25802

  • 22

    DallakyanS.OlsonA. J. (2015). “Small-molecule library screening by docking with PyRx,” in Chemical biology: methods and protocols. Editors HempelJ. E.WilliamsC. H.HongC. C. (New York, NY: Springer New York), 243250.

  • 23

    DeTureM. A.DicksonD. W. (2019). The neuropathological diagnosis of Alzheimer’s disease. Mol. Neurodegener.14 (1), 32. 10.1186/s13024-019-0333-5

  • 24

    DuJ.ShuiH.ChenR.DongY.XiaoC.HuY.et al (2024). Neuraminidase-1 (NEU1): biological roles and therapeutic relevance in human disease. Curr. Issues Mol. Biol.46 (8), 80318052. 10.3390/cimb46080475

  • 25

    DulyA. M. P.KaoF. C. L.TeoW. S.KavallarisM. (2022). Βiii-tubulin gene regulation in health and disease. Front. Cell Dev. Biol.10, 851542. 10.3389/fcell.2022.851542

  • 26

    DzulkoM.PonsM.HenkeA.SchneiderG.KrämerO. H. (2020). The PP2A subunit PR130 is a key regulator of cell development and oncogenic transformation. Biochimica Biophysica Acta (BBA) - Rev. Cancer1874 (2), 188453. 10.1016/j.bbcan.2020.188453

  • 27

    EaswaranV.AlmeleebiaT. M.MantargiM. J. S.KhanN. A.AlshahraniS. M.OrayjK.et al (2023b). Patient safety culture in the southern region of Saudi arabia: a survey among community pharmacies. Healthcare11 (10), 1416. 10.3390/healthcare11101416

  • 28

    EaswaranV.KhanN. A.IqbalM. J.AlshahraniS. M.OrayjK.AlmeleebiaT. M.et al (2023a). The study of healthcare professionals' perspective towards the quality of diabetic care services in Abha. Eur. Rev. Med. Pharmacol. Sci.27 (10), 43284336. 10.26355/eurrev_202305_32437

  • 29

    FangQ.GaoY.ChenM.-f.GuoX.-l.YangX.YangX.et al (2014). Molecular epidemiology and evolution of A(H1N1)pdm09 and H3N2 virus during winter 2012–2013 in Beijing, China. J. Mol. Epidemiol. Evol. Genet. Infect. Dis.26, 228240. 10.1016/j.meegid.2014.05.034

  • 30

    FaridM.JaffarM. (2019). Creation and assessment of buccal mucoadhesive sustained release oral films. Int. J. Pharm. Pharm. Res.16 (2), 3041.

  • 31

    FaridM.JaffarM. (2020). Formulation and evaluation of ofloxacin loaded chitosan nanoparticles. WORLD J. Pharm. Pharm. Sci.9 (10), 113.

  • 32

    FerreiraP. C. L.De BastianiM. A.BellaverB.PovalaG.BrumW. S.RamosV. G.et al (2020). Functional enrichment analysis of differentially expressed genes leads to dysregulation in biological processes networks in Alzheimer’s disease continuum. Alzheimer's & Dementia16 (S2), e039506. 10.1002/alz.039506

  • 33

    ForsythE.RitzlineP. D. (1998). An overview of the etiology, diagnosis, and treatment of alzheimer disease. Phys. Ther.78 (12), 13251331. 10.1093/ptj/78.12.1325

  • 34

    FouadM.ZakiA.ShabanaS.El-HaggarR.HaggarM.JaffarS. (2017). Evaluation of oral hypoglycemic potency of Medicago polymorpha and Zygophyllum simplex: a Drug-Drug interaction study. Medicine6 (6), 648651.

  • 35

    GaoC.JiangJ.TanY.ChenS. (2023). Microglia in neurodegenerative diseases: mechanism and potential therapeutic targets. Signal Transduct. Target. Ther.8 (1), 359. 10.1038/s41392-023-01588-0

  • 36

    GermoushM. O.MantargiM. J. S.SarhanM.AlrashdiB. M.MassoudD.AlzwainS.et al (2024). Molecular docking of eleven snake venom peptides targeting human immunodeficiency virus capsid glycoprotein as inhibitors. Open Veterinary J.14 (11), 29362949. 10.5455/ovj.2024.v14.i11.22

  • 37

    GlanzV. Y.MyasoedovaV. A.GrechkoA. V.OrekhovA. N. (2018). Inhibition of sialidase activity as a therapeutic approach. Drug Des. Devel Ther.12, 34313437. 10.2147/dddt.s176220

  • 38

    GreenM. D.NetteyH.WirtzR. A. (2008). Determination of oseltamivir quality by colorimetric and liquid chromatographic methods. Emerg. Infect. Dis. J.14 (4), 552556. 10.3201/eid1404.061199

  • 39

    GrumatiP.MorozziG.HölperS.MariM.HarwardtM.-L. I. E.YanR.et al (2017). Full length RTN3 regulates turnover of tubular endoplasmic reticulum via selective autophagy. eLife6, e25555. 10.7554/elife.25555

  • 40

    GubarevaL.MohanT. (2022). Antivirals targeting the neuraminidase. Cold Spring Harb. Perspect. Med.12 (1), a038455. 10.1101/cshperspect.a038455

  • 41

    HamaR. (2016). The mechanisms of delayed onset type adverse reactions to oseltamivir. Infect. Dis.48 (9), 651660. 10.1080/23744235.2016.1189592

  • 42

    HataK.KosekiK.YamaguchiK.MoriyaS.SuzukiY.YingsakmongkonS.et al (2008). Limited inhibitory effects of oseltamivir and zanamivir on human sialidases. Antimicrob. Agents Chemother.52 (10), 34843491. 10.1128/aac.00344-08

  • 43

    HoltzmanD. M.MorrisJ. C.GoateA. M. (2011). Alzheimer’s disease: the challenge of the second century. Sci. Transl. Med.3 (77), 77sr1. 10.1126/scitranslmed.3002369

  • 44

    IqbalK.Liu F Fau - GongC. X.Gong Cx Fau - Grundke-IqbalI.Grundke-IqbalI. (2010). Tau in Alzheimer disease and related tauopathies. Curr. Alzheimer Res.7 (8), 656664. 10.2174/156720510793611592

  • 45

    JaffarM.FaridM. (2019). Formulation and in vivo evaluation of risperidone spherical agglomerates prepared for earlier absorption of risperidone. World J. Pharm. Res.8 (13), 14301444.

  • 46

    JheeS. S.YenM.EreshefskyL.LeibowitzM.SchulteM.KaeserB.et al (2008). Low penetration of oseltamivir and its carboxylate into cerebrospinal fluid in healthy Japanese and caucasian volunteers. Antimicrob. Agents Chemother.52 (10), 36873693. 10.1128/aac.00327-08

  • 47

    KalyanG. U.ReddyD. P.ChandrakanthG.PoojaB.AnithaV.VivekD. (2023). “Gene association disease prediction by GEO2R tool,” in 2023 international conference on evolutionary algorithms and soft computing techniques (EASCT).

  • 48

    KeilJ. M.RafnG. R.TuranI. M.AljohaniM. A.Sahebjam-AtabakiR.SunX.-L. (2022). Sialidase inhibitors with different mechanisms. J. Med. Chem.65 (20), 1357413593. 10.1021/acs.jmedchem.2c01258

  • 49

    KelleyB. J.PetersenR. C. (2007). Alzheimer's disease and mild cognitive impairment. Neurol. Clin.25 (3), 577609. 10.1016/j.ncl.2007.03.008

  • 50

    KhatuaB.RoyS.MandalC. (2013). Sialic acids siglec interaction: a unique strategy to circumvent innate immune response by pathogens. Indian J. Med. Res.138 (5), 648662.

  • 51

    KhunnawutmanothamN.SooknualP.BatsomboonP.PloypradithP.ChimnoiN.PatigoA.et al (2024). Synthesis, antiacetylcholinesterase activity, and molecular dynamics simulation of aporphine–benzylpyridinium conjugates. ACS Med. Chem. Lett.15 (1), 132142. 10.1021/acsmedchemlett.3c00467

  • 52

    KimS.ChenJ.ChengT.GindulyteA.HeJ.HeS.et al (2024). PubChem 2025 update. Nucleic Acids Res.53, D1516D1525. 10.1093/nar/gkae1059

  • 53

    KuleshovM. V.JonesM. R.RouillardA. D.FernandezN. F.DuanQ.WangZ.et al (2016). Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res.44 (W1), W90W97. 10.1093/nar/gkw377

  • 54

    LeeS.-M.LeeS. H.JungY.LeeY.YoonJ. H.ChoiJ. Y.et al (2020). FABP3-mediated membrane lipid saturation alters fluidity and induces ER stress in skeletal muscle with aging. Nat. Commun.11 (1), 5661. 10.1038/s41467-020-19501-6

  • 55

    LiangW. S.DunckleyT.BeachT. G.GroverA.MastroeniD.RamseyK.et al (2008b). Altered neuronal gene expression in brain regions differentially affected by Alzheimer's disease: a reference data set. Physiol. Genomics33 (2), 240256. 10.1152/physiolgenomics.00242.2007

  • 56

    LiangW. S.DunckleyT.BeachT. G.GroverA.MastroeniD.WalkerD. G.et al (2007). Gene expression profiles in anatomically and functionally distinct regions of the normal aged human brain. Physiol. Genomics28 (3), 311322. 10.1152/physiolgenomics.00208.2006

  • 57

    LiangW. S.ReimanE. M.VallaJ.DunckleyT.BeachT. G.GroverA.et al (2008a). Alzheimer's disease is associated with reduced expression of energy metabolism genes in posterior cingulate neurons. Proc. Natl. Acad. Sci.105 (11), 44414446. 10.1073/pnas.0709259105

  • 58

    LivingstonG.HuntleyJ.SommerladA.AmesD.BallardC.BanerjeeS.et al (2020). Lancet. Lancet396 (10248), 413446. 10.1016/s0140-6736(20)30367-6

  • 59

    MageshS.SavitaV.MoriyaS.SuzukiT.MiyagiT.IshidaH.et al (2009). Human sialidase inhibitors: design, synthesis, and biological evaluation of 4-acetamido-5-acylamido-2-fluoro benzoic acids. Bioorg. & Med. Chem.17 (13), 45954603. 10.1016/j.bmc.2009.04.065

  • 60

    MageshS.SuzukiT.MiyagiT.IshidaH.KisoM. (2006). Homology modeling of human sialidase enzymes NEU1, NEU3 and NEU4 based on the crystal structure of NEU2: hints for the design of selective NEU3 inhibitors. J. Mol. Graph. Model.25 (2), 196207. 10.1016/j.jmgm.2005.12.006

  • 61

    MaloshR. E.MartinE. T.HeikkinenT.BrooksW. A.WhitleyR. J.MontoA. S. (2017). Efficacy and safety of oseltamivir in children: systematic review and individual patient data meta-analysis of randomized controlled trials. Clin. Infect. Dis.66 (10), 14921500. 10.1093/cid/cix1040

  • 62

    MauriceP.BaudS.BocharovaO. V.BocharovE. V.KuznetsovA. S.KaweckiC.et al (2016). New insights into molecular organization of human neuraminidase-1: transmembrane topology and dimerization ability. Sci. Rep.6 (1), 38363. 10.1038/srep38363

  • 63

    McAuleyJ. L.GilbertsonB. P.TrifkovicS.BrownL. E.McKimm-BreschkinJ. L. (2019). Influenza virus neuraminidase structure and functions. Front. Microbiol.10, 39. 10.3389/fmicb.2019.00039

  • 64

    McKimm-BreschkinJ. L. (2013). Influenza neuraminidase inhibitors: antiviral action and mechanisms of resistance. Influenza Other Respir. Viruses7 (s1), 2536. 10.1111/irv.12047

  • 65

    MillerB. R.IIIMcGeeT. D.Jr.SwailsJ. M.HomeyerN.GohlkeH.RoitbergA. E. (2012). MMPBSA.py: an efficient program for end-state free energy calculations. J. Chem. Theory Comput.8 (9), 33143321. 10.1021/ct300418h

  • 66

    MiyagiT. (2008). “Aberrant expression of sialidase in cancer and diabetes,” in Experimental glycoscience: glycobiology. Editor TaniguchiN., (Tokyo: Springer Japan), 329332.

  • 67

    MuM.SunH.GengS.XuT.SunC.ZhangZ.et al (2024). Neurexin-3 in the paraventricular nucleus of the hypothalamus regulates body weight and glucose homeostasis independently of food intake. Mol. Brain17 (1), 49. 10.1186/s13041-024-01124-3

  • 68

    MuruganN. A.PodobasA.GadioliD.VitaliE.PalermoG.MarkidisS. (2022). A review on parallel virtual screening softwares for high-performance computers. Pharmaceuticals15 (1), 63. 10.3390/ph15010063

  • 69

    NicholsE.SteinmetzJ. D.VollsetS. E.FukutakiK.ChalekJ.Abd-AllahF.et al (2022). Estimation of the global prevalence of dementia in 2019 and forecasted prevalence in 2050: an analysis for the Global Burden of Disease Study 2019. Lancet Public Health7 (2), e105e125. 10.1016/s2468-2667(21)00249-8

  • 70

    O'BoyleN. M.BanckM.JamesC. A.MorleyC.VandermeerschT.HutchisonG. R. (2011). Open Babel: an open chemical toolbox. J. Cheminformatics3 (1), 33. 10.1186/1758-2946-3-33

  • 71

    OlssonB.ZetterbergH.HampelH.BlennowK. (2011). Biomarker-based dissection of neurodegenerative diseases. Prog. Neurobiol.95 (4), 520534. 10.1016/j.pneurobio.2011.04.006

  • 72

    PartonenT. (2023). “Medication effects,” in Encyclopedia of sleep and circadian rhythms. Editor KushidaC. A.Second Edition) (Oxford: Academic Press), 755760.

  • 73

    PasseriE.ElkhouryK.MorsinkM.BroersenK.LinderM.TamayolA.et al (2022). Alzheimer’s disease: treatment strategies and their limitations. Int. J. Mol. Sci.23 (22), 13954. 10.3390/ijms232213954

  • 74

    PoroikovV.FilimonovD.LaguninA.GloriozovaT.ZakharovA. (2007). PASS: identification of probable targets and mechanisms of toxicity†. SAR QSAR Environ. Res.18 (1-2), 101110. 10.1080/10629360601054032

  • 75

    PuigdellívolM.AllendorfD. H.BrownG. C. (2020). Sialylation and galectin-3 in microglia-mediated neuroinflammation and neurodegeneration. Front. Cell. Neurosci.14, 162. 10.3389/fncel.2020.00162

  • 76

    PutraR.ImaniastutiR.NasutionM.KeramiD.TambunanU. (2018). “Molecular docking simulation of neuraminidase influenza A subtype H1N1 with potential inhibitor of disulfide cyclic peptide (DNY, NNY, LRL),” in IOP conference series: materials science and engineering (Bristol, United Kingdom: IOP Publishing).

  • 77

    RathodR.KaleA.JoshiS. (2016). Novel insights into the effect of vitamin B12 and omega-3 fatty acids on brain function. J. Biomed. Sci.23 (1), 17. 10.1186/s12929-016-0241-8

  • 78

    ReadheadB.Haure-MirandeJ.-V.FunkC. C.RichardsM. A.ShannonP.HaroutunianV.et al (2018). Multiscale analysis of independent Alzheimer’s cohorts finds disruption of molecular, genetic, and clinical networks by human herpesvirus. Neuron99 (1), 6482.e7. 10.1016/j.neuron.2018.05.023

  • 79

    ReddyY. P.ChandrasekharK. B.Mohammed JaffarS. (2015). A study of nigella sativa induced growth inhibition of MCF and HepG2 cell lines: an anti-neoplastic study along with its mechanism of action. Pharmacogn. Res.7 (2), 193. 10.4103/0974-8490.150541

  • 80

    SambamurtiK. (2007). “Beta secretase,” in xPharm: the comprehensive pharmacology reference. Editors EnnaS. J.BylundD. B. (New York: Elsevier), 18.

  • 81

    SamsonW. K.ZhangJ. V.Avsian-KretchmerO.CuiK.YostenG. L. C.KleinC.et al (2008). Neuronostatin encoded by the somatostatin gene regulates neuronal, cardiovascular, and metabolic functions. J. Biol. Chem.283 (46), 3194931959. 10.1074/jbc.m804784200

  • 82

    SasumanaJ.KaushikV. (2018). “Molecular dynamics simulation by using NAMD-VMD and gromacs,” in Journal of emerging technologies and innovative research.

  • 83

    SchirmerP.HolodniyM. (2009). Oseltamivir for treatment and prophylaxis of influenza infection. Expert Opin. Drug Saf.8 (3), 357371. 10.1517/14740330902840519

  • 84

    SchützeS.MachleidtT.KrönkeM. (1992). Mechanisms of tumor necrosis factor action. Semin. Oncol.19 (2 Suppl. 4), 1624.

  • 85

    SeoN.LeeH.OhM. J.KimG. H.LeeS. G.AhnJ. K.et al (2021). Isomer-specific monitoring of sialylated N-glycans reveals association of α2,3-linked sialic acid epitope with behcet’s disease. Front. Mol. Biosci.8, 778851. 10.3389/fmolb.2021.778851

  • 86

    SinghM. K.ShinY.JuS.HanS.KimS. S.KangI. (2024). Comprehensive overview of Alzheimer’s disease: etiological insights and degradation strategies. Int. J. Mol. Sci.25 (13), 6901. 10.3390/ijms25136901

  • 87

    SmithJ. R.RaynerC. R.DonnerB.WollenhauptM.KlumppK.DutkowskiR. (2011). Oseltamivir in seasonal, pandemic, and avian influenza: a comprehensive review of 10-years clinical experience. Adv. Ther.28 (11), 927959. 10.1007/s12325-011-0072-7

  • 88

    SrujanaM. P.VigneshwaranE.KumarG. S.JyoshnaK.Jaffar SadiqM. J. (2017). Assessment of quality of life in children with epilepsy in rural settings of South India: a cross sectional study. CHRISMED J. Health Res.4 (2), 110116. 10.4103/2348-3334.201980

  • 89

    TianZ.FengB.WangX.-Q.TianJ. (2022). Focusing on cyclin-dependent kinases 5: a potential target for neurological disorders. Front. Mol. Neurosci.15, 1030639. 10.3389/fnmol.2022.1030639

  • 90

    TreanorJ. J.HaydenF. G.VroomanP. S.BarbarashR.BettisR.RiffD.et al (2000). Efficacy and safety of the oral neuraminidase inhibitor oseltamivir in treating acute InfluenzaA randomized controlled trial. JAMA283 (8), 10161024. 10.1001/jama.283.8.1016

  • 91

    VargheseA. F.ChockalingamM. P.ViswaksenanS.SanthappanJ. S.RavindranA. L.BalacB.et al (2024). Experimental determination of radius of gyration of an object using compound pendulum technique. AIP Conf. Proc.3059 (1), 030005. 10.1063/5.0193845

  • 92

    VesterS. K.BeavilR. L.LynhamS.BeavilA. J.Cunninghame GrahamD. S.McDonnellJ. M.et al (2021). Nucleolin acts as the receptor for C1QTNF4 and supports C1QTNF4-mediated innate immunity modulation. J. Biol. Chem.296, 100513. 10.1016/j.jbc.2021.100513

  • 93

    VolkowN. D.DingY.-S.FowlerJ. S.GatleyS. J. (2001). Imaging brain cholinergic activity with positron emission tomography: its role in the evaluation of cholinergic treatments in Alzheimer’s dementia. Biol. Psychiatry49 (3), 211220. 10.1016/s0006-3223(00)01112-4

  • 94

    Walczak-NowickaŁ. J.HerbetM. (2021). Acetylcholinesterase inhibitors in the treatment of neurodegenerative diseases and the role of acetylcholinesterase in their pathogenesis. Int. J. Mol. Sci.22 (17), 9290. 10.3390/ijms22179290

  • 95

    WangS.QaisarU.YinX.GrammasP. (2012). Gene expression profiling in Alzheimer's disease brain microvessels. J. Alzheimer’s Dis.31 (1), 193205. 10.3233/jad-2012-120454

  • 96

    WißfeldJ.Abou AssaleT.Cuevas-RiosG.LiaoH.NeumannH. (2024). Therapeutic potential to target sialylation and SIGLECs in neurodegenerative and psychiatric diseases. Front. Neurology15, 1330874. 10.3389/fneur.2024.1330874

  • 97

    XiaoM.YaoC.LiuF.XiangW.ZuoY.FengK.et al (2022). Sialic acid ameliorates cognitive deficits by reducing amyloid deposition, nerve fiber production, and neuronal apoptosis in a mice model of Alzheimer’s disease. NeuroSci3 (1), 2840. 10.3390/neurosci3010002

  • 98

    XieZ.BaileyA.KuleshovM. V.ClarkeD. J. B.EvangelistaJ. E.JenkinsS. L.et al (2021). Gene set knowledge discovery with enrichr. Curr. Protoc.1 (3), e90. 10.1002/cpz1.90

  • 99

    YesselmanJ. D.PriceD. J.KnightJ. L.Brooks IiiC. L. (2012). MATCH: an atom-typing toolset for molecular mechanics force fields. J. Comput. Chem.33 (2), 189202. 10.1002/jcc.21963

  • 100

    YuanS.ChanH. C. S.HuZ. (2017). Using PyMOL as a platform for computational drug design. WIREs Comput. Mol. Sci.7 (2), e1298. 10.1002/wcms.1298

  • 101

    ZhangR.YipV. L. Y.WithersS. G. (2010). “8.11 - mechanisms of enzymatic glycosyl transfer,” in Comprehensive natural products II. Editors LiuH.-W.ManderL. (Oxford: Elsevier), 385422.

  • 102

    ZhangW.XiaoD.MaoQ.XiaH. (2023). Role of neuroinflammation in neurodegeneration development. Signal Transduct. Target. Ther.8 (1), 267. 10.1038/s41392-023-01486-5

  • 103

    ZhaoK.-W.NeufeldE. F. (2000). Purification and characterization of recombinant human α-N-acetylglucosaminidase secreted by Chinese hamster ovary cells. Protein Expr. Purif.19 (1), 202211. 10.1006/prep.2000.1230

  • 104

    ZhuW.ZhouY.GuoL.FengS. (2024). Biological function of sialic acid and sialylation in human health and disease. Cell Death Discov.10 (1), 415. 10.1038/s41420-024-02180-3

  • 105

    ZoeteV.CuendetM. A.GrosdidierA.MichielinO. (2011). SwissParam: a fast force field generation tool for small organic molecules. J. Comput. Chem.32 (11), 23592368. 10.1002/jcc.21816

Summary

Keywords

Alzheimer’s disease, oseltamivir, neuraminidase, molecular docking, molecular dynamic simulation, gene expression analysis, gene enrichment analysis, insilico study

Citation

Alzarea SI, Alsaidan OA, Alhassan HH, Alzarea AI, Alsahli TG, Alharbi M, Afzal M and Sadiq Mantargi MJ (2025) Neuraminidase as a novel therapeutic management strategy for Alzheimer’s disease: evidenced through molecular docking, molecular dynamic simulation and gene expression analysis. Front. Chem. 13:1574702. doi: 10.3389/fchem.2025.1574702

Received

11 February 2025

Accepted

01 May 2025

Published

30 May 2025

Volume

13 - 2025

Edited by

Naohiko Yoshikai, Tohoku University, Japan

Reviewed by

Pranab Kishor Mohapatra, C. V. Raman Global University, India

Mamdouh F. A. Mohamed, Sohag University, Egypt

Updates

Copyright

*Correspondence: Muhammad Afzal,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics