Abstract
Galantamine, an acetylcholinesterase (AChE) inhibitor used for symptomatic treatment of Alzheimer's disease (AD), shows substantial inter-individual variability in clinical response. Missense single nucleotide polymorphisms (SNPs) within the AChE active-site gorge may modulate inhibitor recognition. In this computational study, binding residues were defined from human AChE inhibitor co-crystal structures and cross-referenced with dbSNP missense variation, followed by in-silico predictions of variant impact, evolutionary conservation and folding stability, and assessment of ligand engagement by docking and molecular dynamics (MD) with MM/GBSA binding-energy estimation. Using complexes containing galantamine (GNT) and a donepezil-like ligand (E20), 11 of 807 AChE missense variants overlapped binding-site residues, highlighting Phe294 (UniProt Phe326) and His447 (UniProt His479). ConSurf classified His447 as highly conserved, and MUpro predicted decreased folding stability for His447 substitutions. SwissDock docking indicated that His447Gln retains a plausible GNT binding pose and yielded the least favourable docking score among the tested variants, consistent with a potential reduction in binding strength. MD simulations (200 ns) of wild-type and His447Gln AChE-GNT complexes supported preserved global structural integrity of the complex over the simulated timescale, while indicating local remodelling of the GNT binding microenvironment. MM/GBSA estimates from terminal snapshots suggested a modestly less favourable theoretical binding free energy for His447Gln relative to wild-type (approximately 2.0 kcal mol⁻¹). Given that His447 is the catalytic triad histidine, such substitutions may have consequences for catalysis in addition to inhibitor binding; these in-silico findings require experimental validation using site-directed mutagenesis with kinetic and binding assays.