Huperzine A (HUP) plays a crucial role in Alzheimer’s therapy by enhancing cognitive function through increased cholinergic activity as a reversible acetylcholinesterase (AChE) inhibitor. Despite some limitations being seen in AChE inhibitors, ongoing research remains dedicated to finding innovative and more effective treatments for Alzheimer’s disease. To achieve the goal of the discovery of potential HUP analogues with improved physicochemical properties, less toxic properties, and high biological activity, many in silico methods were applied. Based on the acetylcholinesterase–ligand complex, an e-pharmacophore model was developed. Subsequently, a virtual screening involving a collection of 1762 natural compounds, sourced from the PubChem database, was performed. This screening yielded 131 compounds that exhibited compatibility with the established pharmacophoric hypothesis. These selected ligands were then subjected to molecular docking within the active site of the 4EY5 receptor. As a result, we identified four compounds that displayed remarkable docking scores and exhibited low free binding energy to the target. These top four compounds, CID_162895946, CID_44461278, CID_44285285, and CID_81108419, were submitted to ADMET prediction and molecular dynamic simulations, yielding encouraging findings in terms of their pharmacokinetic characteristics and stability. Finally, the molecular dynamic simulation, cross-dynamic correlation matrix, free energy landscape, and MM-PBSA calculations demonstrated that two ligands from the selected ligands formed very resilient complexes with the enzyme acetylcholinesterase, with significant binding affinity. Therefore, these two compounds are recommended for further experimental research as possible (AChE) inhibitors.