Translate this page into:
Investigating the potential of novel thiazolidinone derivatives as anti-alzheimer agents: A comprehensive study through kinetics, in silico, and in vitro approaches
* Corresponding authors: E-mail addresses: arifbiochem@hotmail.com (M. Arif lodhi), petershow@ojlab.ac.cn (P. Shaw)
-
Received: ,
Accepted: ,
Abstract
Alzheimer’s disease (AD) is a progressive neurodegenerative disorder characterized by behavior, memory, and cognitive impairment. Acetylcholinesterase (AChE) is primarily responsible for acetylcholine (ACh) breakdown and is considered a leading cause of AD. The current study was designed to assess thiazolidinone compounds (1-12) as synthetic inhibitors of AChE as potential therapeutic options for AD. For this purpose, inhibitory activities, kinetics analysis, antispasmodic potential, and molecular dynamic (MD) simulations were conducted for 12 newly synthesized thiazolidinone compounds. The results demonstrate significant inhibition of AChE (IC50 ranging from 209.53 ± 1.01 µM to 1656.01 ± 1.60 µM) by all tested compounds. Kinetic analysis revealed that these compounds bind to the allosteric site and reduce AChE proficiency. The non-competitive nature of these compounds is an important finding that suggests that they could effectively reduce ACh hydrolysis even at high substrate concentrations. The obtained molecular docking results suggested that these compounds formed π-π stacking and hydrogen bonding interactions with the anionic sub-site and peripheral anionic site (PAS) of AChE. Their inhibitory efficacy is supported by this structural insight, which also helps in identifying important interactions for future optimization. MD simulation showed that the interaction of selected compounds 12 and 3 with AChE remains stable during simulations. This highlights their potential as effective therapeutic agents. Additionally, the whole series displayed calcium (Ca2+) antagonistic and antispasmodic potential, which could contribute to their therapeutic impacts beyond the inhibition of AChE. Among all tested compounds, compounds 12 and 3 demonstrated the highest proficiency in countering spontaneous and potassium (K+)-induced spasms. Therefore, the anti-AChE activity, Ca2⁺ antagonistic abilities, and safety profile of these compounds position them as potential multi-target drug candidates for AD and its related complications in the future.
Keywords
Acetylcholinesterase
Alzheimer’s therapy
Antispasmodic activity
Kinetics
Molecular dynamic simulation
Thiazolidinone

