Abstract Development of a pharmacophore hypothesis related to small-molecule activity is pivotal to chemical optimization of a series, since it defines features beneficial or detrimental to activity. Although crystal structures may provide detailed 3D interaction information for one molecule with its receptor, docking a different ligand to that model often leads to unreliable results due to protein flexibility. Graham Richards’ lab was one of the first groups to utilize “fuzzy” pattern recognition algorithms taken from the field of image processing to solve problems in protein modeling. Thus, descriptor “fuzziness” was partly able to emulate conformational flexibility of the target while simultaneously enhancing the speed of the search. In this work, we extend these developments to a ligand-based method for describing and aligning molecules in flexible chemical space termed FEature POint PharmacophoreS (FEPOPS), which allows exploration of dynamic biological space. We develop a novel, combinatorial algorithm for molecular comparisons and evaluate it using the WOMBAT dataset. The new approach shows superior retrospective virtual screening performance than earlier shape-based or charge-based algorithms. Additionally, we use target prediction to evaluate how FEPOPS alignments match the molecules biological activity by identifying the atoms and features that make the key contributions to overall chemical similarity. Overall, we find that FEPOPS are sufficiently fuzzy and flexible to find not only new ligand scaffolds, but also challenging molecules that occupy different conformational states of dynamic biological space as from induced fits.