Sphingosine-1-phosphate receptor 1 (S1PR1), a member of the G protein-coupled receptor (GPCR) family, is a crucial therapeutic target for various diseases. Activation of S1PR1 has been recognized as an effective therapeutic strategy for multiple sclerosis (MS), inflammatory bowel disease (IBD), and psoriasis. Natural products (NPs) serve as a rich source of bioactive compounds for drug discovery. Here, we aimed to discover novel S1PR1 agonists from NPs via multi-level virtual screening (VS). Using a validated HipHop pharmacophore model, we screened a database containing 54,642 NPs, followed by molecular docking. Based on binding mode analysis, four candidate S1PR1 agonists (NPC323626, NPC264112, NPC469907, and NPC22192) were selected. Subsequent molecular dynamics (MD) simulations and binding free energy calculations confirmed the stability of the receptor-ligand complexes and their binding affinities. Among the four candidates, NPC469907 exhibited the strongest binding affinity for S1PR1, with a value of -58.08 ± 0.13 kJ/mol. Furthermore, hydrogen bonds formed between NPC469907 and Glu121 of S1PR1 were found to be essential for receptor activation. Quantum mechanical calculations further revealed that the phenyl-ring-attached hydrogen site in NPC469907 could be modified without compromising its ability to activate S1PR1. The analysis of Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADMET) indicated that NPC469907 possessed favorable pharmacokinetic properties and low toxicity. In conclusion, our study identified NPC469907 as a promising natural S1PR1 agonist and established an effective VS strategy for the discovery of novel S1PR1 agonists.