Abstract
We present a new computational approach, named Watermelon, designed for the development of pharmacophore models based on receptor structures. The methodology involves the sampling of potential hotspots for ligand interactions within a protein target's binding site, utilising molecular fragments as probes. By employing docking and molecular dynamics (MD) simulations, the most significant interactions formed by these probes within distinct regions of the binding site are identified. These interactions are subsequently transformed into pharmacophore features that delineates key anchoring sites for potential ligands. The reliability of the approach was experimentally validated using the monoacylglycerol lipase (MAGL) enzyme. The generated pharmacophore model captured features representing ligand-MAGL interactions observed in various X-ray co-crystal structures and was employed to screen a database of commercially available compounds, in combination with consensus docking and MD simulations. The screening successfully identified two new MAGL inhibitors with micromolar potency, thus confirming the reliability of the Watermelon approach.