Casein kinase 1 epsilon (CK1ε) plays a critical role in cancer progression by activating oncogenic signaling pathways, making it a target for cancer therapy. However, no inhibitors are currently available for clinical use, highlighting the need for novel therapeutic candidates. This study aimed to identify potential CK1ε inhibitors. To achieve this, a modified version of a previously reported pharmacophore model was applied to an ultra-large database of over 100 million compounds for virtual screening. Hits were filtered based on drug-likeness and pH-dependent pharmacophore compliance and then grouped according to their structural core. A representative compound from each structural group underwent molecular dynamic (MD) simulations and binding free energy calculations to predict its stability and affinity, allowing extrapolation of the results to the entire set of candidates. Pharmacophore matching initially identified 290 compounds. After energy minimization, and an assessment of drug-likeness and pharmacophore compliance, we selected 29 structurally related candidates. MD simulations showed that most of the compounds representative of structural groups had stable binding modes, favorable intermolecular interactions, and free energies comparable to those of previously reported CK1ε inhibitors. An analysis of additional members of the most promising structural group showed that two 2,4-diaminopyrimidine-based compounds likely inhibit CK1ε. These findings provide structural insights into the design of CK1ε inhibitors, supporting compound optimization and the eventual development of targeted cancer therapeutics.