1. Introduction
Alzheimer’s disease (AD) is a neurodegenerative disease characterized by cognitive decline with memory loss [1]. Current estimates indicate that about 50 million people globally live with dementia, 30-35 million of them suffering from AD [2]. The cortical cholinergic neurotransmission dysfunction is closely related to AD; thus, it was hypothesized that dysfunction of neurons containing acetylcholine (Ach) in the brain contributes significantly to the cognitive decline in those individuals with old age AD. The cholinergic hypothesis is the new therapeutic focus for Alzheimer’s treatment. This hypothesis aims to inhibit acetylcholinesterase (AChE) (the enzymes that hydrolyze ACh) and increase ACh levels in the synaptic cleft [3]. Therefore, the most effective therapies for AD treatment are those that aim to boost the ACh levels, a vital neurotransmitter, the levels of which drop as the disease progresses [4]. The AChE inhibitors, such as galantamine, are efficient, as a result of which the breakdown of ACh is reduced. Still, galantamine’s clinical efficacy is hindered by its adverse effects, including nausea, bradycardia, and loss of appetite. To overcome this limitation, chemists focused on designing and developing novel leads with enhanced safety and therapeutic efficacy to treat AD [5]. In this regard, heterocycles, including thiadiazoles and thiazolidinone, provide a backbone for designing novel drugs. These scaffolds exhibit broad biological activities and can be used in rationally designing drugs with various modifications [6]. Modifying the thiazolidinone structure led to the inhibition of butyrylcholinesterase and AChE [7, 8]. This type of modification creates compounds that inhibit aldose reductase [9], glucosidase [8], urease, and α-amylase activities [10], and also show anticancer, anti-trypanosomal [11-13], COX-2 inhibitory [14], antimicrobial [15], antifungal [16, 17], and anti-inflammatory effects [18].
We designed this study to investigate a novel series of thiazolidinone derivatives (1-12) for their potential to inhibit AChE catalytic activity. We further extended our analysis to explore binding interactions and stability through chemical kinetics, docking, and molecular dynamics (MD) simulations.
2. Materials and Methods
2.1. General procedure for the synthesis of 4-thiazolidinone derivatives based on 2-mercaptobenzimidazol (2-MBI)
Synthesis of starting material hydrazone: 2-MBI derivatives were synthesized via multiple sequence reactions. Initially, 2-MBI was refluxed with bromoethane in EtOH using potassium hydroxide (KOH) as a base for 10 h. Following EtOH evaporation, 2-ethylthio-benzimidazole was produced as bright white needle-like crystals. In the second stage, 2-(2-(ethylthio)benzimidazolyl) acetate was produced in the semisolid form, which involved reacting 2-ethylthio-benzimidazole with ethyl chloroacetate in the presence of anhydrous potassium carbonate using dimethyl formamide (DMF) as a solvent under reflux for 15 h. After refluxing this intermediate with hydrazine hydrate in MeOH for 10 h, 2-((ethylthio)benzimidazolyl)acetohydrazide was formed. This precipitated when added to cold water and was filtered out. The last stage involved dissolving this hydrazide in methanol using a catalytic quantity of glacial acetic acid and then adding an aldehyde. The mixture was refluxed for 5 to 6 h, and when it cooled, the product was obtained as a precipitate. It was then filtered, cleaned, and dried in an open atmosphere.
Synthesis of final derivatives ( 1-12): For synthesis of thiazolidinone derivatives, 20 mL MeOH and 5 mL acetic acid were taken in a 100 mL round bottom flask, to which 0.733 mmol (0.09 g) of zinc chloride in molten form was added and stirred for about 5 min on a magnetic hot plate stirrer. Then 1.00 mmol (0.430 g) of the 2- mercaptobenzimidazole hydrazone derivative was added, and the mixture was stirred until completely dissolved on a hot plate. After 5 min, we added 1.00 mmol (0.14 mL) of mercapto-acetic acid (thioglycolic acid), and the reaction mixture was refluxed for about 8-10 h. We employed Thin Layer Chromatography to monitor the progress of the reaction. After completion of the reaction, the product was precipitated in freezing water, filtered, and washed with Na2CO3 to remove excess thioglycolic acid. The product was dried and preserved. The structures of novel derivatives were expanded upon using spectroscopic techniques like nuclear magnetic resonance (NMR), i.e., 1H NMR and 13C NMR [19].
2.2. General
The current study utilized analytical-grade chemical reagents. AChE (EC; 3.1.1.7) (C3389-2KU), 5-dithiobis 2-nitrobenzoic acid (DTNB) (D8130-10G), Galantamine (Cat. No. 69353-21-5), and ACh (A7000-25G) were used. In addition, Sigma Aldrich sourced ethanol (No. 64-17-5), monobasic sodium phosphate (1342-35-0), and Sodium phosphate dibasic (558-79-4). The thiazolidinone derivatives (1-12) were previously synthesized and described earlier [19].
2.3. Inhibitory evaluation of AChE
We carried out the AChE inhibition assay following the method described earlier [20], where galantamine served as a control. Initially, 140 µL of sodium phosphate buffer (100 mM) pH 8, 10 µL DTNB (0.25 mM), 0.5 mM of thiazolidinone compound dissolved in 20 µL, and 20 µL AChE enzyme were mixed and incubated for 15-20 min at room temperature. Finally, 0.4 mM of substrate ACh (10 µL) was added to the mixture and incubated for a further af 5 min. The reaction started, and ACh was hydrolyzed by AChE enzyme in the presence of DTNB; a yellow-colored 5-thio-2-nitrobenzoate anion was generated, which indicates the reaction completion and was read at 412 nm wavelength. We performed each experiment in biological triplicate and analyzed the data using SoftMax Pro v. 6.3, using the following Eq. (1) to compute the inhibition as a percentage:
2.4. Kinetic study
The amount of compound required for 50% inhibition is known as IC50, which was calculated using EZ-Fit (Perrella Scientific Inc; Amherst, USA) at various concentrations of the test compounds.
In the enzyme-substrate (ES) reaction, the AChE enzyme reacts with substrate ACh to form AChE-ACh complex, forming the product (P). The dissociation constant values were calculated using Lineweaver-Burk, Dixon plots, and their secondary replots [21]. The values of Km, ki, and Vmax were calculated from Lineweaver-Burk and Dixon plots through a non-linear regression equation. The Lineweaver-Burk plot was employed to determine the Ki values; initially, values of 1/V maxapp were found at each line intersection point of each concentration of the test compound on the y-axis. The slope of each line generated by the compound was plotted against various test compound concentrations on the Lineweaver-Burk plot.
2.5. Analysis of statistics
GraFit software [22] was used to evaluate correlation coefficients, slope, intercept, and standard error values.
2.6. AChE and ligand docking
The Molecular Operating Environment (MOE) package (www.chemcomp.com) was used to explore various binding conformations of AChE with different compounds. The 3D structure of Homo sapiens AChE (hAChE) (PDB code: 4EY6, resolution: 2.40 Å) was obtained from PDB (Protein Data Bank) (https://www.rcsb.org).
The online server Play-molecule was used to prepare the target protein, involving chain completion, protonation, and correction of bonds [23]. Crystalized H2O molecules were removed from the retrieved crystal structure using MOE. Subsequently, 3D protonation was applied, and energy minimization was performed using the default parameters of the MOE program. The structures of synthesized thiazolidinone derivatives were built in the MOE program and protonated, and then energy was minimized through MMFP94x force fields with a 0.05 gradient. Each tested compound was docked in the active sites of the target AChE enzyme using the default parameters of the MOE-Dock program. Ten conformations were generated for each compound, and the top-ranked conformations based on docking scores were selected for further in silico studies.
To validate the reliability of the docking procedure, a re-docking experiment was performed using the AChE co-crystalized ligand as a reference (Figure S1). This step ensures the accuracy and consistency of the molecular docking method employed in this study.
2.7. Protein-ligand interaction fingerprint (PLIF)
The accuracy of the molecular docking data was confirmed using PLIF. The PLIF tool was linked with MOE and applied to display the principal interaction efficiency [24]. A docking protocol was employed to create a database. Bonding interactions, including hydrogen bonding, ionic bonding, and surface contacts, were analyzed per residue. The PLIF histogram was generated, which illustrates a database of the complex. The MOE PLIF was utilized to analyze all the molecular docking postures, and the software was executed using the default parameters. All the molecular docking poses were submitted to the MOE PLIF, and the program was run using default parameters.
2.8. MD simulations
A total of two top-ranked compounds (12 and 3) were short-listed based on IC50 and docking scores for further MD simulations (GROMACS 5.1.2) [25]. Initially, a protein ligand complex was generated using the pdb 2 gmx module. Parametrization and topology of protein were generated by the OPLS-AA/L force field [26], and ligand topology was prepared using the Swissparam online web server (http://www.swissparam.ch/). Subsequently, each protein-ligand complex was solvated in a cubic box (10 × 10 × 10) surrounded by approximately 33,439 TiP3P water molecules. In the next step, 9 sodium ions were added to neutralize the complex. Using positional restraint and the steepest descent algorithm, energy minimization was carried out up to a maximum force of 10 KJ/mol [27], eliminating the atomic location mistake and structural variations such as bond length, bond angle, and structural conflicts between the ions, water molecules, and protein complex. The solvent around the protein was equilibrated under the number of atoms, volume, and temperature (NVT) ensemble, and then using the NPT ensemble. NVT and NPT would both perform 50,000 steps for 0.1 ns at 300 K and 1 bar of pressure. Each complex was then put through MD simulations utilizing a leapfrog algorithm for 10,000,000 ps at 300,000 atmospheric NPT ensemble pressure and periodic boundary conditions for 0.002 ps. All hydrogen bonds were contained throughout the procedure using the LINCS algorithm [28]. However, the particle MeshEwald module was functionalized with a 0.16 Å Fourier grid spacing [21].
2.9. Antispasmodic and Ca2⁺ antagonistic activities
The antispasmodic potentials of thiazolidinone compounds (1-12) were evaluated by the examination of spontaneous contracting rabbit jejunum [20]. The local breed rabbits of both sexes (1.5-2.0 kg) used in the current study were generously supplied by Agha Khan University, Pakistan. Animals were allowed free access to water, but their food was stopped 24 h before the experiment started. The animals were sacrificed with a blow to the back of the head, the abdomen was cut open, and a piece of their jejunum was removed. The 2 cm long segments were suspended in Tyrode’s solution (KCl 80 mM, MgCl2 1.05 mM, NaHCO3 11.87mM, NaCl 91.04 mM, NaH2PO4 0.41mM, CaCl2 1.8 mM, and glucose 5.55 mM), which was aerated with a mixture of 95% O2 and 5% CO2, and preserved at 37°C. The spontaneous intestinal movements were observed and recorded isotonically using Harvard transducers and the Harvard Student Oscillograph. Each tissue was permitted to equilibrate for at least 30 min before any drug addition. The jejunum displayed spontaneous rhythmic contractions under these experimental conditions, enabling the study of spasmolytic activity directly without the use of an agonist. In addition, a Ca2⁺ antagonistic activity was conducted to authenticate the capability of the thiazolidinone compound (1-12) to relax the contractions induced by 80 mM K+. [29]
2.10. Cytotoxicity activity
Human neutrophil isolation: Human neutrophils were used to assess the cytotoxic effects of synthesized thiazolidine derivatives. For this purpose, whole blood (heparinized) was obtained from a local blood bank, and the required neutrophils were isolated [30]. The whole blood was then mixed with Ficoll-Paque, and after sedimentation, the unwanted RBCs were layered with the buffy coat on the Ficoll cushion (3 mL). The mixture was centrifuged for 30 min at 1500 rmp, and the supernatant was disposed of. The preparation was then centrifuged, and the neutrophils were cleaned with modified Hanks solution (MHS) and resuspended at a concentration of 1× 107 cells/mL.
Assay procedure: The isolated neutrophils were incubated with thiazolidinone compounds for 30 min, followed by the addition of 0.25 mM WST-1 in a water bath shaker at 37°C [20]. The absorbance change was then measured at 450 nm in a 96-well plate using a Spectra Max, CA, United States. The OD (optical density) is the mean of five experimental replicates.
Percent viability was determined by Eq. (2).
3. Result and Discussion
3.1. Biology
Assessment of thiazolidinone inhibitory efficacy was performed against AChE. All compounds exhibited varying degrees of potential at various concentrations, with IC50 values ranging from 209.53 ± 1.0 µM to 1656.01 ± 1.60 µM, as illustrated in Table 1. The current results suggest that these compounds possess significant potential as inhibitors of AChE, laying the foundation for further investigation. Compound 12 was the most potent inhibitor, with the lowest IC50 (209.53 ± 1.01 µM) among the tested series, but had lower efficacy as compared to standard Galantamine (IC50 = 92.49 ± 0.69 µM). The effectiveness of this compound may be due to the presence of an octyl group in its structure. Compound 3 (IC50 = 249.52 ± 1.25 µM) with a hydroxyl moiety on the para position was the second most potent compound in this category. Compound 11 was the third most potent compound with an IC50 of 345.15 ± 1.47 µM; the effectiveness of this compound is likely attributed to the existence of the 2-methyl furan ring. Compound 4 (IC50 = 372.09 ± 1.04 µM) was ranked fourth in potency against AChE; the efficacy of this compound may be due to the high electronegative -Cl on ortho and para positions of a benzene ring. Compound 8 (IC50 = 528.10 ± 1.45 µM) contains a hydroxyl group at the ortho and an ethoxy group at the meta positions. The reduced activity, compared to compound 3, suggests that the meta-positioned ethoxy group may hinder the inhibitory effect.
| Compound no. | R | IC50 ± SEM (µM) | Docking Score |
|---|---|---|---|
| 1 |
|
1172.64 ± 0.49 | -9.529 |
| 2 |
|
1656.01 ± 1.60 | -9.129 |
| 3 |
|
249.52 ± 1.25 | -12.380 |
| 4 |
|
372.09 ± 1.04 | -11.371 |
| 5 |
|
808.40 ± 0.5 | -10.485 |
| 6 |
|
697.19 ± 0.78 | -10.987 |
| 7 |
|
65.86 ± 0.91 | -10.794 |
| 8 |
|
528.10 ± 1.45 | -11.123 |
| 9 |
|
1227.97 ± 0.86 | -10.321 |
| 10 |
|
801.74 ± 0.89 | -10.699 |
| 11 |
|
345.15 ± 1.47 | -11.842 |
| 12 |
|
209.53 ± 1.01 | -13.511 |
| Galantamine |
|
92.49 ± 0.69 | -14.542 |
Similarly, compared with compound 8, compound 6, which contains the hydroxyl group at the ortho position, exhibited a significant loss in efficacy (IC50 = 697.19 ± 0.78 µM). This was observed when the ethoxy group was replaced with a methoxy group at the meta position of the compound. Compound 7 bearing an electron-withdrawing chloro group at the ortho position was the 7th active compound of the series with IC50 = 765.86 ± 0.91 µM. Compounds 10 and 5 displayed slightly different IC50 values, IC50 = 801.74 ± 0.89 µM, and IC50 = 808.40 ± 0.5 µM. Compound 10, which has hydroxyl groups at the ortho and para positions, demonstrated greater efficacy than compound 5, substituted with an electron-withdrawing nitro group. Compounds 1, 9, and 2 were less potent with IC50 values of 1172.64 ± 0.49, 1227.97 ± 0.86, and 1656.01 ± 1.60, respectively, as shown in Table 1. These results indicate that the activities of the compound may be associated with the electron-withdrawing hydroxyl and chloro groups or the electron-donating alkyl group.
3.2. Kinetics of AChE by thiazolidinone derivatives
Kinetic studies were employed to study the inhibitory efficacy of thiazolinone derivatives (1-12) on AChE. As a result, the tested compounds inhibited our enzyme with Ki values ranging from 41.74 ± 0.8 to 578.21 ± 0.4 µM. Notably, Compound 12 with Ki = 41.74 ± 0.8 µM established strong binding, whereas Compound 3 demonstrated weak interaction (Ki = 578.21 ± 0.4 µM) with AChE. Consequently, the IC50 value for compound 12 (209.53 ± 1.01 µM) is lower and confirms its inhibitory potential compared to compound 3, which demonstrated a weak inhibitory efficacy.
According to kinetic parameters, as shown in Table 2, all compounds of Thiazolidinones acted as noncompetitive inhibitors of AChE. The Km values were consistent, and Vmax decreased, which indicated a noncompetitive inhibition mechanism (Table 2). The steady-state inhibitory data of compound 12 against AChE have been illustrated in Figure 1.
| Compound | Ki (µM) ± SEM | Km (mM) | Kmapp (mM) | Vmax (µmol/min) | Vmaxapp | Inhibition type |
|---|---|---|---|---|---|---|
| 1 | 398.16 ± 0.5 | 0.12 | 0.12 | 5.0 | 2.0 | Noncompetitive |
| 2 | 578.21 ± 0.4 | 0.12 | 0.12 | 5.0 | 2.1 | Noncompetitive |
| 3 | 55.06 ± 1.4 | 0.12 | 0.12 | 5.0 | 2.2 | Noncompetitive |
| 4 | 104.63 ± 1.9 | 0.12 | 0.12 | 5.0 | 2.1 | Noncompetitive |
| 5 | 238.03 ± 0.5 | 0.12 | 5.0 | 5.0 | 2.3 | Noncompetitive |
| 6 | 229.03 ± 1.2 | 0.12 | 0.12 | 5.0 | 2.2 | Noncompetitive |
| 7 | 232.14 ± 0.6 | 0.12 | 0.12 | 5.0 | 2.3 | Noncompetitive |
| 8 | 234.26 ± 0.6 | 0.12 | 0.12 | 5.0 | 2.3 | Noncompetitive |
| 9 | 301.20 ± 1.5 | 0.12 | 0.12 | 5.0 | 2.1 | Noncompetitive |
| 10 | 235.96 ± 0.8 | 0.12 | 0.12 | 5.0 | 2.3 | Noncompetitive |
| 11 | 75.33 ± 0.3 | 0.12 | 0.12 | 5.0 | 2.3 | Noncompetitive |
| 12 | 41.74 ± 0.8 | 0.12 | 0.12 | 5.0 | 2.4 | Noncompetitive |

- Compound 12’s anti-AChE activity is represented in three plots. (a) Is the first plot of, Lineweaver-Burk, which shows the relationship between initial velocities and four constants conc of ACh: 100 µM (□), 150 µM (●), 175 µM (○), in the lack of compound 12 (■). (b) Is the second Dixon plot, which demonstrates the reciprocal relationship between initial velocities and variable concentration of compound 12 at constant conc. of ACh 100 µM (■), 150 µM (□), 175 µM (●) and 200 µM (○). (c) Is the third plot, which displays the correlation between various concentrations of compound 12 and the slopes.
3.3. PLIF analysis
In the present study, all compounds of thiazolidinone were docked in the pocket of 4EY6 to generate barcodes and population plots Figures 2(a) and (b). Among all residues, TYR124 was found to be the most active residue in binding interactions with the target protein (4EY6) active site and followed by some other important interacting residues revealed by fingerprint bits such as TRP86, TYR341, SER125, GLY121, TYR337, TYR72, SER203, ASN74, ASN87, PHE295, HIS447, and TYR449.

- PLIF between thiazolidinone compounds and 4EY6 (a) barcode (b) population.
3.4. Molecular docking analysis
Molecular docking was conducted to validate in vitro results. The results revealed that all the compounds exhibited good binding affinity towards the 4EY6 active site. The top three compounds (12, 3, and 11) among the series with higher docking scores were evaluated for their interactions with the active site of AChE.
Compound 12 was the most potent compound of the series with IC50 = 209.53± + 1.01 µM and docking score = - 13.511. The carbonyl oxygen of this compound acts as an acceptor and forms two hydrogen bonds, with TYR124 and SER125 having bond lengths of 3.1 and 2.8 Å, respectively. Moreover, the amido group of the same compound mediated hydrogen bond interaction with the side chain TYR337 with a bond length of 1.7 Å, as illustrated in Figures 3(a) and (b). Nevertheless, this compound revealed lower effectiveness than Galantamine (IC50 = 92.49 ± 0.69 µM, docking score = ₋ 14.542). Galantamine was found to mediate two interactions with side chains TYR337 and TYR341, with 3.1 Å and 2.9 Å, respectively. The galantamine tetrahydrofuran moiety acted as an acceptor, making two other interactions with TYR124 and SER125 residues of the AChE (Figures 4a and b).

- AChE active site interactions with compound 12 (a, b), Green sticks: amino acid residues; Yellow sticks: ligand; Dashed lines: hydrogen bonds. Red, white, blue, and green colors show oxygen, hydrogen, nitrogen, and carbon moieties respectively.

- AChE active site interactions with galantamine (a, b), Green sticks: amino acid residues; Yellow sticks: ligand; Dashed lines: hydrogen bonds, Red, white, blue and green colors show oxygen, hydrogen, nitrogen and carbon moieties respectively.
Compound 3 was the second most active compound against AChE, having an IC50 = 249.52 ± 1.25 µM and a docking score of ₋12.380. This compound established a strong hydrogen bond with TYR124, having a bond length of 2.0 Å. In addition, the carbonyl oxygen of compound 3 forms a secondary hydrogen bond with SER125 (3.5 Å) and three π-π stacking interactions with TYR337 and TRP86 residues of the anionic subsite (Figures 5a and b). However, compound 3 exhibited lower potency than compound 12 and galantamine.

- AChE active site interactions with compound 3 (a, b), Green sticks: amino acid residues; Yellow sticks: ligand; Dashed lines: hydrogen bonds, Red, white, blue, and green colors show oxygen, hydrogen, nitrogen, and carbon moieties respectively.
Compound 11 (IC50 value of 345.15 ± 1.47 µM, docking score of -11.842) emerged as the third most active compound among the series of thiazolidinone. Like compounds 12 and 3, the carbonyl oxygen of compound 11 established a hydrogen interaction with SER125 (3.5 Å). Additionally, it established a second strong hydrogen interaction with the TYR124 residue of the anionic subsite, with 2.8 Å (Figures 6a and b). Despite these interactions, compound 11 established lower potency than galantamine, compounds 12 and 3.

- AChE active site interactions with compound 11 (a, b), Green sticks: amino acid residues; Yellow sticks: ligand; Dashed lines: hydrogen bonds, Red, white, blue and green colors show oxygen, hydrogen, nitrogen and carbon moieties, respectively.
As a result, our compounds, despite interacting with all four subsites of the AChE active site, precisely interacted with the peripheral anionic site (PAS), anionic subsites, and allosteric subsites by hydrogen bonding, π-π bonding, and hydrophobic interactions. Therefore, docking analysis and kinetic studies are aligned, indicating a non-competitive behavior of our compounds.
Top-ranked compounds (12 and 3) based on high docking scores and lower IC50 values, as well as control (galantamine), were chosen for further 100 ns simulation. The assessment of protein-ligand complex stability was calculated by means of root mean square deviation (RMSD). The RMSD curve of compound 12 was initially unstable and was followed by stabilization from 50-100 ns, shown by values from 0.20 to 0.3 nm. Similarly, the compound 3 complexes demonstrated a comparable trend, with instability up to 50 ns and sustaining stability till 100 ns. Our findings show that the stability and reliability of the compound 12-AChE complex were slightly better than those of the compound 3-AChE complex. The galantamine maintained the lowest value of RMSD, around 0.20-0.25nm, indicating the highest stability with minimal conformational drift when compared to the other two compounds, 12 and 3 (Figure 7(a)).

- (a) RMSD, (b) RMSF (c) Rg, and (d) Hydrogen bonding plots of compound 12 (Black), compound 3 (Red), and Galantamine (Green) 100 ns in MD simulations.
The root mean square fluctuation (RSMF) analysis was performed to calculate the flexibility of each residue in the AChE when complex with compounds 12 and 3 during MD simulations. The current investigation revealed that both complexes exhibited higher fluctuations in the highly flexible regions of the AChE, specifically in the N-terminal, C-terminal, and loop regions. In addition, residues contributing to α, β sheets and active sites establish decreased RMSF, as represented in Figure 7(b). These results are in line with the inherent mobility of these regions, which is known to support the structural flexibility of the enzyme during ligand binding and catalysis [31].
The radius of gyration (Rg) predicts the conformational changes and compactness values of the chosen protein-ligand complex. Rg values for the compound 12-AChE complex were observed in the range of 2.28 nm to 2.32 nm, which displayed the compactness of our protein and remained unchanged upon compound 12 binding. The compound 3-AChE complex exhibited identical values, which may show natural movement or rotations contained by the molecular structure. The galantamine curve at the end showed some conformational change, but overall, the Rg remains stable, suggesting that the structure retains its compactness without major unfolding (Figure 7(c)). Rg results indicate that the overall protein structure remains stable and compact during the MD simulation process, with minor fluctuations in the dimension relating to dynamic flexibility.
We further calculated hydrogen bonds formed between our compounds and AChE; compound 12 established a maximum of 3 H-bonds, and compound 3 formed 2 H-bonds with AChE during 100 ns MD simulations; galantamine exhibited 4 H-bonds, overall the results indicate that compound 12 formed a more stable complex with AChE than compound 3 but less stable complex as compared to galantamine as represented in Figure 7(d). H-bonds play a significant role in the stabilization of protein-ligand complexes by providing specific interactions that improve the binding affinity and structural integrity. The larger number of H-bonds in the compound 12-AChE complex shows stronger and more sustained interactions between the active residues of AChE and the compound. This may help explain its enhanced binding stability and inhibitory effectiveness. On the other hand, the compound 3-AChE complex smaller number of H-bonds, indicating comparatively weaker interactions, which could lead to reduced stability and efficiency.
3.5. Antispasmodic and Ca2⁺ antagonistic activities
Antispasmodic drugs are frequently used to treat spasms in the smooth muscles of the gastrointestinal tract [32]. Among these, the most commonly used drugs are those with anticholinergic properties, which act by blocking the activity of ACh, a neurotransmitter involved in muscle spasms. This mechanism efficiently reduces muscle contraction and eases the discomfort they cause [33]. Therefore, to analyze the anticholinergic characteristics of thiazolidinone derivatives, we also studied the potential of these compounds for putative muscle anti-contractility effect using the animal jejunum. This approach was adopted because the jejunum has a natural contraction ability, exclusive of utilizing any agonist. This helps explicitly investigate new drugs with smooth muscle relaxant activity [34]. Nevertheless, stimulants of smooth muscles can also be used, and their effects are compared with standard spasmogenic agents responsible for spasmogenicity. Current study results of the thiazolidinone series also exhibited suppression of spontaneous jejunal contractility (Table S1). Compounds, specifically 12 and 3 showed good spasmolytic activity with an ED50 of 0.13±0.5 mg/mL (Mean ± SEM; n=3) among all compounds (Figures 8a and b).

- Dose-response graphs illustrating the antispasmodic impact of (a) compound 12 and (b) 3 on both spontaneous/natural and K+ induced contractions in jejunal tissues (Mean ± SEM, n=3).
The spasm of the rabbit jejunum and other smooth muscles is an indication of the increased level of Ca2⁺in the cytoplasm [35], which plays a key role in the interactions of actin and myosin leading to muscle contraction [36]. The initial elevated level of Ca2⁺is the direct consequence of Ca2⁺influx via voltage-gated channels (VDCs) into a cell or its release from the intracellular stores of Ca+2. This depolarization and repolarization, consistent with the influx of Ca2⁺by VDCs, bring forth the intestine’s spontaneous contractile behavior [37]. In the present study, all Thiazolidinone derivatives were able to relax this contractile behavior of rabbit jejunal smooth muscles, probably by the inhibition of Ca2⁺entry via VDCs or by affecting the Ca2⁺release from intracellular Ca2⁺stores (Table S1). In conclusion, among all tested compounds, compounds 12 and 3 exhibited the most potent activities with median ED50 values of 0.25±0.03 µg/mL and 0.26 ± 0.03 µg/mL (Mean ± SEM; n = 3) on K+ (80 mM) (Figures 8a and b).
3.6. Cytotoxicity evaluation of thiazolidinone (1-12)
Thiazolidinone compounds were subjected to cytotoxicity using a human neutrophil viability assay. All compounds of thiazolidinone (1-12) were nontoxic, like standard galantamine (Table S2). Among the whole series, compounds 12 and 3 showed the highest (94.07 ± 1.0 and 92.07 ± 1.4) percent cell viability at conc 200 µg/mL, as displayed in Figure 9.

- A bar chart shows the percentage viability of human neutrophils upon exposure to the thiazolidinone derivatives (1-12). The error bars represent the standard deviation resulting from the triplicate experiments.

- Synthesis of 4-thiazolidinone derivatives. R= Aromatic or Aliphatic aldehyde.
4. Conclusions
In conclusion, this in vitro AChE inhibition and kinetic study demonstrated that synthesized thiazolidinone derivatives (1-12) inhibit AChE non-competitively. The activities of enzymes were explored experimentally (kinetic study) and theoretically (molecular docking and MD simulations) based on molecular interactions between inhibitors and enzymes. However, in vitro and in silico analysis revealed that compounds 12 and 3 showed significant inhibitory capability for AChE activity compared to other series compounds evaluated in this study. At the same time, their efficiency was lower than that of galantamine. The unique structural properties of compounds 12 and 3, such as the octyl group and hydroxyl group (para position), respectively, probably contributed to their inhibitory activity. These characteristics further demonstrate its potential as a promising candidate. Compounds 12 and 3 also demonstrated good spasmolytic, Ca2⁺ antagonistic, and nontoxic effects on neutrophil cells. These findings suggest these compounds are promising drug candidates for AChE inhibition in AD treatment. However, validation through animal model-based studies is required to enhance these findings and pave the way for the advancement of new drugs.
Acknowledgment
The authors thank the Deanship of Scientific Research at King Khalid University for supporting this work through large Groups (Project under grant number GRP/13/45). The study is supported by the Science Foundation of Oujiang Laboratory start up fund (No. Oujiang lab OJQD2022005).
CRediT authorship contribution statement
Farida Begum and Muhammad Arif Lodhi: Designing the study, Farida Begum and Sajid Iqbal: Experimentation, Writing original draft, Formal analysis, Mohammad Y. Alfaifi, Serag Eldin I. Elbehairi, Ali Abdullah Shati and Peter Shaw: Funding acquisition, Methodology, Muhammad Arif Lodhi and Peter Shaw: Supervision, Benedictor Alexander Nguchu, Sajid Iqbal and Haroon Iqbal: Review, editing and visualization. All authors reviewed the manuscript.
Declaration of competing interest
All authors declared no competing intrest.
Declaration of generative AI and AI-assisted technologies in the writing process
The authors confirm that there was no use of artificial intelligence (AI)-assisted technology for assisting in the writing or editing of the manuscript and no images were manipulated using AI.
Supplementary data
Supplementary material to this article can be found online at https://dx.doi.org/10.25259/AJC_95_2025.
References
- Phytomolecules as alzheimer’s therapeutics: A comprehensive review. European Journal of Medicinal Chemistry. 2025;288:117401. https://doi.org/10.1016/j.ejmech.2025.117401
- [Google Scholar]
- The neurobiology and therapeutic potential of multi-targeting β-secretase, glycogen synthase kinase 3β and acetylcholinesterase in Alzheimer’s disease. Ageing Research Reviews. 2023;90:102033. https://doi.org/10.1016/j.arr.2023.102033
- [Google Scholar]
- Innovative approaches in acetylcholinesterase inhibition: a pathway to effective Alzheimer’s disease treatment. Molecular Diversity 2025 https://doi.org/10.1007/s11030-025-11170-1
- [Google Scholar]
- Recent advances in drug development for alzheimer’s disease: A comprehensive review. International Journal of Molecular Sciences. 2025;26:3905. https://doi.org/10.3390/ijms26083905
- [Google Scholar]
- Synthesis and biological evaluation of benzothiazin-4-ones: A possible new class of acetylcholinesterase inhibitors. Journal of Enzyme Inhibition and Medicinal Chemistry. 2019;34:197-203. https://doi.org/10.1080/14756366.2018.1543286
- [Google Scholar]
- Drug-1,3,4-Thiadiazole conjugates as novel mixed-type inhibitors of acetylcholinesterase: Synthesis, molecular docking, pharmacokinetics, and ADMET evaluation. Molecules (Basel, Switzerland). 2019;24:860. https://doi.org/10.3390/molecules24050860
- [Google Scholar]
- One-pot four-component synthesis of thiazolidin-2-imines using CuI/ZnII dual catalysis: A new class of acetylcholinesterase inhibitors. Bioorganic Chemistry. 2019;84:518-528. https://doi.org/10.1016/j.bioorg.2018.12.002
- [Google Scholar]
- Synthesis, characterization and molecular modeling of amino derived thiazolidinones as esterase and glucosidase inhibitors. Journal of Molecular Structure. 2020;1219:128609. https://doi.org/10.1016/j.molstruc.2020.128609
- [Google Scholar]
- In vitro evaluation of 5-arylidene-2-thioxo-4-thiazolidinones active as aldose reductase inhibitors. Bioorganic & Medicinal Chemistry Letters. 2011;21:200-203. https://doi.org/10.1016/j.bmcl.2010.11.041
- [Google Scholar]
- Synthesis of new arylhydrazide bearing Schiff bases/thiazolidinone: α-Amylase, urease activities and their molecular docking studies. Bioorganic Chemistry. 2019;91:103112. https://doi.org/10.1016/j.bioorg.2019.103112
- [Google Scholar]
- Biodistribution and anticancer characteristics of Les-3833, A Novel 4-thiazolidinone-based lead compound. Scientia Pharmaceutica. 2020;88:18. https://doi.org/10.3390/scipharm88020018
- [Google Scholar]
- Design and synthesis of novel 4-thiazolidinone derivatives with promising anti-breast cancer activity: Synthesis, characterization, in vitro and in vivo results. Bioorganic Chemistry. 2020;104:104276. https://doi.org/10.1016/j.bioorg.2020.104276
- [Google Scholar]
- Synthesis of 5-enamine-4-thiazolidinone derivatives with trypanocidal and anticancer activity. Bioorganic Chemistry. 2019;86:126-136. https://doi.org/10.1016/j.bioorg.2019.01.045
- [Google Scholar]
- Molecular modeling, synthesis and screening of some new 4-thiazolidinone derivatives with promising selective COX-2 inhibitory activity. European Journal of Medicinal Chemistry. 2012;57:59-64. https://doi.org/10.1016/j.ejmech.2012.08.046
- [Google Scholar]
- Novel 4-thiazolidinone derivatives as anti-infective agents: Synthesis, characterization, and antimicrobial evaluation. Biochemistry Research International. 2016;2016:8086762. https://doi.org/10.1155/2016/8086762
- [Google Scholar]
- Development of thiazolidinones as fungal carbonic anhydrase inhibitors. International Journal of Molecular Sciences. 2020;21:2960. https://doi.org/10.3390/ijms21082960
- [Google Scholar]
- Design, Synthesis, and biological evaluation of novel thiazolidinone-containing quinoxaline-1,4-di- N-oxides as antimycobacterial and antifungal agents. Frontiers in Chemistry. 2020;8:598. https://doi.org/10.3389/fchem.2020.00598
- [Google Scholar]
- Thiazoles and thiazolidinones as COX/LOX inhibitors. Molecules (Basel, Switzerland). 2018;23:685. https://doi.org/10.3390/molecules23030685
- [Google Scholar]
- Mercaptobenzimidazole based hydrazone derivatives as potential antioxidant and alpha-glucosidase inhibitors. Current Bioactive compounds. 2022;18:77-83. http://doi.org/10.2174/1573407218666220419105140
- [Google Scholar]
- Antibiotic resistance: A current epilogue. Biochemical Pharmacology. 2017;134:139-146. https://doi.org/10.1016/j.bcp.2016.12.005
- [Google Scholar]
- Inhibition of Acetylcholinesterase with novel 1, 3, 4, oxadiazole derivatives: A kinetic, in silico, and in vitro approach. ACS Omega. 2023;8:46816-46829. https://doi.org/10.1021/acsomega.3c06298
- [Google Scholar]
- GraFit Version 4.09, Erithacus Software Ltd, Stains K. St. Louis MO: SYBYL Molecular Modeling Software, Tripos Associated Ltd; 1999.
- playmolecule proteinprepare: A web application for protein preparation for molecular dynamics simulations. Journal of Chemical Information and Modeling. 2017;57:1511-1516. https://doi.org/10.1021/acs.jcim.7b00190
- [Google Scholar]
- Predicting potential antitumor targets of aconitum alkaloids by molecular docking and protein–ligand interaction fingerprint. Medicinal Chemistry Research. 2016;25:1115-1124. https://doi.org/10.1007/s00044-016-1553-7
- [Google Scholar]
- GROMACS: Fast, flexible, and free. Journal of Computational Chemistry. 2005;26:1701-1718. https://doi.org/10.1002/jcc.20291
- [Google Scholar]
- GROMACS: A message-passing parallel molecular dynamics implementation. Computer Physics Communications. 1995;91:43-56. https://doi.org/10.1016/0010-4655(95)00042-e
- [Google Scholar]
- In-Silico drug discovery approach targeting receptor tyrosine kinase-like orphan receptor 1 for cancer treatment. Scientific Reports. 2017;7:1029. https://doi.org/10.1038/s41598-017-01254-w
- [Google Scholar]
- Molecular docking and molecular dynamics (MD) simulation of human anti-complement factor H (CFH) antibody Ab42 and CFH polypeptide. International Journal of Molecular Sciences. 2019;20:2568. https://doi.org/10.3390/ijms20102568
- [Google Scholar]
- Antispasmodic activity of the ethanol extract of Citrullus lanatus seeds: Justifying ethnomedicinal use in Pakistan to treat asthma and diarrhea. Journal of Ethnopharmacology. 2022;295:115314. https://doi.org/10.1016/j.jep.2022.115314
- [Google Scholar]
- Juliflorine: a potent natural peripheral anionic-site-binding inhibitor of acetylcholinesterase with calcium-channel blocking potential, a leading candidate for Alzheimer’s disease therapy. Biochemical and Biophysical Research Communications. 2005;332:1171-1179. http://doi.org/10.1016/j.bbrc.2005.05.068
- [Google Scholar]
- Biological implications of the intrinsic deformability of human acetylcholinesterase induced by diverse compounds: A computational study. Biology. 2024;13:1065. https://doi.org/10.3390/biology13121065
- [Google Scholar]
- Antispasmodic potential of medicinal plants: A comprehensive review. Oxidative Medicine and Cellular Longevity. 2021;2021:4889719. https://doi.org/10.1155/2021/4889719
- [Google Scholar]
- Mechanism of action of cholinergic drugs. In: How Synthetic Drugs Work. Elsevier; p. :27-46. https://doi.org/10.1016/B978-0-323-99855-0.00002-6
- [Google Scholar]
- Sarcorucinine-D inhibits cholinesterases and calcium channels: Molecular dynamics simulation and in vitro mechanistic investigations. Molecules (Basel, Switzerland). 2022;27:3361. https://doi.org/10.3390/molecules27113361
- [Google Scholar]
- Calcium signaling in smooth muscle. Cold Spring Harbor Perspectives in Biology. 2011;3:a004549. https://doi.org/10.1101/cshperspect.a004549
- [Google Scholar]
- Signaling in muscle contraction. Cold Spring Harbor Perspectives in Biology. 2015;7:a006023. https://doi.org/10.1101/cshperspect.a006023
- [Google Scholar]
- How do drugs initiate contraction in smooth muscles? Trends in Pharmacological Sciences. 1981;2:261-265. https://doi.org/10.1016/0165-6147(81)90334-5
- [Google Scholar]
